Spatially distributed infection increases viral load in a computational model of SARS-CoV-2 lung infection.

A key question in SARS-CoV-2 infection is why viral loads and patient outcomes vary dramatically across individuals. Because spatial-temporal dynamics of viral spread and immune response are challenging to study in vivo, we developed Spatial Immune Model of Coronavirus (SIMCoV), a scalable computational model that simulates hundreds of millions of lung cells, including respiratory epithelial cells and T cells. SIMCoV replicates viral growth dynamics observed in patients and shows how spatially dispersed infections can lead to increased viral loads. The model also shows how the timing and strength of the T cell response can affect viral persistence, oscillations, and control. By incorporating spatial interactions, SIMCoV provides a parsimonious explanation for the dramatically different viral load trajectories among patients by varying only the number of initial sites of infection and the magnitude and timing of the T cell immune response. When the branching airway structure of the lung is explicitly represented, we find that virus spreads faster than in a 2D layer of epithelial cells, but much more slowly than in an undifferentiated 3D grid or in a well-mixed differential equation model. These results illustrate how realistic, spatially explicit computational models can improve understanding of within-host dynamics of SARS-CoV-2 infection.

CT Colonography: Improving Interpretive Skill by Avoiding Pitfalls.

An earlier incorrect version of this article appeared online. This article was corrected on December 20, 2019.

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.

More Is Faster: Why Population Size Matters in Biological Search.

Many biological scenarios have multiple cooperating searchers, and the timing of the initial first contact between any one of those searchers and its target is critically important. However, we are unaware of biological models that predict how long it takes for the first of many searchers to discover a target. We present a novel mathematical model that predicts initial first contact times between searchers and targets distributed at random in a volume. We compare this model with the extreme first passage time approach in physics that assumes an infinite number of searchers all initially positioned at the same location. We explore how the number of searchers, the distribution of searchers and targets, and the initial distances between searchers and targets affect initial first contact times. Given a constant density of uniformly distributed searchers and targets, the initial first contact time decreases linearly with both search volume and the number of searchers. However, given only a single target and searchers placed at the same starting location, the relationship between the initial first contact time and the number of searchers shifts from a linear decrease to a logarithmic decrease as the number of searchers grows very large. More generally, we show that initial first contact times can be dramatically faster than the average first contact times and that the initial first contact times decrease with the number of searchers, while the average search times are independent of the number of searchers. We suggest that this is an underappreciated phenomenon in biology and other collective search problems.

Shifts in metabolic scaling, production, and efficiency across major evolutionary transitions of life.

The diversification of life involved enormous increases in size and complexity. The evolutionary transitions from prokaryotes to unicellular eukaryotes to metazoans were accompanied by major innovations in metabolic design. Here we show that the scalings of metabolic rate, population growth rate, and production efficiency with body size have changed across the evolutionary transitions. Metabolic rate scales with body mass superlinearly in prokaryotes, linearly in protists, and sublinearly in metazoans, so Kleiber's 3/4 power scaling law does not apply universally across organisms. The scaling of maximum population growth rate shifts from positive in prokaryotes to negative in protists and metazoans, and the efficiency of production declines across these groups. Major changes in metabolic processes during the early evolution of life overcame existing constraints, exploited new opportunities, and imposed new constraints.

A general basis for quarter-power scaling in animals.

It has been known for decades that the metabolic rate of animals scales with body mass with an exponent that is almost always <1, >2/3, and often very close to 3/4. The 3/4 exponent emerges naturally from two models of resource distribution networks, radial explosion and hierarchically branched, which incorporate a minimum of specific details. Both models show that the exponent is 2/3 if velocity of flow remains constant, but can attain a maximum value of 3/4 if velocity scales with its maximum exponent, 1/12. Quarter-power scaling can arise even when there is no underlying fractality. The canonical "fourth dimension" in biological scaling relations can result from matching the velocity of flow through the network to the linear dimension of the terminal "service volume" where resources are consumed. These models have broad applicability for the optimal design of biological and engineered systems where energy, materials, or information are distributed from a single source.

Quantitative analyses of T cell motion in tissue reveals factors driving T cell search in tissues.

