Simulations of DNA topoisomerase 1B bound to supercoiled DNA reveal changes in the flexibility pattern of the enzyme and a secondary protein–DNA binding site

Human topoisomerase 1B has been simulated covalently bound to a negatively supercoiled DNA minicircle, and its behavior compared to the enzyme bound to a simple linear DNA duplex. The presence of the more realistic supercoiled substrate facilitates the formation of larger number of protein–DNA interactions when compared to a simple linear duplex fragment. The number of protein–DNA hydrogen bonds doubles in proximity to the active site, affecting all of the residues in the catalytic pentad. The clamp over the DNA, characterized by the salt bridge between Lys369 and Glu497, undergoes reduced fluctuations when bound to the supercoiled minicircle. The linker domain of the enzyme, which is implicated in the controlled relaxation of superhelical stress, also displays an increased number of contacts with the minicircle compared to linear DNA. Finally, the more complex topology of the supercoiled DNA minicircle gives rise to a secondary DNA binding site involving four residues located on subdomain III. The simulation trajectories reveal significant changes in the interactions between the enzyme and the DNA for the more complex DNA topology, which are consistent with the experimental observation that the protein has a preference for binding to supercoiled DNA.     

How do site-specific DNA-binding proteins find their targets?

Essentially all the biological functions of DNA depend on site-specific DNA-binding proteins finding their targets, and therefore ‘searching’ through megabases of non-target DNA. In this article, we review current understanding of how this sequence searching is done. We review how simple diffusion through solution may be unable to account for the rapid rates of association observed in experiments on some model systems, primarily the Lac repressor. We then present a simplified version of the ‘facilitated diffusion’ model of Berg, Winter and von Hippel, showing how non-specific DNA–protein interactions may account for accelerated targeting, by permitting the protein to sample many binding sites per DNA encounter. We discuss the 1-dimensional ‘sliding’ motion of protein along non-specific DNA, often proposed to be the mechanism of this multiple site sampling, and we discuss the role of short-range diffusive ‘hopping’ motions. We then derive the optimal range of sliding for a few physical situations, including simple models of chromosomes in vivo, showing that a sliding range of ~100 bp before dissociation optimizes targeting in vivo. Going beyond first-order binding kinetics, we discuss how processivity, the interaction of a protein with two or more targets on the same DNA, can reveal the extent of sliding and we review recent experiments studying processivity using the restriction enzyme EcoRV. Finally, we discuss how single molecule techniques might be used to study the dynamics of DNA site-specific targeting of proteins.     

How Lac repressor finds lac operator in vitro.

Filter-binding and gel mobility shift assays were used to analyse the kinetics of the interaction of Lac repressor with lac operator. A comparison of the two techniques reveals that filter-binding assays with tetrameric Lac repressor have often been misinterpreted. It has been assumed that all complexes of Lac repressor and lac operator DNA bind with equal affinity to nitrocellulose filters. This assumption is wrong. Sandwich or loop complexes where two lac operators bind to one tetrameric Lac repressor are not or are only badly retained on nitrocellulose filters under normal conditions. Taking this into account, dimeric and tetrameric Lac repressor do not show any DNA-length dependence of their association and dissociation rate constants when they bind to DNA fragments smaller than 2455 base-pairs carrying a single symmetric ideal lac operator. A ninefold increased association rate to ideal lac operator on lambda DNA is observed for tetrameric but not dimeric Lac repressor. It is presumably due to intersegment transfer involving lac operator-like sequences.     

From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.     

Inferring primase-DNA specific recognition using a data driven approach

DNA–protein interactions play essential roles in all living cells. Understanding of how features embedded in the DNA sequence affect specific interactions with proteins is both challenging and important, since it may contribute to finding the means to regulate metabolic pathways involving DNA–protein interactions. Using a massive experimental benchmark dataset of binding scores for DNA sequences and a machine learning workflow, we describe the binding to DNA of T7 primase, as a model system for specific DNA–protein interactions. Effective binding of T7 primase to its specific DNA recognition sequences triggers the formation of RNA primers that serve as Okazaki fragment start sites during DNA replication.     

The solution structure of Lac repressor headpiece 62 complexed to a symmetrical lac operator

