Random migration and signal integration promote rapid and robust T cell recruitment

Interpreting T-cell search "strategies" in the light of evolution under constraints

Two decades of in vivo imaging have revealed how diverse T-cell motion patterns can be. Such recordings have sparked the notion of search "strategies": T cells may have evolved ways to search for antigen efficiently depending on the task at hand. Mathematical models have indeed confirmed that several observed T-cell migration patterns resemble a theoretical optimum; for example, frequent turning, stop-and-go motion, or alternating short and long motile runs have all been interpreted as deliberately tuned behaviours, optimising the cell's chance of finding antigen. But the same behaviours could also arise simply because T cells cannot follow a straight, regular path through the tight spaces they navigate. Even if T cells do follow a theoretically optimal pattern, the question remains: which parts of that pattern have truly been evolved for search, and which merely reflect constraints from the cell's migration machinery and surroundings? We here employ an approach from the field of evolutionary biology to examine how cells might evolve search strategies under realistic constraints. Using a cellular Potts model (CPM), where motion arises from intracellular dynamics interacting with cell shape and a constraining environment, we simulate evolutionary optimization of a simple task: explore as much area as possible. We find that our simulated cells indeed evolve their motility patterns. But the evolved behaviors are not shaped solely by what is functionally optimal; importantly, they also reflect mechanistic constraints. Cells in our model evolve several motility characteristics previously attributed to search optimisation-even though these features are not beneficial for the task given here. Our results stress that search patterns may evolve for other reasons than being "optimal". In part, they may be the inevitable side effects of interactions between cell shape, intracellular dynamics, and the diverse environments T cells face in vivo.

Adaptive tracking algorithm for trajectory analysis of cells and layer-by-layer assessment of motility dynamics

Tracking biological objects such as cells or subcellular components imaged with time-lapse microscopy enables us to understand the molecular principles about the dynamics of cell behaviors. However, automatic object detection, segmentation and extracting trajectories remain as a rate-limiting step due to intrinsic challenges of video processing. This paper presents an adaptive tracking algorithm (Adtari) that automatically finds the optimum search radius and cell linkages to determine trajectories in consecutive frames. A critical assumption in most tracking studies is that displacement remains unchanged throughout the movie and cells in a few frames are usually analyzed to determine its magnitude. Tracking errors and inaccurate association of cells may occur if the user does not correctly evaluate the value or prior knowledge is not present on cell movement. The key novelty of our method is that minimum intercellular distance and maximum displacement of cells between frames are dynamically computed and used to determine the threshold distance. Since the space between cells is highly variable in a given frame, our software recursively alters the magnitude to determine all plausible matches in the trajectory analysis. Our method therefore eliminates a major preprocessing step where a constant distance was used to determine the neighbor cells in tracking methods. Cells having multiple overlaps and splitting events were further evaluated by using the shape attributes including perimeter, area, ellipticity and distance. The features were applied to determine the closest matches by minimizing the difference in their magnitudes. Finally, reporting section of our software were used to generate instant maps by overlaying cell features and trajectories. Adtari was validated by using videos with variable signal-to-noise, contrast ratio and cell density. We compared the adaptive tracking with constant distance and other methods to evaluate performance and its efficiency. Our algorithm yields reduced mismatch ratio, increased ratio of whole cell track, higher frame tracking efficiency and allows layer-by-layer assessment of motility to characterize single-cells. Adaptive tracking provides a reliable, accurate, time efficient and user-friendly open source software that is well suited for analysis of 2D fluorescence microscopy video datasets.

Spatiotemporal behavior of T cells in vaccination

Vaccines are the most cost-effective resource to contain and eliminate infectious diseases. Despite decades of research in the field, several pathogens have eluded the effect of conventional vaccines mostly due their failure in inducing strong T cell responses. There is a need for new vaccine technologies that can surpass this problem. Recent advances in imaging techniques have allowed the study of T cell dynamics within their surrounding cellular niches. This information is invaluable to elucidate the main cellular mechanisms to target in order to optimize vaccine efficiency. In this review, we summarize the most recent key discoveries in T cell behavior in the context of vaccination and immunization.

Integrative lymph node-mimicking models created with biomaterials and computational tools to study the immune system

The lymph node (LN) is a vital organ of the lymphatic and immune system that enables timely detection, response, and clearance of harmful substances from the body. Each LN comprises of distinct substructures, which host a plethora of immune cell types working in tandem to coordinate complex innate and adaptive immune responses. An improved understanding of LN biology could facilitate treatment in LN-associated pathologies and immunotherapeutic interventions, yet at present, animal models, which often have poor physiological relevance, are the most popular experimental platforms. Emerging biomaterial engineering offers powerful alternatives, with the potential to circumvent limitations of animal models, for in-depth characterization and engineering of the lymphatic and adaptive immune system. In addition, mathematical and computational approaches, particularly in the current age of big data research, are reliable tools to verify and complement biomaterial works. In this review, we first discuss the importance of lymph node in immunity protection followed by recent advances using biomaterials to create in vitro/vivo LN-mimicking models to recreate the lymphoid tissue microstructure and microenvironment, as well as to describe the related immuno-functionality for biological investigation. We also explore the great potential of mathematical and computational models to serve as in silico supports. Furthermore, we suggest how both in vitro/vivo and in silico approaches can be integrated to strengthen basic patho-biological research, translational drug screening and clinical personalized therapies. We hope that this review will promote synergistic collaborations to accelerate progress of LN-mimicking systems to enhance understanding of immuno-complexity.

Local actin dynamics couple speed and persistence in a cellular Potts model of cell migration

Cell migration is astoundingly diverse. Molecular signatures, cell-cell interactions, and environmental structures each play their part in shaping cell motion, yielding numerous morphologies and migration modes. Nevertheless, in recent years, a simple unifying law was found to describe cell migration across many different cell types and contexts: faster cells turn less frequently. This universal coupling between speed and persistence (UCSP) was explained by retrograde actin flow from front to back, but it remains unclear how this mechanism generalizes to cells with complex shapes and cells migrating in structured environments, which may not have a well-defined front-to-back orientation. Here, we present an in-depth characterization of an existing cellular Potts model, in which cells polarize dynamically from a combination of local actin dynamics (stimulating protrusions) and global membrane tension along the perimeter (inhibiting protrusions). We first show that the UCSP emerges spontaneously in this model through a cross talk of intracellular mechanisms, cell shape, and environmental constraints, resembling the dynamic nature of cell migration in vivo. Importantly, we find that local protrusion dynamics suffice to reproduce the UCSP-even in cases in which no clear global, front-to-back polarity exists. We then harness the spatial nature of the cellular Potts model to show how cell shape dynamics limit both the speed and persistence a cell can reach and how a rigid environment such as the skin can restrict cell motility even further. Our results broaden the range of potential mechanisms underlying the speed-persistence coupling that has emerged as a fundamental property of migrating cells.

Modeling the Dynamics of T-Cell Development in the Thymus

The thymus hosts the development of a specific type of adaptive immune cells called T cells. T cells orchestrate the adaptive immune response through recognition of antigen by the highly variable T-cell receptor (TCR). T-cell development is a tightly coordinated process comprising lineage commitment, somatic recombination of Tcr gene loci and selection for functional, but non-self-reactive TCRs, all interspersed with massive proliferation and cell death. Thus, the thymus produces a pool of T cells throughout life capable of responding to virtually any exogenous attack while preserving the body through self-tolerance. The thymus has been of considerable interest to both immunologists and theoretical biologists due to its multi-scale quantitative properties, bridging molecular binding, population dynamics and polyclonal repertoire specificity. Here, we review experimental strategies aimed at revealing quantitative and dynamic properties of T-cell development and how they have been implemented in mathematical modeling strategies that were reported to help understand the flexible dynamics of the highly dividing and dying thymic cell populations. Furthermore, we summarize the current challenges to estimating in vivo cellular dynamics and to reaching a next-generation multi-scale picture of T-cell development.

The regulation of DOCK family proteins on T and B cells

The dedicator of cytokinesis (DOCK) family proteins consist of 11 members, each of which contains 2 domains, DOCK homology region (DHR)-1 and DHR-2, and as guanine nucleotide exchange factors, they mediate activation of small GTPases. Both DOCK2 and DOCK8 deficiencies in humans can cause severe combined immunodeficiency, but they have different characteristics. DOCK8 defect mainly causes high IgE, allergic disease, refractory skin virus infection, and increased incidence of malignant tumor, whereas DOCK2 defect mainly causes early-onset, invasive infection with less atopy and increased IgE. However, the underlying molecular mechanisms causing the disease remain unclear. This paper discusses the role of DOCK family proteins in regulating B and T cells, including development, survival, migration, activation, immune tolerance, and immune functions. Moreover, related signal pathways or molecule mechanisms are also described in this review. A greater understanding of DOCK family proteins and their regulation of lymphocyte functions may facilitate the development of new therapeutics for immunodeficient patients and improve their prognosis.

Manipulating the TCR signaling network for cellular immunotherapy: Challenges & opportunities

T cells can help confer protective immunity by eliminating infections and tumors or drive immunopathology by damaging host cells. Both outcomes require a series of steps from the activation of naïve T cells to their clonal expansion, differentiation and migration to tissue sites. In addition to specific recognition of the antigen via the T cell receptor (TCR), multiple accessory signals from costimulatory molecules, cytokines and metabolites also influence each step along the progression of the T cell response. Current efforts to modify effector T cell function in many clinical contexts focus on the latter - which encompass antigen-independent and broad, contextual regulators. Not surprisingly, such approaches are often accompanied by adverse events, as they also affect T cells not relevant to the specific treatment. In contrast, fine tuning T cell responses by precisely targeting antigen-specific TCR signals has the potential to radically alter therapeutic strategies in a focused manner. Development of such approaches, however, requires a better understanding of functioning of the TCR and the biochemical signaling network coupled to it. In this article, we review some of the recent advances which highlight important roles of TCR signals throughout the activation and differentiation of T cells during an immune response. We discuss how, an appreciation of specific signaling modalities and variant ligands that influence the function of the TCR has the potential to influence design principles for the next generation of pharmacologic and cellular therapies, especially in the context of tumor immunotherapies involving adoptive cell transfers.

Estimating Residence Times of Lymphocytes in Ovine Lymph Nodes

The ability of lymphocytes to recirculate between blood and secondary lymphoid tissues such as lymph nodes (LNs) and spleen is well established. Sheep have been used as an experimental system to study lymphocyte recirculation for decades and multiple studies document accumulation and loss of intravenously (i.v.) transferred lymphocytes in efferent lymph of various ovine LNs. Yet, surprisingly little work has been done to accurately quantify the dynamics of lymphocyte exit from the LNs and to estimate the average residence times of lymphocytes in ovine LNs. In this work we developed a series of mathematical models based on fundamental principles of lymphocyte recirculation in the body under non-inflammatory (resting) conditions. Our analysis suggested that in sheep, recirculating lymphocytes spend on average 3 h in the spleen and 20 h in skin or gut-draining LNs with a distribution of residence times in LNs following a skewed gamma (lognormal-like) distribution. Our mathematical models also suggested an explanation for a puzzling observation of the long-term persistence of i.v. transferred lymphocytes in the efferent lymph of the prescapular LN (pLN); the model predicted that this is a natural consequence of long-term persistence of the transferred lymphocytes in circulation. We also found that lymphocytes isolated from the skin-draining pLN have a 2-fold increased entry rate into the pLN as opposed to the mesenteric (gut-draining) LN (mLN). Likewise, lymphocytes from mLN had a 3-fold increased entry rate into the mLN as opposed to entry rate into pLN. In contrast, these cannulation data could not be explained by preferential retention of cells in LNs of their origin. Taken together, our work illustrates the power of mathematical modeling in describing the kinetics of lymphocyte migration in sheep and provides quantitative estimates of lymphocyte residence times in ovine LNs.

Sleep Matters: CD4+ T Cell Memory Formation and the Central Nervous System

The mechanisms of CD4⁺ T-cell memory formation in the immune system are debated. With the well-established concept of memory formation in the central nervous system (CNS), we propose that formation of CD4⁺ T-cell memory depends on the interaction of two different cell systems handling two types of stored information. First, information about antigen (event) and challenge (context) is taken up by antigen-presenting cells, as initial storage. Second, event and context information is transferred to CD4⁺ T cells. During activation, two categories of CD4⁺ T cell develop: effector CD4⁺ T cells, carrying event and context information, enabling them to efficiently focus their response to tissues under attack; and persisting CD4⁺ T cells, providing context-independent antigen-specific memories and long-term storage. This novel hypothesis is supported by the observation that mammalian sleep can improve both CNS and CD4⁺ T-cell memory.

Virus dynamics and phyloanatomy: Merging population dynamic and phylogenetic approaches

In evolutionary biology and epidemiology, phylodynamic methods are widely used to infer population biological characteristics, such as the rates of replication, death, migration, or, in the epidemiological context, pathogen spread. More recently, these methods have been used to elucidate the dynamics of viruses within their hosts. Especially the application of phylogeographic approaches has the potential to shed light on anatomical colonization pathways and the exchange of viruses between distinct anatomical compartments. We and others have termed this phyloanatomy. Here, we review the promise and challenges of phyloanatomy, and compare them to more classical virus dynamics and population genetic approaches. We argue that the extremely strong selection pressures that exist within the host may represent the main obstacle to reliable phyloanatomic analysis.

Distributed Adaptive Search in T Cells: Lessons From Ants

