Parallel simulation of apoptotic receptor-clustering on GPGPU many-core architectures

Apoptosis, the programmed cell death, is a physiological process that handles the removal of unwanted or damaged cells in living organisms. The process itself is initiated by signaling through tumor necrosis factor (TNF) receptors and ligands, which form clusters on the cell membrane. The exact function of this process is not yet fully understood and currently subject of basic research. Different mathematical models have been developed to describe and simulate the apoptotic receptor-clustering. In this interdisciplinary work, a previously introduced model of the apoptotic receptor-clustering has been extended by a new receptor type to allow a more precise description and simulation of the signaling process. Due to the high computational requirements of the model, an efficient algorithmic mapping to a modern many-core GPGPU architecture has been developed. Such architectures enable high-performance computing (HPC) simulation tasks on the desktop at low costs. The developed mapping reduces average simulation times from months to days (peak speedup of 256x), allowing the productive use of the model in research.

[1]  Peter Scheurich,et al.  Single chain TNF derivatives with individually mutated receptor binding sites reveal differential stoichiometry of ligand receptor complex formation for TNFR1 and TNFR2. , 2010, Cellular signalling.

[2]  Piet Hut,et al.  A hierarchical O(N log N) force-calculation algorithm , 1986, Nature.

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

[4]  Jeremiah P. Ostriker,et al.  A self-consistent field method for galactic dynamics , 1992 .

[5]  Robert G. Belleman,et al.  High Performance Direct Gravitational N-body Simulations on Graphics Processing Units , 2007, ArXiv.

[6]  S. Swain.  Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences , 1984 .

[7]  Andrea Zappe,et al.  Dual function of cysteine rich domain (CRD) 1 of TNF receptor type 1: conformational stabilization of CRD2 and control of receptor responsiveness. , 2010, Cellular signalling.

[8]  Janet E. Jones.  On the Determination of Molecular Fields. I. From the Variation of the Viscosity of a Gas with Temperature , 1924 .

[9]  C. W. Gardiner,et al.  Handbook of stochastic methods - for physics, chemistry and the natural sciences, Second Edition , 1986, Springer series in synergetics.

[10]  Rüdiger Westermann,et al.  UberFlow: a GPU-based particle engine , 2004, SIGGRAPH '04.

[11]  A. Perelson.  Receptor clustering on a cell surface. I. theory of receptor cross-linking by ligands bearing two chemically identical functional groups , 1980 .

[12]  Peter Scheurich,et al.  Control of Receptor-induced Signaling Complex Formation by the Kinetics of Ligand/Receptor Interaction* , 2002, The Journal of Biological Chemistry.

[13]  H. Levine,et al.  A thermodynamic model for receptor clustering. , 1999, Biophysical journal.

[14]  Martin Falk,et al.  Modeling and visualization of receptor clustering on the cellular membrane , 2011, 2011 IEEE Symposium on Biological Data Visualization (BioVis)..

[15]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .