Strong physical constraints on sequence-specific target location by proteins on DNA molecules

Structural Basis of Enhanced Facilitated Diffusion of DNA-Binding Protein in Crowded Cellular Milieu

Although the fast association between DNA-binding proteins (DBPs) and DNA is explained by a facilitated diffusion mechanism, in which DBPs adopt a weighted combination of three-dimensional diffusion and one-dimensional (1D) sliding and hopping modes of transportation, the role of cellular environment that contains many nonspecifically interacting proteins and other biomolecules is mostly overlooked. By performing large-scale computational simulations with an appropriately tuned model of protein and DNA in the presence of nonspecifically interacting bulk and DNA-bound crowders (genomic crowders), we demonstrate the structural basis of the enhanced facilitated diffusion of DBPs inside a crowded cellular milieu through, to our knowledge, novel 1D scanning mechanisms. In this one-dimensional scanning mode, the protein can float along the DNA under the influence of nonspecific interactions of bulk crowder molecules. The search mode is distinctly different compared to usual 1D sliding and hopping dynamics in which protein diffusion is regulated by the DNA electrostatics. In contrast, the presence of genomic crowders expedites the target search process by transporting the protein over DNA segments through the formation of a transient protein-crowder bridged complex. By analyzing the ruggedness of the associated potential energy landscape, we underpin the molecular origin of the kinetic advantages of these search modes and show that they successfully explain the experimentally observed acceleration of facilitated diffusion of DBPs by molecular crowding agents and crowder-concentration-dependent enzymatic activity of transcription factors. Our findings provide crucial insights into gene regulation kinetics inside the crowded cellular milieu.

Optimal Length of Conformational Transition Region in Protein Search for Targets on DNA

The starting point of many fundamental biological processes is associated with protein molecules finding and recognizing specific sites on DNA. However, despite a large number of experimental and theoretical studies on protein search for targets on DNA, many molecular aspects of underlying mechanisms are still not well understood. Experiments show that proteins bound to DNA can switch between slow recognition and fast search conformations. However, from a theoretical point of view, such conformational transitions should slow down the protein search for specific sites on DNA, in contrast to available experimental observations. In addition, experiments indicate that the nucleotide composition near the target site is more symmetrically homogeneous, leading to stronger effective interactions between proteins and DNA at these locations. However, as has been shown theoretically, this should also make the search less efficient, which is not observed. We propose a possible resolution of these problems by suggesting that conformational transitions occur only within a segment around the target where stronger interactions between proteins and DNA are observed. Two theoretical methods, based on continuum and discrete-state stochastic calculations, are developed, allowing us to obtain a comprehensive dynamic description for the protein search process in this system. The existence of an optimal length of the conformational transition zone with the shortest mean search time is predicted.

Crowding on DNA in Protein Search for Targets

Proteins searching and recognizing specific sites on DNA is required for initiating all major biological processes. While the details of the protein search for targets on DNA in purified in vitro systems are reasonably well understood, the situation in real cells is much less clear. The presence of other types of molecules on DNA should prevent reaching the targets, but experiments show that, surprisingly, the molecular crowding on DNA influences the search dynamics much less than expected. We develop a theoretical method that allowed us to clarify the mechanisms of the protein search on DNA in the presence of crowding. It is found that the dimensionality of the search trajectories specifies whether the crowding will affect the target finding. For 3D search pathways it is minimal, while the strongest effect is for 1D search pathways when the crowding particle can block the search. In addition, for 1D search we determined that the critical parameter is a mobility of crowding agents: highly mobile molecules do not affect the search dynamics, while the slow particles can significantly slow down the process. Physical-chemical explanations of the observed phenomena are presented. Our theoretical predictions thus explain the experimental observations, and they are also supported by extensive numerical simulations.

Biosynthesis and Function of Modified Bases in Bacteria and Their Viruses

