Computational analysis of looping of a large family of highly bent DNA by LacI

Sequence-dependent intrinsic curvature of DNA influences looping by regulatory proteins such as LacI and NtrC. Curvature can enhance stability and control shape, as observed in LacI loops formed with three designed sequences with operators bracketing an A-tract bend. We explore geometric, topological, and energetic effects of curvature with an analysis of a family of highly bent sequences, using the elastic rod model from previous work. A unifying straight-helical-straight representation uses two phasing parameters to describe sequences composed of two straight segments that flank a common helically supercoiled segment. We exercise the rod model over this two-dimensional space of phasing parameters to evaluate looping behaviors. This design space is found to comprise two subspaces that prefer parallel versus anti-parallel binding topologies. The energetic cost of looping varies from 4 to 12 kT. Molecules can be designed to yield distinct binding topologies as well as hyperstable or hypostable loops and potentially loops that can switch conformations. Loop switching could be a mechanism for control of gene expression. Model predictions for linking numbers and sizes of LacI-DNA loops can be tested using multiple experimental approaches, which coupled with theory could address whether proteins or DNA provide the observed flexibility of protein-DNA loops.     

Statistical-mechanical theory of DNA looping

The lack of a rigorous analytical theory for DNA looping has caused many DNA-loop-mediated phenomena to be interpreted using theories describing the related process of DNA cyclization. However, distinctions in the mechanics of DNA looping versus cyclization can have profound quantitative effects on the thermodynamics of loop closure. We have extended a statistical mechanical theory recently developed for DNA cyclization to model DNA looping, taking into account protein flexibility. Notwithstanding the underlying theoretical similarity, we find that the topological constraint of loop closure leads to the coexistence of multiple classes of loops mediated by the same protein structure. These loop topologies are characterized by dramatic differences in twist and writhe; because of the strong coupling of twist and writhe within a loop, DNA looping can exhibit a complex overall helical dependence in terms of amplitude, phase, and deviations from uniform helical periodicity. Moreover, the DNA-length dependence of optimal looping efficiency depends on protein elasticity, protein geometry, and the presence of intrinsic DNA bends. We derive a rigorous theory of loop formation that connects global mechanical and geometric properties of both DNA and protein and demonstrates the importance of protein flexibility in loop-mediated protein-DNA interactions.     

An assessment of three dinucleotide parameters to predict DNA curvature by quantitative comparison with experimental data

Curved DNA fragments are often found near functionally important sites such as promoters and origins of replication, and hence sequence-dependent DNA curvature prediction is of great utility in genomics and bioinformatics. In light of this, an assessment of three different dinucleotide step parameters (based on gel retardation as well as crystal structure data) is carried out. These parameters (BMHT, LB and CS) are evaluated quantitatively for their ability to predict correctly the experimental results of a large set of nucleic acid sequences containing A-tracts as well as GC-rich motifs. This set contained around 40 synthetic as well as natural sequences whose solution properties have been well characterized experimentally. All three models could account reasonably well for curvature in the various DNA sequences. The CS model, where dinucleotide parameters are calculated from crystal structure data, consistently shows slightly better correlation with experimental data. Our simple analysis also indicates that presently available trinucleotide parameters fail to predict curvature in some of the well-characterized sequences. The study shows that the dinucleotide parameters with some further refinement can be used to predict sequence-dependent curvature correctly in genomic sequences.     

Dual role of DNA intrinsic curvature and flexibility in determining nucleosome stability

A statistical mechanistic approach to evaluate the sequence-dependent thermodynamic stability of nucleosomes is proposed. The model is based on the calculation of the DNA intrinsic curvature, obtained by integrating the nucleotide step deviations from the canonical B-DNA structure, and on the evaluation of the first order elastic distortion energy to reach the nucleosomal superstructure. Literature data on the free energy of nucleosome formation as obtained by competitive nucleosome reconstitution of a significant pool of different DNA sequences were compared with the theoretical results, and a satisfactorily good correlation was found. A striking result of the comparison is the emergence of two opposite roles of the DNA intrinsic curvature and flexibility in determining nucleosome stability. Finally, the obtained results suggest that the curvature-dependent DNA hydration should play a relevant role in the sequence-dependent nucleosome stability.     

The loopometer: a quantitative in vivo assay for DNA-looping proteins

