Simulating the decentralized processes of the human immune system in a virtual anatomy model

BACKGROUND: Many physiological processes within the human body can be perceived and modeled as large systems of interacting particles or swarming agents. The complex processes of the human immune system prove to be challenging to capture and illustrate without proper reference to the spacial distribution of immune-related organs and systems. Our work focuses on physical aspects of immune system processes, which we implement through swarms of agents. This is our first prototype for integrating different immune processes into one comprehensive virtual physiology simulation. RESULTS: Using agent-based methodology and a 3-dimensional modeling and visualization environment (LINDSAY Composer), we present an agent-based simulation of the decentralized processes in the human immune system. The agents in our model - such as immune cells, viruses and cytokines - interact through simulated physics in two different, compartmentalized and decentralized 3-dimensional environments namely, (1) within the tissue and (2) inside a lymph node. While the two environments are separated and perform their computations asynchronously, an abstract form of communication is allowed in order to replicate the exchange, transportation and interaction of immune system agents between these sites. The distribution of simulated processes, that can communicate across multiple, local CPUs or through a network of machines, provides a starting point to build decentralized systems that replicate larger-scale processes within the human body, thus creating integrated simulations with other physiological systems, such as the circulatory, endocrine, or nervous system. Ultimately, this system integration across scales is our goal for the LINDSAY Virtual Human project. CONCLUSIONS: Our current immune system simulations extend our previous work on agent-based simulations by introducing advanced visualizations within the context of a virtual human anatomy model. We also demonstrate how to distribute a collection of connected simulations over a network of computers. As a future endeavour, we plan to use parameter tuning techniques on our model to further enhance its biological credibility. We consider these in silico experiments and their associated modeling and optimization techniques as essential components in further enhancing our capabilities of simulating a whole-body, decentralized immune system, to be used both for medical education and research as well as for virtual studies in immunoinformatics

Biological Simulation and Evolutionary Optimization: Modelling the Physiology Behind Influenza A Infection

Using agent-based methodology and a 3-dimensional modelling and visualization environment (LINDSAY Composer), we present an agent-based simulation of the decentralized processes in the human immune system. The agents in our model – such as immune cells, viruses and cytokines – interact through simulated physics in two different, compartmentalized and decentralized 3-dimensional environments namely, (1) within the tissue and (2) inside a lymph node. While the two environments are separated and perform their computations asynchronously, an abstract form of communication is allowed in order to replicate the exchange, transportation and interaction of immune system agents between these sites. The distribution of simulated processes, that can communicate across multiple, local CPUs or through a network of machines, provides a starting point to build decentralized systems that replicate larger-scale processes within the human body, thus creating integrated simulations with other physiological systems, such as the circulatory, endocrine, or nervous system. One of the challenges of modelling biological systems is choosing the parameter values which lend it biological credibility. As a potential solution, we propose a parameter tuning approach using Particle Swarm Optimization. This approach relies on a graphical representation of an expected outcome as the metric for evaluating the feasibility of a particular set of parameters. As part of our experiments, we apply the optimization approach to the parameters of the clonal selection mechanism within the simulated lymph node. The results of the optimization allow us to understand the benefits and limitations of using this approach, as well as predict its applicability to larger, more complex biological simulations

Optimization of Swarm-Based Simulations

In computational swarms, large numbers of reactive agents are simulated. The swarm individuals may coordinate their movements in a “search space” to create efficient routes, to occupy niches, or to find the highest peaks. From a more general perspective though, swarms are a means of representation and computation to bridge the gap between local, individual interactions, and global, emergent phenomena. Computational swarms bear great advantages over other numeric methods, for instance, regarding their extensibility, potential for real-time interaction, dynamic interaction topologies, close translation between natural science theory and the computational model, and the integration of multiscale and multiphysics aspects. However, the more comprehensive a swarm-based model becomes, the more demanding its configuration and the more costly its computation become. In this paper, we present an approach to effectively configure and efficiently compute swarm-based simulations by means of heuristic, population-based optimization techniques. We emphasize the commonalities of several of our recent studies that shed light on top-down model optimization and bottom-up abstraction techniques, culminating in a postulation of a general concept of self-organized optimization in swarm-based simulations.Peer Reviewe

