Deformability and collision-induced reorientation enhance cell topotaxis in dense microenvironments.

In vivo, cells navigate through complex environments filled with obstacles such as other cells and the extracellular matrix. Recently, the term "topotaxis" has been introduced for navigation along topographic cues such as obstacle density gradients. Experimental and mathematical efforts have analyzed topotaxis of single cells in pillared grids with pillar density gradients. A previous model based on active Brownian particles (ABPs) has shown that ABPs perform topotaxis, i.e., drift toward lower pillar densities, due to decreased effective persistence lengths at high pillar densities. The ABP model predicted topotactic drifts of up to 1% of the instantaneous speed, whereas drifts of up to 5% have been observed experimentally. We hypothesized that the discrepancy between the ABP and the experimental observations could be in 1) cell deformability and 2) more complex cell-pillar interactions. Here, we introduce a more detailed model of topotaxis based on the cellular Potts model (CPM). To model persistent cells we use the Act model, which mimics actin-polymerization-driven motility, and a hybrid CPM-ABP model. Model parameters were fitted to simulate the experimentally found motion of Dictyostelium discoideum on a flat surface. For starved D. discoideum, the topotactic drifts predicted by both CPM variants are closer to the experimental results than the previous ABP model due to a larger decrease in persistence length. Furthermore, the Act model outperformed the hybrid model in terms of topotactic efficiency, as it shows a larger reduction in effective persistence time in dense pillar grids. Also pillar adhesion can slow down cells and decrease topotaxis. For slow and less-persistent vegetative D. discoideum cells, both CPMs predicted a similar small topotactic drift. We conclude that deformable cell volume results in higher topotactic drift compared with ABPs, and that feedback of cell-pillar collisions on cell persistence increases drift only in highly persistent cells.

Computational modelling of cell motility modes emerging from cell-matrix adhesion dynamics.

Lymphocytes have been described to perform different motility patterns such as Brownian random walks, persistent random walks, and Lévy walks. Depending on the conditions, such as confinement or the distribution of target cells, either Brownian or Lévy walks lead to more efficient interaction with the targets. The diversity of these motility patterns may be explained by an adaptive response to the surrounding extracellular matrix (ECM). Indeed, depending on the ECM composition, lymphocytes either display a floating motility without attaching to the ECM, or sliding and stepping motility with respectively continuous or discontinuous attachment to the ECM, or pivoting behaviour with sustained attachment to the ECM. Moreover, on the long term, lymphocytes either perform a persistent random walk or a Brownian-like movement depending on the ECM composition. How the ECM affects cell motility is still incompletely understood. Here, we integrate essential mechanistic details of the lymphocyte-matrix adhesions and lymphocyte intrinsic cytoskeletal induced cell propulsion into a Cellular Potts model (CPM). We show that the combination of de novo cell-matrix adhesion formation, adhesion growth and shrinkage, adhesion rupture, and feedback of adhesions onto cell propulsion recapitulates multiple lymphocyte behaviours, for different lymphocyte subsets and various substrates. With an increasing attachment area and increased adhesion strength, the cells' speed and persistence decreases. Additionally, the model predicts random walks with short-term persistent but long-term subdiffusive properties resulting in a pivoting type of motility. For small adhesion areas, the spatial distribution of adhesions emerges as a key factor influencing cell motility. Small adhesions at the front allow for more persistent motility than larger clusters at the back, despite a similar total adhesion area. In conclusion, we present an integrated framework to simulate the effects of ECM proteins on cell-matrix adhesion dynamics. The model reveals a sufficient set of principles explaining the plasticity of lymphocyte motility.

A Novel Function of TLR2 and MyD88 in the Regulation of Leukocyte Cell Migration Behavior During Wounding in Zebrafish Larvae.

Toll-like receptor (TLR) signaling via myeloid differentiation factor 88 protein (MyD88) has been indicated to be involved in the response to wounding. It remains unknown whether the putative role of MyD88 in wounding responses is due to a control of leukocyte cell migration. The aim of this study was to explore in vivo whether TLR2 and MyD88 are involved in modulating neutrophil and macrophage cell migration behavior upon zebrafish larval tail wounding. Live cell imaging of tail-wounded larvae was performed in tlr2 and myd88 mutants and their corresponding wild type siblings. In order to visualize cell migration following tissue damage, we constructed double transgenic lines with fluorescent markers for macrophages and neutrophils in all mutant and sibling zebrafish lines. Three days post fertilization (dpf), tail-wounded larvae were studied using confocal laser scanning microscopy (CLSM) to quantify the number of recruited cells at the wounding area. We found that in both tlr2-/- and myd88-/- groups the recruited neutrophil and macrophage numbers are decreased compared to their wild type sibling controls. Through analyses of neutrophil and macrophage migration patterns, we demonstrated that both tlr2 and myd88 control the migration direction of distant neutrophils upon wounding. Furthermore, in both the tlr2 and the myd88 mutants, macrophages migrated more slowly toward the wound edge. Taken together, our findings show that tlr2 and myd88 are involved in responses to tail wounding by regulating the behavior and speed of leukocyte migration in vivo.

Predicting Metabolism from Gene Expression in an Improved Whole-Genome Metabolic Network Model of Danio rerio.

Zebrafish is a useful modeling organism for the study of vertebrate development, immune response, and metabolism. Metabolic studies can be aided by mathematical reconstructions of the metabolic network of zebrafish. These list the substrates and products of all biochemical reactions that occur in the zebrafish. Mathematical techniques such as flux-balance analysis then make it possible to predict the possible metabolic flux distributions that optimize, for example, the turnover of food into biomass. The only available genome-scale reconstruction of zebrafish metabolism is ZebraGEM. In this study, we present ZebraGEM 2.0, an updated and validated version of ZebraGEM. ZebraGEM 2.0 is extended with gene-protein-reaction associations (GPRs) that are required to integrate genetic data with the metabolic model. To demonstrate the use of these GPRs, we performed an in silico genetic screening for knockouts of metabolic genes and validated the results against published in vivo genetic knockout and knockdown screenings. Among the single knockout simulations, we identified 74 essential genes, whose knockout stopped growth completely. Among these, 11 genes are known have an abnormal knockout or knockdown phenotype in vivo (partial), and 41 have human homologs associated with metabolic diseases. We also added the oxidative phosphorylation pathway, which was unavailable in the published version of ZebraGEM. The updated model performs better than the original model on a predetermined list of metabolic functions. We also determined a minimal feed composition. The oxidative phosphorylation pathways were validated by comparing with published experiments in which key components of the oxidative phosphorylation pathway were pharmacologically inhibited. To test the utility of ZebraGEM2.0 for obtaining new results, we integrated gene expression data from control and Mycobacterium marinum-infected zebrafish larvae. The resulting model predicts impeded growth and altered histidine metabolism in the infected larvae.