Proteins that can bring together separate DNA sites, either on the same or on different DNA molecules, are critical for a variety of DNA-based processes. However, there are no general and technically simple assays to detect proteins capable of DNA looping in vivo nor to quantitate their in vivo looping efficiency. Here, we develop a quantitative in vivo assay for DNA-looping proteins in Escherichia coli that requires only basic DNA cloning techniques and a LacZ assay. The assay is based on loop assistance, where two binding sites for the candidate looping protein are inserted internally to a pair of operators for the E. coli LacI repressor. DNA looping between the sites shortens the effective distance between the lac operators, increasing LacI looping and strengthening its repression of a lacZ reporter gene. Analysis based on a general model for loop assistance enables quantitation of the strength of looping conferred by the protein and its binding sites. We use this ‘loopometer’ assay to measure DNA looping for a variety of bacterial and phage proteins.     

Experimental evaluation of the Liu–Beveridge dinucleotide step model of DNA structure

Methods for predicting DNA curvature have many possible applications. Dinucleotide step models describe DNA shape by characterization of helical twist, deflection angles and the direction of deflection for nearest neighbor base pairs. Liu and Beveridge have extended previous applications of dinucleotide step models with the development and qualitative validation of a predictive method for sequence-dependent DNA curvature (the LB model). We tested whether the LB model accurately predicts experimentally deduced curvature angles and helical repeat parameters for DNA sequences not in its training set, particularly when challenged with quantitative data and subtle sequence phasings. We examined a series of 17 well-characterized DNA sequences to compare electrophoretic and computational results. The LB model is superior to two other models in the prediction of helical repeat parameters. We observed a strong linear correlation between curvature magnitudes predicted using the LB model and those determined by electrophoretic ligation ladder experiments, although the LB model somewhat underestimated apparent curvature. With longer electrophoretic phasing probes the LB model slightly overestimated gel mobility anomalies, with modest deviations in predicted helical repeat parameters. Overall, our analyses suggest that the LB model provides reasonably accurate predictions for the electrophoretic behavior of DNA.     

Gal repressosome contains an antiparallel DNA loop

Gal repressosome assembly and repression of the gal operon in Escherichia coli occurs when two dimeric GalR proteins and the histone-like HU protein bind to cognate sites causing DNA looping. Structure-based genetic analysis defined the GalR surfaces interacting to form a stacked, V-shaped, tetrameric structure. Stereochemical models of the four possible DNA loops compatible with the GalR tetramer configuration were constructed using the sequence-dependent structural parameters of the interoperator DNA and conformation changes caused by GalR and asymmetric HU binding. Evaluation of their DNA elastic energies gave unambiguous preference to a loop structure in which the two gal operators adopt an antiparallel orientation causing undertwisting of DNA.     

A Possible Role of DNA Superstructures in Genome Evolution

The concept of DNA as a simple repository of the gene information has changed in that of a polymorphic macromolecule, which plays a relevant part in the management of the complex biochemical transformations in living matter. As a consequence of the slight stereochemical differences between base pairs, the direction of the DNA double helix axis undergoes deterministic writhing. A useful representation of such sequence-dependent structural distortions is the curvature diagram. Here, it is reported as an evolution simulation obtained by extensive point mutations along a biologically important DNA tract. The curvature changes, consequence of the point mutations. were compared to the related experimental gel electrophoresis mobility. The curvature of most mutants decreases and the mobility increases accordingly, suggesting the curvature of that tract is genetically selected. Moreover, DNA images by scanning force microscopy, show evidence of a sequence-dependent adhesion of curved DNA tracts to inorganic crystal surfaces. In particular, mica shows a large affinity towards the TT-rich dinucleotide sequences. This suggests a possible mechanism of selection of curved DNA regions, characterized by AA.TT dinucleotides in phase with double-helical periodicity, in the very early evolution steps.     

Terminal twist induced continuous writhe of a circular rod with intrinsic curvature

It is well known that a large linking number induces an abrupt writhing of a circular rod with zero intrinsic curvature, i.e., the stress-free state of the rod is straight. We show here that for any rod with a uniform natural curvature, no matter how small the intrinsic curvature is, a twist will induce a continuous writhing from the circular configuration and the abrupt writhing is only the limiting case when the intrinsic curvature is absolutely zero. The implication of this result on elastic models of circular DNA is discussed.     

基于蛋白质-DNA复合物晶体结构的DNA结构动力学特性研究

