Practical hardware for evolvable robots.

The evolutionary robotics field offers the possibility of autonomously generating robots that are adapted to desired tasks by iteratively optimising across successive generations of robots with varying configurations until a high-performing candidate is found. The prohibitive time and cost of actually building this many robots means that most evolutionary robotics work is conducted in simulation, but to apply evolved robots to real-world problems, they must be implemented in hardware, which brings new challenges. This paper explores in detail the design of an example system for realising diverse evolved robot bodies, and specifically how this interacts with the evolutionary process. We discover that every aspect of the hardware implementation introduces constraints that change the evolutionary space, and exploring this interplay between hardware constraints and evolution is the key contribution of this paper. In simulation, any robot that can be defined by a suitable genetic representation can be implemented and evaluated, but in hardware, real-world limitations like manufacturing/assembly constraints and electrical power delivery mean that many of these robots cannot be built, or will malfunction in operation. This presents the novel challenge of how to constrain an evolutionary process within the space of evolvable phenotypes to only those regions that are practically feasible: the viable phenotype space. Methods of phenotype filtering and repair were introduced to address this, and found to degrade the diversity of the robot population and impede traversal of the exploration space. Furthermore, the degrees of freedom permitted by the hardware constraints were found to be poorly matched to the types of morphological variation that would be the most useful in the target environment. Consequently, the ability of the evolutionary process to generate robots with effective adaptations was greatly reduced. The conclusions from this are twofold. 1) Designing a hardware platform for evolving robots requires different thinking, in which all design decisions should be made with reference to their impact on the viable phenotype space. 2) It is insufficient to just evolve robots in simulation without detailed consideration of how they will be implemented in hardware, because the hardware constraints have a profound impact on the evolutionary space.

A study of error diversity in robotic swarms for task partitioning in foraging tasks.

Often in swarm robotics, an assumption is made that all robots in the swarm behave the same and will have a similar (if not the same) error model. However, in reality, this is not the case, and this lack of uniformity in the error model, and other operations, can lead to various emergent behaviors. This paper considers the impact of the error model and compares robots in a swarm that operate using the same error model (uniform error) against each robot in the swarm having a different error model (thus introducing error diversity). Experiments are presented in the context of a foraging task. Simulation and physical experimental results show the importance of the error model and diversity in achieving the expected swarm behavior.

Simulating CXCR5 Dynamics in Complex Tissue Microenvironments.

To effectively navigate complex tissue microenvironments, immune cells sense molecular concentration gradients using G-protein coupled receptors. However, due to the complexity of receptor activity, and the multimodal nature of chemokine gradients in vivo, chemokine receptor activity in situ is poorly understood. To address this issue, we apply a modelling and simulation approach that permits analysis of the spatiotemporal dynamics of CXCR5 expression within an in silico B-follicle with single-cell resolution. Using this approach, we show that that in silico B-cell scanning is robust to changes in receptor numbers and changes in individual kinetic rates of receptor activity, but sensitive to global perturbations where multiple parameters are altered simultaneously. Through multi-objective optimization analysis we find that the rapid modulation of CXCR5 activity through receptor binding, desensitization and recycling is required for optimal antigen scanning rates. From these analyses we predict that chemokine receptor signaling dynamics regulate migration in complex tissue microenvironments to a greater extent than the total numbers of receptors on the cell surface.

B cell zone reticular cell microenvironments shape CXCL13 gradient formation.

Through the formation of concentration gradients, morphogens drive graded responses to extracellular signals, thereby fine-tuning cell behaviors in complex tissues. Here we show that the chemokine CXCL13 forms both soluble and immobilized gradients. Specifically, CXCL13+ follicular reticular cells form a small-world network of guidance structures, with computer simulations and optimization analysis predicting that immobilized gradients created by this network promote B cell trafficking. Consistent with this prediction, imaging analysis show that CXCL13 binds to extracellular matrix components in situ, constraining its diffusion. CXCL13 solubilization requires the protease cathepsin B that cleaves CXCL13 into a stable product. Mice lacking cathepsin B display aberrant follicular architecture, a phenotype associated with effective B cell homing to but not within lymph nodes. Our data thus suggest that reticular cells of the B cell zone generate microenvironments that shape both immobilized and soluble CXCL13 gradients.

Strategies for calibrating models of biology.