BACKGROUND: Lactose repressor protein (Lac) controls the expression of the lactose metabolic genes in Escherichia coli by binding to an operator sequence in the promoter of the lac operon. Binding of inducer molecules to the Lac core domain induces changes in tertiary structure that are propagated to the DNA-binding domain through the connecting hinge region, thereby reducing the affinity for the operator. Protein-protein and protein-DNA interactions involving the hinge region play a crucial role in the allosteric changes occurring upon induction, but have not, as yet, been analyzed in atomic detail. RESULTS: We have used nuclear magnetic resonance (NMR) spectroscopy and restrained molecular dynamics (rMD) to determine the structure of the Lac repressor DNA-binding domain (headpeice 62; HP62) in complex with a symmetrized lac operator. Analysis of the structures reveals specific interactions between Lac repressor and DNA that were not found in previously investigated Lac repressor-DNA complexes. Important differences with the previously reported structures of the HP56-DNA complex were found in the loop following the helix-turn-helix (HTH) motif. The protein-protein and protein-DNA interactions involving the hinge region and the deformations in the DNA structure could be delineated in atomic detail. The structures were also used for comparison with the available crystallographic data on the Lac and Pur repressor-DNA complexes. CONCLUSIONS: The structures of the HP62-DNA complex provide the basis for a better understanding of the specific recognition in the Lac repressor-operator complex. In addition, the structural features of the hinge region provide detailed insight into the protein-protein and protein-DNA interactions responsible for the high affinity of the repressor for operator DNA.     

Structure-specific binding of the proto-oncogene protein DEK to DNA

The ubiquitous proto-oncogene protein DEK has been found to be associated with chromatin during the entire cell cycle. It changes the topology of DNA in chromatin and protein-free DNA through the introduction of positive supercoils. The sequence and structure specificities of DEK–DNA interactions are not completely understood. The binding of DEK to DNA is not sequence specific, but we describe here that DEK has a clear preference for supercoiled and four-way junction DNA. In the presence of topoisomerase II, DEK stimulates intermolecular catenation of circular DNA molecules. DEK also increases the probability of intermolecular ligation of linear DNA molecules by DNA ligase. These binding properties qualify DEK as an architectural protein.     

Aspects of protein–DNA interactions: a review of quantitative thermodynamic theory for modelling synthetic circuits utilising LacI and CI repressors,IPTG and the reporter gene <Emphasis Type="Italic">lacZ</Emphasis>

Protein–DNA interactions are central to the control of gene expression across all forms of life. The development of approaches to rigorously model such interactions has often been hindered both by a lack of quantitative binding data and by the difficulty in accounting for parameters relevant to the intracellular situation, such as DNA looping and thermodynamic non-ideality. Here, we review these considerations by developing a thermodynamically based mathematical model that attempts to simulate the functioning of an Escherichia coli expression system incorporating two of the best characterised prokaryotic DNA binding proteins, Lac repressor and lambda CI repressor. The key aim was to reproduce experimentally observed reporter gene activities arising from the expression of either wild-type CI repressor or one of three positive-control CI mutants. The model considers the role of several potentially important, but sometimes neglected, biochemical features, including DNA looping, macromolecular crowding and non-specific binding, and allowed us to obtain association constants for the binding of CI and its variants to a specific operator sequence.     

Analytical biochemistry of DNA–protein assemblies from crude cell extracts

Purification of specific DNA–protein complexes is a challenging task, as the involved interactions can be both electrostatic/H-bond and hydrophobic. The chromatographic stringency needed to obtain reasonable purifications uses salts and detergents. However, these components elicit the removal of proteins unspecifically bound to the chromatographic support itself, thus contaminating the purification products. In this work, a photocleavable linker connected the target oligonucleotidic sequence to the chromatographic beads so as to allow the irradiation-based release of the purified DNA–protein complexes off the beads. Our bioanalytical conditions were validated by purifying the tetracycline repressor protein onto a specific oligonucleotide. The purification factor was unprecedented, with a single contaminant. The robustness of our method was challenged by applying it to the purification of multiprotein assemblies forming onto DNA damage-mimicking oligonucleotides. The purified components were identified as well-known DNA repair proteins, and were shown to retain their enzymatic activities, as seen by monitoring DNA ligation products. Remarkably, kinase activities, also monitored, were found to be distinct on the beads and on the purified DNA–protein complexes, showing the benefits to uncouple the DNA–protein assemblies from the beads for a proper understanding of biochemical regulatory mechanisms involved in the DNA–protein assemblies.     

