An integrated model of transcription factor diffusion shows the importance of intersegmental transfer and quaternary protein structure for target site finding

Design patterns of biological cells

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

A quantitative modelling approach for DNA repair on a population scale

The great advances of sequencing technologies allow the in vivo measurement of nuclear processes-such as DNA repair after UV exposure-over entire cell populations. However, data sets usually contain only a few samples over several hours, missing possibly important information in between time points. We developed a data-driven approach to analyse CPD repair kinetics over time in Saccharomyces cerevisiae. In contrast to other studies that consider sequencing signals as an average behaviour, we understand them as the superposition of signals from independent cells. By motivating repair as a stochastic process, we derive a minimal model for which the parameters can be conveniently estimated. We correlate repair parameters to a variety of genomic features that are assumed to influence repair, including transcription rate and nucleosome density. The clearest link was found for the transcription unit length, which has been unreported for budding yeast to our knowledge. The framework hence allows a comprehensive analysis of nuclear processes on a population scale.

Modulating binding affinity, specificity and configurations by multivalent interactions

Biological functions of proteins rely on their specific interactions with binding partners. Many proteins contain multiple domains, which can bind to their targets that often have more than one binding site, resulting in multivalent interactions. While it has been shown that multivalent interactions play an crucial role in modulating binding affinity and specificity, other potential effects of multivalent interactions are less explored. Here, we developed a broadly applicable transfer matrix formalism and used it to investigate the binding of two-domain ligands to targets with multiple binding sites. We show that 1) ligands with two specific binding domains can drastically boost both the binding affinity and specificity and down-shift the working concentration range, compared to single-domain ligands, 2) the presence of a positive domain-domain cooperativity or containing a non-specific binding domain can down-shift the working concentration range of ligands by increasing the binding affinity without compromising the binding specificity, 3) the configuration of the bound ligands has a strong concentration dependence, providing important insights into the physical origin of phase-separation processes taking place in living cells. In line with previous studies, our results suggest that multivalent interactions are utilized by cells for highly efficient regulation of target binding involved in a diverse range of cellular processes such as signal transduction, gene transcription, antibody-antigen recognition.

Emergent correlations in gene expression dynamics as footprints of resource competition

Genetic circuits need a cellular environment to operate in, which naturally couples the circuit function with the overall functionality of gene regulatory network. To execute their functions, all gene circuits draw resources in the form of RNA polymerases, ribosomes, and tRNAs. Recent experiments pointed out that the role of resource competition on synthetic circuit outputs could be immense. However, the effect of complexity of the circuit architecture on resource sharing dynamics is yet unexplored. In this paper, we employ mathematical modelling and in-silico experiments to identify the sources of resource trade-off and to quantify its impact on the function of a genetic circuit, keeping our focus on regulation of immediate downstream proteins, which are often used as protein read-outs. We show that estimating gene expression dynamics from readings of downstream protein data might be unreliable when the resource is limited and ribosome affinities are asymmetric. We focus on the impact of mRNA copy number and ribosome binding site (RBS) strength on the nonlinear isocline that emerges with two regimes, prominently separated by a tipping point, and study how correlation and competition dominate each other depending on various circuit parameters. Focusing further on genetic toggle circuit, we have identified major effects of resource competition in this model motif and quantified the observations. The observations are testable in wet-lab experiments, as all the parameters chosen are experimentally relevant.

FOXO4 Transactivation Domain Interaction with Forkhead DNA Binding Domain and Effect on Selective DNA Recognition for Transcription Initiation

Forkhead box O4 (FOXO4) is a human transcription factor (TF) that participates in cell homeostasis. While the structure and DNA binding properties of the conserved forkhead domain (FHD) have been thoroughly investigated, how the transactivation domain (TAD) regulates the DNA binding properties of the protein remains elusive. Here, we investigated the role of TAD in modulating the DNA binding properties of FOXO4 using solution NMR. We found that TAD and FHD form an intramolecular complex mainly governed by hydrophobic interaction. Remarkably, TAD and DNA share the same surface of FHD for binding. While FHD did not differentiate binding to target and non-target DNA, the FHD-TAD complex showed different behaviors depending on the DNA sequence. In the presence of TAD, free and DNA-bound FHD exhibited a slow exchange with target DNA and a fast exchange with non-target DNA. The interaction of the two domains affected the kinetic function of FHD depending on the type of DNA. Based on these findings, we suggest a transcription initiation model by which TAD modulates FOXO4 recognition of its target promoter DNA sequences. This study describes the function of TAD in FOXO4 and provides a new kinetic perspective on target sequence selection by TFs.