Computational and mathematical modelling has become a valuable tool for investigating biological systems. Modelling enables prediction of how biological components interact to deliver system-level properties and extrapolation of biological system performance to contexts and experimental conditions where this is unknown. A model's value hinges on knowing that it faithfully represents the biology under the contexts of use, or clearly ascertaining otherwise and thus motivating further model refinement. These qualities are evaluated through calibration, typically formulated as identifying model parameter values that align model and biological behaviours as measured through a metric applied to both. Calibration is critical to modelling but is often underappreciated. A failure to appropriately calibrate risks unrepresentative models that generate erroneous insights. Here, we review a suite of strategies to more rigorously challenge a model's representation of a biological system. All are motivated by features of biological systems, and illustrative examples are drawn from the modelling literature. We examine the calibration of a model against distributions of biological behaviours or outcomes, not only average values. We argue for calibration even where model parameter values are experimentally ascertained. We explore how single metrics can be non-distinguishing for complex systems, with multiple-component dynamic and interaction configurations giving rise to the same metric output. Under these conditions, calibration is insufficiently constraining and the model non-identifiable: multiple solutions to the calibration problem exist. We draw an analogy to curve fitting and argue that calibrating a biological model against a single experiment or context is akin to curve fitting against a single data point. Though useful for communicating model results, we explore how metrics that quantify heavily emergent properties may not be suitable for use in calibration. Lastly, we consider the role of sensitivity and uncertainty analysis in calibration and the interpretation of model results. Our goal in this manuscript is to encourage a deeper consideration of calibration, and how to increase its capacity to either deliver faithful models or demonstrate them otherwise.

Using Emulation to Engineer and Understand Simulations of Biological Systems.

Modeling and simulation techniques have demonstrated success in studying biological systems. As the drive to better capture biological complexity leads to more sophisticated simulators, it becomes challenging to perform statistical analyses that help translate predictions into increased understanding. These analyses may require repeated executions and extensive sampling of high-dimensional parameter spaces: analyses that may become intractable due to time and resource limitations. Significant reduction in these requirements can be obtained using surrogate models, or emulators, that can rapidly and accurately predict the output of an existing simulator. We apply emulation to evaluate and enrich understanding of a previously published agent-based simulator of lymphoid tissue organogenesis, showing an ensemble of machine learning techniques can reproduce results obtained using a suite of statistical analyses within seconds. This performance improvement permits incorporation of previously intractable analyses, including multi-objective optimization to obtain parameter sets that yield a desired response, and Approximate Bayesian Computation to assess parametric uncertainty. To facilitate exploitation of emulation in simulation-focused studies, we extend our open source statistical package, spartan, to provide a suite of tools for emulator development, validation, and application. Overcoming resource limitations permits enriched evaluation and refinement, easing translation of simulator insights into increased biological understanding.

Application of Modeling Approaches to Explore Vaccine Adjuvant Mode-of-Action.

Novel adjuvant technologies have a key role in the development of next-generation vaccines, due to their capacity to modulate the duration, strength and quality of the immune response. The AS01 adjuvant is used in the malaria vaccine RTS,S/AS01 and in the licensed herpes-zoster vaccine (Shingrix) where the vaccine has proven its ability to generate protective responses with both robust humoral and T-cell responses. For many years, animal models have provided insights into adjuvant mode-of-action (MoA), generally through investigating individual genes or proteins. Furthermore, modeling and simulation techniques can be utilized to integrate a variety of different data types; ranging from serum biomarkers to large scale "omics" datasets. In this perspective we present a framework to create a holistic integration of pre-clinical datasets and immunological literature in order to develop an evidence-based hypothesis of AS01 adjuvant MoA, creating a unified view of multiple experiments. Furthermore, we highlight how holistic systems-knowledge can serve as a basis for the construction of models and simulations supporting exploration of key questions surrounding adjuvant MoA. Using the Systems-Biology-Graphical-Notation, a tool for graphical representation of biological processes, we have captured high-level cellular behaviors and interactions, and cytokine dynamics during the early immune response, which are substantiated by a series of diagrams detailing cellular dynamics. Through explicitly describing AS01 MoA we have built a consensus of understanding across multiple experiments, and so we present a framework to integrate modeling approaches into exploring adjuvant MoA, in order to guide experimental design, interpret results and inform rational design of vaccines.