Naturally occurring modification of the canonical A, G, C, and T bases can be found in the DNA of cellular organisms and viruses from all domains of life. Bacterial viruses (bacteriophages) are a particularly rich but still underexploited source of such modified variant nucleotides. The modifications conserve the coding and base-pairing functions of DNA, but add regulatory and protective functions. In prokaryotes, modified bases appear primarily to be part of an arms race between bacteriophages (and other genomic parasites) and their hosts, although, as in eukaryotes, some modifications have been adapted to convey epigenetic information. The first half of this review catalogs the identification and diversity of DNA modifications found in bacteria and bacteriophages. What is known about the biogenesis, context, and function of these modifications are also described. The second part of the review places these DNA modifications in the context of the arms race between bacteria and bacteriophages. It focuses particularly on the defense and counter-defense strategies that turn on direct recognition of the presence of a modified base. Where modification has been shown to affect other DNA transactions, such as expression and chromosome segregation, that is summarized, with reference to recent reviews.

Searching target sites on DNA by proteins: Role of DNA dynamics under confinement

DNA-binding proteins (DBPs) rapidly search and specifically bind to their target sites on genomic DNA in order to trigger many cellular regulatory processes. It has been suggested that the facilitation of search dynamics is achieved by combining 3D diffusion with one-dimensional sliding and hopping dynamics of interacting proteins. Although, recent studies have advanced the knowledge of molecular determinants that affect one-dimensional search efficiency, the role of DNA molecule is poorly understood. In this study, by using coarse-grained simulations, we propose that dynamics of DNA molecule and its degree of confinement due to cellular crowding concertedly regulate its groove geometry and modulate the inter-communication with DBPs. Under weak confinement, DNA dynamics promotes many short, rotation-decoupled sliding events interspersed by hopping dynamics. While this results in faster 1D diffusion, associated probability of missing targets by jumping over them increases. In contrast, strong confinement favours rotation-coupled sliding to locate targets but lacks structural flexibility to achieve desired specificity. By testing under physiological crowding, our study provides a plausible mechanism on how DNA molecule may help in maintaining an optimal balance between fast hopping and rotation-coupled sliding dynamics, to locate target sites rapidly and form specific complexes precisely.

Method for the analysis of contribution of sliding and hopping to a facilitated diffusion of DNA-binding protein: Application to in vivo data

DNA-binding protein searches for its target, a specific site on DNA, by means of diffusion. The search process consists of many recurrent steps of one-dimensional diffusion (sliding) along the DNA chain and three-dimensional diffusion (hopping) after dissociation of a protein from the DNA chain. Here we propose a computational method that allows extracting the contribution of sliding and hopping to the search process in vivo from the measurements of the kinetics of the target search by the lac repressor in Escherichia coli [P. Hammar et al., Science 336, 1595 (2012)]. The method combines lattice Monte Carlo simulations with the Brownian excursion theory and includes explicitly steric constraints for hopping due to the helical structure of DNA. The simulation results including all experimental data reveal that the in vivo target search is dominated by sliding. The short-range hopping to the same base pair interrupts one-dimensional sliding while long-range hopping does not contribute significantly to the kinetics of the search of the target in vivo.

Real sequence effects on the search dynamics of transcription factors on DNA

Recent experiments show that transcription factors (TFs) indeed use the facilitated diffusion mechanism to locate their target sequences on DNA in living bacteria cells: TFs alternate between sliding motion along DNA and relocation events through the cytoplasm. From simulations and theoretical analysis we study the TF-sliding motion for a large section of the DNA-sequence of a common E. coli strain, based on the two-state TF-model with a fast-sliding search state and a recognition state enabling target detection. For the probability to detect the target before dissociating from DNA the TF-search times self-consistently depend heavily on whether or not an auxiliary operator (an accessible sequence similar to the main operator) is present in the genome section. Importantly, within our model the extent to which the interconversion rates between search and recognition states depend on the underlying nucleotide sequence is varied. A moderate dependence maximises the capability to distinguish between the main operator and similar sequences. Moreover, these auxiliary operators serve as starting points for DNA looping with the main operator, yielding a spectrum of target detection times spanning several orders of magnitude. Auxiliary operators are shown to act as funnels facilitating target detection by TFs.

Obstacles may facilitate and direct DNA search by proteins

