Marcin Tabaka

Bimodal gene expression in noncooperative regulatory systems

Bimodality of gene expression, as a mechanism contributing to phenotypic diversity, enhances the survival of cells in a fluctuating environment. To date, the bimodal response of a gene regulatory system has been attributed to the cooperativity of transcription factor binding or to feedback loops. It has remained unclear whether noncooperative binding of transcription factors can give rise to bimodality in an open-loop system. We study a theoretical model of gene expression in a two-step cascade (a deterministically monostable system) in which the regulatory gene produces transcription factors that have a nonlinear effect on the activity of the target gene. We show that a unimodal distribution of transcription factors over the cell population can generate a bimodal steady-state output without cooperative transcription factor binding. We introduce a simple method of geometric construction that allows one to predict the onset of bimodality. The construction only involves the parameters of bursting of the regulatory gene and the dose–response curve of the target gene. Using this method, we show that the gene expression may switch between unimodal and bimodal as the concentration of inducers or corepressors is varied. These findings may explain the experimentally observed bimodal response of cascades consisting of a fluorescent protein reporter controlled by the tetracycline repressor. The geometric construction provides a useful tool for designing experiments and for interpretation of their results. Our findings may have important implications for understanding the strategies adopted by cell populations to survive in changing environments.

Biologistics--diffusion coefficients for complete proteome of Escherichia coli.

