B Cell Activation Triggered by the Formation of the Small Receptor Cluster: A Computational Study

We proposed a spatially extended model of early events of B cell receptors (BCR) activation, which is based on mutual kinase-receptor interactions that are characteristic for the immune receptors and the Src family kinases. These interactions lead to the positive feedback which, together with two nonlinearities resulting from the double phosphorylation of receptors and Michaelis-Menten dephosphorylation kinetics, are responsible for the system bistability. We demonstrated that B cell can be activated by a formation of a tiny cluster of receptors or displacement of the nucleus. The receptors and Src kinases are activated, first locally, in the locus of the receptor cluster or the region where the cytoplasm is the thinnest. Then the traveling wave of activation propagates until activity spreads over the whole cell membrane. In the models in which we assume that the kinases are free to diffuse in the cytoplasm, we found that the fraction of aggregated receptors, capable to initiate B cell activation decreases with the decreasing thickness of cytoplasm and decreasing kinase diffusion. When kinases are restricted to the cell membrane - which is the case for most of the Src family kinases - even a cluster consisting of a tiny fraction of total receptors becomes activatory. Interestingly, the system remains insensitive to the modest changes of total receptor level. The model provides a plausible mechanism of B cells activation due to the formation of small receptors clusters collocalized by binding of polyvalent antigens or arising during the immune synapse formation

The Morphostatic Limit for a Model of Skeletal Pattern Formation in the Vertebrate Limb

A recently proposed mathematical model of a “core” set of cellular and molecular interactions present in the developing vertebrate limb was shown to exhibit pattern-forming instabilities and limb skeleton-like patterns under certain restrictive conditions, suggesting that it may authentically represent the underlying embryonic process (Hentschel et al., Proc. R. Soc. B 271, 1713–1722,2004). The model, an eight-equation system of partial differential equations, incorporates the behavior of mesenchymal cells as “reactors,” both participating in the generation of morphogen patterns and changing their state and position in response to them. The full system, which has smooth solutions that exist globally in time, is nonetheless highly complex and difficult to handle analytically or numerically. According to a recent classification of developmental mechanisms (Salazar-Ciudad et al., Development 130, 2027–2037, 2003), the limb model of Hentschel et al. is “morphodynamic,” since differentiation of new cell types occurs simultaneously with cell rearrangement. This contrasts with “morphostatic” mechanisms, in which cell identity becomes established independently of cell rearrangement. Under the hypothesis that development of some vertebrate limbs employs the core mechanism in a morphostatic fashion, we derive in an analytically rigorous fashion a pair of equations representing the spatiotemporal evolution of the morphogen fields under the assumption that cell differentiation relaxes faster than the evolution of the overall cell density (i.e., the morphostatic limit of the full system). This simple reaction–diffusion system is unique in having been derived analytically from a substantially more complex system involving multiple morphogens, extracellular matrix deposition, haptotaxis, and cell translocation. We identify regions in the parameter space of the reduced system where Turing-type pattern formation is possible, which we refer to as its “Turing space.” Obtained values of the parameters are used in numerical simulations of the reduced system, using a new Galerkin finite element method, in tissue domains with nonstandard geometry. The reduced system exhibits patterns of spots and stripes like those seen in developing limbs, indicating its potential utility in hybrid continuum-discrete stochastic modeling of limb development. Lastly, we discuss the possible role in limb evolution of selection for increasingly morphostatic developmental mechanisms

On Multiscale Approaches to 3-Dimensional Modeling of Morphogenesis

In this paper we present the foundation of a unified, object-oriented, three-dimensional (3D) biomodeling environment, which allows us to integrate multiple submodels at scales from subcellular to tissues and organs. Our current implementation combines a modified discrete model from statistical mechanics, the Cellular Potts Model (CPM), with a continuum reaction-diffusion (RD) model and a state automaton with well-defined conditions for cell differentiation transitions to model genetic regulation. This environment allows us to rapidly and compactly create computational models of a class of complex developmental phenomena. To illustrate model development, we simulate a simplified version of the formation of the skeletal pattern in a growing embryonic vertebrate limb