Discrete-state stochastic kinetic models for target DNA search by proteins: Theory and experimental applications

To perform their functions, transcription factors and DNA-repair/modifying enzymes randomly search DNA in order to locate their specific targets on DNA. Discrete-state stochastic kinetic models have been developed to explain how the efficiency of the search process is influenced by the molecular properties of proteins and DNA as well as by other factors such as molecular crowding. These theoretical models not only offer explanations on the relation of microscopic processes to macroscopic behavior of proteins, but also facilitate the analysis and interpretation of experimental data. In this review article, we provide an overview on discrete-state stochastic kinetic models and explain how these models can be applied to experimental investigations using stopped-flow, single-molecule, nuclear magnetic resonance (NMR), and other biophysical and biochemical methods.

Transient binding and jumping dynamics of p53 along DNA revealed by sub-millisecond resolved single-molecule fluorescence tracking

Characterization of the target search dynamics of DNA-binding proteins along DNA has been hampered by the time resolution of a standard single-molecule fluorescence microscopy. Here, we achieved the time resolution of 0.5 ms in the fluorescence microscopy measurements by optimizing the fluorescence excitation based on critical angle illumination and by utilizing the time delay integration mode of the electron-multiplying charge coupled device. We characterized the target search dynamics of the tumor suppressor p53 along nonspecific DNA at physiological salt concentrations. We identified a short-lived encounter intermediate before the formation of the long-lived p53–DNA complex. Both the jumps and the one-dimensional diffusion of p53 along DNA were accelerated at higher salt concentrations, suggesting the rotation-uncoupled movement of p53 along DNA grooves and conformational changes in the p53/DNA complex. This method can be used to clarify the unresolved dynamics of DNA-binding proteins previously hidden by time averaging.

Modelling and measuring intracellular competition for finite resources during gene expression

Dissecting the competition between genes for shared expressional resources is of fundamental importance for understanding the interplay between cellular components. Owing to the relationship between gene expression and cellular fitness, genomes are shaped by evolution to improve resource allocation. Whereas experimental approaches to investigate intracellular competition require technical resources and human expertise, computational models and in silico simulations allow vast numbers of experiments to be carried out and controlled easily, and with significantly reduced costs. Thus, modelling competition has a pivotal role in understanding the effects of competition on the biophysics of the cell. In this article, we review various computational models proposed to describe the different types of competition during gene expression. We also present relevant synthetic biology experiments and their biotechnological implications, and discuss the open questions in the field.

How p53 Molecules Solve the Target DNA Search Problem: A Review

Interactions between DNA and DNA-binding proteins play an important role in many essential cellular processes. A key function of the DNA-binding protein p53 is to search for and bind to target sites incorporated in genomic DNA, which triggers transcriptional regulation. How do p53 molecules achieve "rapid" and "accurate" target search in living cells? The search dynamics of p53 were expected to include 3D diffusion in solution, 1D diffusion along DNA, and intersegmental transfer between two different DNA strands. Single-molecule fluorescence microscopy enabled the tracking of p53 molecules on DNA and the characterization of these dynamics quantitatively. Recent intensive single-molecule studies of p53 succeeded in revealing each of these search dynamics. Here, we review these studies and discuss the target search mechanisms of p53.

A Model for the Spatiotemporal Design of Gene Regulatory Circuits †

