Parallelizing the Cellular Potts Model on graphics processing units

Abstract The Cellular Potts Model (CPM) is a lattice based modeling technique used for simulating cellular structures in computational biology. The computational complexity of the model means that current serial implementations restrict the size of simulation to a level well below biological relevance. Parallelization on computing clusters enables scaling the size of the simulation but marginally addresses computational speed due to the limited memory bandwidth between nodes. In this paper we present new data-parallel algorithms and data structures for simulating the Cellular Potts Model on graphics processing units. Our implementations handle most terms in the Hamiltonian, including cell–cell adhesion constraint, cell volume constraint, cell surface area constraint, and cell haptotaxis. We use fine level checkerboards with lock mechanisms using atomic operations to enable consistent updates while maintaining a high level of parallelism. A new data-parallel memory allocation algorithm has been developed to handle cell division. Tests show that our implementation enables simulations of > 10 6 cells with lattice sizes of up to 256 3 on a single graphics card. Benchmarks show that our implementation runs ∼80× faster than serial implementations, and ∼5× faster than previous parallel implementations on computing clusters consisting of 25 nodes. The wide availability and economy of graphics cards mean that our techniques will enable simulation of realistically sized models at a fraction of the time and cost of previous implementations and are expected to greatly broaden the scope of CPM applications.

[1]  Glazier,et al.  Simulation of biological cell sorting using a two-dimensional extended Potts model. , 1992, Physical review letters.

[2]  Yi-Tsann Jiang,et al.  Extended large-Q Potts model simulation of foam drainage , 1996 .

[3]  Benny Pinkas,et al.  Analysis of the Linux random number generator , 2006, 2006 IEEE Symposium on Security and Privacy (S&P'06).

[4]  Nan Chen,et al.  A Parallel Implementation of the Cellular Potts Model for Simulation of Cell-Based Morphogenesis , 2006, ACRI.

[5]  Archana Jagannatam,et al.  Mersenne Twister A Pseudo-Random Number Generator , 2007 .

[6]  Takuji Nishimura,et al.  Mersenne twister: a 623-dimensionally equidistributed uniform pseudo-random number generator , 1998, TOMC.

[7]  Gerson G. H. Cavalheiro,et al.  An Efficient Parallel Algorithm to Evolve Simulations of the Cellular Potts Model , 2005, Parallel Process. Lett..

[8]  Jesús A. Izaguirre,et al.  A framework for three-dimensional simulation of morphogenesis , 2005, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[9]  Wayne Luk,et al.  A comparison of CPUs, GPUs, FPGAs, and massively parallel processor arrays for random number generation , 2009, FPGA '09.

[10]  Y Jiang,et al.  Hysteresis and avalanches in two-dimensional foam rheology simulations. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[11]  S. A. Wright,et al.  Potts-model grain growth simulations: Parallel algorithms and applications , 1997 .

[12]  James A. Glazier,et al.  Foam Drainage: Extended Large-Q Potts Model Simulation , 1996 .

[13]  Barkema,et al.  Parallel simulation of the Ising model. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[14]  Yi Jiang,et al.  Topography of Extracellular Matrix Mediates Vascular Morphogenesis and Migration Speeds in Angiogenesis , 2009, PLoS Comput. Biol..

[15]  Jie Cheng,et al.  Programming Massively Parallel Processors. A Hands-on Approach , 2010, Scalable Comput. Pract. Exp..

[16]  Peter Y. K. Cheung,et al.  Proceedings of the ACM/SIGDA International Symposium on Field Programmable Gate Arrays, FPGA 2001, Monterey, CA, USA, February 11-13, 2001 , 2001, FPGA.

[17]  Philip Heidelberger,et al.  Analysis of parallel replicated simulations under a completion time constraint , 1991, TOMC.

[18]  Brenda M Rubenstein,et al.  The role of extracellular matrix in glioma invasion: a cellular Potts model approach. , 2008, Biophysical journal.