Cell region fingerprints enable highly precise single-cell tracking and lineage reconstruction.

Experimental studies of cell growth, inheritance and their associated processes by microscopy require accurate single-cell observations of sufficient duration to reconstruct the genealogy. However, cell tracking-assigning identical cells on consecutive images to a track-is often challenging, resulting in laborious manual verification. Here, we propose fingerprints to identify problematic assignments rapidly. A fingerprint distance compares the structural information contained in the low frequencies of a Fourier transform to measure the similarity between cells in two consecutive images. We show that fingerprints are broadly applicable across cell types and image modalities, provided the image has sufficient structural information. Our tracker (TracX) uses fingerprints to reject unlikely assignments, thereby increasing tracking performance on published and newly generated long-term data sets. For Saccharomyces cerevisiae, we propose a comprehensive model for cell size control at the single-cell and population level centered on the Whi5 regulator, demonstrating how precise tracking can help uncover previously undescribed single-cell biology.

Model-based inference of neutralizing antibody avidities against influenza virus.

To assess the response to vaccination, quantity (concentration) and quality (avidity) of neutralizing antibodies are the most important parameters. Specifically, an increase in avidity indicates germinal center formation, which is required for establishing long-term protection. For influenza, the classical hemagglutination inhibition (HI) assay, however, quantifies a combination of both, and to separately determine avidity requires high experimental effort. We developed from first principles a biophysical model of hemagglutination inhibition to infer IgG antibody avidities from measured HI titers and IgG concentrations. The model accurately describes the relationship between neutralizing antibody concentration/avidity and HI titer, and explains quantitative aspects of the HI assay, such as robustness to pipetting errors and detection limit. We applied our model to infer avidities against the pandemic 2009 H1N1 influenza virus in vaccinated patients (n = 45) after hematopoietic stem cell transplantation (HSCT) and validated our results with independent avidity measurements using an enzyme-linked immunosorbent assay with urea elution. Avidities inferred by the model correlated with experimentally determined avidities (ρ = 0.54, 95% CI = [0.31, 0.70], P < 10-4). The model predicted that increases in IgG concentration mainly contribute to the observed HI titer increases in HSCT patients and that immunosuppressive treatment is associated with lower baseline avidities. Since our approach requires only easy-to-establish measurements as input, we anticipate that it will help to disentangle causes for poor vaccination outcomes also in larger patient populations. This study demonstrates that biophysical modelling can provide quantitative insights into agglutination assays and complement experimental measurements to refine antibody response analyses.

Assumptions on decision making and environment can yield multiple steady states in microbial community models.

BACKGROUND: Microbial community simulations using genome scale metabolic networks (GSMs) are relevant for many application areas, such as the analysis of the human microbiome. Such simulations rely on assumptions about the culturing environment, affecting if the culture may reach a metabolically stationary state with constant microbial concentrations. They also require assumptions on decision making by the microbes: metabolic strategies can be in the interest of individual community members or of the whole community. However, the impact of such common assumptions on community simulation results has not been investigated systematically. RESULTS: Here, we investigate four combinations of assumptions, elucidate how they are applied in literature, provide novel mathematical formulations for their simulation, and show how the resulting predictions differ qualitatively. Our results stress that different assumption combinations give qualitatively different predictions on microbial coexistence by differential substrate utilization. This fundamental mechanism is critically under explored in the steady state GSM literature with its strong focus on coexistence states due to crossfeeding (division of labor). Furthermore, investigating a realistic synthetic community, where the two involved strains exhibit no growth in isolation, but grow as a community, we predict multiple modes of cooperation, even without an explicit cooperation mechanism. CONCLUSIONS: Steady state GSM modelling of microbial communities relies both on assumed decision making principles and environmental assumptions. In principle, dynamic flux balance analysis addresses both. In practice, our methods that address the steady state directly may be preferable, especially if the community is expected to display multiple steady states.

Efficient design of synthetic gene circuits under cell-to-cell variability.

BACKGROUND: Synthetic biologists use and combine diverse biological parts to build systems such as genetic circuits that perform desirable functions in, for example, biomedical or industrial applications. Computer-aided design methods have been developed to help choose appropriate network structures and biological parts for a given design objective. However, they almost always model the behavior of the network in an average cell, despite pervasive cell-to-cell variability. RESULTS: Here, we present a computational framework and an efficient algorithm to guide the design of synthetic biological circuits while accounting for cell-to-cell variability explicitly. Our design method integrates a Non-linear Mixed-Effects (NLME) framework into a Markov Chain Monte-Carlo (MCMC) algorithm for design based on ordinary differential equation (ODE) models. The analysis of a recently developed transcriptional controller demonstrates first insights into design guidelines when trying to achieve reliable performance under cell-to-cell variability. CONCLUSION: We anticipate that our method not only facilitates the rational design of synthetic networks under cell-to-cell variability, but also enables novel applications by supporting design objectives that specify the desired behavior of cell populations.

