Novel descriptive and model based statistical approaches in immunology and signal transduction

Biological systems are usually complex nonlinear systems of which we only have a limited understanding. Here we show three different aspects of investigating such systems. We present a method to extract detailed knowledge from typical biological trajectory data, which have randomness as a main characteristic. The migration of immune cells, such as leukocytes, are a key example of our study. The application of our methodology leads to the discovery of novel random walk behaviour of leukocyte migration. Furthermore we use the gathered knowledge to construct the under- lying mathematical model that captures the behaviour of leukocytes, or more precisely macrophages and neutrophils, under acute injury. Any model of a biological system has little predictive power if it is not compared to collected data. We present a pipeline of how complex spatio- temporal trajectory data can be used to calibrate our model of leukocyte migration. The pipeline employs approximate methods in a Bayesian framework. Using the same approach we are able to learn additional information about the underlying signalling network, which is not directly apparent in the cell migration data. While these two methods can be seen as data processing and analysis, we show in the last part of this work how to assess the information content of experiments. The choice of an experiment with the highest information content out of a set of possible experiments leads us to the problem of optimal experimental design. We develop and implement an algorithm for simulation based Bayesian experimental design in order to learn parameters of a given model. We validate our algorithm with the help of toy examples and apply it to examples in immunology (Hes1 transcription regulation) and signal transduction (growth factor induced MAPK pathway)

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

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

Accurate Reconstruction of Cell and Particle Tracks from 3D Live Imaging Data

Spatial structures often constrain the 3D movement of cells or particles in vivo, yet this information is obscured when microscopy data are analyzed using standard approaches. Here, we present methods, called unwrapping and Riemannian manifold learning, for mapping particle-tracking data along unseen and irregularly curved surfaces onto appropriate 2D representations. This is conceptually similar to the problem of reconstructing accurate geography from conventional Mercator maps, but our methods do not require prior knowledge of the environments’ physical structure. Unwrapping and Riemannian manifold learning accurately recover the underlying 2D geometry from 3D imaging data without the need for fiducial marks. They outperform standard x-y projections, and unlike standard dimensionality reduction techniques, they also successfully detect both bias and persistence in cell migration modes. We demonstrate these features on simulated data and zebrafish and Drosophila in vivo immune cell trajectory datasets. Software packages that implement unwrapping and Riemannian manifold learning are provided

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

An Unexpected Major Role for Proteasome-Catalyzed Peptide Splicing in Generation of T Cell Epitopes:Is There Relevance for Vaccine Development?

Efficient and safe induction of CD8+ T cell responses is a desired characteristic of vaccines against intracellular pathogens. To achieve this, a new generation of safe vaccines is being developed accommodating single, dominant antigens of pathogens of interest. In particular, the selection of such antigens is challenging, since due to HLA polymorphism the ligand specificities and immunodominance hierarchies of pathogen-specific CD8+ T cell responses differ throughout the human population. A recently discovered mechanism of proteasome-mediated CD8+ T cell epitope generation, i.e., by proteasome-catalyzed peptide splicing (PCPS), expands the pool of peptides and antigens, presented by MHC class I HLA molecules. On the cell surface, one-third of the presented self-peptides are generated by PCPS, which coincides with one-fourth in terms of abundance. Spliced epitopes are targeted by CD8+ T cell responses during infection and, like non-spliced epitopes, can be identified within antigen sequences using a novel in silico strategy. The existence of spliced epitopes, by enlarging the pool of peptides available for presentation by different HLA variants, opens new opportunities for immunotherapies and vaccine design.</p

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

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

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

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

P38 and JNK have opposing effects on persistence of in vivo leukocyte migration in zebrafish

The recruitment and migration of macrophages and neutrophils is an important process during the early stages of the innate immune system in response to acute injury. Transgenic pu.1:EGFP zebrafish permit the acquisition of leukocyte migration trajectories during inflammation. Currently, these high-quality live-imaging data are mainly analysed using general statistics, for example, cell velocity. Here, we present a spatio-temporal analysis of the cell dynamics using transition matrices, which provide information of the type of cell migration. We find evidence that leukocytes exhibit types of migratory behaviour, which differ from previously described random walk processes. Dimethyl sulfoxide treatment decreased the level of persistence at early time points after wounding and ablated temporal dependencies observed in untreated embryos. We then use pharmacological inhibition of p38 and c-Jun N-terminal kinase mitogen-activated protein kinases to determine their effects on in vivo leukocyte migration patterns and discuss how they modify the characteristics of the cell migration process. In particular, we find that their respective inhibition leads to decreased and increased levels of persistent motion in leukocytes following wounding. This example shows the high level of information content, which can be gained from live-imaging data if appropriate statistical tools are used