Luke Czapla

Effects of the nucleoid protein HU on the structure, flexibility, and ring-closure properties of DNA deduced from Monte Carlo simulations.

The histone-like HU (heat unstable) protein plays a key role in the organization and regulation of the Escherichia coli genome. The nonspecific nature of HU binding to DNA complicates analysis of the mechanism by which the protein contributes to the looping of DNA. Conventional models of the looping of HU-bound duplexes attribute the changes in biophysical properties of DNA brought about by the random binding of protein to changes in the effective parameters of an ideal helical wormlike chain. Here, we introduce a novel Monte Carlo approach to study the effects of nonspecific HU binding on the configurational properties of DNA directly. We randomly decorated segments of an ideal double-helical DNA with HU molecules that induce the bends and other structural distortions of the double helix find in currently available X-ray structures. We find that the presence of HU at levels approximating those found in the cell reduces the persistence length by roughly threefold compared with that of naked DNA. The binding of protein has particularly striking effects on the cyclization properties of short duplexes, altering the dependence of ring closure on chain length in a way that cannot be mimicked by a simple wormlike model and accumulating at higher-than-expected levels on successfully closed chains. Moreover, the uptake of protein on small minicircles depends on chain length, taking advantage of the HU-induced deformations of DNA structure to facilitate ligation. Circular duplexes with bound HU show much greater propensity than protein-free DNA to exist as negatively supercoiled topoisomers, suggesting a potential role of HU in organizing the bacterial nucleoid. The local bending and undertwisting of DNA by HU, in combination with the number of bound proteins, provide a structural rationale for the condensation of DNA and the observed expression levels of reporter genes in vivo.

Understanding apparent DNA flexibility enhancement by HU and HMGB architectural proteins.

Understanding and predicting the mechanical properties of protein/DNA complexes are challenging problems in biophysics. Certain architectural proteins bind DNA without sequence specificity and strongly distort the double helix. These proteins rapidly bind and unbind, seemingly enhancing the flexibility of DNA as measured by cyclization kinetics. The ability of architectural proteins to overcome DNA stiffness has important biological consequences, but the detailed mechanism of apparent DNA flexibility enhancement by these proteins has not been clear. Here, we apply a novel Monte Carlo approach that incorporates the precise effects of protein on DNA structure to interpret new experimental data for the bacterial histone-like HU protein and two eukaryotic high-mobility group class B (HMGB) proteins binding to ∼200-bp DNA molecules. These data (experimental measurement of protein-induced increase in DNA cyclization) are compared with simulated cyclization propensities to deduce the global structure and binding characteristics of the closed protein/DNA assemblies. The simulations account for all observed (chain length and concentration dependent) effects of protein on DNA behavior, including how the experimental cyclization maxima, observed at DNA lengths that are not an integral helical repeat, reflect the deformation of DNA by the architectural proteins and how random DNA binding by different proteins enhances DNA cyclization to different levels. This combination of experiment and simulation provides a powerful new approach to resolve a long-standing problem in the biophysics of protein/DNA interactions.

BROMOC-D: Brownian Dynamics/Monte-Carlo Program Suite to Study Ion and DNA Permeation in Nanopores

A theoretical framework is presented to model ion and DNA translocation across a nanopore confinement under an applied electric field. A combined Grand Canonical Monte Carlo Brownian Dynamics (GCMC/BD) algorithm offers a general approach to study ion permeation through wide molecular pores with a direct account of ion–ion and ion–DNA correlations. This work extends previously developed theory by incorporating the recently developed coarse-grain polymer model of DNA by de Pablo and colleagues [Knotts, T. A.; Rathore, N.; Schwartz, D. C.; de Pablo, J. J. J. Chem. Phys. 2007, 126] with explicit ions for simulations of polymer dynamics. Atomistic MD simulations were used to guide model developments. The power of the developed scheme is illustrated with studies of single-stranded DNA (ss-DNA) oligomer translocation in two model cases: a cylindrical pore with a varying radius and a well-studied experimental system, the staphylococcal α-hemolysin channel. The developed model shows good agreement with experimental data for model studies of two homopolymers: ss-poly(dA)n and ss-poly(dC)n. The developed protocol allows for direct evaluation of different factors (charge distribution and pore shape and size) controlling DNA translocation in a variety of nanopores.

DNA topology confers sequence specificity to nonspecific architectural proteins

