  Over the last two decades, scientific literature has been blooming with
various means of simulating epithelial cell colonies. Each of these simulations
can be separated by their respective efficiency (expressed in terms of consumed
computational resources), the amount of cells/the size of tissues that can be
simulated, the time scale of the simulated dynamics and the coarse grained
level of precision. Choosing the right algorithm for the simulation of
epithelial cells and tissues is a compromise between each of these key
elements. Irrespective of the method, each algorithm includes part, or all, of
the following features: short-range membrane-mediated attraction between cells,
soft-core repulsion between cells, cell proliferation, cell death, cell
motility, fluctuations, etc. We will first give a non-exhaustive overview of
commonly used modeling approaches for tissues at a mesoscopic level, giving a
rough idea of the coarse-graining decisions made for every one of them. Then we
will dive into greater detail on how to implement a relaxation procedure
according to the Vertex Model, refreshing aspects of the theoretical
groundwork, describing required data structures and simulation steps and
pointing out details of the simulation that can present pitfalls to a
first-time implementation of the model.