All-in-One Pseudo-MS³ Method for the Analysis of Gas-Phase Cleavable Protein Crosslinking Reactions

Crosslinking mass spectrometry (XL-MS) supports structure analysis of individual proteins and highly complex whole-cell interactomes. The identification of crosslinked peptides from enzymatic digests remains challenging, especially at the cell level. Empirical methods that use gas-phase cleavable crosslinkers can simplify the identification process by enabling an MS3-based strategy that turns crosslink identification into a simpler problem of detecting two separable peptides. However, the method is limited to select instrument platforms and is challenged by duty cycle constraints. Here, we revisit a pseudo-MS3 concept that incorporates in-source fragmentation, where a fast switch between gentle high-transmission source conditions and harsher in-source fragmentation settings liberates peptides for standard MS2-based peptide identification. We present an all-in-one method where retention time matches between the crosslink precursor and the liberated peptides establish linkage, and MS2 sequencing identifies the source-liberated peptides. We demonstrate that DC4, a very labile cleavable crosslinker, generates high-intensity peptides in-source. Crosslinks can be identified from these liberated peptides, as they are chromatographically well-resolved from monolinks. Using bovine serum albumin (BSA) as a crosslinking test case, we detect 27% more crosslinks with pseudo-MS3 over a best-in-class MS3 method. While performance is slightly lower for whole-cell lysates (generating two-thirds of the identifications of a standard method), we find that 60% of these hits are unique, highlighting the complementarity of the method

All-in-One Pseudo-MS³ Method for the Analysis of Gas-Phase Cleavable Protein Crosslinking Reactions

Crosslinking mass spectrometry (XL-MS) supports structure analysis of individual proteins and highly complex whole-cell interactomes. The identification of crosslinked peptides from enzymatic digests remains challenging, especially at the cell level. Empirical methods that use gas-phase cleavable crosslinkers can simplify the identification process by enabling an MS3-based strategy that turns crosslink identification into a simpler problem of detecting two separable peptides. However, the method is limited to select instrument platforms and is challenged by duty cycle constraints. Here, we revisit a pseudo-MS3 concept that incorporates in-source fragmentation, where a fast switch between gentle high-transmission source conditions and harsher in-source fragmentation settings liberates peptides for standard MS2-based peptide identification. We present an all-in-one method where retention time matches between the crosslink precursor and the liberated peptides establish linkage, and MS2 sequencing identifies the source-liberated peptides. We demonstrate that DC4, a very labile cleavable crosslinker, generates high-intensity peptides in-source. Crosslinks can be identified from these liberated peptides, as they are chromatographically well-resolved from monolinks. Using bovine serum albumin (BSA) as a crosslinking test case, we detect 27% more crosslinks with pseudo-MS3 over a best-in-class MS3 method. While performance is slightly lower for whole-cell lysates (generating two-thirds of the identifications of a standard method), we find that 60% of these hits are unique, highlighting the complementarity of the method

Automating data analysis for hydrogen/deuterium exchange mass spectrometry using data-independent acquisition methodology

Abstract We present a hydrogen/deuterium exchange workflow coupled to tandem mass spectrometry (HX-MS2) that supports the acquisition of peptide fragment ions alongside their peptide precursors. The approach enables true auto-curation of HX data by mining a rich set of deuterated fragments, generated by collisional-induced dissociation (CID), to simultaneously confirm the peptide ID and authenticate MS1-based deuteration calculations. The high redundancy provided by the fragments supports a confidence assessment of deuterium calculations using a combinatorial strategy. The approach requires data-independent acquisition (DIA) methods that are available on most MS platforms, making the switch to HX-MS2 straightforward. Importantly, we find that HX-DIA enables a proteomics-grade approach and wide-spread applications. Considerable time is saved through auto-curation and complex samples can now be characterized and at higher throughput. We illustrate these advantages in a drug binding analysis of the ultra-large protein kinase DNA-PKcs, isolated directly from mammalian cells