There are striking similarities between the strategies ant colonies use to forage for food and immune systems use to search for pathogens. Searchers (ants and cells) use the appropriate combination of random and directed motion, direct and indirect agent-agent interactions, and traversal of physical structures to solve search problems in a variety of environments. An effective immune response requires immune cells to search efficiently and effectively for diverse types of pathogens in different tissues and organs, just as different species of ants have evolved diverse search strategies to forage effectively for a variety of resources in a variety of habitats. Successful T cell search is required to initiate the adaptive immune response in lymph nodes and to eradicate pathogens at sites of infection in peripheral tissue. Ant search strategies suggest novel predictions about T cell search. In both systems, the distribution of targets in time and space determines the most effective search strategy. We hypothesize that the ability of searchers to sense and adapt to dynamic targets and environmental conditions enhances search effectiveness through adjustments to movement and communication patterns. We also suggest that random motion is a more important component of search strategies than is generally recognized. The behavior we observe in ants reveals general design principles and constraints that govern distributed adaptive search in a wide variety of complex systems, particularly the immune system.

Spatial Lymphocyte Dynamics in Lymph Nodes Predicts the Cytotoxic T Cell Frequency Needed for HIV Infection Control

The surveillance of host body tissues by immune cells is central for mediating their defense function. In vivo imaging technologies have been used to quantitatively characterize target cell scanning and migration of lymphocytes within lymph nodes (LNs). The translation of these quantitative insights into a predictive understanding of immune system functioning in response to various perturbations critically depends on computational tools linking the individual immune cell properties with the emergent behavior of the immune system. By choosing the Newtonian second law for the governing equations, we developed a broadly applicable mathematical model linking individual and coordinated T-cell behaviors. The spatial cell dynamics is described by a superposition of autonomous locomotion, intercellular interaction, and viscous damping processes. The model is calibrated using in vivo data on T-cell motility metrics in LNs such as the translational speeds, turning angle speeds, and meandering indices. The model is applied to predict the impact of T-cell motility on protection against HIV infection, i.e., to estimate the threshold frequency of HIV-specific cytotoxic T cells (CTLs) that is required to detect productively infected cells before the release of viral particles starts. With this, it provides guidance for HIV vaccine studies allowing for the migration of cells in fibrotic LNs.

Integrative Computational Modeling of the Lymph Node Stromal Cell Landscape

Adaptive immune responses develop in secondary lymphoid organs such as lymph nodes (LNs) in a well-coordinated series of interactions between migrating immune cells and resident stromal cells. Although many processes that occur in LNs are well understood from an immunological point of view, our understanding of the fundamental organization and mechanisms that drive these processes is still incomplete. The aim of systems biology approaches is to unravel the complexity of biological systems and describe emergent properties that arise from interactions between individual constituents of the system. The immune system is greater than the sum of its parts, as is the case with any sufficiently complex system. Here, we review recent work and developments of computational LN models with focus on the structure and organization of the stromal cells. We explore various mathematical studies of intranodal T cell motility and migration, their interactions with the LN-resident stromal cells, and computational models of functional chemokine gradient fields and lymph flow dynamics. Lastly, we discuss briefly the importance of hybrid and multi-scale modeling approaches in immunology and the technical challenges involved.

Deep Sequencing Reveals Transient Segregation of T Cell Repertoires in Splenic T Cell Zones during an Immune Response

Immunological differences between hosts, such as diverse TCR repertoires, are widely credited for reducing the risk of pathogen spread and adaptation in a population. Within-host immunological diversity might likewise be important for robust pathogen control, but to what extent naive TCR repertoires differ across different locations in the same host is unclear. T cell zones (TCZs) in secondary lymphoid organs provide secluded microenvironmental niches. By harboring distinct TCRs, such niches could enhance within-host immunological diversity. In contrast, rapid T cell migration is expected to dilute such diversity. In this study, we combined tissue microdissection and deep sequencing of the TCR β-chain to examine the extent to which TCR repertoires differ between TCZs in murine spleens. In the absence of Ag, we found little evidence for differences between TCZs of the same spleen. Yet, 3 d after immunization with sheep RBCs, we observed a >10-fold rise in the number of clones that appeared to localize to individual zones. Remarkably, these differences largely disappeared at 4 d after immunization, when hallmarks of an ongoing immune response were still observed. These data suggest that in the absence of Ag, any repertoire differences observed between TCZs of the same host can largely be attributed to random clone distribution. Upon Ag challenge, TCR repertoires in TCZs first segregate and then homogenize within days. Such "transient mosaic" dynamics could be an important barrier for pathogen adaptation and spread during an immune response.

Switching between individual and collective motility in B lymphocytes is controlled by cell-matrix adhesion and inter-cellular interactions

