Nunchaku: optimally partitioning data into piece-wise contiguous segments.

MOTIVATION: When analyzing 1D time series, scientists are often interested in identifying regions where one variable depends linearly on the other. Typically, they use an ad hoc and therefore often subjective method to do so. RESULTS: Here, we develop a statistically rigorous, Bayesian approach to infer the optimal partitioning of a dataset not only into contiguous piece-wise linear segments, but also into contiguous segments described by linear combinations of arbitrary basis functions. We therefore present a general solution to the problem of identifying discontinuous change points. Focusing on microbial growth, we use the algorithm to find the range of optical density where this density is linearly proportional to the number of cells and to automatically find the regions of exponential growth for both Escherichia coli and Saccharomyces cerevisiae. For budding yeast, we consequently are able to infer the Monod constant for growth on fructose. Our algorithm lends itself to automation and high throughput studies, increases reproducibility, and should facilitate data analyses for a broad range of scientists. AVAILABILITY AND IMPLEMENTATION: The corresponding Python package, entitled Nunchaku, is available at PyPI: https://pypi.org/project/nunchaku.

Multiple nutrient transporters enable cells to mitigate a rate-affinity tradeoff.

Eukaryotic genomes often encode multiple transporters for the same nutrient. For example, budding yeast has 17 hexose transporters (HXTs), all of which potentially transport glucose. Using mathematical modelling, we show that transporters that use either facilitated diffusion or symport can have a rate-affinity tradeoff, where an increase in the maximal rate of transport decreases the transporter's apparent affinity. These changes affect the import flux non-monotonically, and for a given concentration of extracellular nutrient there is one transporter, characterised by its affinity, that has a higher import flux than any other. Through encoding multiple transporters, cells can therefore mitigate the tradeoff by expressing those transporters with higher affinities in lower concentrations of nutrients. We verify our predictions using fluorescent tagging of seven HXT genes in budding yeast and follow their expression over time in batch culture. Using the known affinities of the corresponding transporters, we show that their regulation in glucose is broadly consistent with a rate-affinity tradeoff: as glucose falls, the levels of the different transporters peak in an order that mostly follows their affinity for glucose. More generally, evolution is constrained by tradeoffs. Our findings indicate that one such tradeoff often occurs in the cellular transport of nutrients.

Analysing and meta-analysing time-series data of microbial growth and gene expression from plate readers.

Responding to change is a fundamental property of life, making time-series data invaluable in biology. For microbes, plate readers are a popular, convenient means to measure growth and also gene expression using fluorescent reporters. Nevertheless, the difficulties of analysing the resulting data can be a bottleneck, particularly when combining measurements from different wells and plates. Here we present omniplate, a Python module that corrects and normalises plate-reader data, estimates growth rates and fluorescence per cell as functions of time, calculates errors, exports in different formats, and enables meta-analysis of multiple plates. The software corrects for autofluorescence, the optical density's non-linear dependence on the number of cells, and the effects of the media. We use omniplate to measure the Monod relationship for the growth of budding yeast in raffinose, showing that raffinose is a convenient carbon source for controlling growth rates. Using fluorescent tagging, we study yeast's glucose transport. Our results are consistent with the regulation of the hexose transporter (HXT) genes being approximately bipartite: the medium and high affinity transporters are predominately regulated by both the high affinity glucose sensor Snf3 and the kinase complex SNF1 via the repressors Mth1, Mig1, and Mig2; the low affinity transporters are predominately regulated by the low affinity sensor Rgt2 via the co-repressor Std1. We thus demonstrate that omniplate is a powerful tool for exploiting the advantages offered by time-series data in revealing biological regulation.

Roadmap on biology in time varying environments.

