Interconvertible lac repressor-DNA loops revealed by single-molecule experiments

Optimizing the random search of a finite-lived target by a Lévy flight

In many random search processes of interest in chemistry, biology, or during rescue operations, an entity must find a specific target site before the latter becomes inactive, no longer available for reaction or lost. We present exact results on a minimal model system, a one-dimensional searcher performing a discrete time random walk, or Lévy flight. In contrast with the case of a permanent target, the capture probability and the conditional mean first passage time can be optimized. The optimal Lévy index takes a nontrivial value, even in the long lifetime limit, and exhibits an abrupt transition as the initial distance to the target is varied. Depending on the target lifetime, this transition is discontinuous or continuous, separated by a nonconventional tricritical point. These results pave the way to the optimization of search processes under time constraints.

A DNA origami-based device for investigating DNA bending proteins by transmission electron microscopy

The DNA origami technique offers precise positioning of nanoscale objects with high accuracy. This has facilitated the development of DNA origami-based functional nanomechanical devices that enable the investigation of DNA-protein interactions at the single particle level. Herein, we used the DNA origami technique to fabricate a nanoscale device for studying DNA bending proteins. For a proof of concept, we used TATA-box binding protein (TBP) to evaluate our approach. Upon binding to the TATA box, TBP causes a bend to DNA of ∼90°. Our device translates this bending into an angular change that is readily observable with a conventional transmission electron microscope (TEM). Furthermore, we investigated the roles of transcription factor II A (TF(II)A) and transcription factor II B (TF(II)B). Our results indicate that TF(II)A introduces additional bending, whereas TF(II)B does not significantly alter the TBP-DNA structure. Our approach can be readily adopted to a wide range of DNA-bending proteins and will aid the development of DNA-origami-based devices tailored for the investigation of DNA-protein interactions.

Designed architectural proteins that tune DNA looping in bacteria

Architectural proteins alter the shape of DNA. Some distort the double helix by introducing sharp kinks. This can serve to relieve strain in tightly-bent DNA structures. Here, we design and test artificial architectural proteins based on a sequence-specific Transcription Activator-like Effector (TALE) protein, either alone or fused to a eukaryotic high mobility group B (HMGB) DNA-bending domain. We hypothesized that TALE protein binding would stiffen DNA to bending and twisting, acting as an architectural protein that antagonizes the formation of small DNA loops. In contrast, fusion to an HMGB domain was hypothesized to generate a targeted DNA-bending architectural protein that facilitates DNA looping. We provide evidence from Escherichia coli Lac repressor gene regulatory loops supporting these hypotheses in living bacteria. Both data fitting to a thermodynamic DNA looping model and sophisticated molecular modeling support the interpretation of these results. We find that TALE protein binding inhibits looping by stiffening DNA to bending and twisting, while the Nhp6A domain enhances looping by bending DNA without introducing twisting flexibility. Our work illustrates artificial approaches to sculpt DNA geometry with functional consequences. Similar approaches may be applicable to tune the stability of small DNA loops in eukaryotes.

Protamine loops DNA in multiple steps

Protamine proteins dramatically condense DNA in sperm to almost crystalline packing levels. Here, we measure the first step in the in vitro pathway, the folding of DNA into a single loop. Current models for DNA loop formation are one-step, all-or-nothing models with a looped state and an unlooped state. However, when we use a Tethered Particle Motion (TPM) assay to measure the dynamic, real-time looping of DNA by protamine, we observe the presence of multiple folded states that are long-lived (∼100 s) and reversible. In addition, we measure folding on DNA molecules that are too short to form loops. This suggests that protamine is using a multi-step process to loop the DNA rather than a one-step process. To visualize the DNA structures, we used an Atomic Force Microscopy (AFM) assay. We see that some folded DNA molecules are loops with a ∼10-nm radius and some of the folded molecules are partial loops—c-shapes or s-shapes—that have a radius of curvature of ∼10 nm. Further analysis of these structures suggest that protamine is bending the DNA to achieve this curvature rather than increasing the flexibility of the DNA. We therefore conclude that protamine loops DNA in multiple steps, bending it into a loop.

Sequence-dependent dynamics of synthetic and endogenous RSSs in V(D)J recombination

Developing lymphocytes of jawed vertebrates cleave and combine distinct gene segments to assemble antigen-receptor genes. This process called V(D)J recombination that involves the RAG recombinase binding and cutting recombination signal sequences (RSSs) composed of conserved heptamer and nonamer sequences flanking less well-conserved 12- or 23-bp spacers. Little quantitative information is known about the contributions of individual RSS positions over the course of the RAG-RSS interaction. We employ a single-molecule method known as tethered particle motion to track the formation, lifetime and cleavage of individual RAG-12RSS-23RSS paired complexes (PCs) for numerous synthetic and endogenous 12RSSs. We reveal that single-bp changes, including in the 12RSS spacer, can significantly and selectively alter PC formation or the probability of RAG-mediated cleavage in the PC. We find that some rarely used endogenous gene segments can be mapped directly to poor RAG binding on their adjacent 12RSSs. Finally, we find that while abrogating RSS nicking with Ca2+ leads to substantially shorter PC lifetimes, analysis of the complete lifetime distributions of any 12RSS even on this reduced system reveals that the process of exiting the PC involves unidentified molecular details whose involvement in RAG-RSS dynamics are crucial to quantitatively capture kinetics in V(D)J recombination.

First Hitting Times to Intermittent Targets

In noisy environments such as the cell, many processes involve target sites that are often hidden or inactive, and thus not always available for reaction with diffusing entities. To understand reaction kinetics in these situations, we study the first hitting time statistics of a one-dimensional Brownian particle searching for a target site that switches stochastically between visible and hidden phases. At high crypticity, an unexpected rate limited power-law regime emerges for the first hitting time density, which markedly differs from the classic t^{-3/2} scaling for steady targets. Our problem admits an asymptotic mapping onto a mixed, or Robin, boundary condition. Similar results are obtained with non-Markov targets and particles diffusing anomalously.