Lymphocytes alternate between phases of individual migration across tissues and phases of clustering during activation and function. The range of lymphocyte motility behaviors and the identity of the factors that govern them remain elusive. To explore this point, we here collected unprecedented statistics pertaining to cell displacements, cell:matrix and cell:cell interactions using a model B cell line as well as primary human B lymphocytes. At low cell density, individual B lymphocytes displayed a high heterogeneity in their speed and diffusivity. Beyond this intrinsic variability, B lymphocytes adapted their motility to the composition of extra-cellular matrix, adopting slow persistent walks over collagen IV and quick Brownian walks over fibronectin. At high cell density, collagen IV favored the self-assembly of B lymphocytes into clusters endowed with collective coordination, while fibronectin stimulated individual motility. We show that this behavioral plasticity is controlled by acto-myosin dependent adhesive and Arp2/3-dependent protrusive actin pools, respectively. Our study reveals the adaptive nature of B lymphocyte motility and group dynamics, which are shaped by an interplay between and cell:matrix and cell:cell interactions.

How Germinal Centers Evolve Broadly Neutralizing Antibodies: the Breadth of the Follicular Helper T Cell Response

Many HIV-1-infected patients evolve broadly neutralizing antibodies (bnAbs). This evolutionary process typically takes several years and is poorly understood as selection taking place in germinal centers occurs on the basis of antibody affinity. B cells with the highest-affinity receptors tend to acquire the most antigen from the follicular dendritic cell (FDC) network and present the highest density of cognate peptides to follicular helper T (Tfh) cells, which provide survival signals to the B cell. bnAbs are therefore expected to evolve only when the B cell lineage evolving breadth is consistently capturing and presenting more peptides to Tfh cells than other lineages of more specific B cells. Here we develop mathematical models of Tfh cells in germinal centers to explicitly define the mechanisms of selection in this complex evolutionary process. Our results suggest that broadly reactive B cells presenting a high density of peptides bound to major histocompatibility complex class II molecules (pMHC) are readily outcompeted by B cells responding to lineages of HIV-1 that transiently dominate the within host viral population. Conversely, if broadly reactive B cells acquire a large variety of several HIV-1 proteins from the FDC network and present a high diversity of several pMHC, they can be rescued by a large fraction of the Tfh cell repertoire in the germinal center. Under such circumstances the evolution of bnAbs is much more consistent. Increasing either the magnitude of the Tfh cell response or the breadth of the Tfh cell repertoire markedly facilitates the evolution of bnAbs. Because both the magnitude and breadth can be increased by vaccination with several HIV-1 proteins, this calls for experimental testing.IMPORTANCE Many HIV-infected patients slowly evolve antibodies that can neutralize a large variety of viruses. Such broadly neutralizing antibodies (bnAbs) could in the future become therapeutic agents. bnAbs appear very late, and patients are typically not protected by them. At the moment, we fail to understand why this takes so long and how the immune system selects for broadly neutralizing capacity. Typically, antibodies are selected based on affinity and not on breadth. We developed mathematical models to study two different mechanisms by which the immune system can select for broadly neutralizing capacity. One of these is based upon the repertoire of different follicular helper T (Tfh) cells in germinal centers. We suggest that broadly reactive B cells may interact with a larger fraction of this repertoire and demonstrate that this would select for bnAbs. Intriguingly, this suggests that broadening the Tfh cell repertoire by vaccination may speed up the evolution of bnAbs.

Insights from mathematical modelling for T cell migration into the central nervous system

The migration of immune cells from peripheral immune organs into the central nervous system (CNS) through the blood-brain barrier (BBB) is a tightly regulated process. The complex interplay between cells of the BBB and immune cells coordinates cell migration as a part of normal immune surveillance while its dysregulation is critically involved in the pathogenesis of various CNS diseases. To develop tools for a deeper understanding of distribution and migratory pattern of immune cells regulated by the BBB, we made use of a mathematical modelling approach derived from Markov chain theory. We present a data-driven model using a derivation of kinetic differential equations from a particle game. According to the theory of gases, these equations allow one to predict the mean behaviour of a large class of cells by modelling cell-cell interactions. We used this model to assess the distribution of naive, central memory and effector memory T lymphocytes in the peripheral blood and cerebrospinal fluid. Our model allows us to evaluate the impact of activation status, migratory capacity and cell death for cell distribution in the peripheral blood and the CNS.

Imaging the pharmacology of nanomaterials by intravital microscopy: Toward understanding their biological behavior

Therapeutic nanoparticles (NPs) can deliver cytotoxic chemotherapeutics and other drugs more safely and efficiently to patients; furthermore, selective delivery to target tissues can theoretically be accomplished actively through coating NPs with molecular ligands, and passively through exploiting physiological "enhanced permeability and retention" features. However, clinical trial results have been mixed in showing improved efficacy with drug nanoencapsulation, largely due to heterogeneous NP accumulation at target sites across patients. Thus, a clear need exists to better understand why many NP strategies fail in vivo and not result in significantly improved tumor uptake or therapeutic response. Multicolor in vivo confocal fluorescence imaging (intravital microscopy; IVM) enables integrated pharmacokinetic and pharmacodynamic (PK/PD) measurement at the single-cell level, and has helped answer key questions regarding the biological mechanisms of in vivo NP behavior. This review summarizes progress to date and also describes useful technical strategies for successful IVM experimentation.