PolyRound: polytope rounding for random sampling in metabolic networks.

SUMMARY: Random flux sampling is a powerful tool for the constraint-based analysis of metabolic networks. The most efficient sampling method relies on a rounding transform of the constraint polytope, but no available rounding implementation can round all relevant models. By removing redundant polytope constraints on the go, PolyRound simplifies the numerical problem and rounds all the 108 models in the BiGG database without parameter tuning, compared to ∼50% for the state-of-the-art implementation. AVAILABILITY AND IMPLEMENTATION: The implementation is available on gitlab: https://gitlab.com/csb.ethz/PolyRound. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cellular and computational models reveal environmental and metabolic interactions in MMUT-type methylmalonic aciduria.

Methylmalonyl-coenzyme A (CoA) mutase (MMUT)-type methylmalonic aciduria is a rare inherited metabolic disease caused by the loss of function of the MMUT enzyme. Patients develop symptoms resembling those of primary mitochondrial disorders, but the underlying causes of mitochondrial dysfunction remain unclear. Here, we examined environmental and genetic interactions in MMUT deficiency using a combination of computational modeling and cellular models to decipher pathways interacting with MMUT. Immortalized fibroblast (hTERT BJ5ta) MMUT-KO (MUTKO) clones displayed a mild mitochondrial impairment in standard glucose-based medium, but they did not to show increased reliance on respiratory metabolism nor reduced growth or viability. Consistently, our modeling predicted MUTKO specific growth phenotypes only for lower extracellular glutamine concentrations. Indeed, two of three MMUT-deficient BJ5ta cell lines showed a reduced viability in glutamine-free medium. Further, growth on 183 different carbon and nitrogen substrates identified increased NADH (nicotinamide adenine dinucleotide) metabolism of BJ5ta and HEK293 MUTKO cells compared with controls on purine- and glutamine-based substrates. With this knowledge, our modeling predicted 13 reactions interacting with MMUT that potentiate an effect on growth, primarily those of secondary oxidation of propionyl-CoA, oxidative phosphorylation and oxygen diffusion. Of these, we validated 3-hydroxyisobutytyl-CoA hydrolase (HIBCH) in the secondary propionyl-CoA oxidation pathway. Altogether, these results suggest compensation for the loss of MMUT function by increasing anaplerosis through glutamine or by diverting flux away from MMUT through the secondary propionyl-CoA oxidation pathway, which may have therapeutic relevance.

Model-based integration of genomics and metabolomics reveals SNP functionality in Mycobacterium tuberculosis.

Human tuberculosis is caused by members of the Mycobacterium tuberculosis complex (MTBC) that vary in virulence and transmissibility. While genome-wide association studies have uncovered several mutations conferring drug resistance, much less is known about the factors underlying other bacterial phenotypes. Variation in the outcome of tuberculosis infection and diseases has been attributed primarily to patient and environmental factors, but recent evidence indicates an additional role for the genetic diversity among MTBC clinical strains. Here, we used metabolomics to unravel the effect of genetic variation on the strain-specific metabolic adaptive capacity and vulnerability. To define the functionality of single-nucleotide polymorphisms (SNPs) systematically, we developed a constraint-based approach that integrates metabolomic and genomic data. Our model-based predictions correctly classify SNP effects in pyruvate kinase and suggest a genetic basis for strain-specific inherent baseline susceptibility to the antibiotic para-aminosalicylic acid. Our method is broadly applicable across microbial life, opening possibilities for the development of more selective treatment strategies.

Association of Host Factors With Antibody Response to Seasonal Influenza Vaccination in Allogeneic Hematopoietic Stem Cell Transplant Patients.