Biological organisms experience constantly changing environments, from sudden changes in physiology brought about by feeding, to the regular rising and setting of the Sun, to ecological changes over evolutionary timescales. Living organisms have evolved to thrive in this changing world but the general principles by which organisms shape and are shaped by time varying environments remain elusive. Our understanding is particularly poor in the intermediate regime with no separation of timescales, where the environment changes on the same timescale as the physiological or evolutionary response. Experiments to systematically characterize the response to dynamic environments are challenging since such environments are inherently high dimensional. This roadmap deals with the unique role played by time varying environments in biological phenomena across scales, from physiology to evolution, seeking to emphasize the commonalities and the challenges faced in this emerging area of research.

Distributed and dynamic intracellular organization of extracellular information.

Although cells respond specifically to environments, how environmental identity is encoded intracellularly is not understood. Here, we study this organization of information in budding yeast by estimating the mutual information between environmental transitions and the dynamics of nuclear translocation for 10 transcription factors. Our method of estimation is general, scalable, and based on decoding from single cells. The dynamics of the transcription factors are necessary to encode the highest amounts of extracellular information, and we show that information is transduced through two channels: Generalists (Msn2/4, Tod6 and Dot6, Maf1, and Sfp1) can encode the nature of multiple stresses, but only if stress is high; specialists (Hog1, Yap1, and Mig1/2) encode one particular stress, but do so more quickly and for a wider range of magnitudes. In particular, Dot6 encodes almost as much information as Msn2, the master regulator of the environmental stress response. Each transcription factor reports differently, and it is only their collective behavior that distinguishes between multiple environmental states. Changes in the dynamics of the localization of transcription factors thus constitute a precise, distributed internal representation of extracellular change. We predict that such multidimensional representations are common in cellular decision-making.

Morphologically constrained and data informed cell segmentation of budding yeast.

Motivation: Although high-content image cytometry is becoming increasingly routine, processing the large amount of data acquired during time-lapse experiments remains a challenge. The majority of approaches for automated single-cell segmentation focus on flat, uniform fields of view covered with a single layer of cells. In the increasingly popular microfluidic devices that trap individual cells for long term imaging, these conditions are not met. Consequently, most techniques for segmentation perform poorly. Although potentially constraining the generalizability of software, incorporating information about the microfluidic features, flow of media and the morphology of the cells can substantially improve performance. Results: Here we present DISCO (Data Informed Segmentation of Cell Objects), a framework for using the physical constraints imposed by microfluidic traps, the shape based morphological constraints of budding yeast and temporal information about cell growth and motion to allow tracking and segmentation of cells in microfluidic devices. Using manually curated datasets, we demonstrate substantial improvements in both tracking and segmentation when compared with existing software. Availability and implementation: The MATLAB code for the algorithm and for measuring performance is available at https://github.com/pswain/segmentation-software and the test images and the curated ground-truth results used for comparing the algorithms are available at http://datashare.is.ed.ac.uk/handle/10283/2002. Contact: mcrane2@uw.edu. Supplementary information: Supplementary data are available at Bioinformatics online.

Mechanistic links between cellular trade-offs, gene expression, and growth.

Intracellular processes rarely work in isolation but continually interact with the rest of the cell. In microbes, for example, we now know that gene expression across the whole genome typically changes with growth rate. The mechanisms driving such global regulation, however, are not well understood. Here we consider three trade-offs that, because of limitations in levels of cellular energy, free ribosomes, and proteins, are faced by all living cells and we construct a mechanistic model that comprises these trade-offs. Our model couples gene expression with growth rate and growth rate with a growing population of cells. We show that the model recovers Monod's law for the growth of microbes and two other empirical relationships connecting growth rate to the mass fraction of ribosomes. Further, we can explain growth-related effects in dosage compensation by paralogs and predict host-circuit interactions in synthetic biology. Simulating competitions between strains, we find that the regulation of metabolic pathways may have evolved not to match expression of enzymes to levels of extracellular substrates in changing environments but rather to balance a trade-off between exploiting one type of nutrient over another. Although coarse-grained, the trade-offs that the model embodies are fundamental, and, as such, our modeling framework has potentially wide application, including in both biotechnology and medicine.

BioJazz: in silico evolution of cellular networks with unbounded complexity using rule-based modeling.