DNA recognition by DNA-binding proteins (DBPs), which is a pivotal event in most gene regulatory processes, is often preceded by an extensive search for the correct site. A facilitated diffusion process in which a DBP combines three-dimensional diffusion in solution with one-dimensional sliding along DNA has been suggested to explain how proteins can locate their target sites on DNA much faster than predicted by three-dimensional diffusion alone. Although experimental and theoretical studies have recently advanced understanding of the biophysical principles underlying the search mechanism, the process under in vivo cellular conditions is poorly understood. In this study, we used various computational approaches to explore how the presence of obstacle proteins on the DNA influences search efficiency. At a low obstacle occupancy (i.e., when few obstacles occupy sites on the DNA), sliding by the searching DBP may be confined, which may impair search efficiency. The obstacles, however, can be bypassed during hopping events, and the number of bypasses is larger for higher obstacle occupancies. Dynamism on the part of the obstacles may even further facilitate search kinetics. Our study shows that the nature and efficiency of the search process may be governed not only by the intrinsic properties of the DBP and the salt concentration of the medium, but also by the in vivo association of DNA with other macromolecular obstacles, their location, and occupancy.

Quantitative influence of macromolecular crowding on gene regulation kinetics

We introduce macromolecular crowding quantitatively into the model for kinetics of gene regulation in Escherichia coli. We analyse and compute the specific-site searching time for 180 known transcription factors (TFs) regulating 1300 operons. The time is between 160 s (e.g. for SoxS M_w = 12.91 kDa) and 1550 s (e.g. for PepA₆ of M_w = 329.28 kDa). Diffusion coefficients for one-dimensional sliding are between for large proteins up to for small monomers or dimers. Three-dimensional diffusion coefficients in the cytoplasm are 2 orders of magnitude larger than 1D sliding coefficients, nevertheless the sliding enhances the binding rates of TF to specific sites by 1–2 orders of magnitude. The latter effect is due to ubiquitous non-specific binding. We compare the model to experimental data for LacI repressor and find that non-specific binding of the protein to DNA is activation- and not diffusion-limited. We show that the target location rate by LacI repressor is optimized with respect to microscopic rate constant for association to non-specific sites on DNA. We analyse the effect of oligomerization of TFs and DNA looping effects on searching kinetics. We show that optimal searching strategy depends on TF abundance.

The effects of transcription factor competition on gene regulation

Transcription factor (TF) molecules translocate by facilitated diffusion (a combination of 3D diffusion around and 1D random walk on the DNA). Despite the attention this mechanism received in the last 40 years, only a few studies investigated the influence of the cellular environment on the facilitated diffusion mechanism and, in particular, the influence of "other" DNA binding proteins competing with the TF molecules for DNA space. Molecular crowding on the DNA is likely to influence the association rate of TFs to their target site and the steady state occupancy of those sites, but it is still not clear how it influences the search in a genome-wide context, when the model includes biologically relevant parameters (such as: TF abundance, TF affinity for DNA and TF dynamics on the DNA). We performed stochastic simulations of TFs performing the facilitated diffusion mechanism, and considered various abundances of cognate and non-cognate TFs. We show that, for both obstacles that move on the DNA and obstacles that are fixed on the DNA, changes in search time are not statistically significant in case of biologically relevant crowding levels on the DNA. In the case of non-cognate proteins that slide on the DNA, molecular crowding on the DNA always leads to statistically significant lower levels of occupancy, which may confer a general mechanism to control gene activity levels globally. When the "other" molecules are immobile on the DNA, we found a completely different behavior, namely: the occupancy of the target site is always increased by higher molecular crowding on the DNA. Finally, we show that crowding on the DNA may increase transcriptional noise through increased variability of the occupancy time of the target sites.

The influence of transcription factor competition on the relationship between occupancy and affinity