T cells are required to clear infection, and T cell motion plays a role in how quickly a T cell finds its target, from initial naive T cell activation by a dendritic cell to interaction with target cells in infected tissue. To better understand how different tissue environments affect T cell motility, we compared multiple features of T cell motion including speed, persistence, turning angle, directionality, and confinement of T cells moving in multiple murine tissues using microscopy. We quantitatively analyzed naive T cell motility within the lymph node and compared motility parameters with activated CD8 T cells moving within the villi of small intestine and lung under different activation conditions. Our motility analysis found that while the speeds and the overall displacement of T cells vary within all tissues analyzed, T cells in all tissues tended to persist at the same speed. Interestingly, we found that T cells in the lung show a marked population of T cells turning at close to 180o, while T cells in lymph nodes and villi do not exhibit this "reversing" movement. T cells in the lung also showed significantly decreased meandering ratios and increased confinement compared to T cells in lymph nodes and villi. These differences in motility patterns led to a decrease in the total volume scanned by T cells in lung compared to T cells in lymph node and villi. These results suggest that the tissue environment in which T cells move can impact the type of motility and ultimately, the efficiency of T cell search for target cells within specialized tissues such as the lung.

Challenges in cybersecurity: Lessons from biological defense systems.

Defending against novel, repeated, or unpredictable attacks, while avoiding attacks on the 'self', are the central problems of both mammalian immune systems and computer systems. Both systems have been studied in great detail, but with little exchange of information across the different disciplines. Here, we present a conceptual framework for structured comparisons across the fields of biological immunity and cybersecurity, by framing the context of defense, considering different (combinations of) defensive strategies, and evaluating defensive performance. Throughout this paper, we pose open questions for further exploration. We hope to spark the interdisciplinary discovery of general principles of optimal defense, which can be understood and applied in biological immunity, cybersecurity, and other defensive realms.

A general model for effects of temperature on ectotherm ontogenetic growth and development.

The temperature size rule (TSR) is the tendency for ectotherms to develop faster but mature at smaller body sizes at higher temperatures. It can be explained by a simple model in which the rate of growth or biomass accumulation and the rate of development have different temperature dependence. The model accounts for both TSR and the less frequently observed reverse-TSR, predicts the fraction of energy allocated to maintenance and synthesis over the course of development, and also predicts that less total energy is expended when developing at warmer temperatures for TSR and vice versa for reverse-TSR. It has important implications for effects of climate change on ectothermic animals.

Population stability, cooperation, and the invasibility of the human species.

The biogeographic expansion of modern humans out of Africa began approximately 50,000 years ago. This expansion resulted in the colonization of most of the land area and habitats throughout the globe and in the replacement of preexisting hominid species. However, such rapid population growth and geographic spread is somewhat unexpected for a large primate with a slow, density-dependent life history. Here, we suggest a mechanism for these outcomes by modifying a simple density-dependent population model to allow varying levels of intraspecific competition for finite resources. Reducing intraspecific competition increases carrying capacities, growth rates, and stability, including persistence times and speed of recovery from perturbations. Our model suggests that the energetic benefits of cooperation in modern humans may have outweighed the slow rate of human population growth, effectively ensuring that once modern humans colonized a region long-term population persistence was near inevitable. Our model also provides insight into the interplay of structural complexity and stability in social species.

Agent-based modeling predicts RAC1 is critical for ovarian cancer metastasis.

Experimental and computational studies pinpoint rate-limiting step(s) in metastasis governed by Rac1. Using ovarian cancer cell and animal models, Rac1 expression was manipulated, and quantitative measurements of cell-cell and cell-substrate adhesion, cell invasion, mesothelial clearance, and peritoneal tumor growth discriminated the tumor behaviors most highly influenced by Rac1. The experimental data were used to parameterize an agent-based computational model simulating peritoneal niche colonization, intravasation, and hematogenous metastasis to distant organs. Increased ovarian cancer cell survival afforded by the more rapid adhesion and intravasation upon Rac1 overexpression is predicted to increase the numbers of and the rates at which tumor cells are disseminated to distant sites. Surprisingly, crowding of cancer cells along the blood vessel was found to decrease the numbers of cells reaching a distant niche irrespective of Rac1 overexpression or knockdown, suggesting that sites for tumor cell intravasation are rate limiting and become accessible if cells intravasate rapidly or are displaced due to diminished viability. Modeling predictions were confirmed through animal studies of Rac1-dependent metastasis to the lung. Collectively, the experimental and modeling approaches identify cell adhesion, rapid intravasation, and survival in the blood as parameters in the ovarian metastatic cascade that are most critically dependent on Rac1.

Liquid brains, solid brains.

Cognitive networks have evolved a broad range of solutions to the problem of gathering, storing and responding to information. Some of these networks are describable as static sets of neurons linked in an adaptive web of connections. These are 'solid' networks, with a well-defined and physically persistent architecture. Other systems are formed by sets of agents that exchange, store and process information but without persistent connections or move relative to each other in physical space. We refer to these networks that lack stable connections and static elements as 'liquid' brains, a category that includes ant and termite colonies, immune systems and some microbiomes and slime moulds. What are the key differences between solid and liquid brains, particularly in their cognitive potential, ability to solve particular problems and environments, and information-processing strategies? To answer this question requires a new, integrative framework. This article is part of the theme issue 'Liquid brains, solid brains: How distributed cognitive architectures process information'.

