Protein surface salt bridges and paths for DNA wrapping

PRC2 direct transfer from G-quadruplex RNA to dsDNA has implications for RNA-binding chromatin modifiers

The chromatin-modifying enzyme, Polycomb Repressive Complex 2 (PRC2), deposits the H3K27me3 epigenetic mark to negatively regulate expression at numerous target genes, and this activity has been implicated in embryonic development, cell differentiation, and various cancers. A biological role for RNA binding in regulating PRC2 histone methyltransferase activity is generally accepted, but the nature and mechanism of this relationship remains an area of active investigation. Notably, many in vitro studies demonstrate that RNA inhibits PRC2 activity on nucleosomes through mutually antagonistic binding, while some in vivo studies indicate that PRC2's RNA-binding activity is critical for facilitating its biological function(s). Here we use biochemical, biophysical, and computational approaches to interrogate PRC2's RNA and DNA-binding kinetics. Our findings demonstrate that PRC2-polynucleotide dissociation rates are dependent on the concentration of free ligand, indicating the potential for direct transfer between nucleic acid ligands without a free-enzyme intermediate. Direct transfer explains the variation in previously reported dissociation kinetics, allows reconciliation of prior in vitro and in vivo studies, and expands the potential mechanisms of RNA-mediated PRC2 regulation. Moreover, simulations indicate that such a direct transfer mechanism could be obligatory for RNA to recruit proteins to chromatin.

Biological Nanopores: Engineering on Demand

Nanopore-based sensing is a powerful technique for the detection of diverse organic and inorganic molecules, long-read sequencing of nucleic acids, and single-molecule analyses of enzymatic reactions. Selected from natural sources, protein-based nanopores enable rapid, label-free detection of analytes. Furthermore, these proteins are easy to produce, form pores with defined sizes, and can be easily manipulated with standard molecular biology techniques. The range of possible analytes can be extended by using externally added adapter molecules. Here, we provide an overview of current nanopore applications with a focus on engineering strategies and solutions.

Designer Histone Complexes: Controlling Protein-DNA Interactions with Protein Charge as an "All-or-None" Digital Switch

An artificial histone is synthesized that functions as a DNA-protein digital switch, where DNA binding is all or none, controlled by a sharp threshold of protein charge. A non-DNA-binding protein, glucose oxidase (GOx), was chemically modified by attaching an increasing number of triethylenetetramine (TETA) side chains to its glutamate/aspartate groups to obtain a small library of covalently modified GOx(n) derivatives. The parameter n denotes the net charge on the protein at pH 7, which was increased from -62 (pristine GOx) to +75 by attaching an increasing number of TETA residues to the protein. All GOx(n) derivatives retained their secondary structure to a good extent, as monitored by UV circular dichroism (CD) spectroscopy, and they also retained oxidase activities to a significant extent. The interaction of the GOx(n) with calf thymus DNA was examined by isothermal titration calorimetry (ITC). Pristine GOx of -62 charge at pH 7 in 10 mM Tris-HCl and 50 mM NaCl buffer had no affinity for the negatively charged DNA helix, and GOx(n) with n < +30 had no affinity for DNA either. However, binding has been turned on abruptly when n ≥ +30 with binding constants (K_b) ranging from (1.5 ± 0.7) × 10⁷ to (7.3 ± 2.8) × 10⁷ M^-1 for n values of +30 and +75, respectively, and this type of "all-or-none" binding based on protein charge is intriguing. Furthermore, thermodynamic analysis of the titration data revealed that binding is entirely entropy-driven with ΔS ranging from 0.09 ± 0.007 to 0.19 ± 0.008 kcal/mol K with enthalpic penalties of 17.0 ± 2.3 and 46.1 ± 2.1 kcal/mol, respectively. The binding had intrinsic propensities (ΔG) ranging from -9.8 ± 0.14 to -10.7 ± 0.25 kcal/mol, independent of n. DNA binding distorted protein-DNA secondary structure, as evidenced by CD spectroscopy, but oxidase activity of GOx(n)/DNA complexes has been unaffected. This is the very first example of an artificial histone (GOx(n)) where the protein charge functioned as a DNA-binding switch; protein charge is in turn under complete chemical control while preserving the biological activity of the protein. The new insight gained here could be useful in the design of novel "on-off" protein switches.

Caenorhabditis elegans NONO-1: Insights into DBHS protein structure, architecture, and function

Members of the Drosophila behavior/human splicing (DBHS) protein family have been characterized in the vertebrates Homo sapiens and Mus musculus, and the invertebrates Drosophila melanogaster and Chironomus tentans. Collectively, both vertebrate and invertebrate DBHS proteins function throughout gene regulation, largely but not always, within the nucleus. In this study, we report a structural and bioinformatic analysis of the DBHS protein family to guide future studies into DBHS protein function. To explore the structural plasticity of the family, we describe the 2.4 Å crystal structure of Caenorhabditis elegans non-POU domain-containing octamer-binding protein 1 (NONO-1). The structure is dimeric, with a domain arrangement consistent with mammalian DBHS proteins. Comparison with the DBHS structures available from H. sapiens reveals that there is inherent domain flexibility within the homologous DBHS region. Mapping amino acid similarity within the family to the NONO-1 dimer highlights the dimer interface, coiled-coil oligomerization motif, and putative RNA binding surfaces. Surprisingly, the interior surface of RNA recognition motif 2 (RRM2) that faces a large internal void is highly variable, but the external β2-β3 loops of RRM2 show remarkable preservation. Overall, the DBHS region is under strong purifying selection, whereas the sequences N- and C-terminal to the DBHS region are less constrained. The findings described in this study provide a molecular basis for further investigation into the mechanistic function of the DBHS protein family in biology.

Investigating the allosteric reverse signalling of PARP inhibitors with microsecond molecular dynamic simulations and fluorescence anisotropy

The inhibition of the poly(ADP-ribose) polymerase (PARP) family members is a strategy pursued for the development of novel therapeutic agents in a range of diseases, including stroke, cardiac ischemia, cancer, inflammation and diabetes. Even though some PARP-1 inhibitors have advanced to clinical setting for cancer therapy, a great deal of attention is being devoted to understand the polypharmacology of current PARP inhibitors. Besides blocking the catalytic activity, recent works have shown that some PARP inhibitors exhibit a poisoning activity, by trapping the enzyme at damaged sites of DNA and forming cytotoxic complexes. In this study we have used microsecond molecular dynamics to study the allosteric reverse signalling that is at the basis of such an effect. We show that Olaparib, but not Veliparib and HYDAMTIQ, is able to induce a specific conformational drift of the WGR domain of PARP-1, which stabilizes PARP-1/DNA complex through the locking of several salt bridge interactions. Fluorescence anisotropy assays support such a mechanism, providing the first experimental evidence that HYDAMTIQ, a potent PARP inhibitor with neuroprotective properties, is less potent than Olaparib to trap PARP-1/DNA complex.

Wrapping/unwrapping transition of double-stranded DNA in DNA-nanosphere complexes induced by multivalent anions