Transcription factors (TFs) are proteins that bind to specific sites on the DNA and regulate gene activity. Identifying where TF molecules bind and how much time they spend on their target sites is key to understanding transcriptional regulation. It is usually assumed that the free energy of binding of a TF to the DNA (the affinity of the site) is highly correlated to the amount of time the TF remains bound (the occupancy of the site). However, knowing the binding energy is not sufficient to infer actual binding site occupancy. This mismatch between the occupancy predicted by the affinity and the observed occupancy may be caused by various factors, such as TF abundance, competition between TFs or the arrangement of the sites on the DNA. We investigated the relationship between the affinity of a TF for a set of binding sites and their occupancy. In particular, we considered the case of the transcription factor lac repressor (lacI) in E.coli, and performed stochastic simulations of the TF dynamics on the DNA for various combinations of lacI abundance and competing TFs that contribute to macromolecular crowding. We also investigated the relationship of site occupancy and the information content of position weight matrices (PWMs) used to represent binding sites. Our results showed that for medium and high affinity sites, TF competition does not play a significant role for genomic occupancy except in cases when the abundance of the TF is significantly increased, or when the PWM displays relatively low information content. Nevertheless, for medium and low affinity sites, an increase in TF abundance (for both cognate and non-cognate molecules) leads to an increase in occupancy at several sites.

Computational models for large-scale simulations of facilitated diffusion

The binding of site-specific transcription factors to their genomic target sites is a key step in gene regulation. While the genome is huge, transcription factors belong to the least abundant protein classes in the cell. It is therefore fascinating how short the time frame is that they require to home in on their target sites. The underlying search mechanism is called facilitated diffusion and assumes a combination of three-dimensional diffusion in the space around the DNA combined with one-dimensional random walk on it. In this review, we present the current understanding of the facilitated diffusion mechanism and identify questions that lack a clear or detailed answer. One way to investigate these questions is through stochastic simulation and, in this manuscript, we support the idea that such simulations are able to address them. Finally, we review which biological parameters need to be included in such computational models in order to obtain a detailed representation of the actual process.

System size reduction in stochastic simulations of the facilitated diffusion mechanism

Background

Site-specific Transcription Factors (TFs) are proteins that bind to specific sites on the DNA and control the activity of a target gene by enhancing or decreasing the rate at which the gene is transcribed by RNA polymerase. The process by which TF molecules locate their target sites is a key component of transcriptional regulation. Therefore it is essential to gain insight into the mechanisms by which TFs search for the target sites.

Research in this area uses experimental and analytical approaches, but also stochastic simulations of the search process. Previous work based on stochastic simulations focussed only on short sequences, primarily for reasons of technical feasibility. Many of these studies had to disregard possible biases introduced by reducing a genome-wide system to a smaller subsystem. In particular, we identified crucial parameters that require adjustment, which were not adequately changed in these previous studies.

Results

We investigated several methods that adequately adapt the parameters of stochastic simulations of the facilitated diffusion, when the full sequence space is reduced to smaller regions of interest. We found two methods that scale the system accordingly: the copy number model and the association rate model. We systematically compared the results produced by simulations of the subsystem with respect to the original system. Our results confirmed that the copy number model is adequate only for high abundance TFs, while for low abundance TFs the association rate model is the only one that reproduces with high accuracy the results of the full system.

Conclusions

We propose a strategy to reduce the size of the system that adequately adapts important parameters to capture the behaviour of the full system. This enables correct simulations of a smaller sequence space (which can be as small as 100 Kbp) and, thus, provides independence from computationally intensive genome-wide simulations of the facilitated diffusion mechanism.

A comprehensive computational model of facilitated diffusion in prokaryotes

MOTIVATION

Gene activity is mediated by site-specific transcription factors (TFs). Their binding to defined regions in the genome determines the rate at which their target genes are transcribed.

RESULTS

We present a comprehensive computational model of the search process of TF for their genomic target site(s). The computational model considers: the DNA sequence, various TF species and the interaction of the individual molecules with the DNA or between themselves. We also demonstrate a systematic approach how to parametrize the system using available experimental data.

Is photocleavage of DNA by YOYO-1 using a synchrotron radiation light source sequence dependent?