Growth of Murine Splenic Tissue Is Suppressed by Lymphotoxin β-Receptor Signaling (LTβR) Originating from Splenic and Non-Splenic Tissues

Development and maintenance of secondary lymphoid organs such as lymph nodes and spleen essentially depend on lymphotoxin β-receptor (LTβR) signaling. It is unclear, however, by which molecular mechanism their size is limited. Here, we investigate whether the LTβR pathway is also growth suppressing. By using splenic tissue transplantation it is possible to analyze a potential contribution of LTβR signaling inside and outside of the implanted tissue. We show that LTβR signaling within the endogenous spleen and within non-splenic tissues both significantly suppressed the regeneration of implanted splenic tissue. The suppressive activity positively correlated with the total number of LTβR expressing cells in the animal (regenerate weights of 115 ± 8 mg in LTβR deficient recipients and of 12 ± 9 mg in wild-type recipients), affected also developed splenic tissue, and was induced but not executed via LTβR signaling. Two-dimensional differential gel electrophoresis and subsequent mass spectrometry of stromal splenic tissue was applied to screen for potential factors mediating the LTβR dependent suppressive activity. Thus, LTβR dependent growth suppression is involved in regulating the size of secondary lymphoid organs, and might be therapeutically used to eradicate tertiary lymphoid tissues during autoimmune diseases.

Large Variations in HIV-1 Viral Load Explained by Shifting-Mosaic Metapopulation Dynamics

The viral population of HIV-1, like many pathogens that cause systemic infection, is structured and differentiated within the body. The dynamics of cellular immune trafficking through the blood and within compartments of the body has also received wide attention. Despite these advances, mathematical models, which are widely used to interpret and predict viral and immune dynamics in infection, typically treat the infected host as a well-mixed homogeneous environment. Here, we present mathematical, analytical, and computational results that demonstrate that consideration of the spatial structure of the viral population within the host radically alters predictions of previous models. We study the dynamics of virus replication and cytotoxic T lymphocytes (CTLs) within a metapopulation of spatially segregated patches, representing T cell areas connected by circulating blood and lymph. The dynamics of the system depend critically on the interaction between CTLs and infected cells at the within-patch level. We show that for a wide range of parameters, the system admits an unexpected outcome called the shifting-mosaic steady state. In this state, the whole body’s viral population is stable over time, but the equilibrium results from an underlying, highly dynamic process of local infection and clearance within T-cell centers. Notably, and in contrast to previous models, this new model can explain the large differences in set-point viral load (SPVL) observed between patients and their distribution, as well as the relatively low proportion of cells infected at any one time, and alters the predicted determinants of viral load variation.

A novel metapopulation model of HIV suggests that within-host infections are characterized by a highly dynamic process of localized infection followed by clearance within T cell centers.

When a person is infected with HIV, the initial peak level of virus in the blood is usually very high before a lower, relatively stable level is reached and maintained for the duration of the chronic infection. This stable level is known as the set-point viral load (SPVL) and is associated with severity of infection. SPVL is also highly variable among patients, ranging from 100 to a million copies of the virus per mL of blood. The replicative capacity of the infecting virus and the strength of the immune response both influence SPVL. However, standard mathematical models show that variation in these two factors cannot easily reproduce the observed distribution of SPVL among patients. Standard models typically treat infected individuals as well-mixed systems, but in reality viral replication is localised in T-cell centres, or patches, found in secondary lymphoid tissue. To account for this population structure, we developed a carefully parameterised metapopulation model. We find the system can reach a steady state at which the viral load in the blood is relatively stable, representing SPVL, but surprisingly, the patches are highly dynamic, characterised by bursts of infection followed by elimination of virus due to localised host immune responses. Significantly, this model can reproduce the wide distribution of SPVLs found among infected individuals for realistic distributions of viral replicative capacity and strength of immune response. Our model can also be used in the future to understand other aspects of chronic HIV infection.

Immune-inspired search strategies for robot swarms

Detection of targets distributed randomly in space is a task common to both robotic and biological systems. Lévy search has previously been used to characterize T cell search in the immune system. We use a robot swarm to evaluate the effectiveness of a Lévy search strategy and map the relationship between search parameters and target configurations. We show that the fractal dimension of the Lévy search which optimizes search efficiency depends strongly on the distribution of targets but only weakly on the number of agents involved in search. Lévy search can therefore be tuned to the target configuration while also being scalable. Implementing search behaviors observed in T cells in a robot swarm provides an effective, adaptable, and scalable swarm robotic search strategy. Additionally, the adaptability and scalability of Lévy search may explain why Lévy-like movement has been observed in T cells in multiple immunological contexts.