Systems biologists aim to decipher the structure and dynamics of signaling and regulatory networks underpinning cellular responses; synthetic biologists can use this insight to alter existing networks or engineer de novo ones. Both tasks will benefit from an understanding of which structural and dynamic features of networks can emerge from evolutionary processes, through which intermediary steps these arise, and whether they embody general design principles. As natural evolution at the level of network dynamics is difficult to study, in silico evolution of network models can provide important insights. However, current tools used for in silico evolution of network dynamics are limited to ad hoc computer simulations and models. Here we introduce BioJazz, an extendable, user-friendly tool for simulating the evolution of dynamic biochemical networks. Unlike previous tools for in silico evolution, BioJazz allows for the evolution of cellular networks with unbounded complexity by combining rule-based modeling with an encoding of networks that is akin to a genome. We show that BioJazz can be used to implement biologically realistic selective pressures and allows exploration of the space of network architectures and dynamics that implement prescribed physiological functions. BioJazz is provided as an open-source tool to facilitate its further development and use. Source code and user manuals are available at: http://oss-lab.github.io/biojazz and http://osslab.lifesci.warwick.ac.uk/BioJazz.aspx.

The scaffold protein Ste5 directly controls a switch-like mating decision in yeast.

Evolution has resulted in numerous innovations that allow organisms to increase their fitness by choosing particular mating partners, including secondary sexual characteristics, behavioural patterns, chemical attractants and corresponding sensory mechanisms. The haploid yeast Saccharomyces cerevisiae selects mating partners by interpreting the concentration gradient of pheromone secreted by potential mates through a network of mitogen-activated protein kinase (MAPK) signalling proteins. The mating decision in yeast is an all-or-none, or switch-like, response that allows cells to filter weak pheromone signals, thus avoiding inappropriate commitment to mating by responding only at or above critical concentrations when a mate is sufficiently close. The molecular mechanisms that govern the switch-like mating decision are poorly understood. Here we show that the switching mechanism arises from competition between the MAPK Fus3 and a phosphatase Ptc1 for control of the phosphorylation state of four sites on the scaffold protein Ste5. This competition results in a switch-like dissociation of Fus3 from Ste5 that is necessary to generate the switch-like mating response. Thus, the decision to mate is made at an early stage in the pheromone pathway and occurs rapidly, perhaps to prevent the loss of the potential mate to competitors. We argue that the architecture of the Fus3-Ste5-Ptc1 circuit generates a novel ultrasensitivity mechanism, which is robust to variations in the concentrations of these proteins. This robustness helps assure that mating can occur despite stochastic or genetic variation between individuals. The role of Ste5 as a direct modulator of a cell-fate decision expands the functional repertoire of scaffold proteins beyond providing specificity and efficiency of information processing. Similar mechanisms may govern cellular decisions in higher organisms and be disrupted in cancer.

A geometric analysis of fast-slow models for stochastic gene expression.

Stochastic models for gene expression frequently exhibit dynamics on several different scales. One potential time-scale separation is caused by significant differences in the lifetimes of mRNA and protein; the ratio of the two degradation rates gives a natural small parameter in the resulting chemical master equation, allowing for the application of perturbation techniques. Here, we develop a framework for the analysis of a family of 'fast-slow' models for gene expression that is based on geometric singular perturbation theory. We illustrate our approach by giving a complete characterisation of a standard two-stage model which assumes transcription, translation, and degradation to be first-order reactions. In particular, we present a systematic expansion procedure for the probability-generating function that can in principle be taken to any order in the perturbation parameter, allowing for an approximation of the corresponding propagator probabilities to that same order. For illustrative purposes, we perform this expansion explicitly to first order, both on the fast and the slow time-scales; then, we combine the resulting asymptotics into a composite fast-slow expansion that is uniformly valid in time. In the process, we extend, and prove rigorously, results previously obtained by Shahrezaei and Swain (Proc Natl Acad Sci USA 105(45):17256-17261, 2008) and Bokes et al. (J Math Biol 64(5):829-854, 2012; J Math Biol 65(3):493-520, 2012). We verify our asymptotics by numerical simulation, and we explore its practical applicability and the effects of a variation in the system parameters and the time-scale separation. Focussing on biologically relevant parameter regimes that induce translational bursting, as well as those in which mRNA is frequently transcribed, we find that the first-order correction can significantly improve the steady-state probability distribution. Similarly, in the time-dependent scenario, inclusion of the first-order fast asymptotics results in a uniform approximation for the propagator probabilities that is superior to the slow dynamics alone. Finally, we discuss the generalisation of our geometric framework to models for regulated gene expression that involve additional stages.