The photocleavage of double-stranded and single-stranded DNA by the fluorescent dye YOYO-1 was investigated in real time by using the synchrotron radiation light source ASTRID (ISA, Denmark) both to initiate the reaction and to monitor its progress using Couette flow linear dichroism (LD) throughout the irradiation period. The dependence of LD signals on DNA sequences and on time in the intense light beam was explored and quantified for single-stranded poly(dA), poly[(dA-dT)(2)], calf thymus DNA (ctDNA) and Micrococcus luteus DNA (mlDNA). The DNA and ligand regions of the spectrum showed different LD kinetic behaviors, and there was significant sequence dependence of the kinetics. However, in contrast to expectations from the literature, we found that poly(dA), mlDNA, low salt ctDNA and low salt poly[(dA-dT)(2)] all had significant populations of groove-bound YOYO. It seems that this mode was predominantly responsible for the catalysis of DNA cleavage. In homopolymeric DNAs, intercalated YOYO was unable to cleave DNA. In mixed-sequence DNAs the data suggest that YOYO in some but not all intercalated binding sites can cause cleavage. It is also likely that cleavage occurs at transient single-stranded regions. The reaction rates for a 100 mA beam current of 0.5-μW power varied from 0.6 h(-1) for single-stranded poly(dA) to essentially zero for low salt poly[(dG-dC)(2)] and high salt poly[(dA-dT)(2)]. At the conclusion of the experiments with each kind of DNA, uncleaved DNA with intercalated YOYO remained.

Dancing on DNA: kinetic aspects of search processes on DNA

Recognition and binding of specific sites on DNA by proteins is central for many cellular functions such as transcription, replication, and recombination. In the search for its target site, the DNA-associated protein is facing both thermodynamic and kinetic difficulties. The thermodynamic challenge lies in recognizing and tightly binding a cognate (specific) site among the billions of other (non-specific) sequences on the DNA. The kinetic difficulty lies in finding a cognate site in mere seconds amidst the crowded cellular environment that is filled with other DNA sequences and proteins. Herein, we discuss the history of the DNA search problem, the theoretical background and the various experimental methods used to study the kinetics of proteins searching for target sites on DNA.

DNA translocation by type III restriction enzymes: a comparison of current models of their operation derived from ensemble and single-molecule measurements

Much insight into the interactions of DNA and enzymes has been obtained using a number of single-molecule techniques. However, recent results generated using two of these techniques-atomic force microscopy (AFM) and magnetic tweezers (MT)-have produced apparently contradictory results when applied to the action of the ATP-dependent type III restriction endonucleases on DNA. The AFM images show extensive looping of the DNA brought about by the existence of multiple DNA binding sites on each enzyme and enzyme dimerisation. The MT experiments show no evidence for looping being a requirement for DNA cleavage, but instead support a diffusive sliding of the enzyme on the DNA until an enzyme-enzyme collision occurs, leading to cleavage. Not only do these two methods appear to disagree, but also the models derived from them have difficulty explaining some ensemble biochemical results on DNA cleavage. In this 'Survey and Summary', we describe several different models put forward for the action of type III restriction enzymes and their inadequacies. We also attempt to reconcile the different models and indicate areas for further experimentation to elucidate the mechanism of these enzymes.

The kinetics of YOYO-1 intercalation into single molecules of double-stranded DNA

The cyanine dye, YOYO-1, has frequently been used in single DNA molecule imaging work to stain double-stranded DNA as it fluoresces strongly when bound. The binding of YOYO-1 lengthens the DNA due to bis-intercalation. We have investigated the kinetics of binding, via this increase in DNA length, for single, hydrodynamically-stretched molecules of lambda DNA observed via Total Internal Reflection Fluorescence (TIRF) microscopy. The rate and degree of lengthening in 40mM NaHCO(3) (pH 8.0) buffer depend upon the free dye concentration with the reaction taking several minutes to reach completion even in relatively high, 40nM, concentrations of YOYO-1. In the absence of overstretching of the DNA molecule, we determine the second order rate constant to be 3.8±0.7×10(5)s(-1)M(-1), the dissociation constant to be 12.1±3.4nM and the maximum DNA molecule extension to be 36±4%. The intercalation time constant (inverse of the pseudo-first order rate constant), τ, decreased from 309 to 62s as YOYO-1 levels increased from 10 to 40nM. The kinetics of binding help with interpretation of the behavior of DNA-YOYO-1 complexes when overstretched and establish defined conditions for the preparation of DNA-YOYO-1 complexes.