Wrapping and unwrapping behaviors of double-stranded DNA around a positively charged nanosphere in solution are studied by using the coarse-grained molecular dynamics (CGMD) simulation method. When monovalent, divalent and trivalent anions are added to the DNA-nanosphere complex solution, double-stranded DNA binds with a nanosphere owing to strong electrostatic attraction. However, when tetravalent anions are added to the DNA-nanosphere complex solution, local charge inversion is observed for a high anion concentration of tetravalent anions and the double-stranded DNA can be unwrapped from the nanosphere because of the local charge inversion near the nanosphere. Moreover, the helical structure of DNA is damaged when double-stranded DNA wraps around the nanosphere and the helical structure can be rebuilt when the double-stranded DNA unwraps from the nanosphere. This study can help us understand how to control the release of DNA in DNA-nanosphere complexes.

A simple but effective modeling strategy for structural properties of non-heme Fe(II) sites in proteins: test of force field models and application to proteins in the AlkB family

To facilitate computational study of proteins in the AlkB family and related α-ketoglutarate/Fe(II)-dependent dioxygenases, we have tested a simple modeling strategy for the non-heme Fe(II) site in which the iron is represented by a simple +2 point charge with Lennard-Jones parameters. Calculations for an AlkB active site model in the gas phase and ∼150 ns molecular dynamics (MD) simulations for two enzyme-dsDNA complexes (E. coli AlkB-dsDNA and ABH2-dsDNA) suggest that this simple modeling strategy provides a satisfactory description of structural properties of the Fe(II) site in AlkB enzymes, provided that care is exercised to control the binding mode of carboxylate (Asp) to the iron. MD simulations using the model for AlkB-dsDNA and ABH2-dsDNA systems find that although the structural features for the latter are overall in good agreement with the crystal structure, the dsDNA, and AlkB-dsDNA interface undergo substantial changes during the MD simulations from the crystal structure. Even for ABH2, new interactions form between a long loop region and dsDNA upon structural relaxation of the loop, supporting the role of this loop in DNA binding despite the lack of interactions between them in the crystal structure. Analysis of DNA backbone torsional distributions helps identify regions that adopt strained conformations. Collectively, the results highlight that crystal packing may have a significant impact on the structure of protein-DNA complexes; the simulations also provide additional insights regarding why AlkB and ABH2 prefer single-strand and double-strand DNA, respectively, as substrate.

DNA conformation and energy in nucleosome core: a theoretical approach

DNA conformation in complex with proteins is far from its canonical B-form. The affinity of complex formation and structure of DNA depend on its attachment configuration and sequence. In this article, we develop a mechanical model to address the problem of DNA structure and energy under deformation. DNA in nucleosome core particle is described as an example. The structure and energy of nucleosomal DNA is calculated based on its sequence and positioning state. The inferred structure has remarkable similarity with X-ray data. Although there is no sequence-specific interaction of bases and the histone core, we found considerable sequence dependency for the nucleosomal DNA positioning. The affinity of nucleosome formation for several sequences is examined and the differences are compatible with observations. We argue that structural energy determines the natural state of nucleosomal DNA and is the main reason for affinity differences in vitro. This theory can be utilized for the DNA structure and energy determination in protein-DNA complexes in general. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:17.

Water-mediated ionic interactions in protein structures

It is well known that water molecules play an indispensable role in the structure and function of biological macromolecules. The water-mediated ionic interactions between the charged residues provide stability and plasticity and in turn address the function of the protein structures. Thus, this study specifically addresses the number of possible water-mediated ionic interactions, their occurrence, distribution and nature found in 90% non-redundant protein chains. Further, it provides a statistical report of different charged residue pairs that are mediated by surface or buried water molecules to form the interactions. Also, it discusses its contributions in stabilizing various secondary structural elements of the protein. Thus, the present study shows the ubiquitous nature of the interactions that imparts plasticity and flexibility to a protein molecule.

GMI1, a structural-maintenance-of-chromosomes-hinge domain-containing protein, is involved in somatic homologous recombination in Arabidopsis

DNA double-strand breaks (DSBs) pose one of the most severe threats to genome integrity, potentially leading to cell death. After detection of a DSB, the DNA damage and repair response is initiated and the DSB is repaired by non-homologous end joining and/or homologous recombination. Many components of these processes are still unknown in Arabidopsis thaliana. In this work, we characterized γ-irradiation and mitomycin C induced 1 (GMI1), a member of the SMC-hinge domain-containing protein family. RT-PCR analysis and promoter-GUS fusion studies showed that γ-irradiation, the radio-mimetic drug bleocin, and the DNA cross-linking agent mitomycin C strongly enhance GMI1 expression particularly in meristematic tissues. The induction of GMI1 by γ-irradiation depends on the signalling kinase Ataxia telangiectasia-mutated (ATM) but not on ATM and Rad3-related (ATR). Epistasis analysis of single and double mutants demonstrated that ATM acts upstream of GMI1 while the atr gmi1-2 double mutant was more sensitive than the respective single mutants. Comet assay revealed a reduced rate of DNA double-strand break repair in gmi1 mutants during the early recovery phase after exposure to bleocin. Moreover, the rate of homologous recombination of a reporter construct was strongly reduced in gmi1 mutant plants upon exposure to bleocin or mitomycin C. GMI1 is the first member of its protein family known to be involved in DNA repair.

Structural and thermodynamic signatures of DNA recognition by Mycobacterium tuberculosis DnaA

An essential protein, DnaA, binds to 9-bp DNA sites within the origin of replication oriC. These binding events are prerequisite to forming an enigmatic nucleoprotein scaffold that initiates replication. The number, sequences, positions, and orientations of these short DNA sites, or DnaA boxes, within the oriCs of different bacteria vary considerably. To investigate features of DnaA boxes that are important for binding Mycobacterium tuberculosis DnaA (MtDnaA), we have determined the crystal structures of the DNA binding domain (DBD) of MtDnaA bound to a cognate MtDnaA-box (at 2.0 Å resolution) and to a consensus Escherichia coli DnaA-box (at 2.3 Å). These structures, complemented by calorimetric equilibrium binding studies of MtDnaA DBD in a series of DnaA-box variants, reveal the main determinants of DNA recognition and establish the [T/C][T/A][G/A]TCCACA sequence as a high-affinity MtDnaA-box. Bioinformatic and calorimetric analyses indicate that DnaA-box sequences in mycobacterial oriCs generally differ from the optimal binding sequence. This sequence variation occurs commonly at the first 2 bp, making an in vivo mycobacterial DnaA-box effectively a 7-mer and not a 9-mer. We demonstrate that the decrease in the affinity of these MtDnaA-box variants for MtDnaA DBD relative to that of the highest-affinity box TTGTCCACA is less than 10-fold. The understanding of DnaA-box recognition by MtDnaA and E. coli DnaA enables one to map DnaA-box sequences in the genomes of M. tuberculosis and other eubacteria.

Understanding the physical basis of the salt dependence of the electrostatic binding free energy of mutated charged ligand-nucleic acid complexes

The predictions of the derivative of the electrostatic binding free energy of a biomolecular complex, ΔG(el), with respect to the logarithm of the 1:1 salt concentration, d(ΔG(el))/d(ln[NaCl]), SK, by the Poisson-Boltzmann equation, PBE, are very similar to those of the simpler Debye-Hückel equation, DHE, because the terms in the PBE's predictions of SK that depend on the details of the dielectric interface are small compared to the contributions from long-range electrostatic interactions. These facts allow one to obtain predictions of SK using a simplified charge model along with the DHE that are highly correlated with both the PBE and experimental binding data. The DHE-based model developed here, which was derived from the generalized Born model, explains the lack of correlation between SK and ΔG(el) in the presence of a dielectric discontinuity, which conflicts with the popular use of this supposed correlation to parse experimental binding free energies into electrostatic and nonelectrostatic components. Moreover, the DHE model also provides a clear justification for the correlations between SK and various empirical quantities, like the number of ion pairs, the ligand charge on the interface, the Coulomb binding free energy, and the product of the charges on the complex's components, but these correlations are weak, questioning their usefulness.