Traveling waves for reaction-diffusion PDE coupled to difference equation with nonlocal dispersal term and time delay

International audienceWe consider a class of biological models represented by a system composed of reactiondiffusion PDE coupled with difference equations (renewal equations) in n -dimensional space, with nonlocal dispersal terms and implicit time delays. The difference equation generally arises, by means of the method of characteristics, from an age-structured partial differential system. Using upper and lower solutions, we study the existence of monotonic planar traveling wave fronts connecting the extinction state to the uniform positive state. The corresponding minimum wave speed is also obtained. In addition, we investigate the effect of the parameters on this minimum wave speed and we give a detailed analysis of its asymptotic behavior

Traveling waves in a coupled reaction-diffusion and difference model of hematopoiesis

International audienceThe formation and development of blood cells is a very complex process, called hematopoiesis. This process involves a small population of cells called hematopoietic stem cells (HSCs). The HSCs are undifferentiated cells, located in the bone marrow before they become mature blood cells and enter the blood stream. They have a unique ability to produce either similar cells (self-renewal), or cells engaged in one of different lineages of blood cells: red blood cells, white cells and platelets (differentiation). The HSCs can be either in a proliferating or in a quiescent phase. In this paper, we distinguish between dividing cells that enter directly to the quiescent phase and dividing cells that return to the proliferating phase to divide again. We propose a mathematical model describing the dynamics of HSC population, taking into account their spatial distribution. The resulting model is an age-structured reaction-diffusion system. The method of characteristics reduces this model to a coupled reaction-diffusion equation and difference equation with delay. We study the existence of traveling wave fronts connecting the zero steady state with the unique positive uniform one. We use a monotone iteration technique coupled with the upper and lower solutions method

Viral Infection Spreading and Mutation in Cell Culture

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY)International audienc

Reaction-Diffusion Model of Atherosclerosis Development.

International audienceAtherosclerosis begins as an inflammation in blood vessel walls (intima). The inflammatory response of the organism leads to the recruitment of monocytes. Trapped in the intima, they differentiate into macrophages and foam cells leading to the production of inflammatory cytokines and further recruitment of white blood cells. This self-accelerating process, strongly influenced by low-density lipoproteins (cholesterol), results in a dramatic increase of the width of blood vessel walls, formation of an atherosclerotic plaque and, possibly, of its rupture. We suggest a 2D mathematical model of the initiation and development of atherosclerosis which takes into account the concentration of blood cells inside the intima and of pro- and anti-inflammatory cytokines. The model represents a reaction-diffusion system in a strip with nonlinear boundary conditions which describe the recruitment of monocytes as a function of the concentration of inflammatory cytokines. We prove the existence of travelling waves described by this system and confirm our previous results which suggest that atherosclerosis develops as a reaction-diffusion wave. The theoretical results are confirmed by the results of numerical simulations

Mathematical modeling of respiratory viral infection and applications to SARS‐CoV‐2 progression

International audienceViral infection in cell culture and tissue is modeled with delay reaction‐diffusion equations. It is shown that progression of viral infection can be characterized by the viral replication number, time‐dependent viral load, and the speed of infection spreading. These three characteristics are determined through the original model parameters including the rates of cell infection and of virus production in the infected cells. The clinical manifestations of viral infection, depending on tissue damage, correlate with the speed of infection spreading, while the infectivity of a respiratory infection depends on the viral load in the upper respiratory tract. Parameter determination from the experiments on Delta and Omicron variants allows the estimation of the infection spreading speed and viral load. Different variants of the SARS‐CoV‐2 infection are compared confirming that Omicron is more infectious and has less severe symptoms than Delta variant. Within the same variant, spreading speed (symptoms) correlates with viral load allowing prognosis of disease progression