How does mobility help distributed systems compute?

Brains are composed of connected neurons that compute by transmitting signals. The neurons are generally fixed in space, but the communication patterns that enable information processing change rapidly. By contrast, other biological systems, such as ant colonies, bacterial colonies, slime moulds and immune systems, process information using agents that communicate locally while moving through physical space. We refer to systems in which agents are strongly connected and immobile as solid, and to systems in which agents are not hardwired to each other and can move freely as liquid. We ask how collective computation depends on agent movement. A liquid cellular automaton (LCA) demonstrates the effect of movement and communication locality on consensus problems. A simple mathematical model predicts how these properties of the LCA affect how quickly information propagates through the system. While solid brains allow complex network structures to move information over long distances, mobility provides an alternative way for agents to transport information when long-range connectivity is expensive or infeasible. Our results show how simple mobile agents solve global information processing tasks more effectively than similar systems that are stationary. This article is part of the theme issue 'Liquid brains, solid brains: How distributed cognitive architectures process information'.

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.

Revisiting a model of ontogenetic growth: estimating model parameters from theory and data.

The ontogenetic growth model (OGM) of West et al. provides a general description of how metabolic energy is allocated between production of new biomass and maintenance of existing biomass during ontogeny. Here, we reexamine the OGM, make some minor modifications and corrections, and further evaluate its ability to account for empirical variation on rates of metabolism and biomass in vertebrates both during ontogeny and across species of varying adult body size. We show that the updated version of the model is internally consistent and is consistent with other predictions of metabolic scaling theory and empirical data. The OGM predicts not only the near universal sigmoidal form of growth curves but also the M(1/4) scaling of the characteristic times of ontogenetic stages in addition to the curvilinear decline in growth efficiency described by Brody. Additionally, the OGM relates the M(3/4) scaling across adults of different species to the scaling of metabolic rate across ontogeny within species. In providing a simple, quantitative description of how energy is allocated to growth, the OGM calls attention to unexplained variation, unanswered questions, and opportunities for future research.

Scaling theory for information networks.

Networks distribute energy, materials and information to the components of a variety of natural and human-engineered systems, including organisms, brains, the Internet and microprocessors. Distribution networks enable the integrated and coordinated functioning of these systems, and they also constrain their design. The similar hierarchical branching networks observed in organisms and microprocessors are striking, given that the structure of organisms has evolved via natural selection, while microprocessors are designed by engineers. Metabolic scaling theory (MST) shows that the rate at which networks deliver energy to an organism is proportional to its mass raised to the 3/4 power. We show that computational systems are also characterized by nonlinear network scaling and use MST principles to characterize how information networks scale, focusing on how MST predicts properties of clock distribution networks in microprocessors. The MST equations are modified to account for variation in the size and density of transistors and terminal wires in microprocessors. Based on the scaling of the clock distribution network, we predict a set of trade-offs and performance properties that scale with chip size and the number of transistors. However, there are systematic deviations between power requirements on microprocessors and predictions derived directly from MST. These deviations are addressed by augmenting the model to account for decentralized flow in some microprocessor networks (e.g. in logic networks). More generally, we hypothesize a set of constraints between the size, power and performance of networked information systems including transistors on chips, hosts on the Internet and neurons in the brain.

Quantitative Measurement of Naïve T Cell Association With Dendritic Cells, FRCs, and Blood Vessels in Lymph Nodes.

T cells play a vital role in eliminating pathogenic infections. To activate, naïve T cells search lymph nodes (LNs) for dendritic cells (DCs). Positioning and movement of T cells in LNs is influenced by chemokines including CCL21 as well as multiple cell types and structures in the LNs. Previous studies have suggested that T cell positioning facilitates DC colocalization leading to T:DC interaction. Despite the influence chemical signals, cells, and structures can have on naïve T cell positioning, relatively few studies have used quantitative measures to directly compare T cell interactions with key cell types. Here, we use Pearson correlation coefficient (PCC) and normalized mutual information (NMI) to quantify the extent to which naïve T cells spatially associate with DCs, fibroblastic reticular cells (FRCs), and blood vessels in LNs. We measure spatial associations in physiologically relevant regions. We find that T cells are more spatially associated with FRCs than with their ultimate targets, DCs. We also investigated the role of a key motility chemokine receptor, CCR7, on T cell colocalization with DCs. We find that CCR7 deficiency does not decrease naïve T cell association with DCs, in fact, CCR7-/- T cells show slightly higher DC association compared with wild type T cells. By revealing these associations, we gain insights into factors that drive T cell localization, potentially affecting the timing of productive T:DC interactions and T cell activation.