依据最新NDB数据库中蛋白质-DNA复合物晶体结构数据,考虑DNA结构的序列依赖特性,对10个二联体和36个约化的四联体计算了DNA动力学结构模型中的关键参数——力常数矩阵,矩阵中的非对角项反映了结构参数间的关联。利用改进的DNA结构动力学模型,可以方便地计算给定序列的结构动力学特性。     

Sequence-dependent DNA curvature and flexibility from scanning force microscopy images

This paper reports a study of the sequence-dependent DNA curvature and flexibility based on scanning force microscopy (SFM) images. We used a palindromic dimer of a 1878-bp pBR322 fragment and collected a large pool of SFM images. The curvature of each imaged chain was measured in modulus and direction. It was found that the ensemble curvature modulus does not allow the separation of static and dynamic contributions to the curvature, whereas the curvature, when its direction in the two dimensions is taken into account, permits the direct separation of the intrinsic curvature contributions static and dynamic contributions. The palindromic symmetry also acted as an internal gauge of the validity of the SFM images statistical analysis. DNA static curvature resulted in good agreement with the predicted sequence-dependent intrinsic curvature. Furthermore, DNA sequence-dependent flexibility was found to correlate with the occurrence of A.T-rich dinucleotide steps along the chain and, in general, with the normalized basepair stacking energy distribution.     

The Effect of Nonspecific Binding of Lambda Repressor on DNA Looping Dynamics

The λ repressor (CI) protein-induced DNA loop maintains stable lysogeny, yet allows efficient switching to lysis. Herein, the kinetics of loop formation and breakdown has been characterized at various concentrations of protein using tethered particle microscopy and a novel, to our knowledge, method of analysis. Our results show that a broad distribution of rate constants and complex kinetics underlie loop formation and breakdown. In addition, comparison of the kinetics of looping in wild-type DNA and DNA with mutated o3 operators showed that these sites may trigger nucleation of nonspecific binding at the closure of the loop. The average activation energy calculated from the rate constant distribution is consistent with a model in which nonspecific binding of CI between the operators shortens their effective separation, thereby lowering the energy barrier for loop formation and broadening the rate constant distribution for looping. Similarly, nonspecific binding affects the kinetics of loop breakdown by increasing the number of loop-securing protein interactions, and broadens the rate constant distribution for this reaction. Therefore, simultaneous increase of the rate constant for loop formation and reduction of that for loop breakdown stabilizes lysogeny. Given these simultaneous changes, the frequency of transitions between the looped and the unlooped state remains nearly constant. Although the loop becomes more stable thermodynamically with increasing CI concentration, it still opens periodically, conferring sensitivity to environmental changes, which may require switching to lytic conditions.     

Studying DNA Looping by Single-Molecule FRET

Bending of double-stranded DNA (dsDNA) is associated with many important biological processes such as DNA-protein recognition and DNA packaging into nucleosomes. Thermodynamics of dsDNA bending has been studied by a method called cyclization which relies on DNA ligase to covalently join short sticky ends of a dsDNA. However, ligation efficiency can be affected by many factors that are not related to dsDNA looping such as the DNA structure surrounding the joined sticky ends, and ligase can also affect the apparent looping rate through mechanisms such as nonspecific binding. Here, we show how to measure dsDNA looping kinetics without ligase by detecting transient DNA loop formation by FRET (Fluorescence Resonance Energy Transfer). dsDNA molecules are constructed using a simple PCR-based protocol with a FRET pair and a biotin linker. The looping probability density known as the J factor is extracted from the looping rate and the annealing rate between two disconnected sticky ends. By testing two dsDNAs with different intrinsic curvatures, we show that the J factor is sensitive to the intrinsic shape of the dsDNA.     

Effects of nucleoid proteins on DNA repression loop formation in Escherichia coli

The intrinsic stiffness of DNA limits its ability to be bent and twisted over short lengths, but such deformations are required for gene regulation. One classic paradigm is DNA looping in the regulation of the Escherichia coli lac operon. Lac repressor protein binds simultaneously to two operator sequences flanking the lac promoter. Analysis of the length dependence of looping-dependent repression of the lac operon provides insight into DNA deformation energetics within cells. The apparent flexibility of DNA is greater in vivo than in vitro, possibly because of host proteins that bind DNA and induce sites of flexure. Here we test DNA looping in bacterial strains lacking the nucleoid proteins HU, IHF or H-NS. We confirm that deletion of HU inhibits looping and that quantitative modeling suggests residual looping in the induced operon. Deletion of IHF has little effect. Remarkably, DNA looping is strongly enhanced in the absence of H-NS, and an explanatory model is proposed. Chloroquine titration, psoralen crosslinking and supercoiling-sensitive reporter assays show that the effects of nucleoid proteins on looping are not correlated with their effects on either total or unrestrained supercoiling. These results suggest that host nucleoid proteins can directly facilitate or inhibit DNA looping in bacteria.     

