Juliane Liepe

ABC-SysBio–approximate Bayesian computation in Python with GPU support

Motivation: The growing field of systems biology has driven demand for flexible tools to model and simulate biological systems. Two established problems in the modeling of biological processes are model selection and the estimation of associated parameters. A number of statistical approaches, both frequentist and Bayesian, have been proposed to answer these questions. Results: Here we present a Python package, ABC-SysBio, that implements parameter inference and model selection for dynamical systems in an approximate Bayesian computation (ABC) framework. ABC-SysBio combines three algorithms: ABC rejection sampler, ABC SMC for parameter inference and ABC SMC for model selection. It is designed to work with models written in Systems Biology Markup Language (SBML). Deterministic and stochastic models can be analyzed in ABC-SysBio. Availability: http://abc-sysbio.sourceforge.net Contact: christopher.barnes@imperial.ac.uk; ttoni@imperial.ac.uk

A large fraction of HLA class I ligands are proteasome-generated spliced peptides

New players in the repertoire Antigen-presenting cells, such as macrophages and dendritic cells, activate immunological T cells by presenting them with antigens bound by major histocompatibility complexes (MHCs). The proteasome typically processes these antigens, which include peptides derived from both self and microbial origins. Liepe et al. now report that, surprisingly, a large fraction of peptides bound to class I MHC on multiple human cell types are spliced together by the proteasome from two different fragments of the same protein. Such merged peptides might turn out to be useful in vaccine or cancer immunotherapy development. Science, this issue p. 354 Spliced peptides make up a major fraction of the epitopes presented by MHC class I on multiple human cell types. The proteasome generates the epitopes presented on human leukocyte antigen (HLA) class I molecules that elicit CD8+ T cell responses. Reports of proteasome-generated spliced epitopes exist, but they have been regarded as rare events. Here, however, we show that the proteasome-generated spliced peptide pool accounts for one-third of the entire HLA class I immunopeptidome in terms of diversity and one-fourth in terms of abundance. This pool also represents a unique set of antigens, possessing particular and distinguishing features. We validated this observation using a range of complementary experimental and bioinformatics approaches, as well as multiple cell types. The widespread appearance and abundance of proteasome-catalyzed peptide splicing events has implications for immunobiology and autoimmunity theories and may provide a previously untapped source of epitopes for use in vaccines and cancer immunotherapy.

Maximizing the Information Content of Experiments in Systems Biology

Our understanding of most biological systems is in its infancy. Learning their structure and intricacies is fraught with challenges, and often side-stepped in favour of studying the function of different gene products in isolation from their physiological context. Constructing and inferring global mathematical models from experimental data is, however, central to systems biology. Different experimental setups provide different insights into such systems. Here we show how we can combine concepts from Bayesian inference and information theory in order to identify experiments that maximize the information content of the resulting data. This approach allows us to incorporate preliminary information; it is global and not constrained to some local neighbourhood in parameter space and it readily yields information on parameter robustness and confidence. Here we develop the theoretical framework and apply it to a range of exemplary problems that highlight how we can improve experimental investigations into the structure and dynamics of biological systems and their behavior.

GPU accelerated biochemical network simulation

Motivation: Mathematical modelling is central to systems and synthetic biology. Using simulations to calculate statistics or to explore parameter space is a common means for analysing these models and can be computationally intensive. However, in many cases, the simulations are easily parallelizable. Graphics processing units (GPUs) are capable of efficiently running highly parallel programs and outperform CPUs in terms of raw computing power. Despite their computational advantages, their adoption by the systems biology community is relatively slow, since differences in hardware architecture between GPUs and CPUs complicate the porting of existing code. Results: We present a Python package, cuda-sim, that provides highly parallelized algorithms for the repeated simulation of biochemical network models on NVIDIA CUDA GPUs. Algorithms are implemented for the three popular types of model formalisms: the LSODA algorithm for ODE integration, the Euler–Maruyama algorithm for SDE simulation and the Gillespie algorithm for MJP simulation. No knowledge of GPU computing is required from the user. Models can be specified in SBML format or provided as CUDA code. For running a large number of simulations in parallel, up to 360-fold decrease in simulation runtime is attained when compared to single CPU implementations. Availability: http://cuda-sim.sourceforge.net/ Contact: christopher.barnes@imperial.ac.uk; m.stumpf@imperial.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.

Proteasome isoforms exhibit only quantitative differences in cleavage and epitope generation

Immunoproteasomes are considered to be optimised to process Ags and to alter the peptide repertoire by generating a qualitatively different set of MHC class I epitopes. Whether the immunoproteasome at the biochemical level, influence the quality rather than the quantity of the immuno‐genic peptide pool is still unclear. Here, we quantified the cleavage‐site usage by human standard‐ and immunoproteasomes, and proteasomes from immuno‐subunit‐deficient mice, as well as the peptides generated from model polypeptides. We show in this study that the different proteasome isoforms can exert significant quantitative differences in the cleavage‐site usage and MHC class I restricted epitope production. However, independent of the proteasome isoform and substrates studied, no evidence was obtained for the abolishment of the specific cleavage‐site usage, or for differences in the quality of the peptides generated. Thus, we conclude that the observed differences in MHC class I restricted Ag presentation between standard‐ and immunoproteasomes are due to quantitative differences in the proteasome‐generated antigenic peptides.

