PI-FLAME: A parallel immune system simulator using the FLAME graphic processing unit environment

Agent-based models (ABMs) are increasingly being used to study population dynamics in complex systems, such as the human immune system. Previously, Folcik et al. (The basic immune simulator: an agent-based model to study the interactions between innate and adaptive immunity. Theor Biol Med Model 2007; 4: 39) developed a Basic Immune Simulator (BIS) and implemented it using the Recursive Porous Agent Simulation Toolkit (RePast) ABM simulation framework. However, frameworks such as RePast are designed to execute serially on central processing units and therefore cannot efficiently handle large model sizes. In this paper, we report on our implementation of the BIS using FLAME GPU, a parallel computing ABM simulator designed to execute on graphics processing units. To benchmark our implementation, we simulate the response of the immune system to a viral infection of generic tissue cells. We compared our results with those obtained from the original RePast implementation for statistical accuracy. We observe that our implementation has a 13× performance advantage over the original RePast implementation.

[1]  J. Harty,et al.  Accelerated CD8+ T-cell memory and prime-boost response after dendritic-cell vaccination , 2005, Nature Medicine.

[2]  David M. Mosser,et al.  T Cell Biasing by Activated Dendritic Cells1 , 2004, The Journal of Immunology.

[3]  Jose L Pruneda-Paz,et al.  Bone marrow plasmacytoid dendritic cells can differentiate into myeloid dendritic cells upon virus infection , 2004, Nature Immunology.

[4]  Mike Holcombe X-machines as a basis for dynamic system specification , 1988, Softw. Eng. J..

[5]  Stefania Gallucci,et al.  Natural adjuvants: Endogenous activators of dendritic cells , 1999, Nature Medicine.

[6]  Simon Coakley,et al.  Cellular Level Agent Based Modelling on the Graphics Processing Unit , 2009, 2009 International Workshop on High Performance Computational Systems Biology.

[7]  P. Matzinger The Danger Model: A Renewed Sense of Self , 2002, Science.

[8]  G. An,et al.  The Basic Immune Simulator: An agent-based model to study the interactions between innate and adaptive immunity , 2007, Theoretical Biology and Medical Modelling.

[9]  A. Khoruts,et al.  In Vivo Detection of Dendritic Cell Antigen Presentation to CD4+ T Cells , 1997, The Journal of experimental medicine.

[10]  Massimo Bernaschi,et al.  OpenMP parallelization of agent-based models , 2005, Parallel Comput..

[11]  M. Cooper,et al.  The evolution of adaptive immunity. , 2006, Annual review of immunology.

[12]  Massimo Bernaschi,et al.  A Parallel Simulator of the Immune Response , 1998, HPCN Europe.

[13]  David Fouchet,et al.  A Population Dynamics Analysis of the Interaction between Adaptive Regulatory T Cells and Antigen Presenting Cells , 2008, PloS one.

[14]  M. Jenkins,et al.  Dendritic cell longevity and T cell persistence is controlled by CD154‐CD40 interactions , 2001, European journal of immunology.

[15]  D. Mosser,et al.  Cutting edge: biasing immune responses by directing antigen to macrophage Fc gamma receptors. , 2002, Journal of immunology.

[16]  R. D'Souza SUGARSCAPE ON STEROIDS : SIMULATING OVER A MILLION AGENTS , 2007 .

[17]  Ferdinando Chiacchio,et al.  Agent-Based Modeling of the Immune System: NetLogo, a Promising Framework , 2014, BioMed research international.

[18]  Marc Parmentier,et al.  Specific Recruitment of Antigen-presenting Cells by Chemerin, a Novel Processed Ligand from Human Inflammatory Fluids , 2003, The Journal of experimental medicine.

[19]  Carl T. Bergstrom,et al.  Models of CD8+ responses: 1. What is the antigen-independent proliferation program. , 2003, Journal of theoretical biology.

[20]  Christopher V. Rao,et al.  A Mathematical Model for Neutrophil Gradient Sensing and Polarization , 2007, PLoS Comput. Biol..

[21]  G. Ricevuti,et al.  Host Tissue Damage by Phagocytes , 1997, Annals of the New York Academy of Sciences.

[22]  Mikola Lysenko,et al.  Data-parallel algorithms for agent-based model simulation of tuberculosis on graphics processing units , 2009, SpringSim '09.

[23]  Uwe Aickelin,et al.  Comparing Stochastic Differential Equations and Agent-Based Modelling and Simulation for Early-Stage Cancer , 2014, PloS one.

[24]  Andreas Beilhack,et al.  Immune traffic: a functional overview. , 2003, Lymphatic research and biology.

[25]  L. Parijs,et al.  A Bcl-2-dependent molecular timer regulates the lifespan and immunogenicity of dendritic cells , 2004, Nature Immunology.

[26]  D. Gillespie Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[27]  C. Lumsden,et al.  Stochastic Simulation of Coupled Reaction-Diffusion Processes , 1996 .

[28]  Keith C. Clarke,et al.  Cellular Automata and Agent-Based Models , 2021, Handbook of Regional Science.

[29]  V. Fadok,et al.  Phosphatidylserine-dependent ingestion of apoptotic cells promotes TGF-beta1 secretion and the resolution of inflammation. , 2002, The Journal of clinical investigation.