The Contribution of Backbone Electrostatic Repulsion to DNA Mechanical Properties is Length-Scale-Dependent

The mechanics of DNA bending is crucially related to many vital biological processes. Recent experiments reported anomalous flexibility for DNA on short length scales, calling into doubt the validity of the harmonic worm-like chain (WLC) model in this region. In the present work, we systematically probed the bending dynamics of DNA at different length scales. In contrast to the remarkable deviation from the WLC description for DNA duplexes of less than three helical turns, our atomistic studies indicate that the neutral "null isomer" behaves in accord with the ideal elastic WLC and exhibits a uniform decay for the directional correlation of local bending. The backbone neutralization weakens the anisotropy in the effective bending preference and the helical periodicity of bend correlation that have previously been observed for normal DNA. The contribution of electrostatic repulsion to stretching cooperativity and the mechanical properties of DNA strands is length-scale-dependent: the phosphate neutralization increases the stiffness of DNA below two helical turns, but it is decreased for longer strands. We find that DNA rigidity is largely determined by base pair stacking, with electrostatic interactions contributing only around 10% of the total persistence length.

Protein-mediated looping of DNA under tension requires supercoiling

Protein-mediated DNA looping is ubiquitous in chromatin organization and gene regulation, but to what extent supercoiling or nucleoid associated proteins promote looping is poorly understood. Using the lac repressor (LacI), a paradigmatic loop-mediating protein, we measured LacI-induced looping as a function of either supercoiling or the concentration of the HU protein, an abundant nucleoid protein in Escherichia coli. Negative supercoiling to physiological levels with magnetic tweezers easily drove the looping probability from 0 to 100% in single DNA molecules under slight tension that likely exists in vivo. In contrast, even saturating (micromolar) concentrations of HU could not raise the looping probability above 30% in similarly stretched DNA or 80% in DNA without tension. Negative supercoiling is required to induce significant looping of DNA under any appreciable tension.

Comparison between Effects of Retroactivity and Resource Competition upon Change in Downstream Reporter Genes of Synthetic Genetic Circuits

Reporter genes have contributed to advancements in molecular biology. Binding of an upstream regulatory protein to a downstream reporter promoter allows quantification of the activity of the upstream protein produced from the corresponding gene. In studies of synthetic biology, analyses of reporter gene activities ensure control of the cell with synthetic genetic circuits, as achieved using a combination of in silico and in vivo experiments. However, unexpected effects of downstream reporter genes on upstream regulatory genes may interfere with in vivo observations. This phenomenon is termed as retroactivity. Using in silico and in vivo experiments, we found that a different copy number of regulatory protein-binding sites in a downstream gene altered the upstream dynamics, suggesting retroactivity of reporters in this synthetic genetic oscillator. Furthermore, by separating the two sources of retroactivity (titration of the component and competition for degradation), we showed that, in the dual-feedback oscillator, the level of the fluorescent protein reporter competing for degradation with the circuits' components is important for the stability of the oscillations. Altogether, our results indicate that the selection of reporter promoters using a combination of in silico and in vivo experiments is essential for the advanced design of genetic circuits.

Phenotypic instability in fungi

Fungi are prone to phenotypic instability, that is, the vegetative phase of these organisms, be they yeasts or molds, undergoes frequent switching between two or more behaviors, often with different morphologies, but also sometime having different physiologies without any obvious morphological outcome. In the context of industrial utilization of fungi, this can have a negative impact on the maintenance of strains and/or on their productivity. Instabilities have been shown to result from various mechanisms, either genetic or epigenetic. This chapter will review different types of instabilities and discuss some lesser-known ones, mostly in filamentous fungi, while it will direct readers to additional literature in the case of well-known phenomena such as the amyloid prions or fungal senescence. It will present in depth the "white/opaque" switch of Candida albicans and the "crippled growth" degeneration of the model fungus Podospora anserina. These are two of the most thoroughly studied epigenetic phenotypic switches. I will also discuss the "sectors" presented by many filamentous ascomycetes, for which a prion-based model exists but is not demonstrated. Finally, I will also describe intriguing examples of phenotypic instability for which an explanation has yet to be provided.

Tethered multifluorophore motion reveals equilibrium transition kinetics of single DNA double helices

Understanding cellular functions and dysfunctions often begins with quantifying the interactions between the binding partners involved in the processes. Learning about the kinetics of the interactions is of particular importance to understand the dynamics of cellular processes. We created a tethered multifluorophore motion assay using DNA origami that enables over 1-hour-long recordings of the statistical binding and unbinding of single pairs of biomolecules directly in equilibrium. The experimental concept is simple and the data interpretation is very direct, which makes the system easy to use for a wide variety of researchers. Due to the modularity and addressability of the DNA origami-based assay, our system may be readily adapted to study various other molecular interactions.

We describe a tethered multifluorophore motion assay based on DNA origami for revealing bimolecular reaction kinetics on the single-molecule level. Molecular binding partners may be placed at user-defined positions and in user-defined stoichiometry; and binding states are read out by tracking the motion of quickly diffusing fluorescent reporter units. Multiple dyes per reporter unit enable singe-particle observation for more than 1 hour. We applied the system to study in equilibrium reversible hybridization and dissociation of complementary DNA single strands as a function of tether length, cation concentration, and sequence. We observed up to hundreds of hybridization and dissociation events per single reactant pair and could produce cumulative statistics with tens of thousands of binding and unbinding events. Because the binding partners per particle do not exchange, we could also detect subtle heterogeneity from molecule to molecule, which enabled separating data reflecting the actual target strand pair binding kinetics from falsifying influences stemming from chemically truncated oligonucleotides. Our data reflected that mainly DNA strand hybridization, but not strand dissociation, is affected by cation concentration, in agreement with previous results from different assays. We studied 8-bp-long DNA duplexes with virtually identical thermodynamic stability, but different sequences, and observed strongly differing hybridization kinetics. Complementary full-atom molecular-dynamics simulations indicated two opposing sequence-dependent phenomena: helical templating in purine-rich single strands and secondary structures. These two effects can increase or decrease, respectively, the fraction of strand collisions leading to successful nucleation events for duplex formation.

The mechanics of DNA loops bridged by proteins unveiled by single-molecule experiments

Protein-induced DNA bridging and looping is a common mechanism for various and essential processes in bacterial chromosomes. This mechanism is preserved despite the very different bacterial conditions and their expected influence on the thermodynamic and kinetic characteristics of the bridge formation and stability. Over the last two decades, single-molecule techniques carried out on in vitro DNA systems have yielded valuable results which, in combination with theoretical works, have clarified the effects of different parameters of nucleoprotein complexes on the protein-induced DNA bridging and looping process. In this review, I will outline the features that can be measured for such processes with various single-molecule techniques in use in the field. I will then describe both the experimental results and the theoretical models that illuminate the contribution of the DNA molecule itself as well as that of the bridging proteins in the DNA looping mechanism at play in the nucleoid of E. coli.

The MalR type regulator AcrC is a transcriptional repressor of acarbose biosynthetic genes in Actinoplanes sp. SE50/110

Background

Acarbose is used in the treatment of diabetes mellitus type II and is produced by Actinoplanes sp. SE50/110. Although the biosynthesis of acarbose has been intensively studied, profound knowledge about transcription factors involved in acarbose biosynthesis and their binding sites has been missing until now. In contrast to acarbose biosynthetic gene clusters in Streptomyces spp., the corresponding gene cluster of Actinoplanes sp. SE50/110 lacks genes for transcriptional regulators.

Results

The acarbose regulator C (AcrC) was identified through an in silico approach by aligning the LacI family regulators of acarbose biosynthetic gene clusters in Streptomyces spp. with the Actinoplanes sp. SE50/110 genome. The gene for acrC, located in a head-to-head arrangement with the maltose/maltodextrin ABC transporter malEFG operon, was deleted by introducing PCR targeting for Actinoplanes sp. SE50/110. Characterization was carried out through cultivation experiments, genome-wide microarray hybridizations, and RT-qPCR as well as electrophoretic mobility shift assays for the elucidation of binding motifs. The results show that AcrC binds to the intergenic region between acbE and acbD in Actinoplanes sp. SE50/110 and acts as a transcriptional repressor on these genes. The transcriptomic profile of the wild type was reconstituted through a complementation of the deleted acrC gene. Additionally, regulatory sequence motifs for the binding of AcrC were identified in the intergenic region of acbE and acbD. It was shown that AcrC expression influences acarbose formation in the early growth phase. Interestingly, AcrC does not regulate the malEFG operon.

Conclusions

This study characterizes the first known transcription factor of the acarbose biosynthetic gene cluster in Actinoplanes sp. SE50/110. It therefore represents an important step for understanding the regulatory network of this organism. Based on this work, rational strain design for improving the biotechnological production of acarbose can now be implemented.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3941-x) contains supplementary material, which is available to authorized users.

Effect of transcription factor resource sharing on gene expression noise

Gene expression is intrinsically a stochastic (noisy) process with important implications for cellular functions. Deciphering the underlying mechanisms of gene expression noise remains one of the key challenges of regulatory biology. Theoretical models of transcription often incorporate the kinetics of how transcription factors (TFs) interact with a single promoter to impact gene expression noise. However, inside single cells multiple identical gene copies as well as additional binding sites can compete for a limiting pool of TFs. Here we develop a simple kinetic model of transcription, which explicitly incorporates this interplay between TF copy number and its binding sites. We show that TF sharing enhances noise in mRNA distribution across an isogenic population of cells. Moreover, when a single gene copy shares it’s TFs with multiple competitor sites, the mRNA variance as a function of the mean remains unaltered by their presence. Hence, all the data for variance as a function of mean expression collapse onto a single master curve independent of the strength and number of competitor sites. However, this result does not hold true when the competition stems from multiple copies of the same gene. Therefore, although previous studies showed that the mean expression follows a universal master curve, our findings suggest that different scenarios of competition bear distinct signatures at the level of variance. Intriguingly, the introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution. These results demonstrate the impact of limited availability of TF resource on the regulation of noise in gene expression.

Genetically identical cells, even when they are exposed to the same environmental conditions, display incredible diversity. Gene expression noise is attributed to be a key source of this phenotypic diversity. Transcriptional dynamics is a dominant source of expression noise. Although scores of theoretical and experimental studies have explored how noise is regulated at the level of transcription, most of them focus on the gene specific, cis regulatory elements, such as the number of transcription factor (TF) binding sites, their binding strength, etc. However, how the global properties of transcription, such as the limited availability of TFs impact noise in gene expression remains rather elusive. Here we build a theoretical model that incorporates the effect of limiting TF pool on gene expression noise. We find that competition between genes for TFs leads to enhanced variability in mRNA copy number across an isogenic population. Moreover, for gene copies sharing TFs with other competitor sites, mRNA variance as a function of the mean shows distinct imprints for one gene copy and multiple gene copies respectively. This stands in sharp contrast to the universal behavior found in mean expression irrespective of the different scenarios of competition. An interesting feature of competition is that introduction of competitor sites can transform a unimodal mRNA distribution into a multimodal distribution, which could lead to phenotypic variability.

Proteins mediating DNA loops effectively block transcription

Loops are ubiquitous topological elements formed when proteins simultaneously bind to two noncontiguous DNA sites. While a loop‐mediating protein may regulate initiation at a promoter, the presence of the protein at the other site may be an obstacle for RNA polymerases (RNAP) transcribing a different gene. To test whether a DNA loop alters the extent to which a protein blocks transcription, the lac repressor (LacI) was used. The outcome of in vitro transcription along templates containing two LacI operators separated by 400 bp in the presence of LacI concentrations that produced both looped and unlooped molecules was visualized with scanning force microscopy (SFM). An analysis of transcription elongation complexes, moving for 60 s at an average of 10 nt/s on unlooped DNA templates, revealed that they more often surpassed LacI bound to the lower affinity O2 operator than to the highest affinity Os operator. However, this difference was abrogated in looped DNA molecules where LacI became a strong roadblock independently of the affinity of the operator. Recordings of transcription elongation complexes, using magnetic tweezers, confirmed that they halted for several minutes upon encountering a LacI bound to a single operator. The average pause lifetime is compatible with RNAP waiting for LacI dissociation, however, the LacI open conformation visualized in the SFM images also suggests that LacI could straddle RNAP to let it pass. Independently of the mechanism by which RNAP bypasses the LacI roadblock, the data indicate that an obstacle with looped topology more effectively interferes with transcription.