Mathematical modeling assists the design of synthetic regulatory networks by providing a detailed mechanistic understanding of biological systems. Models that can predict the performance of a design are fundamental for synthetic biology since they minimize iterations along the design-build-test lifecycle. Such predictability depends crucially on what assumptions (i.e., biological simplifications) the model considers. Here, we challenge a common assumption when it comes to the modeling of bacterial-based gene regulation: considering negligible the effects of intracellular physical space. It is commonly assumed that molecules, such as transcription factors (TF), are homogeneously distributed inside a cell, so there is no need to model their diffusion. We describe a mathematical model that accounts for molecular diffusion and show how simulations of network performance are decisively affected by the distance between its components. Specifically, the model focuses on the search by a TF for its target promoter. The combination of local searches, via one-dimensional sliding along the chromosome, and global searches, via three-dimensional diffusion through the cytoplasm, determine TF-promoter interplay. Previous experimental results with engineered bacteria in which the distance between TF source and target was minimized or enlarged were successfully reproduced by the spatially resolved model we introduce here. This suggests that the spatial specification of the circuit alone can be exploited as a design parameter in synthetic biology to select programmable output levels.

Intrinsically disordered domain of tumor suppressor p53 facilitates target search by ultrafast transfer between different DNA strands

Intersegmental transfer (IST) is an important strategy in the target search used by sequence-specific DNA-binding proteins (DBPs), enabling DBPs to search for targets between multiple DNA strands without dissociation. We examined the IST of the tumor suppressor p53 using ensemble stopped-flow and single-molecule fluorescence measurements. The ensemble measurements demonstrated that p53 exhibits very fast IST, whose rate constant was ∼108 M-1 s-1. To determine the domains of p53 responsible for IST, two mutants with deletions of one of its two DNA binding domains were generated. The mutant lacking the disordered C-terminal (CT) domain (the CoreTet mutant) abolished IST, whereas the mutant lacking the structured core domain (the TetCT mutant) maintained IST, clearly demonstrating the importance of the CT domain. Single-molecule fluorescence measurements further demonstrated the transfer of p53 between two tethered DNA strands. The pseudo-wild-type p53 and the TetCT mutant showed significant transfer efficiencies, whereas the transfer efficiency for the CoreTet mutant was zero. These results suggest that ultrafast IST might be promoted by four copies of the CT domain, by binding to two DNA strands simultaneously. Such ultrafast IST might be important to avoid nearby-bound DBPs during the target search process of p53 in nucleus.

Towards mapping the 3D genome through high speed single-molecule tracking of functional transcription factors in single living cells

How genomic DNA is organized in the nucleus is a long-standing question. We describe a single-molecule bioimaging method utilizing super-localization precision coupled to fully quantitative image analysis tools, towards determining snapshots of parts of the 3D genome architecture of model eukaryote budding yeast Saccharomyces cerevisiae with exceptional millisecond time resolution. We employ astigmatism imaging to enable robust extraction of 3D position data on genomically encoded fluorescent protein reporters that bind to DNA. Our relatively straightforward method enables snippets of 3D architectures of likely single genome conformations to be resolved captured via DNA-sequence specific binding proteins in single functional living cells.

Transcription factors in eukaryotic cells can functionally regulate gene expression by acting in oligomeric assemblies formed from an intrinsically disordered protein phase transition enabled by molecular crowding

High-speed single-molecule fluorescence microscopy in vivo shows that transcription factors in eukaryotes can act in oligomeric clusters mediated by molecular crowding and intrinsically disordered protein. This finding impacts on the longstanding puzzle of how transcription factors find their gene targets so efficiently in the complex, heterogeneous environment of the cell.

Abbreviations CDF - cumulative distribution function; FRAP - fluorescence recovery after photobleaching; GFP - Green fluorescent protein; STORM - stochastic optical reconstruction microscopy; TF - Transcription factor; YFP - Yellow fluorescent protein

Smoldyn: particle-based simulation with rule-based modeling, improved molecular interaction and a library interface

Motivation

Smoldyn is a spatial and stochastic biochemical simulator. It treats each molecule of interest as an individual particle in continuous space, simulating molecular diffusion, molecule-membrane interactions and chemical reactions, all with good accuracy. This article presents several new features.

Results