GABA Regulation of Burst Firing in Hippocampal Astrocyte Neural Circuit: A Biophysical Model.

It is now widely accepted that glia cells and gamma-aminobutyric acidergic (GABA) interneurons dynamically regulate synaptic transmission and neuronal activity in time and space. This paper presents a biophysical model that captures the interaction between an astrocyte cell, a GABA interneuron and pre/postsynaptic neurons. Specifically, GABA released from a GABA interneuron triggers in astrocytes the release of calcium (Ca 2+) from the endoplasmic reticulum via the inositol 1, 4, 5-trisphosphate (IP 3) pathway. This results in gliotransmission which elevates the presynaptic transmission probability rate (PR) causing weight potentiation and a gradual increase in postsynaptic neuronal firing, that eventually stabilizes. However, by capturing the complex interactions between IP 3, generated from both GABA and the 2-arachidonyl glycerol (2-AG) pathway, and PR, this paper shows that this interaction not only gives rise to an initial weight potentiation phase but also this phase is followed by postsynaptic bursting behavior. Moreover, the model will show that there is a presynaptic frequency range over which burst firing can occur. The proposed model offers a novel cellular level mechanism that may underpin both seizure-like activity and neuronal synchrony across different brain regions.

Best Practices to Maximize the Use and Reuse of Quantitative and Systems Pharmacology Models: Recommendations From the United Kingdom Quantitative and Systems Pharmacology Network.

The lack of standardization in the way that quantitative and systems pharmacology (QSP) models are developed, tested, and documented hinders their reproducibility, reusability, and expansion or reduction to alternative contexts. This in turn undermines the potential impact of QSP in academic, industrial, and regulatory frameworks. This article presents a minimum set of recommendations from the UK Quantitative and Systems Pharmacology Network (UK QSP Network) to guide QSP practitioners seeking to maximize their impact, and stakeholders considering the use of QSP models in their environment.

Exploring Self-Repair in a Coupled Spiking Astrocyte Neural Network.

It is now known that astrocytes modulate the activity at the tripartite synapses where indirect signaling via the retrograde messengers, endocannabinoids, leads to a localized self-repairing capability. In this paper, a self-repairing spiking astrocyte neural network (SANN) is proposed to demonstrate a distributed self-repairing capability at the network level. The SANN uses a novel learning rule that combines the spike-timing-dependent plasticity (STDP) and Bienenstock, Cooper, and Munro (BCM) learning rules (hereafter referred to as the BSTDP rule). In this learning rule, the synaptic weight potentiation is not only driven by the temporal difference between the presynaptic and postsynaptic neuron firing times but also by the postsynaptic neuron activity. We will show in this paper that the BSTDP modulates the height of the plasticity window to establish an input-output mapping (in the learning phase) and also maintains this mapping (via self-repair) if synaptic pathways become dysfunctional. It is the functional dependence of postsynaptic neuron firing activity on the height of the plasticity window that underpins how the proposed SANN self-repairs on the fly. The SANN also uses the coupling between the tripartite synapses and γ -GABAergic interneurons. This interaction gives rise to a presynaptic neuron frequency filtering capability that serves to route information, represented as spike trains, to different neurons in the subsequent layers of the SANN. The proposed SANN follows a feedforward architecture with multiple interneuron pathways and astrocytes modulate synaptic activity at the hidden and output neuronal layers. The self-repairing capability will be demonstrated in a robotic obstacle avoidance application, and the simulation results will show that the SANN can maintain learned maneuvers at synaptic fault densities of up to 80% regardless of the fault locations.

Tissue-specific transcriptomic changes associated with AmBisome® treatment of BALB/c mice with experimental visceral leishmaniasis.