Statistical-mechanical lattice models for protein-DNA binding in chromatin

Statistical-mechanical lattice models for protein-DNA binding are well established as a method to describe complex ligand binding equilibria measured in vitro with purified DNA and protein components. Recently, a new field of applications has opened up for this approach since it has become possible to experimentally quantify genome-wide protein occupancies in relation to the DNA sequence. In particular, the organization of the eukaryotic genome by histone proteins into a nucleoprotein complex termed chromatin has been recognized as a key parameter that controls the access of transcription factors to the DNA sequence. New approaches have to be developed to derive statistical-mechanical lattice descriptions of chromatin-associated protein-DNA interactions. Here, we present the theoretical framework for lattice models of histone-DNA interactions in chromatin and investigate the (competitive) DNA binding of other chromosomal proteins and transcription factors. The results have a number of applications for quantitative models for the regulation of gene expression.

Modeling promoter search by E. coli RNA polymerase: one-dimensional diffusion in a sequence-dependent energy landscape

We present a biophysical model of promoter search by Escherichia coli RNA polymerase. We use an unconventional weight matrix derived from promoter strength data to extract the energy landscape common to a large set of known promoters. This exhibits a continuous strengthening of the binding energy when approaching the transcription start site from either side. During promoter search, the RNA polymerase slides along the DNA double helix (one-dimensional diffusion) after randomly binding to it. We discuss the possibility that the sliding has a sequence-dependent component, which implies that the energy landscape influences the movement with respect to speed, direction and efficiency. Based on this assumption, we relate the obtained energy landscape around the promoters to the one-dimensional diffusion of the RNA polymerase. Our analytical results suggest that the sequence-dependent random walk slows down and gets directed upon entering a region of 500 bp around the transcription start site, which significantly increases the efficiency of promoter search. These results may explain how the RNA polymerase is able to find the promoter in biologically relevant times out of a vast excess of non-target sites. Moreover, they provide evidence for a sequence-dependent component of one-dimensional diffusion.

An end to 40 years of mistakes in DNA-protein association kinetics?

Proteins that bind to specific sequences in long DNA molecules have to locate their target sites amid myriad alternative sequences, yet they do so at remarkably rapid rates, sometimes approaching 10(10) M(-1) x s(-1). Hence, it has been asserted widely that binding to specific DNA sites can surpass the maximal rate for 3D (three-dimensional) diffusion through solution and that this could only be accounted for by a reduction in the dimensionality of the search for the target in effect by 1D (one-dimensional) diffusion (or 'sliding') along the DNA contour. It will be shown here that there is, in fact, no known example of a protein binding to a specific DNA site at a rate above the diffusion limit, and that the rapidity of these reactions is due primarily to electrostatic interactions between oppositely charged molecules. It will also be shown that, contrary to popular belief, reduced dimensionality does not, in general, increase the rate of target-site location but instead reduces it. Finally, it will be demonstrated that proteins locate their target sites primarily by multiple dissociation/reassociation events to other (nearby or distant) sites within the same DNA molecule, and that 1D diffusion is limited to local searches covering approximately 50 bp around each landing site.

Sliding and jumping of single EcoRV restriction enzymes on non-cognate DNA