A virtual lymph node model to dissect the requirements for T-cell activation by synapses and kinapses

The initiation of T-cell responses in lymph nodes requires T cells to integrate signals delivered by dendritic cells (DCs) during long-lasting contacts (synapses) or more transient interactions (kinapses). However, it remains extremely challenging to understand how a specific sequence of contacts established by T cells ultimately dictates T-cell fate. Here, we have coupled a computational model of T-cell migration and interactions with DCs with a real-time, flow cytometry-like representation of T-cell activation. In this model, low-affinity peptides trigger T-cell proliferation through kinapses but we show that this process is only effective under conditions of high DC densities and prolonged antigen availability. By contrast, high-affinity peptides favor synapse formation and a vigorous proliferation under a wide range of antigen presentation conditions. In line with the predictions, decreasing the DC density in vivo selectively abolished proliferation induced by the low-affinity peptide. Finally, our results suggest that T cells possess a biochemical memory of previous stimulations of at least 1–2 days. We propose that the stability of T-cell–DC interactions, apart from their signaling potency, profoundly influences the robustness of T-cell activation. By offering the ability to control parameters that are difficult to manipulate experimentally, the virtual lymph node model provides new possibilities to tackle the fundamental mechanisms that regulate T-cell responses elicited by infections or vaccines.

Mathematics in modern immunology

Mathematical and statistical methods enable multidisciplinary approaches that catalyse discovery. Together with experimental methods, they identify key hypotheses, define measurable observables and reconcile disparate results. We collect a representative sample of studies in T-cell biology that illustrate the benefits of modelling–experimental collaborations and that have proven valuable or even groundbreaking. We conclude that it is possible to find excellent examples of synergy between mathematical modelling and experiment in immunology, which have brought significant insight that would not be available without these collaborations, but that much remains to be discovered.

Persistence and Adaptation in Immunity: T Cells Balance the Extent and Thoroughness of Search

Effective search strategies have evolved in many biological systems, including the immune system. T cells are key effectors of the immune response, required for clearance of pathogenic infection. T cell activation requires that T cells encounter antigen-bearing dendritic cells within lymph nodes, thus, T cell search patterns within lymph nodes may be a crucial determinant of how quickly a T cell immune response can be initiated. Previous work suggests that T cell motion in the lymph node is similar to a Brownian random walk, however, no detailed analysis has definitively shown whether T cell movement is consistent with Brownian motion. Here, we provide a precise description of T cell motility in lymph nodes and a computational model that demonstrates how motility impacts T cell search efficiency. We find that both Brownian and Lévy walks fail to capture the complexity of T cell motion. Instead, T cell movement is better described as a correlated random walk with a heavy-tailed distribution of step lengths. Using computer simulations, we identify three distinct factors that contribute to increasing T cell search efficiency: 1) a lognormal distribution of step lengths, 2) motion that is directionally persistent over short time scales, and 3) heterogeneity in movement patterns. Furthermore, we show that T cells move differently in specific frequently visited locations that we call "hotspots" within lymph nodes, suggesting that T cells change their movement in response to the lymph node environment. Our results show that like foraging animals, T cells adapt to environmental cues, suggesting that adaption is a fundamental feature of biological search.

A spatial model of the efficiency of T cell search in the influenza-infected lung

Emerging strains of influenza, such as avian H5N1 and 2009 pandemic H1N1, are more virulent than seasonal H1N1 influenza, yet the underlying mechanisms for these differences are not well understood. Subtle differences in how a given strain interacts with the immune system are likely a key factor in determining virulence. One aspect of the interaction is the ability of T cells to locate the foci of the infection in time to prevent uncontrolled expansion. Here, we develop an agent based spatial model to focus on T cell migration from lymph nodes through the vascular system to sites of infection. We use our model to investigate whether different strains of influenza modulate this process. We calibrate the model using viral and chemokine secretion rates we measure in vitro together with values taken from literature. The spatial nature of the model reveals unique challenges for T cell recruitment that are not apparent in standard differential equation models. In this model comparing three influenza viruses, plaque expansion is governed primarily by the replication rate of the virus strain, and the efficiency of the T cell search-and-kill is limited by the density of infected epithelial cells in each plaque. Thus for each virus there is a different threshold of T cell search time above which recruited T cells are unable to control further expansion. Future models could use this relationship to more accurately predict control of the infection.

T cell and reticular network co-dependence in HIV infection