Background: Liposomal amphotericin B (AmBisome®) as a treatment modality for visceral leishmaniasis (VL) has had significant impact on patient care in some but not all regions where VL is endemic.  As the mode of action of AmBisome® in vivo is poorly understood, we compared the tissue-specific transcriptome in drug-treated vs untreated mice with experimental VL.    Methods:  BALB/c mice infected with L. donovani were treated with 8mg/kg AmBisome®, resulting in parasite elimination from liver and spleen over a 7-day period. At day 1 and day 7 post treatment (R x+1 and R x+7), transcriptomic profiling was performed on spleen and liver tissue from treated and untreated mice and uninfected mice.  BALB/c mice infected with M. bovis BCG (an organism resistant to amphotericin B) were analysed to distinguish between direct effects of AmBisome® and those secondary to parasite death.   Results: AmBisome® treatment lead to rapid parasitological clearance.  At R x+1, spleen and liver displayed only 46 and 88 differentially expressed (DE) genes (P<0.05; 2-fold change) respectively. In liver, significant enrichment was seen for pathways associated with TNF, fatty acids and sterol biosynthesis.  At R x+7, the number of DE genes was increased (spleen, 113; liver 400).  In spleen, these included many immune related genes known to be involved in anti-leishmanial immunity. In liver, changes in transcriptome were largely accounted for by loss of granulomas.   PCA analysis indicated that treatment only partially restored homeostasis.  Analysis of BCG-infected mice treated with AmBisome® revealed a pattern of immune modulation mainly targeting macrophage function.   Conclusions: Our data indicate that the tissue response to AmBisome® treatment varies between target organs and that full restoration of homeostasis is not achieved at parasitological cure.  The pathways required to restore homeostasis deserve fuller attention, to understand mechanisms associated with treatment failure and relapse and to promote more rapid restoration of immune competence.

Macrophage Transactivation for Chemokine Production Identified as a Negative Regulator of Granulomatous Inflammation Using Agent-Based Modeling.

Cellular activation in trans by interferons, cytokines, and chemokines is a commonly recognized mechanism to amplify immune effector function and limit pathogen spread. However, an optimal host response also requires that collateral damage associated with inflammation is limited. This may be particularly so in the case of granulomatous inflammation, where an excessive number and/or excessively florid granulomas can have significant pathological consequences. Here, we have combined transcriptomics, agent-based modeling, and in vivo experimental approaches to study constraints on hepatic granuloma formation in a murine model of experimental leishmaniasis. We demonstrate that chemokine production by non-infected Kupffer cells in the Leishmania donovani-infected liver promotes competition with infected KCs for available iNKT cells, ultimately inhibiting the extent of granulomatous inflammation. We propose trans-activation for chemokine production as a novel broadly applicable mechanism that may operate early in infection to limit excessive focal inflammation.

Livestock Disease Management for Trading Across Different Regulatory Regimes.

The maintenance of livestock health depends on the combined actions of many different actors, both within and across different regulatory frameworks. Prior work recognised that private risk management choices have the ability to reduce the spread of infection to trading partners. We evaluate the efficiency of farmers' alternative biosecurity choices in terms of their own-benefits from unilateral strategies and quantify the impact they may have in filtering the disease externality of trade. We use bovine viral diarrhoea (BVD) in England and Scotland as a case study, since this provides an example of a situation where contrasting strategies for BVD management occur between selling and purchasing farms. We use an agent-based bioeconomic model to assess the payoff dependence of farmers connected by trade but using different BVD management strategies. We compare three disease management actions: test-cull, test-cull with vaccination and vaccination alone. For a two-farm trading situation, all actions carried out by the selling farm provide substantial benefits to the purchasing farm in terms of disease avoided, with the greatest benefit resulting from test-culling with vaccination on the selling farm. Likewise, unilateral disease strategies by purchasers can be effective in reducing disease risks created through trade. We conclude that regulation needs to balance the trade-off between private gains from those bearing the disease management costs and the positive spillover effects on others.

The Need for Combining Implicit and Explicit Communication in Cooperative Robotic Systems.

As the number of robots used in warehouses and manufacturing increases, so too does the need for robots to be able to manipulate objects, not only independently, but also in collaboration with humans and other robots. Our ability to effectively coordinate our actions with fellow humans encompasses several behaviours that are collectively referred to as joint action, and has inspired advances in human-robot interaction by leveraging our natural ability to interpret implicit cues. However, our capacity to efficiently coordinate on object manipulation tasks remains an advantageous process that is yet to be fully exploited in robotic applications. Humans achieve this form of coordination by combining implicit communication (where information is inferred) and explicit communication (direct communication through an established channel) in varying degrees according to the task at hand. Although these two forms of communication have previously been implemented in robotic systems, no system exists that integrates the two in a task-dependent adaptive manner. In this paper, we review existing work on joint action in human-robot interaction, and analyse the state-of-the-art in robot-robot interaction that could act as a foundation for future cooperative object manipulation approaches. We identify key mechanisms that must be developed in order for robots to collaborate more effectively, with other robots and humans, on object manipulation tasks in shared autonomy spaces.

