Computer Simulation of Cellular Patterning Within the Drosophila Pupal Eye

We present a computer simulation and associated experimental validation of assembly of glial-like support cells into the interweaving hexagonal lattice that spans the Drosophila pupal eye. This process of cell movements organizes the ommatidial array into a functional pattern. Unlike earlier simulations that focused on the arrangements of cells within individual ommatidia, here we examine the local movements that lead to large-scale organization of the emerging eye field. Simulations based on our experimental observations of cell adhesion, cell death, and cell movement successfully patterned a tracing of an emerging wild-type pupal eye. Surprisingly, altering cell adhesion had only a mild effect on patterning, contradicting our previous hypothesis that the patterning was primarily the result of preferential adhesion between IRM-class surface proteins. Instead, our simulations highlighted the importance of programmed cell death (PCD) as well as a previously unappreciated variable: the expansion of cells' apical surface areas, which promoted rearrangement of neighboring cells. We tested this prediction experimentally by preventing expansion in the apical area of individual cells: patterning was disrupted in a manner predicted by our simulations. Our work demonstrates the value of combining computer simulation with in vivo experiments to uncover novel mechanisms that are perpetuated throughout the eye field. It also demonstrates the utility of the Glazier–Graner–Hogeweg model (GGH) for modeling the links between local cellular interactions and emergent properties of developing epithelia as well as predicting unanticipated results in vivo

Conjectures on exact solution of three - dimensional (3D) simple orthorhombic Ising lattices

We report the conjectures on the three-dimensional (3D) Ising model on simple orthorhombic lattices, together with the details of calculations for a putative exact solution. Two conjectures, an additional rotation in the fourth curled-up dimension and the weight factors on the eigenvectors, are proposed to serve as a boundary condition to deal with the topologic problem of the 3D Ising model. The partition function of the 3D simple orthorhombic Ising model is evaluated by spinor analysis, by employing these conjectures. Based on the validity of the conjectures, the critical temperature of the simple orthorhombic Ising lattices could be determined by the relation of KK* = KK' + KK'' + K'K'' or sinh 2K sinh 2(K' + K'' + K'K''/K) = 1. For a simple cubic Ising lattice, the critical point is putatively determined to locate exactly at the golden ratio xc = exp(-2Kc) = (sq(5) - 1)/2, as derived from K* = 3K or sinh 2K sinh 6K = 1. If the conjectures would be true, the specific heat of the simple orthorhombic Ising system would show a logarithmic singularity at the critical point of the phase transition. The spontaneous magnetization and the spin correlation functions of the simple orthorhombic Ising ferromagnet are derived explicitly. The putative critical exponents derived explicitly for the simple orthorhombic Ising lattices are alpha = 0, beta = 3/8, gamma = 5/4, delta = 13/3, eta = 1/8 and nu = 2/3, showing the universality behavior and satisfying the scaling laws. The cooperative phenomena near the critical point are studied and the results obtained based on the conjectures are compared with those of the approximation methods and the experimental findings. The 3D to 2D crossover phenomenon differs with the 2D to 1D crossover phenomenon and there is a gradual crossover of the exponents from the 3D values to the 2D ones.Comment: 176 pages, 4 figure

Active contraction of cardiac cells: a reduced model for sarcomere dynamics with cooperative interactions

We propose a reduced ODE model for the mechanical activation of cardiac myofilaments, which is based on explicit spatial representation of nearest-neighbour interactions. Our model is derived from the cooperative Markov Chain model of Washio et al. (Cell Mol Bioeng 5(1):113–126, 2012), under the assumption of conditional independence of specific sets of events. This physically motivated assumption allows to drastically reduce the number of degrees of freedom, thus resulting in a significantly large computational saving. Indeed, the original Markov Chain model involves a huge number of degrees of freedom (order of (Formula presented.)) and is solved by means of the Monte Carlo method, which notoriously reaches statistical convergence in a slow fashion. With our reduced model, instead, numerical simulations can be carried out by solving a system of ODEs, reducing the computational time by more than 10, 000 times. Moreover, the reduced model is accurate with respect to the original Markov Chain model. We show that the reduced model is capable of reproducing physiological steady-state force–calcium and force–length relationships with the observed asymmetry in apparent cooperativity near the calcium level producing half activation. Finally, we also report good qualitative and quantitative agreement with experimental measurements under dynamic conditions