Identifying sources of variation and the flow of information in biochemical networks.

To understand how cells control and exploit biochemical fluctuations, we must identify the sources of stochasticity, quantify their effects, and distinguish informative variation from confounding "noise." We present an analysis that allows fluctuations of biochemical networks to be decomposed into multiple components, gives conditions for the design of experimental reporters to measure all components, and provides a technique to predict the magnitude of these components from models. Further, we identify a particular component of variation that can be used to quantify the efficacy of information flow through a biochemical network. By applying our approach to osmosensing in yeast, we can predict the probability of the different osmotic conditions experienced by wild-type yeast and show that the majority of variation can be informational if we include variation generated in response to the cellular environment. Our results are fundamental to quantifying sources of variation and thus are a means to understand biological "design."

Stochastic branching-diffusion models for gene expression.

A challenge to both understanding and modeling biochemical networks is integrating the effects of diffusion and stochasticity. Here, we use the theory of branching processes to give exact analytical expressions for the mean and variance of protein numbers as a function of time and position in a spatial version of an established model of gene expression. We show that both the mean and the magnitude of fluctuations are determined by the protein's Kuramoto length--the typical distance a protein diffuses over its lifetime--and find that the covariance between local concentrations of proteins often increases if there are substantial bursts of synthesis during translation. Using high-throughput data, we estimate that the Kuramoto length of cytoplasmic proteins in budding yeast to be an order of magnitude larger than the cell diameter, implying that many such proteins should have an approximately uniform concentration. For constitutively expressed proteins that live substantially longer than their mRNA, we give an exact expression for the deviation of their local fluctuations from Poisson fluctuations. If the Kuramoto length of mRNA is sufficiently small, we predict that such local fluctuations become approximately Poisson in bacteria in much of the cell, unless translational bursting is exceptionally strong. Our results therefore demonstrate that diffusion can act to both increase and decrease the complexity of fluctuations in biochemical networks.

Unmixing of fluorescence spectra to resolve quantitative time-series measurements of gene expression in plate readers.

BACKGROUND: To connect gene expression with cellular physiology, we need to follow levels of proteins over time. Experiments typically use variants of Green Fluorescent Protein (GFP), and time-series measurements require specialist expertise if single cells are to be followed. Fluorescence plate readers, however, a standard in many laboratories, can in principle provide similar data, albeit at a mean, population level. Nevertheless, extracting the average fluorescence per cell is challenging because autofluorescence can be substantial. RESULTS: Here we propose a general method for correcting plate reader measurements of fluorescent proteins that uses spectral unmixing and determines both the fluorescence per cell and the errors on that fluorescence. Combined with strain collections, such as the GFP fusion collection for budding yeast, our methodology allows quantitative measurements of protein levels of up to hundreds of genes and therefore provides complementary data to high throughput studies of transcription. We illustrate the method by following the induction of the GAL genes in Saccharomyces cerevisiae for over 20 hours in different sugars and argue that the order of appearance of the Leloir enzymes may be to reduce build-up of the toxic intermediate galactose-1-phosphate. Further, we quantify protein levels of over 40 genes, again over 20 hours, after cells experience a change in carbon source (from glycerol to glucose). CONCLUSIONS: Our methodology is sensitive, scalable, and should be applicable to other organisms. By allowing quantitative measurements on a per cell basis over tens of hours and over hundreds of genes, it should increase our understanding of the dynamic changes that drive cellular behaviour.