Geometrical correlations in the nucleosomal DNA conformation and the role of the covalent bonds rigidity

We develop a simple elastic model to study the conformation of DNA in the nucleosome core particle. In this model, the changes in the energy of the covalent bonds that connect the base pairs of each strand of the DNA double helix, as well as the lateral displacements and the rotation of adjacent base pairs are considered. We show that because of the rigidity of the covalent bonds in the sugar-phosphate backbones, the base pair parameters are highly correlated, especially, strong twist-roll-slide correlation in the conformation of the nucleosomal DNA is vividly observed in the calculated results. This simple model succeeds to account for the detailed features of the structure of the nucleosomal DNA, particularly, its more important base pair parameters, roll and slide, in good agreement with the experimental results.

Disruption and formation of surface salt bridges are coupled to DNA binding by the integration host factor: a computational analysis

Revealing the thermodynamic driving force of protein-DNA interactions is crucial to the understanding of factors that dictate the properties and function of protein-DNA complexes. For the binding of DNA to DNA-wrapping proteins, such as the integration host factor (IHF), Record and co-workers proposed that the disruption of a large number of preexisting salt bridges is coupled with the binding process [Holbrook, J. A., et al. (2001) J. Mol. Biol. 310, 379]. To test this proposal, we have conducted explicit solvent MD simulations (multiple ∼25-50 ns trajectories for each salt concentration) to examine the behavior of charged residues in IHF, especially concerning their ability to form salt bridges at different salt concentrations. Of the 17 cationic residues noted by Record and co-workers, most are engaged in salt bridge interactions for a significant portion of the trajectories, especially in the absence of salt. This observation suggests that, from a structural point of view, their proposal is plausible. However, the complex behaviors of charged residues observed in the MD simulations also suggest that the unusual thermodynamic characteristics of IHF-DNA binding likely arise from the interplay between complex dynamics of charged residues both in and beyond the DNA binding site. Moreover, a comparison of MD simulations at different salt concentrations suggests that the strong dependence of the IHF-DNA binding enthalpy on salt concentration may not be due to a significant decrease in the number of stable salt bridges in apo IHF at high salt concentrations. In addition to the Hofmeister effects quantified in more recent studies of IHF-DNA binding, we recommend consideration of the variation of the enthalpy change of salt bridge disruption at different salt concentrations. Finally, the simulation study presented here explicitly highlights the fact that the electrostatic properties of DNA-binding proteins can be rather different in the apo and DNA-bound states, which has important implications for the design of robust methods for predicting DNA binding sites in proteins.

Charge effects in the selection of NPF motifs by the EH domain of EHD1

The Eps15 homology (EH) domain is found in proteins associated with endocytosis and vesicle trafficking. EH domains bind to their target proteins through an asparagine-proline-phenylalanine (NPF) motif. We have measured the interaction energetics of the EH domain from EHD1 with peptides derived from two of its binding partners: Rabenosyn-5 (Ac-GPSLNPFDEED-NH(2)) and Rab11-Fip2 (Ac-YESTNPFTAK-NH(2)). Heteronuclear single quantum coherence (HSQC) spectroscopy shows that both peptides bind in the canonical binding pocket of EHD1 EH and induce identical structural changes, yet the affinity of the negatively charged Ac-GPSLNPFDEED-NH(2) (K(a) = 8 x 10(5) M(-1)) is tighter by 2 orders of magnitude. The thermodynamic profiles (DeltaG, DeltaH, DeltaS) were measured for both peptides as a function of temperature. The enthalpies of binding are essentially identical, and the difference in affinity is a consequence of the difference in entropic cost. Ac-GPSLNPFDEED-NH(2) binding is salt-dependent, demonstrating an electrostatic component to the interaction, whereas Ac-YESTNPFTAK-NH(2) binding is independent of salt. Successive replacement of acidic residues in Ac-GPSLNPFDEED-NH(2) with neutral residues showed that all are important. Lysine side chains in EHD1 EH create a region of strong positive surface potential near the NPF binding pocket. Contributions by lysine epsilon-amino groups to complex formation with Ac-GPSLNPFDEED-NH(2) was shown using direct-observe (15)N NMR spectroscopy. These experiments have enabled us to define a new extended interaction motif for EHD proteins, N-P-F-[DE]-[DE]-[DE], which we have used to predict new interaction partners and hence broaden the range of cellular activities involving the EHD proteins.

Probing DNA binding, DNA opening, and assembly of a downstream clamp/jaw in Escherichia coli RNA polymerase-lambdaP(R) promoter complexes using salt and the physiological anion glutamate

Transcription by all RNA polymerases (RNAPs) requires a series of large-scale conformational changes to form the transcriptionally competent open complex RP(o). At the lambdaP(R) promoter, Escherichia coli sigma(70) RNAP first forms a wrapped, closed 100 bp complex I(1). The subsequent step opens the entire DNA bubble, creating the relatively unstable (open) complex I(2). Additional conformational changes convert I(2) to the stable RP(o). Here we probe these events by dissecting the effects of Na(+) salts of Glu(-), F(-), and Cl(-) on each step in this critical process. Rapid mixing and nitrocellulose filter binding reveal that the binding constant for I(1) at 25 degrees C is approximately 30-fold larger in Glu(-) than in Cl(-) at the same Na(+) concentration, with the same log-log salt concentration dependence for both anions. In contrast, both the rate constant and equilibrium constant for DNA opening (I(1) to I(2)) are only weakly dependent on salt concentration, and the opening rate constant is insensitive to replacement of Cl(-) with Glu(-). These very small effects of salt concentration on a process (DNA opening) that is strongly dependent on salt concentration in solution may indicate that the backbones of both DNA strands interact with polymerase throughout the process and/or that compensation is present between ion uptake and release. Replacement of Cl(-) with Glu(-) or F(-) at 25 degrees C greatly increases the lifetime of RP(o) and greatly reduces its salt concentration dependence. By analogy to Hofmeister salt effects on protein folding, we propose that the excluded anions Glu(-) and F(-) drive the folding and assembly of the RNAP clamp/jaw domains in the conversion of I(2) to RP(o), while Cl(-) does not. Because the Hofmeister effect of Glu(-) or F(-) largely compensates for the destabilizing Coulombic effect of any salt on the binding of this assembly to downstream promoter DNA, RP(o) remains long-lived even at 0.5 M Na(+) in Glu(-) or F(-) salts. The observation that Esigma(70) RP(o) complexes are exceedingly long-lived at moderate to high Glu(-) concentrations argues that Esigma(70) RNAP does not dissociate from strong promoters in vivo when the cytoplasmic glutamate concentration increases during osmotic stress.

One-step DNA melting in the RNA polymerase cleft opens the initiation bubble to form an unstable open complex