Adaptive Online Fault Diagnosis in Autonomous Robot Swarms.

Previous work has shown that robot swarms are not always tolerant to the failure of individual robots, particularly those that have only partially failed and continue to contribute to collective behaviors. A case has been made for an active approach to fault tolerance in swarm robotic systems, whereby the swarm can identify and resolve faults that occur during operation. Existing approaches to active fault tolerance in swarms have so far omitted fault diagnosis, however we propose that diagnosis is a feature of active fault tolerance that is necessary if swarms are to obtain long-term autonomy. This paper presents a novel method for fault diagnosis that attempts to imitate some of the observed functions of natural immune system. The results of our simulated experiments show that our system is flexible, scalable, and improves swarm tolerance to various electro-mechanical faults in the cases examined.

Tissue and host species-specific transcriptional changes in models of experimental visceral leishmaniasis.

Background: Human visceral leishmaniasis, caused by infection with Leishmania donovani or L. infantum, is a potentially fatal disease affecting 50,000-90,000 people yearly in 75 disease endemic countries, with more than 20,000 deaths reported. Experimental models of infection play a major role in understanding parasite biology, host-pathogen interaction, disease pathogenesis, and parasite transmission. In addition, they have an essential role in the identification and pre-clinical evaluation of new drugs and vaccines. However, our understanding of these models remains fragmentary. Although the immune response to Leishmania donovani infection in mice has been extensively characterized, transcriptomic analysis capturing the tissue-specific evolution of disease has yet to be reported. Methods: We provide an analysis of the transcriptome of spleen, liver and peripheral blood of BALB/c mice infected with L. donovani. Where possible, we compare our data in murine experimental visceral leishmaniasis with transcriptomic data in the public domain obtained from the study of L. donovani-infected hamsters and patients with human visceral leishmaniasis. Digitised whole slide images showing the histopathology in spleen and liver are made available via a dedicated website, www.leishpathnet.org. Results: Our analysis confirms marked tissue-specific alterations in the transcriptome of infected mice over time and identifies previously unrecognized parallels and differences between murine, hamster and human responses to infection. We show commonality of interferon-regulated genes whilst confirming a greater activation of type 2 immune pathways in infected hamsters compared to mice. Cytokine genes and genes encoding immune checkpoints were markedly tissue specific and dynamic in their expression, and pathways focused on non-immune cells reflected tissue specific immunopathology. Our data also addresses the value of measuring peripheral blood transcriptomics as a potential window into underlying systemic disease.  Conclusions: Our transcriptomic data, coupled with histopathologic analysis of the tissue response, provide an additional resource to underpin future mechanistic studies and to guide clinical research.

Generic, scalable and decentralized fault detection for robot swarms.

Robot swarms are large-scale multirobot systems with decentralized control which means that each robot acts based only on local perception and on local coordination with neighboring robots. The decentralized approach to control confers number of potential benefits. In particular, inherent scalability and robustness are often highlighted as key distinguishing features of robot swarms compared with systems that rely on traditional approaches to multirobot coordination. It has, however, been shown that swarm robotics systems are not always fault tolerant. To realize the robustness potential of robot swarms, it is thus essential to give systems the capacity to actively detect and accommodate faults. In this paper, we present a generic fault-detection system for robot swarms. We show how robots with limited and imperfect sensing capabilities are able to observe and classify the behavior of one another. In order to achieve this, the underlying classifier is an immune system-inspired algorithm that learns to distinguish between normal behavior and abnormal behavior online. Through a series of experiments, we systematically assess the performance of our approach in a detailed simulation environment. In particular, we analyze our system's capacity to correctly detect robots with faults, false positive rates, performance in a foraging task in which each robot exhibits a composite behavior, and performance under perturbations of the task environment. Results show that our generic fault-detection system is robust, that it is able to detect faults in a timely manner, and that it achieves a low false positive rate. The developed fault-detection system has the potential to enable long-term autonomy for robust multirobot systems, thus increasing the usefulness of robots for a diverse repertoire of upcoming applications in the area of distributed intelligent automation.

ASPASIA: A toolkit for evaluating the effects of biological interventions on SBML model behaviour.