BACKGROUND: Influenza vaccination efficacy is reduced after hematopoietic stem cell transplantation (HSCT) and patient factors determining vaccination outcomes are still poorly understood. METHODS: We investigated the antibody response to seasonal influenza vaccination in 135 HSCT patients and 69 healthy volunteers (HVs) in a prospective observational multicenter cohort study. We identified patient factors associated with hemagglutination inhibition titers against A/California/2009/H1N1, A/Texas/2012/H3N2, and B/Massachusetts/2012 by multivariable regression on the observed titer levels and on seroconversion/seroprotection categories for comparison. RESULTS: Both regression approaches yielded consistent results but regression on titers estimated associations with higher precision. HSCT patients required 2 vaccine doses to achieve average responses comparable to a single dose in HVs. Prevaccination titers were positively associated with time after transplantation, confirming that HSCT patients can elicit potent antibody responses. However, an unrelated donor, absolute lymphocyte counts below the normal range, and treatment with calcineurin inhibitors lowered the odds of responding. CONCLUSIONS: HSCT patients show a highly heterogeneous vaccine response but, overall, patients benefited from the booster shot and can acquire seroprotective antibodies over the years after transplantation. Several common patient factors lower the odds of responding, urging identification of additional preventive strategies in the poorly responding groups. CLINICAL TRIALS REGISTRATION: NCT03467074.

Probabilistic thermodynamic analysis of metabolic networks.

MOTIVATION: Random sampling of metabolic fluxes can provide a comprehensive description of the capabilities of a metabolic network. However, current sampling approaches do not model thermodynamics explicitly, leading to inaccurate predictions of an organism's potential or actual metabolic operations. RESULTS: We present a probabilistic framework combining thermodynamic quantities with steady-state flux constraints to analyze the properties of a metabolic network. It includes methods for probabilistic metabolic optimization and for joint sampling of thermodynamic and flux spaces. Applied to a model of Escherichia coli, we use the methods to reveal known and novel mechanisms of substrate channeling, and to accurately predict reaction directions and metabolite concentrations. Interestingly, predicted flux distributions are multimodal, leading to discrete hypotheses on E.coli's metabolic capabilities. AVAILABILITY AND IMPLEMENTATION: Python and MATLAB packages available at https://gitlab.com/csb.ethz/pta. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A synthetic multifunctional mammalian pH sensor and CO2 transgene-control device.

All metabolic activities operate within a narrow pH range that is controlled by the CO2-bicarbonate buffering system. We hypothesized that pH could serve as surrogate signal to monitor and respond to the physiological state. By functionally rewiring the human proton-activated cell-surface receptor TDAG8 to chimeric promoters, we created a synthetic signaling cascade that precisely monitors extracellular pH within the physiological range. The synthetic pH sensor could be adjusted by organic acids as well as gaseous CO2 that shifts the CO2-bicarbonate balance toward hydrogen ions. This enabled the design of gas-programmable logic gates, provided remote control of cellular behavior inside microfluidic devices, and allowed for CO2-triggered production of biopharmaceuticals in standard bioreactors. When implanting cells containing the synthetic pH sensor linked to production of insulin into type 1 diabetic mice developing diabetic ketoacidosis, the prosthetic network automatically scored acidic pH and coordinated an insulin expression response that corrected ketoacidosis.

Inflammatory signals directly instruct PU.1 in HSCs via TNF.

The molecular mechanisms governing the transition from hematopoietic stem cells (HSCs) to lineage-committed progenitors remain poorly understood. Transcription factors (TFs) are powerful cell intrinsic regulators of differentiation and lineage commitment, while cytokine signaling has been shown to instruct the fate of progenitor cells. However, the direct regulation of differentiation-inducing hematopoietic TFs by cell extrinsic signals remains surprisingly difficult to establish. PU.1 is a master regulator of hematopoiesis and promotes myeloid differentiation. Here we report that tumor necrosis factor (TNF) can directly and rapidly upregulate PU.1 protein in HSCs in vitro and in vivo. We demonstrate that in vivo, niche-derived TNF is the principal PU.1 inducing signal in HSCs and is both sufficient and required to relay signals from inflammatory challenges to HSCs.

TopoFilter: a MATLAB package for mechanistic model identification in systems biology.