Smoldyn now supports two types of rule-based modeling. These are a wildcard method, which is very convenient, and the BioNetGen package with extensions for spatial simulation, which is better for complicated models. Smoldyn also includes new algorithms for simulating the diffusion of surface-bound molecules and molecules with excluded volume. Both are exact in the limit of short time steps and reasonably good with longer steps. In addition, Smoldyn supports single-molecule tracking simulations. Finally, the Smoldyn source code can be accessed through a C/C ++ language library interface.

Availability and Implementation

Smoldyn software, documentation, code, and examples are at http://www.smoldyn.org .

Contact

steven.s.andrews@gmail.com.

Processive searching ability varies among members of the gap-filling DNA polymerase X family

DNA repair proteins must locate rare damaged sites within the genome. DNA polymerase β (Pol β), a member of the DNA polymerase X family that is involved in base excision repair, uses a processive hopping search mechanism to locate substrates. This effectively enhances its search footprint on DNA, increasing the probability of locating damaged sites. Processive searching has been reported or proposed for many DNA-binding proteins, raising the question of how widespread or specific to certain enzymes the ability to perform this function is. To provide insight into this question, we compared the ability of three homologous DNA Pol X family members to perform a processive search for 1-nucleotide gaps in DNA using a previously developed biochemical assay. We found that at near-predicted physiological ionic strengths, the intramolecular searching ability of Pol β is at least 4-fold higher than that of Pol μ and ∼2-fold higher than that of Pol λ. Pol β also was able to perform intersegmental transfer with the intersegmental searching ability of Pol β being at least 6- and ∼2-fold higher than that of Pols μ and λ, respectively. Mutational analysis suggested that differences in the N-terminal domains of these polymerases are responsible for the varying degrees of searching competence. Of note, the differences in processive searching ability observed among the DNA Pol X family members correlated with their proposed biological functions in base excision repair and nonhomologous end joining.

Transcription factor clusters regulate genes in eukaryotic cells

Transcription is regulated through binding factors to gene promoters to activate or repress expression, however, the mechanisms by which factors find targets remain unclear. Using single-molecule fluorescence microscopy, we determined in vivo stoichiometry and spatiotemporal dynamics of a GFP tagged repressor, Mig1, from a paradigm signaling pathway of Saccharomyces cerevisiae. We find the repressor operates in clusters, which upon extracellular signal detection, translocate from the cytoplasm, bind to nuclear targets and turnover. Simulations of Mig1 configuration within a 3D yeast genome model combined with a promoter-specific, fluorescent translation reporter confirmed clusters are the functional unit of gene regulation. In vitro and structural analysis on reconstituted Mig1 suggests that clusters are stabilized by depletion forces between intrinsically disordered sequences. We observed similar clusters of a co-regulatory activator from a different pathway, supporting a generalized cluster model for transcription factors that reduces promoter search times through intersegment transfer while stabilizing gene expression.

Tetrameric structure of the restriction DNA glycosylase R.PabI in complex with nonspecific double-stranded DNA

R.PabI is a type II restriction enzyme that recognizes the 5′-GTAC-3′ sequence and belongs to the HALFPIPE superfamily. Although most restriction enzymes cleave phosphodiester bonds at specific sites by hydrolysis, R.PabI flips the guanine and adenine bases of the recognition sequence out of the DNA helix and hydrolyzes the N-glycosidic bond of the flipped adenine in a similar manner to DNA glycosylases. In this study, we determined the structure of R.PabI in complex with double-stranded DNA without the R.PabI recognition sequence by X-ray crystallography. The 1.9 Å resolution structure of the complex showed that R.PabI forms a tetrameric structure to sandwich the double-stranded DNA and the tetrameric structure is stabilized by four salt bridges. DNA binding and DNA glycosylase assays of the R.PabI mutants showed that the residues that form the salt bridges (R70 and D71) are essential for R.PabI to find the recognition sequence from the sea of nonspecific sequences. R.PabI is predicted to utilize the tetrameric structure to bind nonspecific double-stranded DNA weakly and slide along it to find the recognition sequence.