Though opening of the start site (+1) region of promoter DNA is required for transcription by RNA polymerase (RNAP), surprisingly little is known about how and when this occurs in the mechanism. Early events at the lambdaP(R) promoter load this region of duplex DNA into the active site cleft of Escherichia coli RNAP, forming the closed, permanganate-unreactive intermediate I(1). Conversion to the subsequent intermediate I(2) overcomes a large enthalpic barrier. Is I(2) open? Here we create a burst of I(2) by rapidly destabilizing open complexes (RP(o)) with 1.1 M NaCl. Fast footprinting reveals that thymines at positions from -11 to +2 in I(2) are permanganate-reactive, demonstrating that RNAP opens the entire initiation bubble in the cleft in a single step. Rates of decay of all observed thymine reactivities are the same as the I(2) to I(1) conversion rate determined by filter binding. In I(2), permanganate reactivity of the +1 thymine on the template (t) strand is the same as the RP(o) control, whereas nontemplate (nt) thymines are significantly less reactive than in RP(o). We propose that: (i) the +1(t) thymine is in the active site in I(2); (ii) conversion of I(2) to RP(o) repositions the nt strand in the cleft; and (iii) movements of the nt strand are coupled to the assembly and DNA binding of the downstream clamp and jaw that occurs after DNA opening and stabilizes RP(o). We hypothesize that unstable open intermediates at the lambdaP(R) promoter resemble the unstable, transcriptionally competent open complexes formed at ribosomal promoters.

Electrostatic interactions guide the active site face of a structure-specific ribonuclease to its RNA substrate

Restrictocin, a member of the α-sarcin family of site-specific endoribonucleases, uses electrostatic interactions to bind to the ribosome and to RNA oligonucleotides, including the minimal specific substrate, the sarcin/ricin loop (SRL) of 23S−28S rRNA. Restrictocin binds to the SRL by forming a ground-state E:S complex that is stabilized predominantly by Coulomb interactions and depends on neither the sequence nor structure of the RNA, suggesting a nonspecific complex. The 22 cationic residues of restrictocin are dispersed throughout this protein surface, complicating a priori identification of a Coulomb interacting surface. Structural studies have identified an enzyme−substrate interface, which is expected to overlap with the electrostatic E:S interface. Here, we identified restrictocin residues that contribute to binding in the E:S complex by determining the salt dependence [∂ log(k₂/K_1/2)/∂ log[KCl]] of cleavage of the minimal SRL substrate for eight point mutants within the protein designed to disrupt contacts in the crystallographically defined interface. Relative to the wild-type salt dependence of −4.1, a subset of the mutants clustering near the active site shows significant changes in salt dependence, with differences of magnitude being ≥0.4. This same subset was identified using calculated salt dependencies for each mutant derived from solutions to the nonlinear Poisson−Boltzmann equation. Our findings support a mechanism in which specific residues on the active site face of restrictocin (primarily K110, K111, and K113) contribute to formation of the E:S complex, thereby positioning the SRL substrate for site-specific cleavage. The same restrictocin residues are expected to facilitate targeting of the SRL on the surface of the ribosome.

Cooperative DNA binding and communication across the dimer interface in the TREX2 3' --> 5'-exonuclease

The activity of human TREX2-catalyzed 3' --> 5'-deoxyribonuclease has been analyzed in steady-state and single turnover kinetic assays and in equilibrium DNA binding studies. These kinetic data provide evidence for cooperative DNA binding within TREX2 and for coordinated catalysis between the TREX2 active sites supporting a model for communication between the protomers of a TREX2 dimer. Mobile loops positioned adjacent to the active sites provide the major DNA binding contribution and facilitate subsequent binding into the active sites. Mutations of three arginine residues on these loops cause decreased TREX2 activities by up to 60-fold. Steady-state kinetic assays of these arginine to alanine TREX2 variants result in increased K(m) values for DNA substrate with no effect on k(cat) values indicating contributions exclusively to DNA binding by all three of the loop arginines. TREX2 heterodimers were prepared to determine whether exonuclease activity in one protomer is communicated to the opposing protomer. Evidence for communication across the dimer interface is provided by the 7-fold lower catalytic activity measured in the TREX2(WT/H188A) heterodimer compared with the TREX2(WT) homodimer, contrasting the 2-fold lower activity measured in the TREX2(WT/R163A,R165A,R167A) heterodimer. The measured activity in TREX2(WT/H188A) heterodimer indicates that defective catalysis in one protomer reduces activity in the opposing protomer. A DNA binding analysis of TREX2 and the heterodimers indicates a cooperative binding effect within the TREX2 protomer. Finally, single turnover kinetic assays identify DNA binding as the rate-limiting step in TREX2 catalysis.

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.

The salt dependence of the interferon regulatory factor 1 DNA binding domain binding to DNA reveals ions are localized around protein and DNA

The equilibrium dissociation constant of the DNA binding domain of interferon regulatory factor 1 (IRF1 DBD) for its DNA binding site depends strongly on salt concentration and salt type. These dependencies are consistent with IRF1 DBD binding to DNA, resulting in the release of cations from the DNA and both release of anions from the protein and uptake of a cation by the protein. We demonstrated this by utilizing the fact that the release of fluoride from protein upon complex formation does not contribute to the salt concentration dependence of binding and by studying mutants in which charged residues in IRF1 DBD that form salt bridges with DNA phosphates are changed to alanine. The salt concentration dependencies of the dissociation constants of wild-type IRF1 DBD and the mutants R64A, D73A, K75A, and D73A/K75A were measured in buffer containing NaF, NaCl, or NaBr. The salt concentration and type dependencies of the mutants relative to wild-type IRF1 DBD provide evidence of charge neutralization by solution ions for R64 and by a salt bridge between D73 and K75 in buffer containing chloride or bromide salts. These data also allowed us to determine the number, type, and localization of condensed ions around both IRF1 DBD and its DNA binding site.

Molecular beacon-equilibrium cyclization detection of DNA-protein complexes

Molecular beacon detection of equilibrium cyclization (MBEC) is a novel, high sensitivity technique that can allow DNA-protein complex formation to be studied under diverse conditions in a cost effective and rapid manner that can be adapted to high throughput screening. To demonstrate the ease and utility of applying MBEC to the investigation of the K(D) values of protein-DNA complexes, the sequence-specific Escherichia coli integration host factor (IHF) protein has been used as a test system. Competition between a labeled MBEC DNA construct and unlabeled duplex DNA for IHF binding allows the determination of K(D) values as a function of the DNA duplex sequence. This allows sequence specificity to be monitored while using only a single molecular beacon-labeled DNA. The robustness of MBEC for monitoring protein-DNA complex formation has been further demonstrated by determining the K(D) values as a function of salt concentration to investigate the net number of salt bridges formed in sequence-specific and -nonspecific IHF-DNA complexes. These MBEC results have been compared with those from other approaches.

Real-time footprinting of DNA in the first kinetically significant intermediate in open complex formation by Escherichia coli RNA polymerase