Fourteen Ways to Reroute Cooperative Communication in the Lactose Repressor: Engineering Regulatory Proteins with Alternate Repressive Functions

The lactose repressor (LacI) is a classic genetic switch that has been used as a fundamental component in a host of synthetic genetic networks. To expand the function of LacI for use in the development of novel networks and other biotechnological applications, we engineered alternate communication in the LacI scaffold via laboratory evolution. Here we produced 14 new regulatory elements based on the LacI topology that are responsive to isopropyl β-d-1-thiogalactopyranoside (IPTG) with variation in repression strengths and ligand sensitivities-on solid media. The new variants exhibit repressive as well as antilac (i.e., inverse-repression + IPTG) functions and variations in the control of gene output upon exposure to different concentrations of IPTG. In addition, examination of this collection of variants in solution results in the controlled output of a canonical florescent reporter, demonstrating the utility of this collection of new regulatory proteins under standard conditions.

Macromolecular Crowding Regulates the Gene Expression Profile by Limiting Diffusion

We seek to elucidate the role of macromolecular crowding in transcription and translation. It is well known that stochasticity in gene expression can lead to differential gene expression and heterogeneity in a cell population. Recent experimental observations by Tan et al. have improved our understanding of the functional role of macromolecular crowding. It can be inferred from their observations that macromolecular crowding can lead to robustness in gene expression, resulting in a more homogeneous cell population. We introduce a spatial stochastic model to provide insight into this process. Our results show that macromolecular crowding reduces noise (as measured by the kurtosis of the mRNA distribution) in a cell population by limiting the diffusion of transcription factors (i.e. removing the unstable intermediate states), and that crowding by large molecules reduces noise more efficiently than crowding by small molecules. Finally, our simulation results provide evidence that the local variation in chromatin density as well as the total volume exclusion of the chromatin in the nucleus can induce a homogenous cell population.

The cellular nucleus is packed with macromolecules such as DNAs and proteins, which leaves limited space for other molecules to move around. Recent experimental results by C. Tan et al. have shown that macromolecular crowding can regulate gene expression, resulting in a more homogenous cell population. We introduce a computational model to uncover the mechanism by which macromolecular crowding functions. Our results suggest that macromolecular crowding limits the diffusion of the transcription factors and attenuates the transcriptional bursting, which leads to a more homogenous cell population. Regulation of gene expression noise by macromolecules depends on the size of the crowders, i.e. larger macromolecules can reduce the noise more effectively than smaller macromolecules. We also demonstrate that local variation of chromatin density can affect the noise of gene expression. This shows the importance of the chromatin structure in gene expression regulation.

Supercoiling biases the formation of loops involved in gene regulation

The function of DNA as a repository of genetic information is well-known. The post-genomic effort is to understand how this information-containing filament is chaperoned to manage its compaction and topological states. Indeed, the activities of enzymes that transcribe, replicate, or repair DNA are regulated to a large degree by access. Proteins that act at a distance along the filament by binding at one site and contacting another site, perhaps as part of a bigger complex, create loops that constitute topological domains and influence regulation. DNA loops and plectonemes are not necessarily spontaneous, especially large loops under tension for which high energy is required to bring their ends together, or small loops that require accessory proteins to facilitate DNA bending. However, the torsion in stiff filaments such as DNA dramatically modulates the topology, driving it from extended and genetically accessible to more looped and compact, genetically secured forms. Furthermore, there are accessory factors that bias the response of the DNA filament to supercoiling. For example, small molecules like polyamines, which neutralize the negative charge repulsions along the phosphate backbone, enhance flexibility and promote writhe over twist in response to torsion. Such increased flexibility likely pushes the topological equilibrium from twist toward writhe at tensions thought to exist in vivo. A predictable corollary is that stiffening DNA antagonizes looping and bending. Certain sequences are known to be more or less flexible or to exhibit curvature, and this may affect interactions with binding proteins. In vivo all of these factors operate simultaneously on DNA that is generally negatively supercoiled to some degree. Therefore, in order to better understand gene regulation that involves protein-mediated DNA loops, it is critical to understand the thermodynamics and kinetics of looping in DNA that is under tension, negatively supercoiled, and perhaps exposed to molecules that alter elasticity. Recent experiments quantitatively reveal how much negatively supercoiling DNA lowers the free energy of looping, possibly biasing the operation of genetic switches.

Force distribution in a semiflexible loop

Loops undergoing thermal fluctuations are prevalent in nature. Ringlike or cross-linked polymers, cyclic macromolecules, and protein-mediated DNA loops all belong to this category. Stability of these molecules are generally described in terms of free energy, an average quantity, but it may also be impacted by local fluctuating forces acting within these systems. The full distribution of these forces can thus give us insights into mechanochemistry beyond the predictive capability of thermodynamics. In this paper, we study the force exerted by an inextensible semiflexible polymer constrained in a looped state. By using a simulation method termed "phase-space sampling," we generate the equilibrium distribution of chain conformations in both position and momentum space. We compute the constraint forces between the two ends of the loop in this chain ensemble using Lagrangian mechanics, and show that the mean of these forces is equal to the thermodynamic force. By analyzing kinetic and potential contributions to the forces, we find that the mean force acts in the direction of increasing extension not because of bending stress, but in spite of it. Furthermore, we obtain a distribution of constraint forces as a function of chain length, extension, and stiffness. Notably, increasing contour length decreases the average force, but the additional freedom allows fluctuations in the constraint force to increase. The force distribution is asymmetric and falls off less sharply than a Gaussian distribution. Our work exemplifies a system where large-amplitude fluctuations occur in a way unforeseen by a purely thermodynamic framework, and offers computational tools useful for efficient, unbiased simulation of a constrained system.