[30]  N. Romani,et al.  High Level IL-12 Product ion by Murine Dendrit ic Cells : Upregulation via MHC Class II and CD 40 Molecules and Downregulat ion by IL-4 and IL-10 By , 2003 .

[31]  A. Segal,et al.  How neutrophils kill microbes. , 2005, Annual review of immunology.

[32]  Denis Noble,et al.  The rise of computational biology , 2002, Nature Reviews Molecular Cell Biology.

[33]  A. M. Turing,et al.  The chemical basis of morphogenesis , 1952, Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences.

[34]  Manuel Rubio,et al.  Human Monocyte–Derived Dendritic Cells Induce Naive T Cell Differentiation into T Helper Cell Type 2 (Th2) or Th1/Th2 Effectors , 2000, The Journal of experimental medicine.

[35]  James E. Evans,et al.  Molecular identification of a danger signal that alerts the immune system to dying cells , 2003, Nature.

[36]  Daniela M. Romano,et al.  High performance cellular level agent-based simulation with FLAME for the GPU , 2010, Briefings Bioinform..

[37]  S J Merrill,et al.  A model of the role of natural killer cells in immune surveillance — I , 1981, Journal of mathematical biology.

[38]  Eric Bonabeau,et al.  Agent-based modeling: Methods and techniques for simulating human systems , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[39]  P E Seiden,et al.  A model for simulating cognate recognition and response in the immune system. , 1992, Journal of theoretical biology.

[40]  Franco Celada,et al.  IMMSIM, a flexible model for in machina experiments on immune system responses , 2002, Future Gener. Comput. Syst..

[41]  Michael J. North,et al.  Parallel agent-based simulation with Repast for High Performance Computing , 2013, Simul..

[42]  L. D. de Pillis,et al.  A cellular automata model of tumor-immune system interactions. , 2006, Journal of theoretical biology.

[43]  P. Kalinski,et al.  Development of Th1-Inducing Capacity in Myeloid Dendritic Cells Requires Environmental Instruction1 , 2000, The Journal of Immunology.

[44]  Jens H. Krüger,et al.  A Survey of General‐Purpose Computation on Graphics Hardware , 2007, Eurographics.

[45]  P. Ward,et al.  Role of C5a in inflammatory responses. , 2005, Annual review of immunology.

[46]  A. V. Grimstone Molecular biology of the cell (3rd edn) , 1995 .

[47]  Christophe Caux,et al.  Dendritic Cells Enhance Growth and Differentiation of CD40-activated B Lymphocytes , 1997, The Journal of experimental medicine.

[48]  Michael Garland,et al.  Designing efficient sorting algorithms for manycore GPUs , 2009, 2009 IEEE International Symposium on Parallel & Distributed Processing.

[49]  D. Green,et al.  Activation‐induced cell death in T cells , 2003, Immunological reviews.

[50]  Mitchell Kronenberg,et al.  Constitutive Cytokine mRNAs Mark Natural Killer (NK) and NK T Cells Poised for Rapid Effector Function , 2003, The Journal of experimental medicine.

[51]  Michael J. North,et al.  Experiences creating three implementations of the repast agent modeling toolkit , 2006, TOMC.

[52]  L. Allen,et al.  Stochastic models for virus and immune system dynamics. , 2011, Mathematical biosciences.

[53]  Massimo Bernaschi,et al.  C-ImmSim : playing with the immune response , 2004 .

[54]  Shalin H. Naik,et al.  Steady-state and inflammatory dendritic-cell development , 2007, Nature Reviews Immunology.

[55]  P. Srivastava Roles of heat-shock proteins in innate and adaptive immunity , 2002, Nature Reviews Immunology.

[56]  Samuel Eilenberg,et al.  Automata, languages, and machines. A , 1974, Pure and applied mathematics.

[57]  Marc Thilo Figge,et al.  Optimization of Immunoglobulin Substitution Therapy by a Stochastic Immune Response Model , 2009, PloS one.

[58]  H. Van Dyke Parunak,et al.  Agent-Based Modeling vs. Equation-Based Modeling: A Case Study and Users' Guide , 1998, MABS.

[59]  D. Hart,et al.  Dendritic cells: unique leukocyte populations which control the primary immune response. , 1997, Blood.

[60]  David M. Mosser,et al.  Cutting Edge: Biasing Immune Responses by Directing Antigen to Macrophage Fcγ Receptors1 , 2002, The Journal of Immunology.

[61]  Daniela M. Romano,et al.  Template-Driven Agent-Based Modeling and Simulation with CUDA , 2011 .

[62]  F. Varela,et al.  Dynamics of a class of immune networks. I. Global stability of idiotype interactions. , 1990, Journal of theoretical biology.

[63]  Richard A Flavell,et al.  Instruction of Distinct CD4 T Helper Cell Fates by Different Notch Ligands on Antigen-Presenting Cells , 2004, Cell.

[64]  Arturo Casadevall,et al.  Antibody-mediated regulation of cellular immunity and the inflammatory response. , 2003, Trends in immunology.

[65]  Dietmar Tamandl,et al.  Balance between NF-kappaB and JNK/AP-1 activity controls dendritic cell life and death. , 2005, Blood.