The architecture of cellular RNA polymerases (RNAPs) dictates that transcription can begin only after promoter DNA bends into a deep channel and the start site nucleotide (+1) binds in the active site located on the channel floor. Formation of this transcriptionally competent "open" complex (RP(o)) by Escherichia coli RNAP at the lambdaP(R) promoter is greatly accelerated by DNA upstream of base pair -47 (with respect to +1). Here we report real-time hydroxyl radical (*OH) and potassium permanganate (KMnO4) footprints obtained under conditions selected for optimal characterization of the first kinetically significant intermediate (I(1)) in RP(o) formation. .OH footprints reveal that the DNA backbone from -71 to -81 is engulfed by RNAP in I(1) but not in RP(o); downstream protection extends to approximately +20 in both complexes. KMnO4 footprinting detects solvent-accessible thymine bases in RP(o), but not in I(1). We conclude that upstream DNA wraps more extensively on RNAP in I(1) than in RP(o) and that downstream DNA (-11 to +20) occupies the active-site channel in I(1) but is not yet melted. Mapping of the footprinting data onto available x-ray structures provides a detailed model of a kinetic intermediate in bacterial transcription initiation and suggests how transient contacts with upstream DNA in I(1) might rearrange the channel to favor entry of downstream duplex DNA.

Pressure dissociation of integration host factor-DNA complexes reveals flexibility-dependent structural variation at the protein-DNA interface

E. coli Integration host factor (IHF) condenses the bacterial nucleoid by wrapping DNA. Previously, we showed that DNA flexibility compensates for structural characteristics of the four consensus recognition elements associated with specific binding (Aeling et al., J. Biol. Chem. 281, 39236-39248, 2006). If elements are missing, high-affinity binding occurs only if DNA deformation energy is low. In contrast, if all elements are present, net binding energy is unaffected by deformation energy. We tested two hypotheses for this observation: in complexes containing all elements, (1) stiff DNA sequences are less bent upon binding IHF than flexible ones; or (2) DNA sequences with differing flexibility have interactions with IHF that compensate for unfavorable deformation energy. Time-resolved Förster resonance energy transfer (FRET) shows that global topologies are indistinguishable for three complexes with oligonucleotides of different flexibility. However, pressure perturbation shows that the volume change upon binding is smaller with increasing flexibility. We interpret these results in the context of Record and coworker's model for IHF binding (J. Mol. Biol. 310, 379-401, 2001). We propose that the volume changes reflect differences in hydration that arise from structural variation at IHF-DNA interfaces while the resulting energetic compensation maintains the same net binding energy.

Physicochemical analysis of electrostatic foundation for DNA-protein interactions in chromatin transformations

Electrostatic interactions between DNA and DNA-packaging proteins, the histones, contribute substantially to stability of eukaryotic chromatin on all levels of its organization and are particularly important in formation of its elementary structural unit, the nucleosome. The release of DNA from the histones is an unavoidable stage in reading the DNA code. In the present review, we discuss the disassembly/assembly process of the nucleosome from a thermodynamic standpoint by considering it as a competition between an excess of polyanions (DNA and acidic/phosphorylated domains of the nuclear proteins) for binding to a limited pool of polycations (the histones). Results obtained in model systems are used to discuss conditions for the electrostatic component of DNA-protein interactions contributing to chromatin statics and dynamics. We propose a simple set of "electrostatic conditions" for the disassembly/assembly of nucleosome/chromatin and apply these to put forward a number of new interpretations for the observations reported in literature on chromatin. The approach sheds light on the functions of acidic domains in the nuclear proteins (nucleoplasmin and other histone chaperones, HMG proteins, the activation domains in transcriptional activators). It results in a putative explanation for the molecular mechanisms behind epigenetic regulation through histone acetylation, phosphorylation, and other alterations ("the language of covalent histone modification"). We also propose a new explanation for the role of phosphorylation of C-terminal domain of RNA polymerase II for regulation of the DNA transcription. Several other examples from literature on chromatin are discussed to support applicability of electrostatic rules for description of chromatin structure and dynamics.

The temperature dependence of salt-protein association is sequence specific

Molecular dynamics (MD) simulations are used to probe the origin of the unexpected temperature dependence of salt accumulation in the C-terminal region of the protein human lymphotactin. As in previous MD simulations, sodium ions accumulate in an enhanced manner near the C-terminal helix at the lower temperature, while the temperature dependence of chloride accumulation is much weaker and slightly positive. In a designed mutant in which all positively charged residues in the C-terminal helix are replaced with neutral polar groups (Ser), the unexpected temperature dependence of the sodium ions is no longer observed. Therefore, these simulations convincingly verified the previous hypothesis that the temperature dependence of ion-protein association is sensitive to the local sequence. This is explained qualitatively in terms of the entropy of association between charged species in solution. These findings have general implications for the interpretation of thermodynamic quantities associated with binding events where ion release is important, such as protein-DNA interactions.

Influence of the nature of quantum dot surface cations on interactions with DNA

Quantum dots are semiconductor nanoparticles that are approximately 1-10nm in diameter, similar to small proteins, and their photoluminescence is sensitive to the presence and nature of adsorbates. We have deployed these nanomaterials as luminescent probes of DNA structure. Sequence dependent conformational flexibility of DNA is of great interest due to its implications for drug-DNA and DNA-protein interactions. The counterion atmosphere surrounding DNA plays an important role in its structure, dynamics, and packaging. In this paper, we investigate the effect that various monovalent and divalent cations have on the binding of 4.5 nm CdS quantum dots to oligonucleotides that have sequence-directed intrinsic structure.

Indirect recognition in sequence-specific DNA binding by Escherichia coli integration host factor: the role of DNA deformation energy

Integration host factor (IHF) is a bacterial histone-like protein whose primary biological role is to condense the bacterial nucleoid and to constrain DNA supercoils. It does so by binding in a sequence-independent manner throughout the genome. However, unlike other structurally related bacterial histone-like proteins, IHF has evolved a sequence-dependent, high affinity DNA-binding motif. The high affinity binding sites are important for the regulation of a wide range of cellular processes. A remarkable feature of IHF is that it employs an indirect readout mechanism to bind and wrap DNA at both the nonspecific and high affinity (sequence-dependent) DNA sites. In this study we assessed the contributions of pre-formed and protein-induced DNA conformations to the energetics of IHF binding. Binding energies determined experimentally were compared with energies predicted for the IHF-induced deformation of the DNA helix (DNA deformation energy) in the IHF-DNA complex. Combinatorial sets of de novo DNA sequences were designed to systematically evaluate the influence of sequence-dependent structural characteristics of the conserved IHF recognition elements of the consensus DNA sequence. We show that IHF recognizes pre-formed conformational characteristics of the consensus DNA sequence at high affinity sites, whereas at all other sites relative affinity is determined by the deformational energy required for nearest-neighbor base pairs to adopt the DNA structure of the bound DNA-IHF complex.

Effects of temperature and salt concentration on the structural stability of human lymphotactin: insights from molecular simulations

Extensive molecular dynamics (MD) simulations ( approximately 70 ns total) with explicit solvent molecules and salt ions are carried out to probe the effects of temperature and salt concentration on the structural stability of the human Lymphotactin (hLtn). The distribution of ions near the protein surface and the stability of various structural motifs are observed to exhibit interesting dependence on the local sequence and structure. Whereas chloride association to the protein is overall enhanced as the temperature increases, the sodium distribution in the C-terminal helical region and, to a smaller degree, the chloride distribution in the same region are found higher at the lower temperature. The similar trend is also observed in nonlinear Poisson-Boltzmann calculations with a temperature-dependent water dielectric constant, once conformational averaging over a series of MD snapshots is done. The unexpected temperature dependence in the ion distribution is explained on the basis of the cancellation of association entropy for ion-side chain pairs of opposite-charge and like-charge characters, which have positive and negative contributions, respectively. The C-terminal helix is observed to partially melt whereas a short beta strand forms at the higher temperature with little salt dependence. The N-terminal region, by contrast, develops partial helical structure at a higher salt concentration. These observed behaviors are consistent with solvent and salt screening playing an important role in stabilizing the canonical chemokine fold of hLtn.