Sequence Affects the Cyclization of DNA Minicircles

Understanding how the sequence of a DNA molecule affects its dynamic properties is a central problem affecting biochemistry and biotechnology. The process of cyclizing short DNA, as a critical step in molecular cloning, lacks a comprehensive picture of the kinetic process containing sequence information. We have elucidated this process by using coarse-grained simulations, enhanced sampling methods, and recent theoretical advances. We are able to identify the types and positions of structural defects during the looping process at a base-pair level. Correlations along a DNA molecule dictate critical sequence positions that can affect the looping rate. Structural defects change the bending elasticity of the DNA molecule from a harmonic to subharmonic potential with respect to bending angles. We explore the subelastic chain as a possible model in loop formation kinetics. A sequence-dependent model is developed to qualitatively predict the relative loop formation time as a function of DNA sequence.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.

A highly sensitive, direct and label-free technique for Hg(2+) detection using Kelvin probe force microscopy

For several decades, various nanomaterials have been used in a wide range of industrial fields, research areas, and commercial products. Among many nanomaterials, nano-sized mercury materials are one of the most widely used nanomaterials in real life. However, due to the high toxicity of Hg(2+), it is imperative to develop an effective and practical detection method for Hg(2+) to protect human health and environment. In this study, a highly sensitive, label-free method of detecting Hg(2+) that requires only a single drop of solution was developed. The detection mechanism is based on the different surface potential arising from Hg(2+) binding to mismatched thymine-thymine sequences, creating a very stable base pair. The surface potential is measured with Kelvin probe force microscopy (KPFM) to a molecular resolution. The developed method is capable of detecting 2 fmol of Hg(2+), which is 500 times more sensitive than previously reported techniques. Moreover, our method can selectively detect Hg(2+) and can also be applied to tap water and river water. This KPFM-based Hg(2+) detection method can be used as an early detection technique for practical applications.

Dynamic Release of Bending Stress in Short dsDNA by Formation of a Kink and Forks

Bending with high curvature is one of the major mechanical properties of double-stranded DNA (dsDNA) that is essential for its biological functions. The emergence of a kink arising from local melting in the middle of dsDNA has been suggested as a mechanism of releasing the energy cost of bending. Herein, we report that strong bending induces two types of short dsDNA deformations, induced by two types of local melting, namely, a kink in the middle and forks at the ends, which we demonstrate using D-shaped DNA nanostructures. The two types of deformed dsDNA structures dynamically interconvert on a millisecond timescale. The transition from a fork to a kink is dominated by entropic contribution (anti-Arrhenius behavior), while the transition from a kink to a fork is dominated by enthalpic contributions. The presence of mismatches in dsDNA accelerates kink formation, and the transition from a kink to a fork is removed when the mismatch size is three base pairs.

AFM of self-assembled lambda DNA-histone networks

Atomic force microscopy (AFM) was used to investigate the self-assembly behavior of λ-DNA and histones at varying histone:DNA ratios. Without histones and at the lowest histone:DNA ratio (less than one histone per 1000 base pairs of DNA), the DNA appeared as individual (uncomplexed), double-stranded DNA molecules. At increasing histone concentrations (one histone per 500, 250 and 167 base pairs of DNA), the DNA molecules started to form extensive polygonal networks of mostly pentagons and hexagons. The observed networks might be one of the naturally occurring, stable DNA-histone structures. The condensing effects of the divalent cations Mg(2+) and Ca(2+) on the DNA-histone complexes were also investigated. The networks persisted at high Mg(2+) concentration (20mM) and the highest histone concentration. At high Ca(2+) concentration and the highest histone concentration, the polygonal network disappeared and, instead, individual, tightly condensed aggregates were formed.

Species differential regulation of COX2 can be described by an NFκB-dependent logic AND gate

Cyclooxygenase 2 (COX2), a key regulatory enzyme of the prostaglandin/eicosanoid pathway, is an important target for anti-inflammatory therapy. It is highly induced by pro-inflammatory cytokines in a Nuclear factor kappa B (NFκB)-dependent manner. However, the mechanisms determining the amplitude and dynamics of this important pro-inflammatory event are poorly understood. Furthermore, there is significant difference between human and mouse COX2 expression in response to the inflammatory stimulus tumor necrosis factor alpha (TNFα). Here, we report the presence of a molecular logic AND gate composed of two NFκB response elements (NREs) which controls the expression of human COX2 in a switch-like manner. Combining quantitative kinetic modeling and thermostatistical analysis followed by experimental validation in iterative cycles, we show that the human COX2 expression machinery regulated by NFκB displays features of a logic AND gate. We propose that this provides a digital, noise-filtering mechanism for a tighter control of expression in response to TNFα, such that a threshold level of NFκB activation is required before the promoter becomes active and initiates transcription. This NFκB-regulated AND gate is absent in the mouse COX2 promoter, most likely contributing to its differential graded response in promoter activity and protein expression to TNFα. Our data suggest that the NFκB-regulated AND gate acts as a novel mechanism for controlling the expression of human COX2 to TNFα, and its absence in the mouse COX2 provides the foundation for further studies on understanding species-specific differential gene regulation.

Promoter architecture dictates cell-to-cell variability in gene expression

Variability in gene expression among genetically identical cells has emerged as a central preoccupation in the study of gene regulation; however, a divide exists between the predictions of molecular models of prokaryotic transcriptional regulation and genome-wide experimental studies suggesting that this variability is indifferent to the underlying regulatory architecture. We constructed a set of promoters in Escherichia coli in which promoter strength, transcription factor binding strength, and transcription factor copy numbers are systematically varied, and used messenger RNA (mRNA) fluorescence in situ hybridization to observe how these changes affected variability in gene expression. Our parameter-free models predicted the observed variability; hence, the molecular details of transcription dictate variability in mRNA expression, and transcriptional noise is specifically tunable and thus represents an evolutionarily accessible phenotypic parameter.

Biochemical characterization of RecA variants that contribute to extreme resistance to ionizing radiation