The fidelity of dynamic signaling by noisy biomolecular networks.

Cells live in changing, dynamic environments. To understand cellular decision-making, we must therefore understand how fluctuating inputs are processed by noisy biomolecular networks. Here we present a general methodology for analyzing the fidelity with which different statistics of a fluctuating input are represented, or encoded, in the output of a signaling system over time. We identify two orthogonal sources of error that corrupt perfect representation of the signal: dynamical error, which occurs when the network responds on average to other features of the input trajectory as well as to the signal of interest, and mechanistic error, which occurs because biochemical reactions comprising the signaling mechanism are stochastic. Trade-offs between these two errors can determine the system's fidelity. By developing mathematical approaches to derive dynamics conditional on input trajectories we can show, for example, that increased biochemical noise (mechanistic error) can improve fidelity and that both negative and positive feedback degrade fidelity, for standard models of genetic autoregulation. For a group of cells, the fidelity of the collective output exceeds that of an individual cell and negative feedback then typically becomes beneficial. We can also predict the dynamic signal for which a given system has highest fidelity and, conversely, how to modify the network design to maximize fidelity for a given dynamic signal. Our approach is general, has applications to both systems and synthetic biology, and will help underpin studies of cellular behavior in natural, dynamic environments.

Ultrasensitivity in phosphorylation-dephosphorylation cycles with little substrate.

Cellular decision-making is driven by dynamic behaviours, such as the preparations for sunrise enabled by circadian rhythms and the choice of cell fates enabled by positive feedback. Such behaviours are often built upon ultrasensitive responses where a linear change in input generates a sigmoidal change in output. Phosphorylation-dephosphorylation cycles are one means to generate ultrasensitivity. Using bioinformatics, we show that in vivo levels of kinases and phosphatases frequently exceed the levels of their corresponding substrates in budding yeast. This result is in contrast to the conditions often required by zero-order ultrasensitivity, perhaps the most well known means for how such cycles become ultrasensitive. We therefore introduce a mechanism to generate ultrasensitivity when numbers of enzymes are higher than numbers of substrates. Our model combines distributive and non-distributive actions of the enzymes with two-stage binding and concerted allosteric transitions of the substrate. We use analytical and numerical methods to calculate the Hill number of the response. For a substrate with [Formula: see text] phosphosites, we find an upper bound of the Hill number of [Formula: see text], and so even systems with a single phosphosite can be ultrasensitive. Two-stage binding, where an enzyme must first bind to a binding site on the substrate before it can access the substrate's phosphosites, allows the enzymes to sequester the substrate. Such sequestration combined with competition for each phosphosite provides an intuitive explanation for the sigmoidal shifts in levels of phosphorylated substrate. Additionally, we find cases for which the response is not monotonic, but shows instead a peak at intermediate levels of input. Given its generality, we expect the mechanism described by our model to often underlay decision-making circuits in eukaryotic cells.

Determining growth rates from bright-field images of budding cells through identifying overlaps.

Much of biochemical regulation ultimately controls growth rate, particularly in microbes. Although time-lapse microscopy visualises cells, determining their growth rates is challenging, particularly for those that divide asymmetrically, like Saccharomyces cerevisiae, because cells often overlap in images. Here, we present the Birth Annotator for Budding Yeast (BABY), an algorithm to determine single-cell growth rates from label-free images. Using a convolutional neural network, BABY resolves overlaps through separating cells by size and assigns buds to mothers by identifying bud necks. BABY uses machine learning to track cells and determine lineages and estimates growth rates as the rates of change of volumes. Using BABY and a microfluidic device, we show that bud growth is likely first sizer- then timer-controlled, that the nuclear concentration of Sfp1, a regulator of ribosome biogenesis, varies before the growth rate does, and that growth rate can be used for real-time control. By estimating single-cell growth rates and so fitness, BABY should generate much biological insight.

Trade-offs and constraints in allosteric sensing.