Compaction of DNA on nanoscale three-dimensional templates

There exist several important in vivo examples, where a DNA chain is compacted on interacting with nanoscale objects such as proteins, thereby forming complexes with a well defined molecular architecture. One of the well known manifestations of such a natural organization of a semi-flexible DNA chain on nanoscale objects is hierarchical DNA molecule assembly into a chromosome, which is mediated by cationic histone proteins at the first stages of compaction. The biological importance of this and other natural nanostructural organizations of the DNA molecule has inspired many theoretical and numerical studies to gain physical insight into this problem. On the other hand, the experimental model systems containing DNA and nanoobjects, which are important to extend our knowledge beyond natural systems, were almost unavailable until the last decade. Accelerating progress in nanoscale chemistry and materials science has brought about various nanoscale three-dimensional structures such as dendrimers, nanoparticles, and nanotubes, and thus has provided a basis for the next important step in creating novel DNA-containing nanostructures, modelling of natural DNA compaction, and verification of accumulated theoretical predictions on the interaction between DNA and nanoscale templates. This review is written to highlight this early stage of nano-inspired progress and it is focused on physico-chemical and biophysical experimental investigations as well as theoretical and numerical studies dedicated to the compaction of DNA on nanoscale three-dimensional templates.

Cation binding linked to a sequence-specific CAP-DNA interaction

The equilibrium association constant observed for many DNA-protein interactions in vitro (K(obs)) is strongly dependent on the salt concentration of the reaction buffer ([MX]). This dependence is often used to estimate the number of ionic contacts between protein and DNA by assuming that release of cations from the DNA is the dominant involvement of ions in the binding reaction. With this assumption, the graph of logK(obs) versus log[MX] is predicted to have a constant slope proportional to the number of ions released from the DNA upon protein binding. However, experimental data often deviate from log-linearity at low salt concentrations. Here we show that for the sequence-specific interaction of CAP with its primary site in the lactose promoter, ionic stoichiometries depend strongly on cation identity and weakly on anion identity. This outcome is consistent with a simple linkage model in which cation binding by the protein accompanies its association with DNA. The order of ion affinities deduced from analysis of DNA binding is the same as that inferred from urea-denaturation experiments performed in the absence of DNA, suggesting that ion binding to free CAP contributes significantly to the ionic stoichiometry of DNA binding. In living cells, the coupling of ion-uptake and DNA binding mechanisms could reduce the sensitivity of gene-regulatory interactions to changes in environmental salt concentration.

Effects of monovalent anions on a temperature-dependent heat capacity change for Escherichia coli SSB tetramer binding to single-stranded DNA

We have previously shown that the linkage of temperature-dependent protonation and DNA base unstacking equilibria contribute significantly to both the negative enthalpy change (DeltaH(obs)) and the negative heat capacity change (DeltaC(p,obs)) for Escherichia coli SSB homotetramer binding to single-stranded (ss) DNA. Using isothermal titration calorimetry we have now examined DeltaH(obs) over a much wider temperature range (5-60 degrees C) and as a function of monovalent salt concentration and type for SSB binding to (dT)(70) under solution conditions that favor the fully wrapped (SSB)(65) complex (monovalent salt concentration >or=0.20 M). Over this wider temperature range we observe a strongly temperature-dependent DeltaC(p,obs). The DeltaH(obs) decreases as temperature increases from 5 to 35 degrees C (DeltaC(p,obs) <0) but then increases at higher temperatures up to 60 degrees C (DeltaC(p,obs) >0). Both salt concentration and anion type have large effects on DeltaH(obs) and DeltaC(p,obs). These observations can be explained by a model in which SSB protein can undergo a temperature- and salt-dependent conformational transition (below 35 degrees C), the midpoint of which shifts to higher temperature (above 35 degrees C) for SSB bound to ssDNA. Anions bind weakly to free SSB, with the preference Br(-) > Cl(-) > F(-), and these anions are then released upon binding ssDNA, affecting both DeltaH(obs) and DeltaC(p,obs). We conclude that the experimentally measured values of DeltaC(p,obs) for SSB binding to ssDNA cannot be explained solely on the basis of changes in accessible surface area (ASA) upon complex formation but rather result from a series of temperature-dependent equilibria (ion binding, protonation, and protein conformational changes) that are coupled to the SSB-ssDNA binding equilibrium. This is also likely true for many other protein-nucleic acid interactions.

Solute probes of conformational changes in open complex (RPo) formation by Escherichia coli RNA polymerase at the lambdaPR promoter: evidence for unmasking of the active site in the isomerization step and for large-scale coupled folding in the subsequent conversion to RPo

Transcription initiation is a multistep process involving a series of requisite conformational changes in RNA polymerase (R) and promoter DNA (P) that create the open complex (RP(o)). Here, we use the small solutes urea and glycine betaine (GB) to probe the extent and type of surface area changes in the formation of RP(o) between Esigma(70) RNA polymerase and lambdaP(R) promoter DNA. Effects of urea quantitatively reflect changes in amide surface and are particularly well-suited to detect coupled protein folding events. GB provides a qualitative probe for the exposure or burial of anionic surface. Kinetics of formation and dissociation of RP(o) reveal strikingly large effects of the solutes on the final steps of RP(o) formation: urea dramatically increases the dissociation rate constant k(d), whereas GB decreases the rate of dissociation. Formation of the first kinetically significant intermediate I(1) is disfavored in urea, and moderately favored by GB. GB slows the rate-determining step that converts I(1) to the second kinetically significant intermediate I(2); urea has no effect on this step. The most direct interpretation of these data is that recognition of promoter DNA in I(1) involves only limited conformational changes. Notably, the data support the following hypotheses: (1) the negatively charged N-terminal domain of sigma(70) remains bound in the "jaws" of polymerase in I(1); (2) the subsequent rate-determining isomerization step involves ejecting this domain from the jaws, thereby unmasking the active site; and (3) final conversion to RP(o) involves coupled folding of the mobile downstream clamp of polymerase.

Common and variable contributions of Fis residues to high-affinity binding at different DNA sequences

Fis is a nucleoid-associated protein that interacts with poorly related DNA sequences with a high degree of specificity. A difference of more than 3 orders of magnitude in apparent Kd values was observed between specific (Kd, approximately 1 to 4 nM) and nonspecific (Kd, approximately 4 microM) DNA binding. To examine the contributions of Fis residues to the high-affinity binding at different DNA sequences, 13 alanine substitutions were generated in or near the Fis helix-turn-helix DNA binding motif, and the resulting proteins were purified. In vitro binding assays at three different Fis sites (fis P II, hin distal, and lambda attR) revealed that R85, T87, R89, K90, and K91 played major roles in high-affinity DNA binding and that R85, T87, and K90 were consistently vital for binding to all three sites. Other residues made variable contributions to binding, depending on the binding site. N84 was required only for binding to the lambda attR Fis site, and the role of R89 was dramatically altered by the lambda attR DNA flanking sequence. The effects of Fis mutations on fis P II or hin distal site binding in vitro generally correlated with their abilities to mediate fis P repression or DNA inversion in vivo, demonstrating that the in vitro DNA-binding effects are relevant in vivo. The results suggest that while Fis is able to recognize a minimal common set of DNA sequence determinants at different binding sites, it is also equipped with a number of residues that contribute to the binding strength, some of which play variable roles.

