Gene Regulatory Networks from Single Cell Data for Exploring Cell Fate Decisions.

Single cell experimental techniques now allow us to quantify gene expression in up to thousands of individual cells. These data reveal the changes in transcriptional state that occur as cells progress through development and adopt specialized cell fates. In this chapter we describe in detail how to use our network inference algorithm (PIDC)-and the associated software package NetworkInference.jl-to infer functional interactions between genes from the observed gene expression patterns. We exploit the large sample sizes and inherent variability of single cell data to detect statistical dependencies between genes that indicate putative (co-)regulatory relationships, using multivariate information measures that can capture complex statistical relationships. We provide guidelines on how best to combine this analysis with other complementary methods designed to explore single cell data, and how to interpret the resulting gene regulatory network models to gain insight into the processes regulating cell differentiation.

Gene Regulatory Network Inference from Single-Cell Data Using Multivariate Information Measures.

While single-cell gene expression experiments present new challenges for data processing, the cell-to-cell variability observed also reveals statistical relationships that can be used by information theory. Here, we use multivariate information theory to explore the statistical dependencies between triplets of genes in single-cell gene expression datasets. We develop PIDC, a fast, efficient algorithm that uses partial information decomposition (PID) to identify regulatory relationships between genes. We thoroughly evaluate the performance of our algorithm and demonstrate that the higher-order information captured by PIDC allows it to outperform pairwise mutual information-based algorithms when recovering true relationships present in simulated data. We also infer gene regulatory networks from three experimental single-cell datasets and illustrate how network context, choices made during analysis, and sources of variability affect network inference. PIDC tutorials and open-source software for estimating PID are available. PIDC should facilitate the identification of putative functional relationships and mechanistic hypotheses from single-cell transcriptomic data.

How to deal with parameters for whole-cell modelling.

Dynamical systems describing whole cells are on the verge of becoming a reality. But as models of reality, they are only useful if we have realistic parameters for the molecular reaction rates and cell physiological processes. There is currently no suitable framework to reliably estimate hundreds, let alone thousands, of reaction rate parameters. Here, we map out the relative weaknesses and promises of different approaches aimed at redressing this issue. While suitable procedures for estimation or inference of the whole (vast) set of parameters will, in all likelihood, remain elusive, some hope can be drawn from the fact that much of the cellular behaviour may be explained in terms of smaller sets of parameters. Identifying such parameter sets and assessing their behaviour is now becoming possible even for very large systems of equations, and we expect such methods to become central tools in the development and analysis of whole-cell models.

Model of Host-Pathogen Interaction Dynamics Links In Vivo Optical Imaging and Immune Responses.

Tracking disease progression in vivo is essential for the development of treatments against bacterial infection. Optical imaging has become a central tool for in vivo tracking of bacterial population development and therapeutic response. For a precise understanding of in vivo imaging results in terms of disease mechanisms derived from detailed postmortem observations, however, a link between the two is needed. Here, we develop a model that provides that link for the investigation of Citrobacter rodentium infection, a mouse model for enteropathogenic Escherichia coli (EPEC). We connect in vivo disease progression of C57BL/6 mice infected with bioluminescent bacteria, imaged using optical tomography and X-ray computed tomography, to postmortem measurements of colonic immune cell infiltration. We use the model to explore changes to both the host immune response and the bacteria and to evaluate the response to antibiotic treatment. The developed model serves as a novel tool for the identification and development of new therapeutic interventions.

MEANS: python package for Moment Expansion Approximation, iNference and Simulation.

MOTIVATION: Many biochemical systems require stochastic descriptions. Unfortunately these can only be solved for the simplest cases and their direct simulation can become prohibitively expensive, precluding thorough analysis. As an alternative, moment closure approximation methods generate equations for the time-evolution of the system's moments and apply a closure ansatz to obtain a closed set of differential equations; that can become the basis for the deterministic analysis of the moments of the outputs of stochastic systems. RESULTS: We present a free, user-friendly tool implementing an efficient moment expansion approximation with parametric closures that integrates well with the IPython interactive environment. Our package enables the analysis of complex stochastic systems without any constraints on the number of species and moments studied and the type of rate laws in the system. In addition to the approximation method our package provides numerous tools to help non-expert users in stochastic analysis. AVAILABILITY AND IMPLEMENTATION: https://github.com/theosysbio/means CONTACTS: m.stumpf@imperial.ac.uk or e.lakatos13@imperial.ac.uk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Topological sensitivity analysis for systems biology.

Mathematical models of natural systems are abstractions of much more complicated processes. Developing informative and realistic models of such systems typically involves suitable statistical inference methods, domain expertise, and a modicum of luck. Except for cases where physical principles provide sufficient guidance, it will also be generally possible to come up with a large number of potential models that are compatible with a given natural system and any finite amount of data generated from experiments on that system. Here we develop a computational framework to systematically evaluate potentially vast sets of candidate differential equation models in light of experimental and prior knowledge about biological systems. This topological sensitivity analysis enables us to evaluate quantitatively the dependence of model inferences and predictions on the assumed model structures. Failure to consider the impact of structural uncertainty introduces biases into the analysis and potentially gives rise to misleading conclusions.

Kinetic and computational evidence for an intermediate in the hydrolysis of sulfonate esters.

The hydrolytic reactions of sulfonate esters have previously been considered to occur by concerted mechanisms. We now report the observation of a break in a Brønsted correlation for the alkaline hydrolysis of aryl benzenesulfonates. On either side of a break-point, ß(leaving group) values of -0.27 (pK(a) < 8.5) and -0.97 (pK(a) > 8.5) are measured. These data are consistent with a two-step mechanism involving a pentavalent intermediate that is also supported by QM/MM calculations. The emerging scenario can be explained by the combined effect of a strong nucleophile with a poor leaving group that compel a usually concerted reaction to favour a stepwise process.

An efficient, multiply promiscuous hydrolase in the alkaline phosphatase superfamily.

We report a catalytically promiscuous enzyme able to efficiently promote the hydrolysis of six different substrate classes. Originally assigned as a phosphonate monoester hydrolase (PMH) this enzyme exhibits substantial second-order rate accelerations ((k(cat)/K(M))/k(w)), ranging from 10(7) to as high as 10(19), for the hydrolyses of phosphate mono-, di-, and triesters, phosphonate monoesters, sulfate monoesters, and sulfonate monoesters. This substrate collection encompasses a range of substrate charges between 0 and -2, transition states of a different nature, and involves attack at two different reaction centers (P and S). Intrinsic reactivities (half-lives) range from 200 days to 10(5) years under near neutrality. The substantial rate accelerations for a set of relatively difficult reactions suggest that efficient catalysis is not necessarily limited to efficient stabilization of just one transition state. The crystal structure of PMH identifies it as a member of the alkaline phosphatase superfamily. PMH encompasses four of the native activities previously observed in this superfamily and extends its repertoire by two further activities, one of which, sulfonate monoesterase, has not been observed previously for a natural enzyme. PMH is thus one of the most promiscuous hydrolases described to date. The functional links between superfamily activities can be presumed to have played a role in functional evolution by gene duplication.

Efficient catalytic promiscuity for chemically distinct reactions.

High catalytic proficiencies observed for the native and promiscuous reaction of the Pseudomonas aeruginosa arylsulfatase (PAS; the picture shows transition states of the two substrates with corresponding binding constants K(tx)) suggest that the trade-off between high activity and tight specificity can be substantially relaxed.