Fibroblastic reticular cells (FRC) are arranged on a network in the T cell zone of lymph nodes, forming a scaffold for T cell migration, and providing survival factors, especially interleukin-7 (IL-7). Conversely, CD4(+) T cells are the major producers of lymphotoxin-β (LT-β), necessary for the construction and maintenance of the FRC network. This interdependence creates the possibility of a vicious cycle, perpetuating loss of both FRC and T cells. Furthermore, evidence that HIV infection is responsible for collagenation of the network suggests that long term loss of network function might be responsible for the attenuated recovery in T cell count seen in HIV patients undergoing antiretroviral therapy (ART). We present computational and mathematical models of this interaction mechanism and subsequent naive CD4(+) T-cell depletion in which (1) collagen deposition impedes access of naive T cells to IL-7 on the FRC and loss of IL-7 production by loss of FRC network itself, leading to the depletion of naive T cells through increased apoptosis; and (2) depletion of naive T cells as the source of LT-β on which the FRC depend for survival leads to loss of the network, thereby amplifying and perpetuating the cycle of depletion of both naive T cells and stromal cells. Our computational model explicitly includes an FRC network and its cytokine exchange with a heterogeneous T-cell population. We also derive lumped models, in terms of partial differential equations and reduced to ordinary differential equations, that provide additional insight into the mechanisms at work. The central conclusions are that (1) damage to the reticular network, caused by HIV infection is a plausible mechanism for attenuated recovery post-ART; (2) within this, the production of T cell survival factors by FRCs may be the key rate-limiting step; and (3) the methods of model reduction and analysis presented are useful for both immunological studies and other contexts in which agent-based models are severely limited by computational cost.

Dynamics of HIV infection in lymphoid tissue network

Human immunodeficiency virus (HIV) is a fast replicating ribonucleic acid virus, which can easily mutate in order to escape the effects of drug administration. Hence, understanding the basic mechanisms underlying HIV persistence in the body is essential in the development of new therapies that could eradicate HIV infection. Lymphoid tissues are the primary sites of HIV infection. Despite the recent progress in real-time monitoring technology, HIV infection dynamics in a whole body is unknown. Mathematical modeling and simulations provide speculations on global behavior of HIV infection in the lymphatic system. We propose a new mathematical model that describes the spread of HIV infection throughout the lymphoid tissue network. In order to represent the volume difference between lymphoid tissues, we propose the proportionality of several kinetic parameters to the lymphoid tissues' volume distribution. Under this assumption, we perform extensive numerical computations in order to simulate the spread of HIV infection in the lymphoid tissue network. Numerical computations simulate single drug treatments of an HIV infection. One of the important biological speculations derived from this study is a drug saturation effect generated by lymphoid network connection. This implies that a portion of reservoir lymphoid tissues to which drug is not sufficiently delivered would inhibit HIV eradication despite of extensive drug injection.

Age-Dependent Cell Trafficking Defects in Draining Lymph Nodes Impair Adaptive Immunity and Control of West Nile Virus Infection

Impaired immune responses in the elderly lead to reduced vaccine efficacy and increased susceptibility to viral infections. Although several groups have documented age-dependent defects in adaptive immune priming, the deficits that occur prior to antigen encounter remain largely unexplored. Herein, we identify novel mechanisms for compromised adaptive immunity that occurs with aging in the context of infection with West Nile virus (WNV), an encephalitic flavivirus that preferentially causes disease in the elderly. An impaired IgM and IgG response and enhanced vulnerability to WNV infection during aging was linked to delayed germinal center formation in the draining lymph node (DLN). Adoptive transfer studies and two-photon intravital microscopy revealed a decreased trafficking capacity of donor naïve CD4+ T cells from old mice, which manifested as impaired T cell diapedesis at high endothelial venules and reduced cell motility within DLN prior to antigen encounter. Furthermore, leukocyte accumulation in the DLN within the first few days of WNV infection or antigen-adjuvant administration was diminished more generally in old mice and associated with a second aging-related defect in local cytokine and chemokine production. Thus, age-dependent cell-intrinsic and environmental defects in the DLN result in delayed immune cell recruitment and antigen recognition. These deficits compromise priming of early adaptive immune responses and likely contribute to the susceptibility of old animals to acute WNV infection.

Visualization and dynamic analysis of host-pathogen interactions

To contain invading microbes, the immune system must efficiently recognize the presence of the invader, mobilize cells to the site of infection, and deploy effector function. Rare antigen-specific T cells must find small numbers of antigen-presenting cells, proliferate and differentiate in secondary lymphoid tissues, then traffic to the infected site and be activated by antigen again to contribute to host defense. Our understanding of the dynamic processes involved has benefited enormously from tools that enable the visualization of cell location and behavior in complex tissue environments. Here we summarize recent insights into T cell trafficking and migration through secondary lymphoid organs and at peripheral infection sites, highlighting cell-intrinsic and extrinsic factors optimizing antigen surveillance at steady-state and delivery of an effector response during infection.