The electrostatic character of the ribosomal surface enables extraordinarily rapid target location by ribotoxins

Alpha-sarcin ribotoxins comprise a unique family of ribonucleases that cripple the ribosome by catalyzing endoribonucleolytic cleavage of ribosomal RNA at a specific location in the sarcin/ricin loop (SRL). The SRL structure alone is cleaved site-specifically by the ribotoxin, but the ribosomal context enhances the reaction rate by several orders of magnitude. We show that, for the alpha-sarcin-like ribotoxin restrictocin, this catalytic advantage arises from favorable electrostatic interactions with the ribosome. Restrictocin binds at many sites on the ribosomal surface and under certain conditions cleaves the SRL with a second-order rate constant of 1.7 x 10(10) M(-1) s(-1), a value that matches the predicted frequency of random restrictocin-ribosome encounters. The results suggest a mechanism of target location whereby restrictocin encounters ribosomes randomly and diffuses within the ribosomal electrostatic field to the SRL. These studies show a role for electrostatics in protein-ribosome recognition.

Simulation of non-specific protein-mRNA interactions

Protein–nucleic acid interactions exhibit varying degrees of specificity. Relatively high affinity, sequence-specific interactions, can be studied with structure determination, but lower affinity, non-specific interactions are also of biological importance. We report simulations that predict the population of nucleic acid paths around protein surfaces, and give binding constant differences for changes in the protein scaffold. The method is applied to the non-specific component of interactions between eIF4Es and messenger RNAs that are bound tightly at the cap site. Adding a fragment of eIF4G to the system changes both the population of mRNA paths and the protein–mRNA binding affinity, suggesting a potential role for non-specific interactions in modulating translational properties. Generally, the free energy simulation technique could work in harness with characterized tethering points to extend analysis of nucleic acid conformation, and its modulation by protein scaffolds.

Surface salt bridges modulate the DNA site size of bacterial histone-like HU proteins

Bacterial histone-like DNA-binding proteins are best known for their role in compacting the genomic DNA. Of these proteins, HU is ubiquitous and highly conserved across the eubacterial kingdom. Using the HBsu (Bacillus subtilis-encoded HU homologue) as a model, we explore here the molecular basis for the ability of some HU homologues to engage a longer approx. 35 bp DNA site as opposed to the much shorter sites reported for other homologues. Using electrophoretic mobility-shift assays, we show that the DNA site size for HBsu is approx. 10-13 bp and that a specific surface salt bridge limits the DNA site size for HBsu. Surface exposure of the highly conserved Lys3, achieved by substitution of its salt-bridging partner Asp26 with Ala, leads to enhanced DNA compaction by HBsu-D26A (where D26A stands for the mutant Asp26-->Ala), consistent with the interaction of Lys3 with the ends of a 25 bp duplex. Both HBsu and HBsu-D26A bend DNA, as demonstrated by their equivalent ability to promote ligase-mediated DNA cyclization, indicating that residues involved in mediating DNA kinks are unaltered in the mutant protein. We suggest that Lys3 is important for DNA wrapping due to its position at a distance from the DNA kinks where it can exert optimal leverage on flanking DNA and that participation of Lys3 in a surface salt bridge competes for its interaction with DNA phosphates, thereby reducing the occluded site size.

Single-molecule spectroscopic determination of lac repressor-DNA loop conformation

The Escherichia coli lactose repressor protein (LacI) provides a classic model for understanding protein-induced DNA looping. LacI has a C-terminal four-helix bundle tetramerization domain that may act as a flexible hinge. In previous work, several DNA constructs, each containing two lac operators bracketing a sequence-induced bend, were designed to stabilize different possible looping geometries. The resulting hyperstable LacI-DNA loops exist as both a compact "closed" form with a V-shaped repressor and also a more "open" form with an extended hinge. The "9C14" construct was of particular interest because footprinting, electrophoretic mobility shift, and ring closure experiments suggested that it forms both geometries. Previous fluorescence resonance energy transfer (FRET) measurements gave an efficiency of energy transfer (ET) of 70%, confirming the existence of a closed form. These measurements could not determine whether open form or intermediate geometries are populated or the timescale of interconversion. We have now applied single-molecule FRET to Cy3, Cy5 double-labeled LacI-DNA loops diffusing freely in solution. By using multiple excitation wavelengths and by carefully examining the behavior of the zero-ET peak during titration with LacI, we show that the LacI-9C14 loop exists exclusively in a single closed form exhibiting essentially 100% ET.

Mg2+-modulated KMnO4 reactivity of thymines in the open transcription complex reflects variation in the negative electrostatic potential along the separated DNA strands. Footprinting of Escherichia coli RNA polymerase complex at the lambdaP(R) promoter revisited

There is still a controversy over the mechanism of promoter DNA strand separation upon open transcription complex (RPo) formation by Escherichia coli RNA polymerase: is it a single or a stepwise process controlled by Mg2+ ions and temperature? To resolve this question, the kinetics of pseudo-first-order oxidation of thymine residues by KMnO4 in the -11 ... +2 DNA region of RPo at the lambdaP(R) promoter was examined under single-hit conditions as a function of temperature (13-37 degrees C) in the absence or presence of 10 mm MgCl2. The reaction was also studied with respect to thymidine and its nucleotides (TMP, TTP and TpT) as a function of temperature and [MgCl2]. The kinetic parameters, (ox)k and (ox)E(a), and Mg-induced enhancement of (ox)k proved to be of the same order of magnitude for RPo-lambdaP(R) and the nucleotides. Unlike the complex, (ox)E(a) for the nucleotides was found to be Mg-independent. The isothermal increase in (ox)k with increasing [Mg2+] was thus interpreted in terms of a simple model of screening of the negative charges on phosphate groups by Mg2+ ions, lowering the electrostatic barrier to the diffusion of MnO4- anions to the reactive double bond of thymine. Similar screening isotherms were determined for the oxidation of two groups of thymines in RPo at a consensus-like Pa promoter, differing in the magnitude of the Mg effect. Together, the findings show that: (a) the two DNA strands in the -11...+2 region of RPo-lambdaP(R) are completely separated over the whole range of temperatures investigated (13-37 degrees C) in the absence of Mg2+ (b) Mg2+ ions induce an increase in the rate of the oxidation reaction by screening negatively charged phosphate and carboxylate groups; and (c) the observed thymine reactivity and the magnitude of the Mg effect reflect variation in the strength of the electrostatic potential along the separated DNA strands, in agreement with the current structural model of RPo.

Elastic correlations in nucleosomal DNA structure