BACKGROUND: To develop mechanistic dynamic models in systems biology, one often needs to identify all (or minimal) representations of the biological processes that are consistent with experimental data, out of a potentially large set of hypothetical mechanisms. However, a simple enumeration of all alternatives becomes quickly intractable when the number of model parameters grows. Selecting appropriate dynamic models out of a large ensemble of models, taking the uncertainty in our biological knowledge and in the experimental data into account, is therefore a key current problem in systems biology. RESULTS: The TopoFilter package addresses this problem in a heuristic and automated fashion by implementing the previously described topological filtering method for Bayesian model selection. It includes a core heuristic for searching the space of submodels of a parametrized model, coupled with a sampling-based exploration of the parameter space. Recent developments of the method allow to balance exhaustiveness and speed of the model space search, to efficiently re-sample parameters, to parallelize the search, and to use custom scoring functions. We use a theoretical example to motivate these features and then demonstrate TopoFilter's applicability for a yeast signaling network with more than 250'000 possible model structures. CONCLUSIONS: TopoFilter is a flexible software framework that makes Bayesian model selection and reduction efficient and scalable to network models of a complexity that represents contemporary problems in, for example, cell signaling. TopoFilter is open-source, available under the GPL-3.0 license at https://gitlab.com/csb.ethz/TopoFilter. It includes installation instructions, a quickstart guide, a description of all package options, and multiple examples.

Controlling cell-to-cell variability with synthetic gene circuits.

Cell-to-cell variability originating, for example, from the intrinsic stochasticity of gene expression, presents challenges for designing synthetic gene circuits that perform robustly. Conversely, synthetic biology approaches are instrumental in uncovering mechanisms underlying variability in natural systems. With a focus on reducing noise in individual genes, the field has established a broad synthetic toolset. This includes noise control by engineering of transcription and translation mechanisms either individually, or in combination to achieve independent regulation of mean expression and its variability. Synthetic feedback circuits use these components to establish more robust operation in closed-loop, either by drawing on, but also by extending traditional engineering concepts. In this perspective, we argue that major conceptual advances will require new theory of control adapted to biology, extensions from single genes to networks, more systematic considerations of origins of variability other than intrinsic noise, and an exploration of how noise shaping, instead of noise reduction, could establish new synthetic functions or help understanding natural functions.

Steady-State Differential Dose Response in Biological Systems.

In pharmacology and systems biology, it is a fundamental problem to determine how biological systems change their dose-response behavior upon perturbations. In particular, it is unclear how topologies, reactions, and parameters (differentially) affect the dose response. Because parameters are often unknown, systematic approaches should directly relate network structure and function. Here, we outline a procedure to compare general non-monotone dose-response curves and subsequently develop a comprehensive theory for differential dose responses of biochemical networks captured by non-equilibrium steady-state linear framework models. Although these models are amenable to analytical derivations of non-equilibrium steady states in principle, their size frequently increases (super) exponentially with model size. We extract general principles of differential responses based on a model's graph structure and thereby alleviate the combinatorial explosion. This allows us, for example, to determine reactions that affect differential responses, to identify classes of networks with equivalent differential, and to reject hypothetical models reliably without needing to know parameter values. We exemplify such applications for models of insulin signaling.

Inferring causal metabolic signals that regulate the dynamic TORC1-dependent transcriptome.

Cells react to nutritional cues in changing environments via the integrated action of signaling, transcriptional, and metabolic networks. Mechanistic insight into signaling processes is often complicated because ubiquitous feedback loops obscure causal relationships. Consequently, the endogenous inputs of many nutrient signaling pathways remain unknown. Recent advances for system-wide experimental data generation have facilitated the quantification of signaling systems, but the integration of multi-level dynamic data remains challenging. Here, we co-designed dynamic experiments and a probabilistic, model-based method to infer causal relationships between metabolism, signaling, and gene regulation. We analyzed the dynamic regulation of nitrogen metabolism by the target of rapamycin complex 1 (TORC1) pathway in budding yeast. Dynamic transcriptomic, proteomic, and metabolomic measurements along shifts in nitrogen quality yielded a consistent dataset that demonstrated extensive re-wiring of cellular networks during adaptation. Our inference method identified putative downstream targets of TORC1 and putative metabolic inputs of TORC1, including the hypothesized glutamine signal. The work provides a basis for further mechanistic studies of nitrogen metabolism and a general computational framework to study cellular processes.

Efficient Characterization of Parametric Uncertainty of Complex (Bio)chemical Networks.

Parametric uncertainty is a particularly challenging and relevant aspect of systems analysis in domains such as systems biology where, both for inference and for assessing prediction uncertainties, it is essential to characterize the system behavior globally in the parameter space. However, current methods based on local approximations or on Monte-Carlo sampling cope only insufficiently with high-dimensional parameter spaces associated with complex network models. Here, we propose an alternative deterministic methodology that relies on sparse polynomial approximations. We propose a deterministic computational interpolation scheme which identifies most significant expansion coefficients adaptively. We present its performance in kinetic model equations from computational systems biology with several hundred parameters and state variables, leading to numerical approximations of the parametric solution on the entire parameter space. The scheme is based on adaptive Smolyak interpolation of the parametric solution at judiciously and adaptively chosen points in parameter space. As Monte-Carlo sampling, it is "non-intrusive" and well-suited for massively parallel implementation, but affords higher convergence rates. This opens up new avenues for large-scale dynamic network analysis by enabling scaling for many applications, including parameter estimation, uncertainty quantification, and systems design.