The restriction endonuclease EcoRV can rapidly locate a short recognition site within long non-cognate DNA using 'facilitated diffusion'. This process has long been attributed to a sliding mechanism, in which the enzyme first binds to the DNA via nonspecific interaction and then moves along the DNA by 1D diffusion. Recent studies, however, provided evidence that 3D translocations (hopping/jumping) also help EcoRV to locate its target site. Here we report the first direct observation of sliding and jumping of individual EcoRV molecules along nonspecific DNA. Using fluorescence microscopy, we could distinguish between a slow 1D diffusion of the enzyme and a fast translocation mechanism that was demonstrated to stem from 3D jumps. Salt effects on both sliding and jumping were investigated, and we developed numerical simulations to account for both the jump frequency and the jump length distribution. We deduced from our study the 1D diffusion coefficient of EcoRV, and we estimated the number of jumps occurring during an interaction event with nonspecific DNA. Our results substantiate that sliding alternates with hopping/jumping during the facilitated diffusion of EcoRV and, furthermore, set up a framework for the investigation of target site location by other DNA-binding proteins.

Spatial effects on the speed and reliability of protein-DNA search

Strong experimental and theoretical evidence shows that transcription factors (TFs) and other specific DNA-binding proteins find their sites using a two-mode search: alternating between three-dimensional (3D) diffusion through the cell and one-dimensional (1D) sliding along the DNA. We show that, due to the 1D component of the search process, the search time of a TF can depend on the initial position of the TF. We formalize this effect by discriminating between two types of searches: global and local. Using analytical calculations and simulations, we estimate how close a TF and binding site need to be to make a local search likely. We then use our model to interpret the wide range of experimental measurements of this parameter. We also show that local and global searches differ significantly in average search time and the variability of search time. These results lead to a number of biological implications, including suggestions of how prokaryotes achieve rapid gene regulation and the relationship between the search mechanism and noise in gene expression. Lastly, we propose a number of experiments to verify the existence and quantify the extent of spatial effects on the TF search process in prokaryotes.

Thermodynamic characterization of specific interactions between the human Lon protease and G-quartet DNA

Lon is an ATP-powered protease that binds DNA. However, the function of DNA binding by Lon remains elusive. Studies suggest that human Lon (hLon) binds preferentially to a G-rich single-stranded DNA (ssDNA) sequence overlapping the light strand promoter of mitochondrial DNA. This sequence is contained within a 24-base oligonucleotide referred to as LSPas. Here, we use biochemical and biophysical approaches to elucidate the structural properties of ssDNAs bound by hLon, as well as the thermodynamics of DNA binding by hLon. Electrophoretic mobility shift assay and circular dichroism show that ssDNAs with a propensity for forming parallel G-quartets are specifically bound by hLon. Isothermal titration calorimetry demonstrates that hLon binding to LSPas is primarily driven by enthalpy change associated with a significant reduction in heat capacity. Differential scanning calorimetry pinpoints an excess heat capacity upon hLon binding to LSPas. By contrast, hLon binding to an 8-base G-rich core sequence is entropically driven with a relatively negligible change in heat capacity. A considerable enhancement of thermal stability accompanies hLon binding to LSPas as compared to the G-rich core. Taken together, these data support the notion that hLon binds G-quartets through rigid-body binding and that binding to LSPas is coupled with structural adaptation.

General transfer matrix formalism to calculate DNA-protein-drug binding in gene regulation: application to OR operator of phage lambda

The transfer matrix methodology is proposed as a systematic tool for the statistical-mechanical description of DNA-protein-drug binding involved in gene regulation. We show that a genetic system of several cis-regulatory modules is calculable using this method, considering explicitly the site-overlapping, competitive, cooperative binding of regulatory proteins, their multilayer assembly and DNA looping. In the methodological section, the matrix models are solved for the basic types of short- and long-range interactions between DNA-bound proteins, drugs and nucleosomes. We apply the matrix method to gene regulation at the O(R) operator of phage lambda. The transfer matrix formalism allowed the description of the lambda-switch at a single-nucleotide resolution, taking into account the effects of a range of inter-protein distances. Our calculations confirm previously established roles of the contact CI-Cro-RNAP interactions. Concerning long-range interactions, we show that while the DNA loop between the O(R) and O(L) operators is important at the lysogenic CI concentrations, the interference between the adjacent promoters P(R) and P(RM) becomes more important at small CI concentrations. A large change in the expression pattern may arise in this regime due to anticooperative interactions between DNA-bound RNA polymerases. The applicability of the matrix method to more complex systems is discussed.