Dissecting protein-induced DNA looping dynamics in real time

Many proteins that interact with DNA perform or enhance their specific functions by binding simultaneously to multiple target sites, thereby inducing a loop in the DNA. The dynamics and energies involved in this loop formation influence the reaction mechanism. Tethered particle motion has proven a powerful technique to study in real time protein-induced DNA looping dynamics while minimally perturbing the DNA–protein interactions. In addition, it permits many single-molecule experiments to be performed in parallel. Using as a model system the tetrameric Type II restriction enzyme SfiI, that binds two copies of its recognition site, we show here that we can determine the DNA–protein association and dissociation steps as well as the actual process of protein-induced loop capture and release on a single DNA molecule. The result of these experiments is a quantitative reaction scheme for DNA looping by SfiI that is rigorously compared to detailed biochemical studies of SfiI looping dynamics. We also present novel methods for data analysis and compare and discuss these with existing methods. The general applicability of the introduced techniques will further enhance tethered particle motion as a tool to follow DNA–protein dynamics in real time.     

DNA-binding properties of a lac repressor mutant incapable of forming tetramers

The interaction of proteins bound to sites widely separated on the genome is a recurrent motif in both prokaryotic and eukaryotic regulatory systems. Lac repressor mediates the formation of "DNA loops" by the simultaneous interaction of a single protein tetramer with two DNA-binding sites. The DNA-binding properties of a Lac repressor mutant (LacIadi) deficient in the association of protein dimers to tetramers was investigated. The results of quantitative footprint and gel mobility-shift titrations suggest that the wild-type Lac repressor (LacI+) binds cooperatively to two operator sites separated by 11 helical turns on a linear DNA restriction fragment by the formation of a "looped complex." LacIadi binds to this two-site operator non-cooperatively and without formation of a looped complex. These results demonstrate that the dimer-tetramer association of LacI+ is directly responsible for its cooperative binding and its ability to mediate formation of a looped complex. The Iadi mutation disrupts the monomer-dimer as well as eliminating the dimer-tetramer association equilibria while the DNA binding affinity of LacIadi to a single site is unchanged relative to the wild-type protein. These results suggest that DNA binding and dimer-tetramer association are functionally unlinked. The similarity of the DNA-binding properties of LacIadi and Gal repressor, a protein believed to function by mediating the formation of a looped complex, are discussed.     

Influence of supercoiling and sequence context on operator DNA binding with lac repressor

The dissociation of the repressor-operator complex from a series of negatively supercoiled plasmid DNAs was examined as a function of the sequence context, orientation, and spacing. The plasmids were grouped into four classes, each with common sequence context. The highest dissociation rate constants were observed for the plasmids containing only a single operator (or pseudooperator) sequence, while approximately 10-fold lower rate constants were measured for plasmids with the I gene pseudooperator in conjunction with either the Z gene pseudooperator or the primary operator. Comparison of the behavior of these two classes of plasmids demonstrated the importance of two operator sequences and supported a model of DNA loop formation to stabilize the repressor-operator complex (Whitson, P. A., and Matthews, K. S. (1986) Biochemistry 25, 3845-3852; Whitson, P. A., Olson, J. S., and Matthews, K. S. (1986) Biochemistry 25, 3852-3858; Whitson, P. A., Hsieh, W. T., Wells, R. D., and Matthews, K. S. (1987) J. Biol. Chem. 262, 4943-4946; Kr?mer, H., Niem?ller, M., Amouyal, M., Revet, B., von Wilcken-Bergmann, B., and Müller-Hill, B. (1987) EMBO J. 6, 1481-1491). The third class, with intermediate dissociation rate constants, was comprised of plasmids which contained the primary operator and the higher affinity pseudooperator normally located in the Z gene. Neither the additional presence of the I gene pseudooperator nor the orientation of the primary operator relative to the Z gene pseudooperator significantly affected the dissociation rate constants. The binding characteristics of this group of plasmids demonstrated the essential role of the Z gene pseudooperator in the formation of intramolecular ternary complex and suggested an in vivo function for this pseudooperator. Plasmids containing two primary operator sequences were the class with lowest dissociation rate constants from lac repressor, and minimal effects of salt or spacing on dissociation of this class were observed. These data are consistent with formation of an intramolecular complex with a looped DNA segment stabilized by the combination of increased local concentration of binding sites and torsional stresses on the DNA which favor binding in supercoiled DNA.     

