Search

Dendritic cells (DCs) are the most effective professional antigen-presenting cell. They ferry antigen from the extremities to T cells and are essential for the initiation of an adaptive immune response. Despite interest in how DCs respond to chemical stimuli, there have been few attempts to model DC migration. In this paper, we simulate the motility of DCs by modeling the generation of forces by filopodia and a force balance on the cell. The direction of fliopodial extension is coupled to differential occupancy of cognate chemokine receptors across the cell. Our model simulates chemokinesis and chemotaxis in a variety of chemical and mechanical environments. Simulated DCs undergoing chemokinesis were measured to have a speed of 5.1 ± 0.07 μm·min-1 and a persistence time of 3.2 ± 0.46 min, consistent with experiment. Cells undergoing chemotaxis exhibited a stronger chemotactic response when exposed to lower average chemokine concentrations, also consistent with experiment. We predicted that when placed in two opposing gradients, cells will cluster in a line, which we call the “line of equistimulation;” this clustering has also been observed. We calculated the effect of varying gradient steepness on the line of equistimulation, with steeper gradients resulting in tighter clustering. Moreover, gradients are found to be most potent when cells are in a gradient of chemokine whose mean concentration is close to the binding of the Kd to the receptor, and least potent when the mean concentration is 0.1Kd. Comparing our simulations to experiment, we can give a quantitative measure of the strength of certain chemokines relative to others. Assigning the signal of CCL19 binding CCR7 a baseline strength of 1, we found CCL21 binding CCR7 had a strength of 0.28, and CXCL12 binding CXCR4 had a strength of 0.30. These differences emerge despite both chemokines having virtually the same Kd, suggesting a mechanism of signal amplification in DCs requiring further study. [ABSTRACT FROM AUTHOR]

The Notch pathway plays a vital role in determining whether cells in the intestinal epithelium adopt a secretory or an absorptive phenotype. Cell fate specification is coordinated via Notch’s interaction with the canonical Wnt pathway. Here, we propose a new mathematical model of the Notch and Wnt pathways, in which the Hes1 promoter acts as a hub for pathway crosstalk. Computational simulations of the model can assist in understanding how healthy intestinal tissue is maintained, and predict the likely consequences of biochemical knockouts upon cell fate selection processes. Chemical reaction network theory (CRNT) is a powerful, generalised framework which assesses the capacity of our model for monostability or multistability, by analysing properties of the underlying network structure without recourse to specific parameter values or functional forms for reaction rates. CRNT highlights the role of β-catenin in stabilising the Notch pathway and damping oscillations, demonstrating that Wnt-mediated actions on the Hes1 promoter can induce dynamic transitions in the Notch system, from multistability to monostability. Time-dependent model simulations of cell pairs reveal the stabilising influence of Wnt upon the Notch pathway, in which β-catenin- and Dsh-mediated action on the Hes1 promoter are key in shaping the subcellular dynamics. Where Notch-mediated transcription of Hes1 dominates, there is Notch oscillation and maintenance of fate flexibility; Wnt-mediated transcription of Hes1 favours bistability akin to cell fate selection. Cells could therefore regulate the proportion of Wnt- and Notch-mediated control of the Hes1 promoter to coordinate the timing of cell fate selection as they migrate through the intestinal epithelium and are subject to reduced Wnt stimuli. Furthermore, mutant cells characterised by hyperstimulation of the Wnt pathway may, through coupling with Notch, invert cell fate in neighbouring healthy cells, enabling an aberrant cell to maintain its neighbours in mitotically active states. [ABSTRACT FROM AUTHOR]