DNA modeling reveals an extended lac repressor conformation in classic in vitro binding assays

Protein-mediated DNA looping, such as that induced by the lactose repressor (LacI) of Escherichia coli, is a well-known gene regulation mechanism. Although researchers have given considerable attention to DNA looping by LacI, many unanswered questions about this mechanism, including the role of protein flexibility, remain. Recent single-molecule observations suggest that the two DNA-binding domains of LacI are capable of splaying open about the tetramerization domain into an extended conformation. We hypothesized that if recent experiments were able to reveal the extended conformation, it is possible that such structures occurred in previous studies as well. In this study, we tested our hypothesis by reevaluating two classic in vitro binding assays using a computational rod model of DNA. The experiments and computations evaluate the looping of both linear DNA and supercoiled DNA minicircles over a broad range of DNA interoperator lengths. The computed energetic minima align well with the experimentally observed interoperator length for optimal loop stability. Of equal importance, the model reveals that the most stable loops for linear DNA occur when LacI adopts the extended conformation. In contrast, for DNA minicircles, optimal stability may arise from either the closed or the extended protein conformation depending on the degree of supercoiling and the interoperator length.     

Fluorescence Resonance Energy Transfer over ∼130 Basepairs in Hyperstable Lac Repressor-DNA Loops

Lac repressor (LacI) binds two operator DNA sites, looping the intervening DNA. DNA molecules containing two lac operators bracketing a sequence-directed bend were previously shown to form hyperstable LacI-looped complexes. Biochemical studies suggested that orienting the operators outward relative to the bend direction (in construct 9C14) stabilizes a positively supercoiled closed form, with a V-shaped LacI, but that the most stable loop construct (11C12) is a more open form. Here, fluorescence resonance energy transfer (FRET) is measured on DNA loops, between fluorescein and TAMRA attached near the two operators, ~130 basepairs apart. For 9C14, efficient LacI-induced energy transfer (~74% based on donor quenching) confirms that the designed DNA shape can force the looped complex into a closed form. From enhanced acceptor emission, correcting for observed donor-dependent quenching of acceptor fluorescence, ~52% transfer was observed. Time-resolved FRET suggests that this complex exists in both closed- and open form populations. Less efficient transfer, ~10%, was detected for DNA-LacI sandwiches and 11C12-LacI, consistent with an open form loop. This demonstration of long-range FRET in large DNA loops confirms that appropriate DNA design can control loop geometry. LacI flexibility may allow it to maintain looping with other proteins bound or under different intracellular conditions.     

Sequence periodicity in nucleosomal DNA and intrinsic curvature

Background

Most eukaryotic DNA contained in the nucleus is packaged by wrapping DNA around histone octamers. Histones are ubiquitous and bind most regions of chromosomal DNA. In order to achieve smooth wrapping of the DNA around the histone octamer, the DNA duplex should be able to deform and should possess intrinsic curvature. The deformability of DNA is a result of the non-parallelness of base pair stacks. The stacking interaction between base pairs is sequence dependent. The higher the stacking energy the more rigid the DNA helix, thus it is natural to expect that sequences that are involved in wrapping around the histone octamer should be unstacked and possess intrinsic curvature. Intrinsic curvature has been shown to be dictated by the periodic recurrence of certain dinucleotides. Several genome-wide studies directed towards mapping of nucleosome positions have revealed periodicity associated with certain stretches of sequences. In the current study, these sequences have been analyzed with a view to understand their sequence-dependent structures.

Results

Higher order DNA structures and the distribution of molecular bend loci associated with 146 base nucleosome core DNA sequence from C. elegans and chicken have been analyzed using the theoretical model for DNA curvature. The curvature dispersion calculated by cyclically permuting the sequences revealed that the molecular bend loci were delocalized throughout the nucleosome core region and had varying degrees of intrinsic curvature.

Conclusions

The higher order structures associated with nucleosomes of C.elegans and chicken calculated from the sequences revealed heterogeneity with respect to the deviation of the DNA axis. The results points to the possibility of context dependent curvature of varying degrees to be associated with nucleosomal DNA.