The structure of DNA in the nucleosome core particle is studied using an elastic model that incorporates anisotropy in the bending energetics and twist-bend coupling. Using the experimentally determined structure of nucleosomal DNA [T. J. Richmond and C. A. Davey, Nature (London) 423, 145 (2003)], it is shown that elastic correlations exist between twist, roll, tilt, and stretching of DNA, as well as the distance between phosphate groups. The twist-bend coupling term is shown to be able to capture these correlations to a large extent, and a fit to the experimental data yields a new estimate of G = 25 nm for the value of the twist-bend coupling constant.

Induced fit and the entropy of structural adaptation in the complexation of CAP and lambda-repressor with cognate DNA sequences

Molecular dynamics (MD) simulations of 5 ns on protein-DNA complexes of catabolite-activator protein (CAP), lambda-repressor, and their corresponding uncomplexed protein and DNA, are reported. These cases represent two extremes of DNA bending, with CAP DNA bent severely and the lambda-operator nearly straight when complexed with protein. The calculations were performed using the AMBER suite of programs and the parm94 force field, validated for these studies by good agreement with experimental nuclear magnetic resonance data on DNA. An explicit computational model of structural adaptation and computation of the quasiharmonic entropy of association were obtained from the MD. The results indicate that, with respect to canonical B-form DNA, the extreme bending of the DNA in the complex with CAP is approximately 60% protein-induced and 40% intrinsic to the sequence-dependent structure of the free oligomer. The DNA in the complex is an energetically strained form, and the MD results are consistent with a conformational-capture mechanism. The calculated quasiharmonic entropy change accounts for the entropy difference between the two cases. The calculated entropy was decomposed into contributions from protein adaptation, DNA adaptation, and protein-DNA structural correlations. The origin of the entropy difference between CAP and lambda-repressor complexation arises more from the additional protein adaptation in the case of lambda, than to DNA bending and entropy contribution from DNA bending. The entropy arising from protein DNA cross-correlations, a contribution not previously discussed, is surprisingly large.

IHF and HU: flexible architects of bent DNA

The energetic cost of bending short segments of DNA is very high. This bending is critical for the packaging of DNA and is exploited to regulate many cellular processes. In prokaryotes, IHF and HU are key architectural proteins present at high concentrations. New protein-DNA co-crystal structures, and the adaptation of advanced biophysical and biochemical techniques have led to an improved understanding of how these proteins interact with DNA. These techniques include time-resolved synchrotron X-ray footprinting, differential scanning calorimetry, isothermal titration calorimetry and single-molecule experiments.

Downstream DNA selectively affects a paused conformation of human RNA polymerase II

Transcriptional pausing by human RNA polymerase II (RNAPII) in the HIV-1 LTR is caused principally by a weak RNA:DNA hybrid that allows rearrangement of reactive or catalytic groups in the enzyme's active site. This rearrangement creates a transiently paused state called the unactivated intermediate that can backtrack into a more long-lived paused species. We report that three different regions of the not-yet-transcribed DNA just downstream of the pause site affect the duration of the HIV-1 pause, and also can influence pause formation. Downstream DNA in at least one region, a T-tract from +5 to +8, increases pause duration by specifically affecting the unactivated intermediate, without corresponding effects on the active or backtracked states. We suggest this effect depends on RNAPII-modulated DNA plasticity and speculate it is mediated by the "trigger loop" thought to participate in RNAP's catalytic cycle. These findings provide a new framework for understanding downstream DNA effects on RNAP.

Unwrapping of DNA-protein complexes under external stretching

A DNA-protein complex modeled by a semiflexible chain and an attractive spherical core is studied in the situation when an external stretching force is acting on one end monomer of the chain while the other end monomer is kept fixed in space. Without a stretching force, the chain is wrapped around the core. By applying an external stretching force, unwrapping of the complex is induced. We study the statics and dynamics of the unwrapping process by computer simulations and simple phenomenological theory. We find two different scenarios depending on the chain stiffness: For a flexible chain, the extension of the complex scales linearly with the external force applied. The sphere-chain complex is disordered; i.e., there is no clear winding of the chain around the sphere. For a stiff chain, on the other hand, the complex structure is ordered, which is reminiscent of nucleosome. There is a clear winding number, and the unwrapping process under external stretching is discontinuous with jumps of the distance-force curve. This is associated with discrete unwinding processes of the complex. Our predictions are of relevance for experiments, which measure force-extension curves of DNA-protein complexes, such as nucleosome, using optical tweezers.

Adsorption of DNA to mica mediated by divalent counterions: a theoretical and experimental study

The adsorption of DNA molecules onto a flat mica surface is a necessary step to perform atomic force microscopy studies of DNA conformation and observe DNA-protein interactions in physiological environment. However, the phenomenon that pulls DNA molecules onto the surface is still not understood. This is a crucial issue because the DNA/surface interactions could affect the DNA biological functions. In this paper we develop a model that can explain the mechanism of the DNA adsorption onto mica. This model suggests that DNA attraction is due to the sharing of the DNA and mica counterions. The correlations between divalent counterions on both the negatively charged DNA and the mica surface can generate a net attraction force whereas the correlations between monovalent counterions are ineffective in the DNA attraction. DNA binding is then dependent on the fractional surface densities of the divalent and monovalent cations, which can compete for the mica surface and DNA neutralizations. In addition, the attraction can be enhanced when the mica has been pretreated by transition metal cations (Ni(2+), Zn(2+)). Mica pretreatment simultaneously enhances the DNA attraction and reduces the repulsive contribution due to the electrical double-layer force. We also perform end-to-end distance measurement of DNA chains to study the binding strength. The DNA binding strength appears to be constant for a fixed fractional surface density of the divalent cations at low ionic strength (I < 0.1 M) as predicted by the model. However, at higher ionic strength, the binding is weakened by the screening effect of the ions. Then, some equations were derived to describe the binding of a polyelectrolyte onto a charged surface. The electrostatic attraction due to the sharing of counterions is particularly effective if the polyelectrolyte and the surface have nearly the same surface charge density. This characteristic of the attraction force can explain the success of mica for performing single DNA molecule observation by AFM. In addition, we explain how a reversible binding of the DNA molecules can be obtained with a pretreated mica surface.

DNA spools under tension

DNA spools, structures in which DNA is wrapped and helically coiled onto itself or onto a protein core, are ubiquitous in nature. We develop a general theory describing the nonequilibrium behavior of DNA spools under linear tension. Two puzzling and seemingly unrelated recent experimental findings, the sudden quantized unwrapping of nucleosomes and that of DNA toroidal condensates under tension, are theoretically explained and shown to be of the same origin. The study provides new insights into nucleosome and chromatin fiber stability and dynamics.

Osmo-regulation of bacterial transcription via poised RNA polymerase

Adaptation to high-salt environments is critical for the survival of a wide range of cells, especially for pathogenic bacteria that colonize the animal gut and urinary tract. The adaptation strategy involves production of the salt potassium glutamate, which induces a specific gene expression program that produces electro-neutral osmolytes while inhibiting general sigma(70) transcription. These data show that in Escherichia coli potassium glutamate stimulates transcription by disengaging inhibitory polymerase interactions at a sigma(38) promoter. These occur in an upstream region that is marked by an osmotic shock promoter DNA consensus sequence. The disruption activates a poised RNA polymerase to transcribe. This transcription program leads to the production of osmolytes that are shown to have only minor effects on transcription and therefore help to restore normal cell function. An osmotic shock gene expression cycle is discussed.