Significance The organization of long genomes in the confined spaces of a cell requires special facilitating mechanisms. A variety of architectural proteins play key roles in these processes. The bacterial heat-unstable (HU) protein helps to condense DNA by introducing sharp turns along its pathway. The protein binds in a sequence-neutral fashion, randomly distorting linear DNA when introduced in computer-simulated structures at levels comparable to those found in the cell. The natural resistance of DNA to severe deformation, however, restricts the nonspecific protein to specific loci when the molecule is covalently closed or looped by a protein. The interplay of DNA topology and protein-induced bending provides insights into ways in which gene fragments may be organized and linked to biological function. Topological constraints placed on short fragments of DNA change the disorder found in chain molecules randomly decorated by nonspecific, architectural proteins into tightly organized 3D structures. The bacterial heat-unstable (HU) protein builds up, counter to expectations, in greater quantities and at particular sites along simulated DNA minicircles and loops. Moreover, the placement of HU along loops with the “wild-type” spacing found in the Escherichia coli lactose (lac) and galactose (gal) operons precludes access to key recognition elements on DNA. The HU protein introduces a unique spatial pathway in the DNA upon closure. The many ways in which the protein induces nearly the same closed circular configuration point to the statistical advantage of its nonspecificity. The rotational settings imposed on DNA by the repressor proteins, by contrast, introduce sequential specificity in HU placement, with the nonspecific protein accumulating at particular loci on the constrained duplex. Thus, an architectural protein with no discernible DNA sequence-recognizing features becomes site-specific and potentially assumes a functional role upon loop formation. The locations of HU on the closed DNA reflect long-range mechanical correlations. The protein responds to DNA shape and deformability—the stiff, naturally straight double-helical structure—rather than to the unique features of the constituent base pairs. The structures of the simulated loops suggest that HU architecture, like nucleosomal architecture, which modulates the ability of regulatory proteins to recognize their binding sites in the context of chromatin, may influence repressor–operator interactions in the context of the bacterial nucleoid.

Interplay of Protein and DNA Structure Revealed in Simulations of the lac Operon

The E. coli Lac repressor is the classic textbook example of a protein that attaches to widely spaced sites along a genome and forces the intervening DNA into a loop. The short loops implicated in the regulation of the lac operon suggest the involvement of factors other than DNA and repressor in gene control. The molecular simulations presented here examine two likely structural contributions to the in-vivo looping of bacterial DNA: the distortions of the double helix introduced upon association of the highly abundant, nonspecific nucleoid protein HU and the large-scale deformations of the repressor detected in low-resolution experiments. The computations take account of the three-dimensional arrangements of nucleotides and amino acids found in crystal structures of DNA with the two proteins, the natural rest state and deformational properties of protein-free DNA, and the constraints on looping imposed by the conformation of the repressor and the orientation of bound DNA. The predicted looping propensities capture the complex, chain-length-dependent variation in repression efficacy extracted from gene expression studies and in vitro experiments and reveal unexpected chain-length-dependent variations in the uptake of HU, the deformation of repressor, and the folding of DNA. Both the opening of repressor and the presence of HU, at levels approximating those found in vivo, enhance the probability of loop formation. HU affects the global organization of the repressor and the opening of repressor influences the levels of HU binding to DNA. The length of the loop determines whether the DNA adopts antiparallel or parallel orientations on the repressor, whether the repressor is opened or closed, and how many HU molecules bind to the loop. The collective behavior of proteins and DNA is greater than the sum of the parts and hints of ways in which multiple proteins may coordinate the packaging and processing of genetic information.

The Role of Architectural Proteins in Lac-Mediated DNA Looping

The widely abundant nucleoid protein HU contributes to both the spatial organization and biological processing of bacterial DNA. One such contribution of HU to biological processing is tied to gene regulation mediated by the looping of DNA induced by proteins, such as the tetrameric Lac Repressor protein, which simultaneously binds two distant operator DNA sites in the bacterial genome. Current understanding of how HU might contribute to the transcription of the lac genes derives from indirect theoretical and computational analysis of the effects of chain length on gene expression in E. coli cells with and without HU. We take a more direct approach to the question of how architectural proteins mediate the structure and looping of DNA in vivo, incorporating the structural effects of both HU and Lac Repressor on DNA in Monte Carlo simulations, taking advantage of new methods that we have developed to understand the properties of protein-bound DNA. We present our predictions from models for understanding Lac-mediated gene expression in vivo, and also demonstrate state-of-the-art calculations that incorporate our structure-based models in order to obtain more detailed insights into the phenomena observed in ring-closure experiments of DNA in the presence of architectural proteins.