Motivation: Biologistics provides data for quantitative analysis of transport (diffusion) processes and their spatio-temporal correlations in cells. Mobility of proteins is one of the few parameters necessary to describe reaction rates for gene regulation. Although understanding of diffusion-limited biochemical reactions in vivo requires mobility data for the largest possible number of proteins in their native forms, currently, there is no database that would contain the complete information about the diffusion coefficients (DCs) of proteins in a given cell type. Results: We demonstrate a method for the determination of in vivo DCs for any molecule—regardless of its molecular weight, size and structure—in any type of cell. We exemplify the method with the database of in vivo DC for all proteins (4302 records) from the proteome of K12 strain of Escherichia coli, together with examples of DC of amino acids, sugars, RNA and DNA. The database follows from the scale-dependent viscosity reference curve (sdVRC). Construction of sdVRC for prokaryotic or eukaryotic cell requires ~20 in vivo measurements using techniques such as fluorescence correlation spectroscopy (FCS), fluorescence recovery after photobleaching (FRAP), nuclear magnetic resonance (NMR) or particle tracking. The shape of the sdVRC would be different for each organism, but the mathematical form of the curve remains the same. The presented method has a high predictive power, as the measurements of DCs of several inert, properly chosen probes in a single cell type allows to determine the DCs of thousands of proteins. Additionally, obtained mobility data allow quantitative study of biochemical interactions in vivo. Contact: rholyst@ichf.edu.pl Supplementary information: Supplementary data are available at Bioinformatics Online.

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 Mw = 12.91 kDa) and 1550 s (e.g. for PepA6 of Mw = 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.

Motion of nanoprobes in complex liquids within the framework of the length-scale dependent viscosity model.

This paper deals with the recent phenomenological model of the motion of nanoscopic objects (colloidal particles, proteins, nanoparticles, molecules) in complex liquids. We analysed motion in polymer, micellar, colloidal and protein solutions and the cytoplasm of living cells using the length-scale dependent viscosity model. Viscosity monotonically approaches macroscopic viscosity as the size of the object increases and thus gives a single, coherent picture of motion at the nano and macro scale. The model includes interparticle interactions (solvent-solute), temperature and the internal structure of a complex liquid. The depletion layer ubiquitously occurring in complex liquids is also incorporated into the model. We also discuss the biological aspects of crowding in terms of the length-scale dependent viscosity model.

Transcriptional leakage versus noise: A simple mechanism of conversion between binary and graded response in autoregulated genes

We study the response of an autoregulated gene to a range of concentrations of signal molecules. We show that transcriptional leakage and noise due to translational bursting have the opposite effects. In a positively autoregulated gene, increasing the noise converts the response from graded to binary, while increasing the leakage converts the response from binary to graded. Our findings support the hypothesis that, being a common phenomenon, leaky expression may be a relatively easy way for evolutionary tuning of the type of gene response without changing the type of regulation from positive to negative.

The effect of macromolecular crowding on mobility of biomolecules, association kinetics, and gene expression in living cells

We discuss a quantitative influence of macromolecular crowding on biological processes: motion, bimolecular reactions, and gene expression in prokaryotic and eukaryotic cells. We present scaling laws relating diffusion coefficient of an object moving in a cytoplasm of cells to a size of this object and degree of crowding. Such description leads to the notion of the length scale dependent viscosity characteristic for all living cells. We present an application of the length-scale dependent viscosity model to the description of motion in the cytoplasm of both eukaryotic and prokaryotic living cells. We compare the model with all recent data on diffusion of nanoscopic objects in HeLa, and E. coli cells. Additionally a description of the mobility of molecules in cell nucleus is presented. Finally we discuss the influence of crowding on the bimolecular association rates and gene expression in living cells.

Characterization of Caulobacter crescentus FtsZ Protein Using Dynamic Light Scattering

Background: Self-assembly of the tubulin-homologue FtsZ is critical in bacterial cell division. Results: Dynamic light scattering (DLS) measurements provide insight into the kinetics and stable length of Caulobacter crescentus FtsZ in vitro. Conclusion: C. crescentus FtsZ forms short linear polymers in solution with the assembly rate depending on the concentrations of GTP and GDP. Significance: DLS is a valuable technique for studying the polymerization of cytoskeletal proteins. The self-assembly of the tubulin homologue FtsZ at the mid-cell is a critical step in bacterial cell division. We introduce dynamic light scattering (DLS) spectroscopy as a new method to study the polymerization kinetics of FtsZ in solution. Analysis of the DLS data indicates that the FtsZ polymers are remarkably monodisperse in length, independent of the concentrations of GTP, GDP, and FtsZ monomers. Measurements of the diffusion coefficient of the polymers demonstrate that their length is remarkably stable until the free GTP is consumed. We estimated the mean size of the FtsZ polymers within this interval of stable length to be between 9 and 18 monomers. The rates of FtsZ polymerization and depolymerization are likely influenced by the concentration of GDP, as the repeated addition of GTP to FtsZ increased the rate of polymerization and slowed down depolymerization. Increasing the FtsZ concentration did not change the size of FtsZ polymers; however, it increased the rate of the depolymerization reaction by depleting free GTP. Using transmission electron microscopy we observed that FtsZ forms linear polymers in solutions which rapidly convert to large bundles upon contact with surfaces at time scales as short as several seconds. Finally, the best studied small molecule that binds to FtsZ, PC190723, had no stabilizing effect on Caulobacter crescentus FtsZ filaments in vitro, which complements previous studies with Escherichia coli FtsZ and confirms that this class of small molecules binds Gram-negative FtsZ weakly.

Taylor dispersion analysis in coiled capillaries at high flow rates.

Taylor Dispersion Analysis (TDA) has been performed for analytes moving at high flow rates in long, coiled capillaries. A thin injection zone of the analyte is stretched by the flow and final distribution of concentration of the analyte at the end of the capillary has the gaussian shape. The high flow rates in coiled capillary generate vortices. They convectively mix the analyte across the capillary. This mixing reduces the width of the gaussian distribution several times in comparison to the width obtained in a straight capillary in standard TDA. We have determined an empirical, scaling equation for the width as a function of the flow rate, molecular diffusion coefficient of the analyte, viscosity of the carrier phase, internal radius of the cylindrical capillary, and external radius of the coiled capillary. This equation can be used for different sizes of capillaries in a wide range of parameters without an additional calibration procedure. Our experimental results of flow in the coiled capillary could not be explained by current models based on approximate solutions of the Navier-Stokes equation. We applied the technique to determine the diffusion coefficients of the following analytes: salts, drugs, single amino acids, peptides (from dipeptides to hexapeptides), and proteins.

Small Crowders Slow Down Kinesin-1 Stepping by Hindering Motor Domain Diffusion.

The dimeric motor protein kinesin-1 moves processively along microtubules against forces of up to 7 pN. However, the mechanism of force generation is still debated. Here, we point to the crucial importance of diffusion of the tethered motor domain for the stepping of kinesin-1: small crowders stop the motor at a viscosity of 5 mPa·s-corresponding to a hydrodynamic load in the sub-fN (~10^{-4} pN) range-whereas large crowders have no impact even at viscosities above 100 mPa·s. This indicates that the scale-dependent, effective viscosity experienced by the tethered motor domain is a key factor determining kinesins functionality. Our results emphasize the role of diffusion in the kinesin-1 stepping mechanism and the general importance of the viscosity scaling paradigm in nanomechanics.

Implications of macromolecular crowding for protein–protein association kinetics in the cytoplasm of living cells

We present a method of prediction of in vivo association rate constants between proteins from in vitro rate constants. This model accounts for the changes in length-scale dependent viscosity of the cytoplasm of a living cell and interaction potential between proteins. We explain the unexpected results obtained for protein–protein association rate constants in the cytoplasm of HeLa cells determined by Phillip et al. Proc. Natl. Acad. Sci. U. S. A., 2012, 109, 1461.