Individual-based model for the movement of cells with a cell group in response to applied electric fields

Code written by Calina Copos (c.copos@northeastern.edu)

Part of manuscript titled "Galvanotactic directionality of cell groups depends on group size"
C.Copos et al. Proc Nat Acad Sci 2025

To run the code for polarity establishment in a pair of cells, open MATLAB and enter the command >> galvano_2d

The code in its current state will generate a small cell cluster (N=28) in the presence of the PI3K inhibitor, pre- and post-electric field application. Changing lines 18-21 to:
pos0        = load('DiskMeshConfig/UnitCirclePointsN128.txt');
bdd         = load('DiskMeshConfig/UnitCircleBoundaryN128.txt');
neighbors   = load('DiskMeshConfig/UnitCircleVertexNeighborsN128.txt');
T           = load('DiskMeshConfig/UnitCircleTrianglesN128.txt');
will run a new simulation for a PI3K inhibited large cell group (N=128).

Any questions can be addressed to c.copos@northeastern.edu.