A tunable synthetic mammalian oscillator.

Autonomous and self-sustained oscillator circuits mediating the periodic induction of specific target genes are minimal genetic time-keeping devices found in the central and peripheral circadian clocks. They have attracted significant attention because of their intriguing dynamics and their importance in controlling critical repair, metabolic and signalling pathways. The precise molecular mechanism and expression dynamics of this mammalian circadian clock are still not fully understood. Here we describe a synthetic mammalian oscillator based on an auto-regulated sense-antisense transcription control circuit encoding a positive and a time-delayed negative feedback loop, enabling autonomous, self-sustained and tunable oscillatory gene expression. After detailed systems design with experimental analyses and mathematical modelling, we monitored oscillating concentrations of green fluorescent protein with tunable frequency and amplitude by time-lapse microscopy in real time in individual Chinese hamster ovary cells. The synthetic mammalian clock may provide an insight into the dynamics of natural periodic processes and foster advances in the design of prosthetic networks in future gene and cell therapies.

Cutting the wires: modularization of cellular networks for experimental design.

Understanding naturally evolved cellular networks requires the consecutive identification and revision of the interactions between relevant molecular species. In this process, initially often simplified and incomplete networks are extended by integrating new reactions or whole subnetworks to increase consistency between model predictions and new measurement data. However, increased consistency with experimental data alone is not sufficient to show the existence of biomolecular interactions, because the interplay of different potential extensions might lead to overall similar dynamics. Here, we present a graph-based modularization approach to facilitate the design of experiments targeted at independently validating the existence of several potential network extensions. Our method is based on selecting the outputs to measure during an experiment, such that each potential network extension becomes virtually insulated from all others during data analysis. Each output defines a module that only depends on one hypothetical network extension, and all other outputs act as virtual inputs to achieve insulation. Given appropriate experimental time-series measurements of the outputs, our modules can be analyzed, simulated, and compared to the experimental data separately. Our approach exemplifies the close relationship between structural systems identification and modularization, an interplay that promises development of related approaches in the future.

Bridging the gaps in systems biology.

Systems biology aims at creating mathematical models, i.e., computational reconstructions of biological systems and processes that will result in a new level of understanding-the elucidation of the basic and presumably conserved "design" and "engineering" principles of biomolecular systems. Thus, systems biology will move biology from a phenomenological to a predictive science. Mathematical modeling of biological networks and processes has already greatly improved our understanding of many cellular processes. However, given the massive amount of qualitative and quantitative data currently produced and number of burning questions in health care and biotechnology needed to be solved is still in its early phases. The field requires novel approaches for abstraction, for modeling bioprocesses that follow different biochemical and biophysical rules, and for combining different modules into larger models that still allow realistic simulation with the computational power available today. We have identified and discussed currently most prominent problems in systems biology: (1) how to bridge different scales of modeling abstraction, (2) how to bridge the gap between topological and mechanistic modeling, and (3) how to bridge the wet and dry laboratory gap. The future success of systems biology largely depends on bridging the recognized gaps.

Near-optimal experimental design for model selection in systems biology.

MOTIVATION: Biological systems are understood through iterations of modeling and experimentation. Not all experiments, however, are equally valuable for predictive modeling. This study introduces an efficient method for experimental design aimed at selecting dynamical models from data. Motivated by biological applications, the method enables the design of crucial experiments: it determines a highly informative selection of measurement readouts and time points. RESULTS: We demonstrate formal guarantees of design efficiency on the basis of previous results. By reducing our task to the setting of graphical models, we prove that the method finds a near-optimal design selection with a polynomial number of evaluations. Moreover, the method exhibits the best polynomial-complexity constant approximation factor, unless P = NP. We measure the performance of the method in comparison with established alternatives, such as ensemble non-centrality, on example models of different complexity. Efficient design accelerates the loop between modeling and experimentation: it enables the inference of complex mechanisms, such as those controlling central metabolic operation. AVAILABILITY: Toolbox 'NearOED' available with source code under GPL on the Machine Learning Open Source Software Web site (mloss.org).