Sensing extracellular changes initiates signal transduction and is the first stage of cellular decision-making. Yet relatively little is known about why one form of sensing biochemistry has been selected over another. To gain insight into this question, we studied the sensing characteristics of one of the biochemically simplest of sensors: the allosteric transcription factor. Such proteins, common in microbes, directly transduce the detection of a sensed molecule to changes in gene regulation. Using the Monod-Wyman-Changeux model, we determined six sensing characteristics--the dynamic range, the Hill number, the intrinsic noise, the information transfer capacity, the static gain, and the mean response time--as a function of the biochemical parameters of individual sensors and of the number of sensors. We found that specifying one characteristic strongly constrains others. For example, a high dynamic range implies a high Hill number and a high capacity, and vice versa. Perhaps surprisingly, these constraints are so strong that most of the space of characteristics is inaccessible given biophysically plausible ranges of parameter values. Within our approximations, we can calculate the probability distribution of the numbers of input molecules that maximizes information transfer and show that a population of one hundred allosteric transcription factors can in principle distinguish between more than four bands of input concentrations. Our results imply that allosteric sensors are unlikely to have been selected for high performance in one sensing characteristic but for a compromise in the performance of many.

Scalable rule-based modelling of allosteric proteins and biochemical networks.

Much of the complexity of biochemical networks comes from the information-processing abilities of allosteric proteins, be they receptors, ion-channels, signalling molecules or transcription factors. An allosteric protein can be uniquely regulated by each combination of input molecules that it binds. This "regulatory complexity" causes a combinatorial increase in the number of parameters required to fit experimental data as the number of protein interactions increases. It therefore challenges the creation, updating, and re-use of biochemical models. Here, we propose a rule-based modelling framework that exploits the intrinsic modularity of protein structure to address regulatory complexity. Rather than treating proteins as "black boxes", we model their hierarchical structure and, as conformational changes, internal dynamics. By modelling the regulation of allosteric proteins through these conformational changes, we often decrease the number of parameters required to fit data, and so reduce over-fitting and improve the predictive power of a model. Our method is thermodynamically grounded, imposes detailed balance, and also includes molecular cross-talk and the background activity of enzymes. We use our Allosteric Network Compiler to examine how allostery can facilitate macromolecular assembly and how competitive ligands can change the observed cooperativity of an allosteric protein. We also develop a parsimonious model of G protein-coupled receptors that explains functional selectivity and can predict the rank order of potency of agonists acting through a receptor. Our methodology should provide a basis for scalable, modular and executable modelling of biochemical networks in systems and synthetic biology.

Analytical distributions for stochastic gene expression.

Gene expression is significantly stochastic making modeling of genetic networks challenging. We present an approximation that allows the calculation of not only the mean and variance, but also the distribution of protein numbers. We assume that proteins decay substantially more slowly than their mRNA and confirm that many genes satisfy this relation by using high-throughput data from budding yeast. For a two-stage model of gene expression, with transcription and translation as first-order reactions, we calculate the protein distribution for all times greater than several mRNA lifetimes and thus qualitatively predict the distribution of times for protein levels to first cross an arbitrary threshold. If in addition the fluctuates between inactive and active states, we can find the steady-state protein distribution, which can be bimodal if fluctuations of the promoter are slow. We show that our assumptions imply that protein synthesis occurs in geometrically distributed bursts and allows mRNA to be eliminated from a master equation description. In general, we find that protein distributions are asymmetric and may be poorly characterized by their mean and variance. Through maximum likelihood methods, our expressions should therefore allow more quantitative comparisons with experimental data. More generally, we introduce a technique to derive a simpler, effective dynamics for a stochastic system by eliminating a fast variable.

Tracing the sources of cellular variation.

As the adage says, variety is the spice of life, and despite our best attempts, cells, even those with the same genome, never seem to behave the same. By combining mathematical and experimental analyses, Colman-Lerner and colleagues propose, in a recent issue of Nature, a method to delicately unravel the sources of this variation. Applying their technique to the pheromone response in budding yeast, they show that much of the observed variation originates from cell cycle effects and is dependent on levels of pathway input.