The 20S Proteasome Splicing Activity Discovered by SpliceMet

The identification of proteasome-generated spliced peptides (PSP) revealed a new unpredicted activity of the major cellular protease. However, so far characterization of PSP was entirely dependent on the availability of patient-derived cytotoxic CD8+ T lymphocytes (CTL) thus preventing a systematic investigation of proteasome-catalyzed peptide splicing (PCPS). For an unrestricted PSP identification we here developed SpliceMet, combining the computer-based algorithm ProteaJ with in vitro proteasomal degradation assays and mass spectrometry. By applying SpliceMet for the analysis of proteasomal processing products of four different substrate polypeptides, derived from human tumor as well as viral antigens, we identified fifteen new spliced peptides generated by PCPS either by cis or from two separate substrate molecules, i.e., by trans splicing. Our data suggest that 20S proteasomes represent a molecular machine that, due to its catalytic and structural properties, facilitates the generation of spliced peptides, thereby providing a pool of qualitatively new peptides from which functionally relevant products may be selected.

Quantitative time-resolved analysis reveals intricate, differential regulation of standard- and immuno-proteasomes

Proteasomal protein degradation is a key determinant of protein half-life and hence of cellular processes ranging from basic metabolism to a host of immunological processes. Despite its importance the mechanisms regulating proteasome activity are only incompletely understood. Here we use an iterative and tightly integrated experimental and modelling approach to develop, explore and validate mechanistic models of proteasomal peptide-hydrolysis dynamics. The 20S proteasome is a dynamic enzyme and its activity varies over time because of interactions between substrates and products and the proteolytic and regulatory sites; the locations of these sites and the interactions between them are predicted by the model, and experimentally supported. The analysis suggests that the rate-limiting step of hydrolysis is the transport of the substrates into the proteasome. The transport efficiency varies between human standard- and immuno-proteasomes thereby impinging upon total degradation rate and substrate cleavage-site usage. DOI: http://dx.doi.org/10.7554/eLife.07545.001

Systems Analysis of the Dynamic Inflammatory Response to Tissue Damage Reveals Spatiotemporal Properties of the Wound Attractant Gradient

Summary In the acute inflammatory phase following tissue damage, cells of the innate immune system are rapidly recruited to sites of injury by pro-inflammatory mediators released at the wound site. Although advances in live imaging allow us to directly visualize this process in vivo, the precise identity and properties of the primary immune damage attractants remain unclear, as it is currently impossible to directly observe and accurately measure these signals in tissues. Here, we demonstrate that detailed information about the attractant signals can be extracted directly from the in vivo behavior of the responding immune cells. By applying inference-based computational approaches to analyze the in vivo dynamics of the Drosophila inflammatory response, we gain new detailed insight into the spatiotemporal properties of the attractant gradient. In particular, we show that the wound attractant is released by wound margin cells, rather than by the wounded tissue per se, and that it diffuses away from this source at rates far slower than those of previously implicated signals such as H2O2 and ATP, ruling out these fast mediators as the primary chemoattractant. We then predict, and experimentally test, how competing attractant signals might interact in space and time to regulate multi-step cell navigation in the complex environment of a healing wound, revealing a period of receptor desensitization after initial exposure to the damage attractant. Extending our analysis to model much larger wounds, we uncover a dynamic behavioral change in the responding immune cells in vivo that is prognostic of whether a wound will subsequently heal or not. Video Abstract

Information theory and signal transduction systems: from molecular information processing to network inference.

Sensing and responding to the environment are two essential functions that all biological organisms need to master for survival and successful reproduction. Developmental processes are marshalled by a diverse set of signalling and control systems, ranging from systems with simple chemical inputs and outputs to complex molecular and cellular networks with non-linear dynamics. Information theory provides a powerful and convenient framework in which such systems can be studied; but it also provides the means to reconstruct the structure and dynamics of molecular interaction networks underlying physiological and developmental processes. Here we supply a brief description of its basic concepts and introduce some useful tools for systems and developmental biologists. Along with a brief but thorough theoretical primer, we demonstrate the wide applicability and biological application-specific nuances by way of different illustrative vignettes. In particular, we focus on the characterisation of biological information processing efficiency, examining cell-fate decision making processes, gene regulatory network reconstruction, and efficient signal transduction experimental design.

Inference of random walk models to describe leukocyte migration.

While the majority of cells in an organism are static and remain relatively immobile in their tissue, migrating cells occur commonly during developmental processes and are crucial for a functioning immune response. The mode of migration has been described in terms of various types of random walks. To understand the details of the migratory behaviour we rely on mathematical models and their calibration to experimental data. Here we propose an approximate Bayesian inference scheme to calibrate a class of random walk models characterized by a specific, parametric particle re-orientation mechanism to observed trajectory data. We elaborate the concept of transition matrices (TMs) to detect random walk patterns and determine a statistic to quantify these TM to make them applicable for inference schemes. We apply the developed pipeline to in vivo trajectory data of macrophages and neutrophils, extracted from zebrafish that had undergone tail transection. We find that macrophage and neutrophils exhibit very distinct biased persistent random walk patterns, where the strengths of the persistence and bias are spatio-temporally regulated. Furthermore, the movement of macrophages is far less persistent than that of neutrophils in response to wounding.