Among strains of Escherichia coli that have evolved to survive extreme exposure to ionizing radiation, mutations in the recA gene are prominent and contribute substantially to the acquired phenotype. Changes at amino acid residue 276, D276A and D276N, occur repeatedly and in separate evolved populations. RecA D276A and RecA D276N exhibit unique adaptations to an environment that can require the repair of hundreds of double strand breaks. These two RecA protein variants (a) exhibit a faster rate of filament nucleation on DNA, as well as a slower extension under at least some conditions, leading potentially to a distribution of the protein among a higher number of shorter filaments, (b) promote DNA strand exchange more efficiently in the context of a shorter filament, and (c) are markedly less inhibited by ADP. These adaptations potentially allow RecA protein to address larger numbers of double strand DNA breaks in an environment where ADP concentrations are higher due to a compromised cellular metabolism.

Tethered particle motion reveals that LacI·DNA loops coexist with a competitor-resistant but apparently unlooped conformation

The lac repressor protein (LacI) efficiently represses transcription of the lac operon in Escherichia coli by binding to two distant operator sites on the bacterial DNA and causing the intervening DNA to form a loop. We employed single-molecule tethered particle motion to observe LacI-mediated loop formation and breakdown in DNA constructs that incorporate optimized operator binding sites and intrinsic curvature favorable to loop formation. Previous bulk competition assays indirectly measured the loop lifetimes in these optimized DNA constructs as being on the order of days; however, we measured these same lifetimes to be on the order of minutes for both looped and unlooped states. In a range of single-molecule DNA competition experiments, we found that the resistance of the LacI-DNA complex to competitive binding is a function of both the operator strength and the interoperator sequence. To explain these findings, we present what we believe to be a new kinetic model of loop formation and DNA competition. In this proposed new model, we hypothesize a new unlooped state in which the unbound DNA-binding domain of the LacI protein interacts nonspecifically with nonoperator DNA adjacent to the operator site at which the second LacI DNA-binding domain is bound.

Arabidopsis SEPALLATA proteins differ in cooperative DNA-binding during the formation of floral quartet-like complexes

The SEPALLATA (SEP) genes of Arabidopsis thaliana encode MADS-domain transcription factors that specify the identity of all floral organs. The four Arabidopsis SEP genes function in a largely yet not completely redundant manner. Here, we analysed interactions of the SEP proteins with DNA. All of the proteins were capable of forming tetrameric quartet-like complexes on DNA fragments carrying two sequence elements termed CArG-boxes. Distances between the CArG-boxes for strong cooperative DNA-binding were in the range of 4-6 helical turns. However, SEP1 also bound strongly to CArG-box pairs separated by smaller or larger distances, whereas SEP2 preferred large and SEP4 preferred small inter-site distances for binding. Cooperative binding of SEP3 was comparatively weak for most of the inter-site distances tested. All SEP proteins constituted floral quartet-like complexes together with the floral homeotic proteins APETALA3 (AP3) and PISTILLATA (PI) on the target genes AP3 and SEP3. Our results suggest an important part of an explanation for why the different SEP proteins have largely, but not completely redundant functions in determining floral organ identity: they may bind to largely overlapping, but not identical sets of target genes that differ in the arrangement and spacing of the CArG-boxes in their cis-regulatory regions.

Multiple LacI-mediated loops revealed by Bayesian statistics and tethered particle motion

The bacterial transcription factor LacI loops DNA by binding to two separate locations on the DNA simultaneously. Despite being one of the best-studied model systems for transcriptional regulation, the number and conformations of loop structures accessible to LacI remain unclear, though the importance of multiple coexisting loops has been implicated in interactions between LacI and other cellular regulators of gene expression. To probe this issue, we have developed a new analysis method for tethered particle motion, a versatile and commonly used in vitro single-molecule technique. Our method, vbTPM, performs variational Bayesian inference in hidden Markov models. It learns the number of distinct states (i.e. DNA–protein conformations) directly from tethered particle motion data with better resolution than existing methods, while easily correcting for common experimental artifacts. Studying short (roughly 100 bp) LacI-mediated loops, we provide evidence for three distinct loop structures, more than previously reported in single-molecule studies. Moreover, our results confirm that changes in LacI conformation and DNA-binding topology both contribute to the repertoire of LacI-mediated loops formed in vitro, and provide qualitatively new input for models of looping and transcriptional regulation. We expect vbTPM to be broadly useful for probing complex protein–nucleic acid interactions.

High-throughput single-molecule studies of protein-DNA interactions

Fluorescence and force-based single-molecule studies of protein-nucleic acid interactions continue to shed critical insights into many aspects of DNA and RNA processing. As single-molecule assays are inherently low-throughput, obtaining statistically relevant datasets remains a major challenge. Additionally, most fluorescence-based single-molecule particle-tracking assays are limited to observing fluorescent proteins that are in the low-nanomolar range, as spurious background signals predominate at higher fluorophore concentrations. These technical limitations have traditionally limited the types of questions that could be addressed via single-molecule methods. In this review, we describe new approaches for high-throughput and high-concentration single-molecule biochemical studies. We conclude with a discussion of outstanding challenges for the single-molecule biologist and how these challenges can be tackled to further approach the biochemical complexity of the cell.

Poly(dA:dT)-rich DNAs are highly flexible in the context of DNA looping