Trapping DNA-protein binding reactions with neutral osmolytes for the analysis by gel mobility shift and self-cleavage assays

We take advantage of our previous observation that neutral osmolytes can strongly slow down the rate of DNA–protein complex dissociation to develop a method that uses osmotic stress to ‘freeze’ mixtures of DNA–protein complexes and prevent further reaction enabling analysis of the products. We apply this approach to the gel mobility shift assay and use it to modify a self-cleavage assay that uses the nuclease activity of the restriction endonucleases to measure sensitively their specific binding to DNA. At sufficiently high concentrations of neutral osmolytes the cleavage reaction can be triggered at only those DNA fragments with initially bound enzyme. The self-cleavage assay allows measurement of binding equilibrium and kinetics directly in solution avoiding the intrinsic problems of gel mobility shift and filter binding assays while providing the same sensitivity level. Here we compare the self-cleavage and gel mobility shift assays applied to the DNA binding of EcoRI and BamHI restriction endonucleases. Initial results indicate that BamHI dissociation from its specific DNA sequence is strongly linked to water activity with the half-life time of the specific complex increasing ~20-fold from 0 to 1 osmolal betaine.     

Crystallographic analysis of Lac repressor bound to natural operator O1

Previous structures of Lac repressor bound to DNA used a fully symmetric "ideal" operator sequence that is missing the central G-C base-pair present in the three natural operator sequences. Here we have determined the X-ray crystal structure of a dimeric Lac repressor bound to a 22 base-pair DNA with the natural operator O1 sequence and the anti-inducer ONPF, at 4.0 A resolution. The natural operator is bent in the same way as the symmetric sequence, due to the binding of the hinge helices of the repressor to the minor groove at the central GCGG sequence of O1. Comparison of the structures of the repressor bound to the natural and symmetric operators shows very similar overall structures, with only slight rearrangements of the headpiece domains of the repressor. Analysis of crystals with iodinated DNA shows that the operator is uniquely positioned and allows for the sequence registration of the DNA relative to the repressor to be determined. The kink in the operator is centered between the left half-site and the central G-C base-pair of O1. Our results are most consistent with a previously proposed model in which, relative to the complex with the symmetric operator, the repressor accommodates binding to the natural operator sequence by shifting the position of the right headpiece by one base-pair step towards the center of O1.     

A closer view of the conformation of the Lac repressor bound to operator

Crystal structures of the Lac repressor, with and without isopropyithiogalactoside (IPTG), and the repressor bound to operator have provided a model for how the binding of the inducer reduces the affinity of the repressor for the operator. However, because of the low resolution of the operator-bound structure (4.8 A), the model for the allosteric transition was presented in terms of structural elements rather than in terms of side chain interactions. Here we have constructed a dimeric Lac repressor and determined its structure at 2.6 A resolution in complex with a symmetric operator and the anti-inducer orthonitrophenylfucoside (ONPF). The structure enables the induced (IPTG-bound) and repressed (operator-bound) conformations of the repressor to be compared in atomic detail. An extensive network of interactions between the DNA-binding and core domains of the repressor suggests a possible mechanism for the allosteric transition.     

Interaction of effecting ligands with lac repressor and repressor-operator complex.

The equilibrium association constants for the binding of a wide variety of effecting ligands of the lac repressor were measured by equilibrium dialysis. Also, detailed investigations of the apparent rate of dissociation of repressor-operator comples as a function of ligand concentration were carried out for several inducers and anti-inducers. The affinity of repressor-ligand comples for operator DNA was evaluated from the specific rate constants at saturating concentrations of effecting ligand. By fitting the experimental data depicting the functional dependence of the rate of dissociation upon ligand concentrations to calculated curves, assuming simple models of the induction mechanism, the equilibrium association constant for the binding of effecting ligand to repressor-operator comples was determined. Inducers reduce the affinity of lac repressor for operator DNA by a factor of approximately 1000 under standard conditions; the extent of destabilization depends on Mg2+ ion concentration. Anti-inducers increase the affinity of repressor for operator at most a factor of five. Only one neutral ligand, which binds to repressor without altering the stability of repressor-operator comples, was found. No homotropic or heterotropic interactions in the binding of effecting ligands either to repressor or to repressor-operator complex are evident.