Phase-amplitude coupling (PAC), a type of cross-frequency coupling (CFC) where the phase of a low-frequency rhythm modulates the amplitude of a higher frequency, is becoming an important indicator of information transmission in the brain. However, the neurobiological mechanisms underlying its generation remain undetermined. A realistic, yet tractable computational model of the phenomenon is thus needed. Here we analyze a neural mass model of a cortical column, comprising fourteen neuronal populations distributed across four layers (L2/3, L4, L5 and L6). A control analysis showed that the conditional transfer entropy (cTE) measure is able to correctly estimate the flow of information between neuronal populations. Then, we computed cTE from the phases to the amplitudes of the oscillations generated in the cortical column. This approach provides information regarding directionality by distinguishing PAC from APC (amplitude-phase coupling), i.e. the information transfer from amplitudes to phases, and can be used to estimate other types of CFC such as amplitude-amplitude coupling (AAC) and phase-phase coupling (PPC). While experiments often only focus on one or two PAC combinations (e.g., theta-gamma or alpha-gamma), we found that a cortical column can simultaneously generate almost all possible PAC combinations, depending on connectivity parameters, time constants, and external inputs. PAC interactions with and without an anatomical equivalent (direct and indirect interactions, respectively) were analyzed. We found that the strength of PAC between two populations was strongly correlated with the strength of the effective connections between the populations and, on average, did not depend on whether the PAC connection was direct or indirect. When considering a cortical column circuit as a complex network, we found that neuronal populations making indirect PAC connections had, on average, higher local clustering coefficient, efficiency, and betweenness centrality than populations making direct connections and populations not involved in PAC connections. This suggests that their interactions were more effective when transmitting information. Since approximately 60% of the obtained interactions represented indirect connections, our results highlight the importance of the topology of cortical circuits for the generation of the PAC phenomenon. Finally, our results demonstrated that indirect PAC interactions can be explained by a cascade of direct CFC and same-frequency band interactions, suggesting that PAC analysis of experimental data should be accompanied by the estimation of other types of frequency interactions for an integrative understanding of the phenomenon. [ABSTRACT FROM AUTHOR]

Sequence variation can affect the physiological state of the immune system. Major experimental efforts targeted at understanding the genetic control of the abundance of immune cell subpopulations. However, these studies are typically focused on a limited number of immune cell types, mainly due to the use of relatively low throughput cell-sorting technologies. Here we present an algorithm that can reveal the genetic basis of inter-individual variation in the abundance of immune cell types using only gene expression and genotyping measurements as input. Our algorithm predicts the abundance of immune cell subpopulations based on the RNA levels of informative marker genes within a complex tissue, and then provides the genetic control on these predicted immune traits as output. A key feature of the approach is the integration of predictions from various sets of marker genes and refinement of these sets to avoid spurious signals. Our evaluation of both synthetic and real biological data shows the significant benefits of the new approach. Our method, VoCAL, is implemented in the freely available R package ComICS. [ABSTRACT FROM AUTHOR]

To understand the process of innate immune fungal recognition, we developed computational tools for the rigorous quantification and comparison of receptor recruitment and distribution at cell-cell contact sites. We used these tools to quantify pattern recognition receptor spatiotemporal distributions in contacts between primary human dendritic cells and the fungal pathogens C. albicans, C. parapsilosis and the environmental yeast S. cerevisiae, imaged using 3D multichannel laser scanning confocal microscopy. The detailed quantitative analysis of contact sites shows that, despite considerable biochemical similarity in the composition and structure of these species' cell walls, the receptor spatiotemporal distribution in host-microbe contact sites varies significantly between these yeasts. Our findings suggest a model where innate immune cells discriminate fungal microorganisms based on differential mobilization and coordination of receptor networks. Our analysis methods are also broadly applicable to a range of cell-cell interactions central to many biological problems. [ABSTRACT FROM AUTHOR]

Johne's disease (JD), a persistent and slow progressing infection of ruminants such as cows and sheep, is caused by slow replicating bacilli Mycobacterium avium subspecies paratuberculosis (MAP) infecting macrophages in the gut. Infected animals initially mount a cell-mediated CD4 T cell response against MAP which is characterized by the production of interferon (Th1 response). Over time, Th1 response diminishes in most animals and antibody response to MAP antigens becomes dominant (Th2 response). The switch from Th1 to Th2 response occurs concomitantly with disease progression and shedding of the bacteria in feces. Mechanisms controlling this Th1/Th2 switch remain poorly understood. Because Th1 and Th2 responses are known to cross-inhibit each other, it is unclear why initially strong Th1 response is lost over time. Using a novel mathematical model of the immune response to MAP infection we show that the ability of extracellular bacteria to persist outside of macrophages naturally leads to switch of the cellular response to antibody production. Several additional mechanisms may also contribute to the timing of the Th1/Th2 switch including the rate of proliferation of Th1/Th2 responses at the site of infection, efficiency at which immune responses cross-inhibit each other, and the rate at which Th1 response becomes exhausted over time. Our basic model reasonably well explains four different kinetic patterns of the Th1/Th2 responses in MAP-infected sheep by variability in the initial bacterial dose and the efficiency of the MAP-specific T cell responses. Taken together, our novel mathematical model identifies factors of bacterial and host origin that drive kinetics of the immune response to MAP and provides the basis for testing the impact of vaccination or early treatment on the duration of infection. [ABSTRACT FROM AUTHOR]