Large-scale DNA deformation is ubiquitous in transcriptional regulation in prokaryotes and eukaryotes alike. Though much is known about how transcription factors and constellations of binding sites dictate where and how gene regulation will occur, less is known about the role played by the intervening DNA. In this work we explore the effect of sequence flexibility on transcription factor-mediated DNA looping, by drawing on sequences identified in nucleosome formation and ligase-mediated cyclization assays as being especially favorable for or resistant to large deformations. We examine a poly(dA:dT)-rich, nucleosome-repelling sequence that is often thought to belong to a class of highly inflexible DNAs; two strong nucleosome positioning sequences that share a set of particular sequence features common to nucleosome-preferring DNAs; and a CG-rich sequence representative of high G+C-content genomic regions that correlate with high nucleosome occupancy in vivo. To measure the flexibility of these sequences in the context of DNA looping, we combine the in vitro single-molecule tethered particle motion assay, a canonical looping protein, and a statistical mechanical model that allows us to quantitatively relate the looping probability to the looping free energy. We show that, in contrast to the case of nucleosome occupancy, G+C content does not positively correlate with looping probability, and that despite sharing sequence features that are thought to determine nucleosome affinity, the two strong nucleosome positioning sequences behave markedly dissimilarly in the context of looping. Most surprisingly, the poly(dA:dT)-rich DNA that is often characterized as highly inflexible in fact exhibits one of the highest propensities for looping that we have measured. These results argue for a need to revisit our understanding of the mechanical properties of DNA in a way that will provide a basis for understanding DNA deformation over the entire range of biologically relevant scenarios that are impacted by DNA deformability.

Modern biophysical approaches probe transcription-factor-induced DNA bending and looping

The genetic information of every living organism is stored in its genomic DNA that is perceived as a chemically stable and robust macromolecule. But at the same time, to fulfil its functions properly, it also needs to be highly dynamic and flexible. This includes partial melting of the double helix or compaction and bending of the DNA often brought about by protein factors that are able to interact with DNA stretches in a specific and non-specific manner. The conformational changes in the DNA need to be understood in order to describe biological systems in detail. As these events play out on the nanometre scale, new biophysical approaches have been employed to monitor conformational changes in this regime at the single-molecule level. Focusing on transcription factor action on promoter DNA, we discuss how current biophysical techniques are able to quantitatively describe this molecular process.

Optical Methods to Study Protein-DNA Interactions in Vitro and in Living Cells at the Single-Molecule Level

The maintenance of intact genetic information, as well as the deployment of transcription for specific sets of genes, critically rely on a family of proteins interacting with DNA and recognizing specific sequences or features. The mechanisms by which these proteins search for target DNA are the subject of intense investigations employing a variety of methods in biology. A large interest in these processes stems from the faster-than-diffusion association rates, explained in current models by a combination of 3D and 1D diffusion. Here, we present a review of the single-molecule approaches at the forefront of the study of protein-DNA interaction dynamics and target search in vitro and in vivo. Flow stretch, optical and magnetic manipulation, single fluorophore detection and localization as well as combinations of different methods are described and the results obtained with these techniques are discussed in the framework of the current facilitated diffusion model.

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.

LacI-DNA-IPTG loops: equilibria among conformations by single-molecule FRET

The E. coli Lac repressor (LacI) tetramer binds simultaneously to a promoter-proximal DNA binding site (operator) and an auxiliary operator, resulting in a DNA loop, which increases repression efficiency. Induction of the lac operon by allolactose reduces the affinity of LacI for DNA, but induction does not completely prevent looping in vivo. Our previous work on the conformations of LacI loops used a hyperstable model DNA construct, 9C14, that contains a sequence directed bend flanked by operators. Single-molecule fluorescence resonance energy transfer (SM-FRET) on a dual fluorophore-labeled LacI-9C14 loop showed that it adopts a single, stable, high-FRET V-shaped LacI conformation. Ligand-induced changes in loop geometry can affect loop stability, and the current work assesses loop population distributions for LacI-9C14 complexes containing the synthetic inducer IPTG. SM-FRET confirms that the high-FRET LacI-9C14 loop is only partially destabilized by saturating IPTG. LacI titration experiments and FRET fluctuation analysis suggest that the addition of IPTG induces loop conformational dynamics and re-equilibration between loop population distributions that include a mixture of looped states that do not exhibit high-efficiency FRET. The results show that repression by looping even at saturating IPTG should be considered in models for regulation of the operon. We propose that persistent DNA loops near the operator function biologically to accelerate rerepression upon exhaustion of inducer.

Extreme bendability of DNA less than 100 base pairs long revealed by single-molecule cyclization

The classical view of DNA posits that DNA must be stiff below the persistence length [<150 base pairs (bp)], but recent studies addressing this have yielded contradictory results. We developed a fluorescence-based, protein-free assay for studying the cyclization of single DNA molecules in real time. The assay samples the equilibrium population of a sharply bent, transient species that is entirely suppressed in single-molecule mechanical measurements and is biologically more relevant than the annealed species sampled in the traditional ligase-based assay. The looping rate has a weak length dependence between 67 and 106 bp that cannot be described by the worm-like chain model. Many biologically important protein-DNA interactions that involve looping and bending of DNA below 100 bp likely use this intrinsic bendability of DNA.

Operator sequence alters gene expression independently of transcription factor occupancy in bacteria

A canonical quantitative view of transcriptional regulation holds that the only role of operator sequence is to set the probability of transcription factor binding, with operator occupancy determining the level of gene expression. In this work, we test this idea by characterizing repression in vivo and the binding of RNA polymerase in vitro in experiments where operators of various sequences were placed either upstream or downstream from the promoter in Escherichia coli. Surprisingly, we find that operators with a weaker binding affinity can yield higher repression levels than stronger operators. Repressor bound to upstream operators modulates promoter escape, and the magnitude of this modulation is not correlated with the repressor-operator binding affinity. This suggests that operator sequences may modulate transcription by altering the nature of the interaction of the bound transcription factor with the transcriptional machinery, implying a new layer of sequence dependence that must be confronted in the quantitative understanding of gene expression.

Sequence dependence of transcription factor-mediated DNA looping

DNA is subject to large deformations in a wide range of biological processes. Two key examples illustrate how such deformations influence the readout of the genetic information: the sequestering of eukaryotic genes by nucleosomes and DNA looping in transcriptional regulation in both prokaryotes and eukaryotes. These kinds of regulatory problems are now becoming amenable to systematic quantitative dissection with a powerful dialogue between theory and experiment. Here, we use a single-molecule experiment in conjunction with a statistical mechanical model to test quantitative predictions for the behavior of DNA looping at short length scales and to determine how DNA sequence affects looping at these lengths. We calculate and measure how such looping depends upon four key biological parameters: the strength of the transcription factor binding sites, the concentration of the transcription factor, and the length and sequence of the DNA loop. Our studies lead to the surprising insight that sequences that are thought to be especially favorable for nucleosome formation because of high flexibility lead to no systematically detectable effect of sequence on looping, and begin to provide a picture of the distinctions between the short length scale mechanics of nucleosome formation and looping.