A calibrated computational model reflects behaviours that are expected or observed in a complex system, providing a baseline upon which sensitivity analysis techniques can be used to analyse pathways that may impact model responses. However, calibration of a model where a behaviour depends on an intervention introduced after a defined time point is difficult, as model responses may be dependent on the conditions at the time the intervention is applied. We present ASPASIA (Automated Simulation Parameter Alteration and SensItivity Analysis), a cross-platform, open-source Java toolkit that addresses a key deficiency in software tools for understanding the impact an intervention has on system behaviour for models specified in Systems Biology Markup Language (SBML). ASPASIA can generate and modify models using SBML solver output as an initial parameter set, allowing interventions to be applied once a steady state has been reached. Additionally, multiple SBML models can be generated where a subset of parameter values are perturbed using local and global sensitivity analysis techniques, revealing the model's sensitivity to the intervention. To illustrate the capabilities of ASPASIA, we demonstrate how this tool has generated novel hypotheses regarding the mechanisms by which Th17-cell plasticity may be controlled in vivo. By using ASPASIA in conjunction with an SBML model of Th17-cell polarisation, we predict that promotion of the Th1-associated transcription factor T-bet, rather than inhibition of the Th17-associated transcription factor RORγt, is sufficient to drive switching of Th17 cells towards an IFN-γ-producing phenotype. Our approach can be applied to all SBML-encoded models to predict the effect that intervention strategies have on system behaviour. ASPASIA, released under the Artistic License (2.0), can be downloaded from http://www.york.ac.uk/ycil/software.

Extending and Applying Spartan to Perform Temporal Sensitivity Analyses for Predicting Changes in Influential Biological Pathways in Computational Models.

Through integrating real time imaging, computational modelling, and statistical analysis approaches, previous work has suggested that the induction of and response to cell adhesion factors is the key initiating pathway in early lymphoid tissue development, in contrast to the previously accepted view that the process is triggered by chemokine mediated cell recruitment. These model derived hypotheses were developed using spartan, an open-source sensitivity analysis toolkit designed to establish and understand the relationship between a computational model and the biological system that model captures. Here, we extend the functionality available in spartan to permit the production of statistical analyses that contrast the behavior exhibited by a computational model at various simulated time-points, enabling a temporal analysis that could suggest whether the influence of biological mechanisms changes over time. We exemplify this extended functionality by using the computational model of lymphoid tissue development as a time-lapse tool. By generating results at twelve- hour intervals, we show how the extensions to spartan have been used to suggest that lymphoid tissue development could be biphasic, and predict the time-point when a switch in the influence of biological mechanisms might occur.

Leukocyte Motility Models Assessed through Simulation and Multi-objective Optimization-Based Model Selection.

The advent of two-photon microscopy now reveals unprecedented, detailed spatio-temporal data on cellular motility and interactions in vivo. Understanding cellular motility patterns is key to gaining insight into the development and possible manipulation of the immune response. Computational simulation has become an established technique for understanding immune processes and evaluating hypotheses in the context of experimental data, and there is clear scope to integrate microscopy-informed motility dynamics. However, determining which motility model best reflects in vivo motility is non-trivial: 3D motility is an intricate process requiring several metrics to characterize. This complicates model selection and parameterization, which must be performed against several metrics simultaneously. Here we evaluate Brownian motion, Lévy walk and several correlated random walks (CRWs) against the motility dynamics of neutrophils and lymph node T cells under inflammatory conditions by simultaneously considering cellular translational and turn speeds, and meandering indices. Heterogeneous cells exhibiting a continuum of inherent translational speeds and directionalities comprise both datasets, a feature significantly improving capture of in vivo motility when simulated as a CRW. Furthermore, translational and turn speeds are inversely correlated, and the corresponding CRW simulation again improves capture of our in vivo data, albeit to a lesser extent. In contrast, Brownian motion poorly reflects our data. Lévy walk is competitive in capturing some aspects of neutrophil motility, but T cell directional persistence only, therein highlighting the importance of evaluating models against several motility metrics simultaneously. This we achieve through novel application of multi-objective optimization, wherein each model is independently implemented and then parameterized to identify optimal trade-offs in performance against each metric. The resultant Pareto fronts of optimal solutions are directly contrasted to identify models best capturing in vivo dynamics, a technique that can aid model selection more generally. Our technique robustly determines our cell populations' motility strategies, and paves the way for simulations that incorporate accurate immune cell motility dynamics.