The physics of protein-DNA interaction networks in the control of gene expression

Protein-DNA interaction networks play a central role in many fundamental cellular processes. In gene regulation, physical interactions and reactions among the molecular components together with the physical properties of DNA control how genes are turned on and off. A key player in all these processes is the inherent flexibility of DNA, which provides an avenue for long-range interactions between distal DNA elements through DNA looping. Such versatility enables multiple interactions and results in additional complexity that is remarkably difficult to address with traditional approaches. This topical review considers recent advances in statistical physics methods to study the assembly of protein-DNA complexes with loops, their effects in the control of gene expression, and their explicit application to the prototypical lac operon genetic system of the E. coli bacterium. In the last decade, it has been shown that the underlying physical properties of DNA looping can actively control transcriptional noise, cell-to-cell variability, and other properties of gene regulation, including the balance between robustness and sensitivity of the induction process. These physical properties are largely dependent on the free energy of DNA looping, which accounts for DNA bending and twisting effects. These new physical methods have also been used in reverse to uncover the actual in vivo free energy of looping double-stranded DNA in living cells, which was not possible with existing experimental techniques. The results obtained for DNA looping by the lac repressor inside the E. coli bacterium showed a more malleable DNA than expected as a result of the interplay of the simultaneous presence of two distinct conformations of looped DNA.

Towards a molecular view of transcriptional control

The accumulation of experimental data over recent years has fueled theoretical work on how transcription factors (TFs) search for and recognise their DNA target sites, how they interact with one another, or with other DNA-binding proteins, and how they cope with the compaction of DNA within bacterial nucleoids or within eukaryotic chromatin. Many models have been built to study the kinetic, thermodynamic and mechanistic aspects of these questions. In some cases they have resulted in a relatively clear consensus view, but a number of questions remain controversial. We present an overview of recent work, with an emphasis on models that provide, or can inspire, a better understanding of transcriptional control at a detailed molecular level.

Probing DNA topology using tethered particle motion

Transcription factors mediate the formation of nucleoprotein complexes that are critical for efficient regulation of epigenetic switches. In these complexes, DNA is frequently bent or looped by the protein; other times, strong interactions lead the DNA to fully wrap the regulatory protein(s). The equilibrium between the bending, looping, full and partial wrapping of DNA governs the level of transcriptional regulation and is tuned by biophysical parameters. Characterization of the structure, kinetics, and thermodynamics of formation of such nucleoprotein complexes is fundamental to the understanding of the molecular mechanisms that underlie the operation of the genetic switches controlled by them. Here, we describe in detail how to perform tethered particle motion experiments aimed at understanding how protein-DNA interactions influence the formation and breakdown of these regulatory complexes.

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.

Single-molecule recognition of biomolecular interaction via Kelvin probe force microscopy

We report the scanning probe microscope (SPM)-based single-molecule recognition of biomolecular interactions between protein kinase and small ligands (i.e., ATP and Imatinib). In general, it is difficult to sense and detect the small ligands bound to protein kinase (at single-molecule resolution) using a conventional atomic force microscope (AFM) due to the limited resolution of conventional AFM for detecting the miniscule changes in molecular size driven by ligand binding. In this study, we have demonstrated that Kelvin probe force microscopy (KPFM) is able to articulate the surface potential of biomolecules interacting with ligands (i.e., the protein kinase-ATP interactions and inhibition phenomena induced by antagonistic molecules) in a label-free manner. Furthermore, measured surface potentials for biomolecular interactions enable quantitative descriptions on the ability of protein kinase to interact with small ligands such as ATP or antagonistic molecules. Our study sheds light on KPFM that allows the precise recognition of single-molecule interactions, which opens a new avenue for the design and development of novel molecular therapeutics.

Investigating Deinococcus radiodurans RecA protein filament formation on double-stranded DNA by a real-time single-molecule approach

With the aid of an efficient, precise, and almost error-free DNA repair system, Deinococcus radiodurans can survive hundreds of double-strand breaks inflicted by high doses of irradiation or desiccation. RecA of D. radiodurans (DrRecA) plays a central role both in the early phase of repair by an extended synthesis-dependent strand annealing process and in the later more general homologous recombination phase. Both roles likely require DrRecA filament formation on duplex DNA. We have developed single-molecule tethered particle motion experiments to study the assembly dynamics of RecA proteins on individual duplex DNA molecules by observing changes in DNA tether length resulting from RecA binding. We demonstrate that DrRecA nucleation on double-stranded DNA is much faster than that of Escherichia coli RecA protein (EcRecA), but the extension is slower. This combination of attributes would tend to increase the number and decrease the length of DrRecA filaments relative to those of EcRecA, a feature that may reflect the requirement to repair hundreds of genomic double-strand breaks concurrently in irradiated Deinococcus cells.

Effect of promoter architecture on the cell-to-cell variability in gene expression

According to recent experimental evidence, promoter architecture, defined by the number, strength and regulatory role of the operators that control transcription, plays a major role in determining the level of cell-to-cell variability in gene expression. These quantitative experiments call for a corresponding modeling effort that addresses the question of how changes in promoter architecture affect variability in gene expression in a systematic rather than case-by-case fashion. In this article we make such a systematic investigation, based on a microscopic model of gene regulation that incorporates stochastic effects. In particular, we show how operator strength and operator multiplicity affect this variability. We examine different modes of transcription factor binding to complex promoters (cooperative, independent, simultaneous) and how each of these affects the level of variability in transcriptional output from cell-to-cell. We propose that direct comparison between in vivo single-cell experiments and theoretical predictions for the moments of the probability distribution of mRNA number per cell can be used to test kinetic models of gene regulation. The emphasis of the discussion is on prokaryotic gene regulation, but our analysis can be extended to eukaryotic cells as well.