Next Article in Journal
Heterogeneous Tau Oligomers as Molecular Targets for Alzheimer’s Disease and Related Tauopathies
Next Article in Special Issue
Adherent Moving of Polymers in Spherical Confined Binary Semiflexible Ring Polymer Mixtures
Previous Article in Journal / Special Issue
The Modification of Polyvinyl Alcohol for Ice Nucleation Based upon the Structures of Antifreeze Glycoproteins Found in Antarctic Fish
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RBL-2H3 Mast Cell Receptor Dynamics in the Immunological Synapse

by
Ming Chih Tsai
1,2,* and
Kathrin Spendier
2,3
1
Department of Physics and Energy Science, University of Colorado Colorado Springs, Colorado Springs, CO 80918, USA
2
BioFrontiers Center, University of Colorado Colorado Springs, Colorado Springs, CO 80918, USA
3
Quantinuum, Broomfield, CO 80021, USA
*
Author to whom correspondence should be addressed.
Biophysica 2022, 2(4), 428-439; https://doi.org/10.3390/biophysica2040038
Submission received: 21 September 2022 / Revised: 26 October 2022 / Accepted: 3 November 2022 / Published: 7 November 2022
(This article belongs to the Collection Feature Papers in Biophysics)

Abstract

:
The RBL-2H3 mast cell immunological synapse dynamics is often simulated with reaction–diffusion and Fokker–Planck equations. The equations focus on how the cell synapse captures receptors following an immune response, where the receptor capture at the immunological site appears to be a delayed process. This article investigates the physical nature and mathematics behind such time-dependent delays. Using signal processing methods, convolution and cross-correlation-type delay capture simulations give a χ-squared range of 22 to 60, in good agreement with experimental results. The cell polarization event is offered as a possible explanation for these capture delays, where polarizing rates measure how fast the cell polarization event occurs. In the case of RBL-2H3 mast cells, polarization appears to be associated with cytoskeletal rearrangement; thus, both cytoskeletal and diffusional components are considered. From these simulations, a maximum polarizing rate ranging from 0.0057 s−2 to 0.031 s−2 is obtained. These results indicate that RBL-2H3 mast cells possess both temporal and spatial memory, and cell polarization is possibly linked to a Turing-type pattern formation.

1. Introduction

The RBL-2H3 mast cell signaling pathway plays an important role in human inflammatory and allergic responses. Mast cells signal by an immunological synapse that forms upon receptor contact with antigens [1,2]. In a synapse, receptors that have been activated by binding usually form clusters all over the cell, but eventually aggregate in a central patch at the point of first contact. The tracking of single receptors within these clusters reveals what appears outwardly to be a diffusion-mediated reaction process [3]. To quantify this process mathematically, various forms of the reaction–diffusion equation, i.e., the Smoluchowski equation, were solved both analytically and numerically by several authors [4,5,6]. Each form of the Smoluchowski equation employed a slightly different capture mechanism for the central patch aggregate, yet the results of their simulations had significant differences. Capture mechanisms used include those of infinite (perfect) absorption, finite absorption, and enhanced diffusion at the capture site. Infinite absorption resulted in a standard Gaussian probability distribution that moves from initial location (Dirac Delta initial condition) toward the capture site, while finite absorption produced a slower overall concentration decrease that can be controlled by a capture rate constant, making it ideal for fitting to real experimental data. Capture rate constants can also be used in enhanced diffusion, where a capture potential gathers the probability density distribution into an energy minimum at the capture site, rather than absorbing them out of the system. In a harmonic potential, probability density is Gaussian as usual, but its dynamics differ from that of finite absorption. With a linear V-potential, probability density is no longer Gaussian. The boundary conditions used to solve these differential equations also contribute to differences in the results. Thus, by choosing the appropriate capture term, boundary condition, and initial condition, these diffusion–reaction equations could have a wide range of applications in molecular systems such as excitons, photosynthetic units, crystal aggregates, and receptor clusters [4].
In the case of receptor clusters, very specific capture conditions are required. In a recent experiment [7], receptors in an immunological synapse on mast cells were observed. Receptors form initial clusters spread over the area of contact, followed by a slower movement towards the point of first contact to form the central patch. While the mechanism for RBL-2H3 mast cell synapse cluster formation has not been directly found, evidence [1] supports the theory of cytoskeletal rearrangement that occurs during cell contact with antigens. This cluster formation causes a delay in the aggregate of the central patch, since the central patch must compete with other clusters for free receptors. However, cytoskeletal rearrangement favors the point of initial contact, and all receptors in the clusters eventually diffuse into the central patch. Therefore, according to these observations, the initial probability distribution is uniform throughout the cell, Neumann boundary conditions are given by the boundary of the contact area, and the capture process at the central patch is delayed by peripheral clusters, where the radius of the central patch increases over time.
A study of mast cell receptor dynamics using some of these conditions is given in reference [7]. The authors improve upon previous theories by incorporating a dynamic radius, R(t), for the central patch, which increases its area (radius) in proportion to the amount of free receptors remaining outside the central patch,
R 2 ( t ) R 2 ( 0 ) R 2 ( ) R 2 ( 0 ) = Q ( 0 ) Q ( t ) Q ( 0 ) Q ( )
where R ( 0 ) is the initial central patch radius, R ( ) the final radius; Q ( 0 ) the initial fraction of total receptors remaining outside of the central patch, Q ( ) the final fraction. Q ( t ) is calculated by integrating (in polar coordinates) the probability density function, P ( r , θ , t ) , of receptors over the area outside the central patch, thus giving an equation for R ( t ) once the probability density is found.
Using Laplace transform, P ( r , θ , t ) is then calculated from the solution to the Smoluchowski equation with a specific capture term,
d P ( r , θ , t ) d t = D   2 P ( r , θ , t ) 2   π r   δ ( r R )   0 t d s   C ( t s ) P ( r , θ , s ) ,
where the Laplacian term represents the diffusion process with diffusion coefficient D, and the Dirac Delta term represents the reaction process, or the capture, with discrete point traps at radius R from the origin. In polar coordinates, the delta function forms the boundary of the central patch, and the integral from time zero to t is the convolution-type memory function which will be discussed below. Since the delta function represents infinite capture and is thus unphysical in this case (the Dirac Delta representation is used in quantum systems, which is not the case here), the authors replace the infinite delta function solution with a finite reaction rate obtained from their previous work [5]. The result is an analytical solution Q ( ε , r 0 , R ( 0 ) ) , where ε is the Laplace variable, and r 0 refers to an initial condition of receptors laying along a ring of radius r 0 . Then, approximating the peripheral clusters by circular rings outside the central patch and integrating Q ( ε , r 0 , R ( 0 ) ) over the area outside the central patch, the authors obtain a solution for Q ( ε ) subject to absorption at a stationary central patch boundary, with an initial condition of concentric circular receptor rings. The combined effect of the central patch capture and peripheral cluster distribution on the quantity Q ( ε ) is described by
Q ( ε ) = N R ( 0 ) r 0 Q ( ε , r 0 , R ( 0 ) )   r 0   d r 0 ,
where N is the normalization factor that normalizes the maximum of Q ( ε ) to one, r 0 the entire contact area radius, and Q ( ε , r 0 , R ( 0 ) ) is an analytical solution that not only depends on the Laplace variable ε but also on the initial central patch radius R ( 0 ) .
However, we note that the central patch boundary is a dynamic quantity, with the radius increase described by R ( t ) in Equation (1). Increasing the central patch size over time is equivalent to dynamically increasing the overall absorption rate. Replacing R ( 0 ) in Equation (3) with R ( t ) drastically increases the complexity of the integral, making it less attainable by analytical methods. Therefore, we suggest the use of numerical approaches to solve such problems.
Additionally, the convolution-type function in the capture term reflects and quantifies the physical situation of the delayed central patch aggregation caused by peripheral clusters competing for free receptors. This convolution-type function has deep physical meanings not only in signal processing but also in quantitative studies of memory and learning [8,9]. In reference [7], the authors refer to these types of functions as memory functions to approach this problem from a cell recollection point of view. With this approach, a cell is able to “learn” from what happened to it in the past, in this case its past history of free receptors. This memory term also makes it possible to study signals that cause the cluster formation and, therefore, illuminate possible underlying mechanisms. However, using memory capture, the theory was fitted to experimental data and a polarizing time constant of 10−12 s−1 was obtained. The inverse of this number, 1012 s, far exceeds the time duration of most cellular processes.
Nevertheless, in a related field of study [10], it is shown that cells are able to communicate through local polarization and depolarization upon contact, providing a new insight into the inner workings of the mast cell immunological synapse. In fact, polarization plays an important role in the immunological synapse of many cells [11,12]. It is also known that polarization in different types of cells is stimuli-specific. Most cell types respond to chemical stimuli, but some will respond to mechanical and other types of contact. Moreover, different cell types also react to polarization stimuli in different ways: some cell types rearrange their cytoskeletal components, while others alter their cytosolic signaling molecules [13]. In the case of RBL-2H3 mast cell, polarization appears to be initiated by antigen binding and is mediated by actin and myosin, resulting in cytoskeletal rearrangement [14]. Thus, in the present study, we continue to build into this theory a capture term that consists of both convolution-type and cross-correlation-type signals, which corresponds to the depolarization and polarization event, respectively. Additionally, to transition from polarization to depolarization, we build into the theory a time-delayed Fokker–Planck potential term. Using these theories, we then obtain polarizing time constants ranging from 0.022 s−1 to 0.031 s−1, and maximum polarization is reached approximately 112 s after polarization initiation. This is consistent with experimental observations in the literature [13], where cell polarization can take seconds to minutes, depending on the cell type.

2. Materials and Methods

Due to high symmetry, the 2-D Smoluchowski equation, Equation (2), can be reduced to 1-D when the center of the capturing trap is located at the origin. The only asymmetric component in this problem would be the initial distribution of receptor clusters around the central patch. When these clusters are approximated as circular rings, as was carried out in the previous study, the problem can still be reduced to 1-D. Reducing the problem to 1-D or to 2-D circular rings saves both computational memory and time, although doing so does not allow details of the 2-D solution such as movement speed of independent clusters toward the central patch to be simulated. With 1-D capture, such features are then hidden in the capture term’s delay-type memory function, i.e., C ( t ) in Equation (2).
The memory function in reference [7] is an important improvement over previous diffusion reaction capture terms, without the use of such functions the central patch aggregation delay could not be properly simulated in 1-D. Illustrations in Figure 1 and Figure 2 show two major types of memory functions used in this study.
Figure 1 shows convolution-type capture delay. The convolution-type capture memory function, C ( t ) = c 2 α   E x p ( α t ) is shown in Figure 1a, while an idealized rectangular probability density function for free receptors is shown in Figure 1b. Variable c 2 is the capture strength, α the polarizing time constant, and x 0 represents an arbitrary point in the 1-D probability density function, P ( x 0 , s ) , in terms of convolution time s . The rectangular shape is used here for illustrative purposes only and does not depict the real density profile of free receptors. Figure 2 follows the same format for cross-correlation-type capture delay, with C ( t ) replaced by a piecewise function of c 2 α   E x p ( α t ) . Using the rectangular probability density function, differences in convolution and cross-correlation results can be seen clearly in Figure 1c and Figure 2c. Assuming the upward increase in capture rate (Figure 1c) is caused by mast cell polarization upon contact with antigens and the downward curve by depolarization, the convolution-type capture gives fast initial polarization and depolarization. The same reasoning in Figure 2c gives cross-correlation-type capture slow initial polarization and depolarization.
Convolution and cross-correlation provide a physical way to distinguish between polarization and depolarization events, but the Dirac Delta function in front of the capture term in Equation (2) is still unphysical because it represents an infinite capture rate. Thus, replacing the Dirac Delta function with a rescaled normal distribution, n ( t ) = E x p ( r 2 ( R ( t ) / 2 ) 2 ) , we obtain a capture trap size of height one, with the radius R ( t ) set to 2√2 standard deviations away from the center of trap (see Figure 3).
Lastly, we note that, in Figure 1b, P ( x 0 , t ) drops to zero at t = 400 s, yet the capture process in Figure 1c does not stop until t = 800 s. This effectively gives a negative P ( x 0 , t ) for times greater than t = 400 s and is again unphysical. The same is true for Figure 2. Therefore, we rescale the capture term with the probability function itself at time t , i.e., multiply the capture term in Equation (2) by P ( x 0 , t ) / P ( x 0 , 0 ) , to obtain the physical capture results for convolution in Figure 4a and cross-correlation in Figure 4b.
Altogether, the 1-D Smoluchowski equation now reads:
d P ( r , t ) d t = D   2 P ( r , t ) E x p ( r 2 ( R ( t ) 2 ) 2 )   P ( r , t ) P ( r , 0 )   0 t d s   C ( t s ) P ( r , s ) ,
with initial condition P ( r , 0 ) = 1 μm−1 uniform throughout the surface of the cell (the initial peripheral clusters are now taken into account by the delay memory functions, instead of being built into the initial condition as 2-D circular rings) and Neumann conditions at the cell contact boundary r = ±   r 0 .
Moreover, evidence shows that the RBL-2H3 mast cell synapse receptor dynamics is affected by their distance from the central patch, where the peripheral receptor cluster average 2-D speeds are faster when they are further away from the central patch [6]. Assuming the radial 1-D speeds of the peripheral clusters are also faster than the central region and the clusters’ overall movement is toward the center, this results in an enhanced absorption where peripheral clusters are pulled into the central patch. Such movements appear to be due to an underlying cytoskeleton rearrangement [14] associated with cell polarization. Representing this cytoskeletal rearrangement by a specific 1-D diffusional potential and linking it to the transition from polarization to depolarization, we add a Fokker–Planck potential term to Equation (4), obtaining
d P ( r , t ) d t = D   2 P ( r , t ) d d r ( P ( r , t ) d d r u ( r ) ) P ( r , t ) P ( r , 0 )   0 t d s   Π ( t s ) P ( r , s ) E x p ( r 2 ( R ( t ) 2 ) 2 )   P ( r , t ) P ( r , 0 )   0 t d s   C ( t s ) P ( r , s ) ,
where Π ( t ) is equal to Γ α   E x p ( α t ) for convolution, and a piecewise function of Γ α   E x p ( α t ) for cross-correlation. Π ( t ) has the same time delay dynamics (same polarizing time constant α ) as C ( t ) , and, similar to how c 2 describes the capture strength, Γ describes the strength of the potential. Thus, the capture term would describe the external receptor aggregation of the central patch, while the potential term describes the internal cytoskeleton rearrangement of the cell. Here, the V-potential is used, u ( r ) = | r | . This V-potential is chosen not only because it creates a richer motion than the harmonic potential [6], but also it is the median of 3 possible potential configurations (see Figure 5). Note that, due to the delays caused by C ( t ) and Π ( t ) , the effect of both capture and potential is most prominent around the time of maximum capture, t m a x (see Figure 1c and Figure 2c). This t m a x also marks the transition from polarization to depolarization, so the memory function is in fact a mathematical representation of such transitions.
To determine Q ( t ) , which represents the fraction of total free receptors remaining at time t , we integrate the fraction P ( r , t ) / P ( r , 0 ) over the cell contact area,
Q ( t ) = 0 r 0 d r   P ( r , t ) P ( r , 0 ) .
Equations (4) and (5), together with their initial and boundary conditions, are solved numerically by Wolfram Mathematica version 13. Memory-function-type equations such as these are called delay partial differential equations (DPDE) [15]. The quantity R ( t ) is an equation that depends on Q ( t ) , which is a function that depends on the solution P ( r , t ) itself. More importantly, P ( r , s ) is a memory-linked function that depends on the history of P ( r , t ) . These features make it difficult to solve the equations analytically, so numerical approximation methods are usually used [16]. One effective strategy, documented in the Python language [17], uses what is known as “history functions”. The history function represents a function that holds the approximate value of the solution P ( r , t ) while the differential solver slowly corrects and pushes it towards the true solution. This strategy employs solving an initial partial differential equation (PDE) that is similar to the task at hand, using that solution as an initial history function for P ( r , s ) in the DPDE, and then repeatedly solving the DPDE in a loop, each time updating the history function P ( r , s ) , until the history function becomes almost identical to the solution P ( r , t ) . The Mathematica NDSolve function currently does not support non-constant delays [18], so users must implement these loops manually (see accompanying code). Here, we have chosen Mathematica for this study due to our familiarity and experience with the language.

3. Results

To compare experiments with Mathematica simulations, observed mean data on the central patch radius in the RBL-2H3 mast cell synapse is taken from reference [7] and fitted with an arctan curve that resembles the experimental data points (Figure 6). This is carried out to provide visual comparisons with simulation curves and also to obtain better estimates for the initial central patch radius R ( 0 ) and the final central patch radius R ( ) . This estimate is necessary because the initial and final mean radius data points do not appear to lie on the same curve as the rest of the data. If the initial and final mean data points were used directly as R ( 0 ) and R ( ) , simulation results on central patch radius would deviate greatly from a large number of the data points. The source of this error might be due to the difficulty in estimating polarization initiation and end times when these measurements were made.
The shifted and scaled arctan curve thus obtained achieves a high-accuracy fit with an χ-Squared value of 61.2542. The choice of this arctan fit is arbitrary, as any sigmoidal curve with a sufficiently high R 2 and low χ 2 values can be used. Thus, the arctan function is not the only sigmoidal curve that can be used to fit the data, but we note that other sigmoidal curves do not produce any χ 2 values that are significantly different. We then compare the arctan curve with our simulation results using the same initial and final central patch radii as the arctan curve in our simulations.
From the arctan curve, R a r c ( t ) , we determine an initial central patch radius of 1.31 μm and a final central patch radius (by taking the arctan limit to ∞) of 3.15 μm. We also calculate the approximate time when maximum capture rate is achieved, t m a x , from R a r c ( t ) , by setting its second derivative to zero. This gives a t m a x value of 112 s, and, as we shall see below, this value is important for estimating cross-correlation-type capture delay.
Then, taking experimental values from [7], we use an initial Q ( 0 ) of 1, a final Q ( ) of 0.26, a cell contact radius of 5.75 μm, and a diffusion constant of 0.011 μm2 s−1 to numerically solve variations in Equation (4) for P ( r , t ) and determine Q ( t ) from Equation (6). Finally, we plot the corresponding R ( t ) using Equation (1) for several different cases of capture in the subsequent results.
To see the effect of un-delayed capture, Equation (4) is modified using an un-delayed constant capture rate c 0 = 0.4 μm−1 s−1, as in Equation (7). The equation is then solved numerically, and the resulting central patch radius, R u c ( t ) , is plotted in Figure 7.
d P ( r , t ) d t = D   2 P ( r , t ) E x p ( r 2 ( R ( t ) 2 ) 2 )   c 0   P ( r , t ) P ( r , 0 )
To see the effects of delayed capture, C ( t ) in Equation (4) needs to be defined for both cross-correlation and convolution. From the shape of the delayed captures, Figure 1c and Figure 2c, one could see that the experimental curve in Figure 6 must be due to a cross-correlation-type polarization before t m a x , and a convolution-type depolarization afterward. Additionally, unlike convolution, cross correlation memory function must be constructed in pieces. Thus, we define a piecewise function for cross-correlation (Equation (8)) and a non-piecewise function for convolution (Equation (9)),
C ( t ) = {   c 2   α 2   E x p ( α   t ) , 0 t t m a x 0 , e l s e w h e r e
and
C ( t ) = c 2   α   E x p ( α   t )
where the extra dimensionless α factor in cross-correlation is due to its integration form (see Figure 2c). Then, by setting α   E x p ( α   t m a x ) to c 2 , adjusting the value of c 2 and solving for α , we obtain c 2 = 0.3 s−1 and α = 0.0229 s−1 for the cross-correlation type. Convolution values for c 2 and α are similarly obtained, with values c 2 = 1 s−1 and α = 0.031 s−1. The cross-correlation and convolution values are then used to numerically solve Equation (4) using Equation (8) and Equation (9), respectively, and the results of their central patch radii, R c c ( t ) for cross-correlation and R c o ( t ) for convolution, are plotted in Figure 8. Both R c c ( t ) and R c o ( t ) are more accurate (lower χ2) than R a r c ( t ) in their respective time intervals ( t < 112 for cross correlation and t > 160 for convolution).
Finally, to see the effect of a time-dependent V-potential, Equation (5) is numerically solved for the case of cross correlation, where Π ( t ) is described by C ( t ) in Equation (8), with c 2 replaced by Γ . Using the values Γ = 0.185 μm2s−1, c 2 = 0.26 s−1, and α = 0.022 s−1, the resulting radius, R c c v ( t ) , is plotted in Figure 9, where R c c v ( t ) is, again, more accurate than R a r c ( t ) .
The resulting time and strength constants from the above simulations are summarized in Table 1, where c 2 is the capture strength in units of s−1, α the polarizing time constant in s−1, c 2   α the maximum polarizing rate in s−2, and Γ the potential strength in μm2 s−1.

4. Discussion

As expected, the un-delayed capture radius does not match experimental data points (Figure 7). However, the latter part of R u c ( t ) in Figure 7 does coincide with the data points and the arctan curve. Physically, this is due to the receptors being absorbed out of the system (into the central patch) naturally, which has nothing to do with depolarization. This reflects the fact that, unlike polarization, depolarization is simply not observed in many systems due to the natural decline of the primary cellular process. In Figure 8, the central patch radius increase is well-described by cross correlation for t < 112 s, and by convolution for t > 160 s. This corresponds to a slow initial polarization and a fast initial depolarization in the RBL-2H3 mast cell, although it remains unclear whether the decrease in the capture rate for times t > 160 s is due to depolarization or the natural decrease in the number of free receptors. The maximum polarizing rate, c 2   α , is 0.0069 s−2 for the cross-correlation type, and 0.031 s−2 for the convolution type. The fact that these values are different for polarization and depolarization suggest that these processes are controlled by separate events in the RBL-2H3 mast cell.
Additionally, the transition from polarization to depolarization, from time t = 112 s to t = 160 s, is not well-described by either the cross-correlation or the convolution-type capture delay, so we suspect some kind of phase change has occurred in the mast cell during this time. This phase change can be simulated by adding a diffusional potential delay, i.e., by considering cytoskeleton rearrangement. The result of this is shown in Figure 9, where experimental data are well-described by using a cross-correlation-type delay in both the capture and potential terms in our simulations (Equation (5)). Thus, the receptor clusters are dynamically attracted to the center of the trap, which is also consistent with experimental observations [6]. The maximum polarizing rate obtained with the inclusion of the V-potential is 0.0057 s−2.
The key features of Equation (5) are the memory function, C ( t ) , and the V-potential, u ( r ) . We note that if the V-potential is replaced by the Y-potential (Figure 5c), slower radial diffusion speeds would bring the latter part of the cross-correlation curve, R c c v ( t ) in Figure 9, slightly closer to a plateau, which could make the curve even more accurate. The success of this equation has deep physical implications. It means that the cell has not only a temporal memory from C ( t ) , but also a spatial memory from u ( r ) , especially in the case of the Y-potential. The cell can “remember” the exact location at r = 0 , the point of initial contact, throughout the entire process of immunological synapse polarization. The application of the potential and memory functions here is especially crucial, since it could be the first of many types of cell memory that can be represented mathematically.
The mechanism behind the cluster formations that caused these capture and potential delays has to do with polarization and cytoskeleton rearrangement, and perhaps some insights can be gathered from the famous Turing patterns [20]. The Turing patterns are formed by a system of two reaction–diffusion equations, each describing the 2-D diffusion of a different molecule, one inhibiting the other. The pattern thus formed describes the color pigment formation of spots on animal coats. We suspect something similar occurs in the cluster formation on mast cells, where the stability of the diffusion field that forms the peripheral clusters is broken and disturbed by the increase in the central patch. This breaking point then marks the beginning of the transition from polarization to depolarization. If this was the case, then a molecule that inhibits polarization as well as a molecule that causes it would have to be simulated at the same time in 2-D. Thus, what appears to be delayed captures and potentials in 1-D would then be equivalent to the Turing pattern formation in 2-D with a central patch that increases in time but captures receptors at a constant rate at its boundary.
The idea of using 2-D pattern formation to describe polarization is not new, and much work has been conducted in this area [13,21]. The Turing-type diffusion reaction equations are widely used to form polarization patterns in these works, which consist of both simulations and experiments [22,23]. Specifically, polarization activators are represented by the active form of either cytosolic or cytoskeletal proteins, while polarization inhibitors by the inactive form [24,25]. Thus, it is the interplay between the active and inactive form that shapes the stable regions in the Turing-type diffusion field. Mathematical theories on polarization go beyond Turing-type simulations when intracellular dynamics is involved, however. Monte Carlo and phase fields are both effective tools when intracellular regions greatly affect the overall polarization [13].
In summary, the choice between delay function and Turing-type simulation has to do with simulation time and biological background considerations. The disadvantage of simulating delayed functions is that it requires an extremely long run time, even for 1-D models. With most common hardware specifications, 1-D simulations take approximately 1 h, while 2-D simulations take significantly longer, anywhere from days to weeks depending on how efficiently the code is programmed and what language is used. With Turing-type simulations, the delay functions are not required since the delay emergence can be built into the Turing-type diffusion field dynamics. Without the delay functions, the 2-D Turing-type reaction–diffusion equations take mere minutes to simulate.
The advantage of simulating delayed functions in 1-D is that it gives numerical values for capture strength and polarization time constants that would not be so readily acquired in Turing-type simulations. These constants can provide a basis for comparison between two different types of cells. For example, one can compare the capture strength of normal cells with that of cancer cells, if a difference in polarization is suspected. On the other hand, Turing-type simulations provide more information on cluster formation, but it does so at the expense of making the assumption that the polarization event is based on a two-component reaction–diffusion system (inactive and active proteins). Making such assumptions requires the researcher to have in-depth knowledge on the polarization mechanism of the interested cell type. Such assumptions are not necessary for delay function simulations. Thus, the choice of 1-D delay function vs. 2-D Turing-type simulation depends on how much experimental biological information the researcher has on the interested cell and whether polarization capture strength and time constants are required for comparison purposes.

5. Conclusions

Using convolution and cross-correlation from signal processing, we can characterize the polarization and depolarization speed of a mast cell synapse formation in 1-D. We found that polarization follows a cross-correlation-type signature, while depolarization follows a convolution-type signature, giving the mast cell a slow initial polarization and a fast initial depolarization. When transition effects are simulated by a time dependent potential term, cross-correlational delay alone is enough to describe the data points to a high degree of accuracy. Two-dimensional simulations of these events could reveal much more interesting details about the nature of cluster formation and possible cellular phase changes, and that might provide a way to link these events to molecular changes in RBL-2H3 mast cells.

Author Contributions

Conceptualization, M.C.T. and K.S.; methodology, M.C.T. and K.S.; software, M.C.T.; validation, M.C.T. and K.S.; formal analysis, M.C.T.; investigation, M.C.T. and K.S.; resources, K.S.; data curation, K.S.; writing—original draft preparation, M.C.T. and K.S.; writing—review and editing, M.C.T. and K.S.; visualization, M.C.T.; supervision, K.S. All authors have read and agreed to the published version of the manuscript.

Funding

Research reported in this publication was supported by the National Institute of General Medical Sciences of the National Institutes of Health under award number 2R15GM128166-02.

Data Availability Statement

The data presented in this study are openly available in Mendeley at http://dx.doi.org/10.17632/6KVZV95W7R.1, accessed on 5 August 2022, reference number [19].

Acknowledgments

The research reported in this publication was supported by the University of Colorado Colorado Springs (UCCS) Center for the University of Colorado BioFrontiers Institute.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Carroll-Portillo, A.; Spendier, K.; Pfeiffer, J.; Griffiths, G.; Li, H.; Lidke, K.A.; Oliver, J.M.; Lidke, D.S.; Thomas, J.L.; Wilson, B.S.; et al. Formation of a Mast Cell Synapse: FcεRI Membrane Dynamics upon Binding Mobile or Immobilized Ligands on Surfaces. J. Immunol. 2010, 184, 1328–1338. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Spendler, K.; Carroll-Portillo, A.; Lidke, K.A.; Wilson, B.S.; Timlin, J.A.; Thomas, J.L. Distribution and dynamics of rat basophilic leukemia immunoglobulin e receptors (fcεri) on planar liqand-presentinq surfaces. Biophys. J. 2010, 99, 388–397. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Spendier, K.; Lidke, K.A.; Lidke, D.S.; Thomas, J.L. Single-particle tracking of immunoglobulin e receptors (FcRI) in micron-sized clusters and receptor patches. FEBS Lett. 2012, 586, 416–421. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Spendier, K.; Sugaya, S.; Kenkre, V.M. Reaction-diffusion theory in the presence of an attractive harmonic potential. Phys. Rev. E Stat. Nonlinear Soft Matter Phys. 2013, 88, 1062142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Spendier, K.; Kenkre, V.M. Analytic solutions for some reaction-diffusion scenarios. J. Phys. Chem. B 2013, 117, 15639–15650. [Google Scholar] [CrossRef] [PubMed]
  6. Chase, M.; Spendier, K.; Kenkre, V.M. Analysis of Confined Random Walkers with Applications to Processes Occurring in Molecular Aggregates and Immunological Systems. J. Phys. Chem. B 2022, 120, 3072–3080. [Google Scholar] [CrossRef] [PubMed]
  7. Kenkre, V.M.; Spendier, K. A theory of coalescence of signaling receptor clusters in immune cells. Phys. A Stat. Mech. Its Appl. 2022, 602, 127650. [Google Scholar] [CrossRef] [PubMed]
  8. Perspectives on Learning and Memory—Google Books. Available online: https://books.google.com.hk/books?hl=en&lr=&id=ex0BAwAAQBAJ&oi=fnd&pg=PA105&dq=convolution+vs+cross+correlation&ots=218okfKUH4&sig=Ff0pb-yNjV_BiMqFKWY9FijwhVk&redir_esc=y#v=onepage&q=convolution%20vs%20cross%20correlation&f=false (accessed on 20 August 2022).
  9. A Comprehensive Introduction to Different Types of Convolutions in Deep Learning | by Kunlun Bai | Towards Data Science. Available online: https://towardsdatascience.com/a-comprehensive-introduction-to-different-types-of-convolutions-in-deep-learning-669281e58215 (accessed on 20 August 2022).
  10. Alert, R.; Trepat, X. Physical Models of Collective Cell Migration. Annu. Rev. Condens. Matter Phys. 2020, 11, 77–101. [Google Scholar] [CrossRef] [Green Version]
  11. Barcia, C.; Sanderson, N.S.; Barrett, R.J.; Wawrowsky, K.; Kroeger, K.M.; Puntel, M.; Liu, C.; Castro, M.G.; Lowenstein, P.R. T Cells’ Immunological Synapses Induce Polarization of Brain Astrocytes In Vivo and In Vitro: A Novel Astrocyte Response Mechanism to Cellular Injury. PLoS ONE 2008, 3, e2977. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Angus, K.L.; Griffiths, G.M. Cell polarisation and the immunological synapse. Curr. Opin. Cell Biol. 2013, 25, 85–91. [Google Scholar] [CrossRef] [PubMed]
  13. Rappel, W.J.; Edelstein-Keshet, L. Mechanisms of cell polarization. Curr. Opin. Syst. Biol. 2017, 3, 43–53. [Google Scholar] [CrossRef] [PubMed]
  14. Dráber, P.; Sulimenko, V.; Dráberová, E. Cytoskeleton in Mast Cell Signaling. Front. Immunol. 2012, 3, 130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Wu, J. Theory and Applications of Partial Functional Differential Equations; Springer: New York, NY, USA, 1996; Volume 119. [Google Scholar]
  16. Bellen, A.; Zennaro, M. Numerical Methods for Delay Differential Equations; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  17. Messina, E. Solving Delay Differential Equations Using Numerical Methods in Python. Available online: https://computationalmindset.com/en/mathematics/solving-delay-differential-equations-in-python-using-numerical-methods.html (accessed on 5 August 2022).
  18. Delay Differential Equations—Wolfram Language Documentation. Available online: https://reference.wolfram.com/language/tutorial/NDSolveDelayDifferentialEquations.html (accessed on 5 August 2022).
  19. Spendier, K. TIRF Microscopy Image Sequences of FcϵRI-Centric Synapse Formation in RBL-2H3 Cells Dataset, Mendeley Data V1. 2020. Available online: https://data.mendeley.com/datasets/6kvzv95w7r/1 (accessed on 5 August 2022).
  20. IPython Cookbook—12.4. Simulating a Partial Differential Equation—Reaction-Diffusion Systems and Turing Patterns. [Online]. Available online: https://ipython-books.github.io/124-simulating-a-partial-differential-equation-reaction-diffusion-systems-and-turing-patterns/ (accessed on 8 August 2022).
  21. Gierer, A.; Meinhardt, H. A theory of biological pattern formation. Kybernetik 1972, 12, 30–39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Xu, J.; Wang, F.; Van Keymeulen, A.; Herzmark, P.; Straight, A.; Kelly, K.; Takuwa, Y.; Sugimoto, N.; Mitchison, T.; Bourne, H.R. Divergent signals and cytoskeletal assemblies regulate self-organizing polarity in neutrophils. Cell 2003, 114, 201–214. [Google Scholar] [CrossRef]
  23. Pablo, M.; Ramirez, S.A.; Elston, T.C. Particle-based simulations of polarity establishment reveal stochastic promotion of Turing pattern formation. PLoS Comput. Biol. 2018, 14, e1006016. [Google Scholar] [CrossRef] [PubMed]
  24. Halatek, J.; Brauns, F.; Frey, E. Self-organization principles of intracellular pattern formation. Philos. Trans. R. Soc. B Biol. Sci. 2018, 373. [Google Scholar] [CrossRef] [PubMed]
  25. Meinhardt, H.; Gierer, A. Pattern formation by local self-activation and lateral inhibition. BioEssays 2000, 22, 753–760. [Google Scholar] [CrossRef]
Figure 1. Convolution-Type Signal Delay: (a) Using idealized values α = 0.01 s−1, c 2 = 1 s−1, and t m a x = 400 s, the memory function, C ( t ) = c 2 α   E x p ( α t ) , is replaced by C ( t s ) = c 2 α   E x p ( α ( t s ) ) in convolution time s . This has the effect of reflecting C ( t ) across s = 0 s and shifting it to the right as time t progresses. (b) the idealized rectangular receptor probability density function has the value of one at s = 0 s and drops to zero instantaneously at s = 400 s. (c) the integral 0 t C ( t s ) P ( x 0 , s ) d s gives the convolution of C ( t ) with P ( x 0 , t ) . This convolution integral describes the capture rate of the central patch. Notice the way capture is delayed: Starting from zero at time t = 0 s, the capture rate increases quickly initially, slows near tmax, and drops back down to zero the same way. Both real time t and convolution time s are measured in seconds.
Figure 1. Convolution-Type Signal Delay: (a) Using idealized values α = 0.01 s−1, c 2 = 1 s−1, and t m a x = 400 s, the memory function, C ( t ) = c 2 α   E x p ( α t ) , is replaced by C ( t s ) = c 2 α   E x p ( α ( t s ) ) in convolution time s . This has the effect of reflecting C ( t ) across s = 0 s and shifting it to the right as time t progresses. (b) the idealized rectangular receptor probability density function has the value of one at s = 0 s and drops to zero instantaneously at s = 400 s. (c) the integral 0 t C ( t s ) P ( x 0 , s ) d s gives the convolution of C ( t ) with P ( x 0 , t ) . This convolution integral describes the capture rate of the central patch. Notice the way capture is delayed: Starting from zero at time t = 0 s, the capture rate increases quickly initially, slows near tmax, and drops back down to zero the same way. Both real time t and convolution time s are measured in seconds.
Biophysica 02 00038 g001
Figure 2. Cross-Correlation-Type Signal Delay: (a) Using idealized values α = 0.0112 s−1, c 2 = 1 s−1, and t m a x = 400 s, the piecewise memory function, C ( t ) = c 2 α   E x p ( α t ) for 0 t t m a x and 0 elsewhere, is replaced by C ( t s ) = c 2 α   E x p ( α ( t s ) ) for 0 ( t s ) t m a x and 0 elsewhere in cross-correlation time s . This has the effect of reflecting C ( t ) across s = 0 s and shifting it to the right as time t progresses, while achieving maximum C ( t s ) at s = t m a x . (b) the idealized rectangular receptor probability density function has the value of one at s = 0 s, and drops to zero instantaneously at s = 400 s. (c) the integral α 0 t C ( t s ) P ( x 0 , s ) d s gives the cross-correlation of C ( t ) with P ( x 0 , t ) . The extra factor of α makes sure the integral is normalized to the maximum of one. This normalization is not required for convolution due to the form of convolution function C ( t s ) naturally decreasing as time t progresses (see Figure 1a). The cross-correlational integral describes the capture rate of the central patch. Notice the way capture is delayed: starting from zero at time t = 0 s, the capture rate increases slowly initially, quickens near tmax, and drops back down to zero the same way. Both real time t and cross-correlation time s are measured in seconds.
Figure 2. Cross-Correlation-Type Signal Delay: (a) Using idealized values α = 0.0112 s−1, c 2 = 1 s−1, and t m a x = 400 s, the piecewise memory function, C ( t ) = c 2 α   E x p ( α t ) for 0 t t m a x and 0 elsewhere, is replaced by C ( t s ) = c 2 α   E x p ( α ( t s ) ) for 0 ( t s ) t m a x and 0 elsewhere in cross-correlation time s . This has the effect of reflecting C ( t ) across s = 0 s and shifting it to the right as time t progresses, while achieving maximum C ( t s ) at s = t m a x . (b) the idealized rectangular receptor probability density function has the value of one at s = 0 s, and drops to zero instantaneously at s = 400 s. (c) the integral α 0 t C ( t s ) P ( x 0 , s ) d s gives the cross-correlation of C ( t ) with P ( x 0 , t ) . The extra factor of α makes sure the integral is normalized to the maximum of one. This normalization is not required for convolution due to the form of convolution function C ( t s ) naturally decreasing as time t progresses (see Figure 1a). The cross-correlational integral describes the capture rate of the central patch. Notice the way capture is delayed: starting from zero at time t = 0 s, the capture rate increases slowly initially, quickens near tmax, and drops back down to zero the same way. Both real time t and cross-correlation time s are measured in seconds.
Biophysica 02 00038 g002
Figure 3. Increase in capture trap size from time t 1 to a later time t 2 : Using n ( t ) = E x p ( r 2 ( R ( t ) / 2 ) 2 ) at time t = t 1 and t 2 , (a) gives the capture trap size, n ( t 1 ) , at time t 1 when the central patch radius, R ( t 1 ) , is at 1 μm, while (b) gives the trap size, R ( t 2 ) , at a later time t 2 when the central patch radius, R ( t 2 ) , has increased to 2 μm.
Figure 3. Increase in capture trap size from time t 1 to a later time t 2 : Using n ( t ) = E x p ( r 2 ( R ( t ) / 2 ) 2 ) at time t = t 1 and t 2 , (a) gives the capture trap size, n ( t 1 ) , at time t 1 when the central patch radius, R ( t 1 ) , is at 1 μm, while (b) gives the trap size, R ( t 2 ) , at a later time t 2 when the central patch radius, R ( t 2 ) , has increased to 2 μm.
Biophysica 02 00038 g003
Figure 4. Physical capture results: (a) gives the scaled convolutional capture term as a function of real time t (in seconds), while (b) gives the same for cross-correlational capture. Unlike the unphysical captures in Figure 1c and Figure 2c, the capture processes here stop when there is nothing more to capture at the trap site.
Figure 4. Physical capture results: (a) gives the scaled convolutional capture term as a function of real time t (in seconds), while (b) gives the same for cross-correlational capture. Unlike the unphysical captures in Figure 1c and Figure 2c, the capture processes here stop when there is nothing more to capture at the trap site.
Biophysica 02 00038 g004
Figure 5. The three possible potential configurations are depicted here, with potential strengths Γ = 1 μm2s−1. (a) u ( r ) = r 2 is the harmonic, or U-potential, that causes receptors to move with high speeds into the potential minimum at r = 0 . (b) u ( r ) = | r | is the V-potential with medium-speed receptors, and (c) u ( r ) = | r 1 / 2 | describes the Y-potential with low-speed receptors.
Figure 5. The three possible potential configurations are depicted here, with potential strengths Γ = 1 μm2s−1. (a) u ( r ) = r 2 is the harmonic, or U-potential, that causes receptors to move with high speeds into the potential minimum at r = 0 . (b) u ( r ) = | r | is the V-potential with medium-speed receptors, and (c) u ( r ) = | r 1 / 2 | describes the Y-potential with low-speed receptors.
Biophysica 02 00038 g005
Figure 6. Experimental Curve Fitting: R a r c ( t ) is a shifted and scaled arctan curve that resembles experimental data points taken from reference [19]. These data points represent the observed mean, but each point in the data also records an observed standard deviation error (not shown in the above figure). Since data points appear sporadic, the arctan curve is used to estimate the initial and final central patch radii. This curve is also used to estimate the time of maximum capture, t m a x . The curve fit equation is shown with its χ-squared and R-squared values. The χ-squared value is given by the formula χ 2 = i = 1 258 ( R a r c ( i ) μ ( i ) σ ( i ) ) 2 , where μ ( i ) is the observed data mean radius at time i with standard deviation σ ( i ) . The R-squared value is given by Mathematica version 13.0 as a result of fitting the curve.
Figure 6. Experimental Curve Fitting: R a r c ( t ) is a shifted and scaled arctan curve that resembles experimental data points taken from reference [19]. These data points represent the observed mean, but each point in the data also records an observed standard deviation error (not shown in the above figure). Since data points appear sporadic, the arctan curve is used to estimate the initial and final central patch radii. This curve is also used to estimate the time of maximum capture, t m a x . The curve fit equation is shown with its χ-squared and R-squared values. The χ-squared value is given by the formula χ 2 = i = 1 258 ( R a r c ( i ) μ ( i ) σ ( i ) ) 2 , where μ ( i ) is the observed data mean radius at time i with standard deviation σ ( i ) . The R-squared value is given by Mathematica version 13.0 as a result of fitting the curve.
Biophysica 02 00038 g006
Figure 7. Un-delayed capture radius: The central patch radius resulting from un-delayed capture, Ruc(t), is plotted against the arctan curve and the data points. Their χ-squared values are given by the formula χ 2 = i = 1 258 ( R ( i ) μ ( i ) σ ( i ) ) 2 , where R ( i ) is the corresponding radius at time i .
Figure 7. Un-delayed capture radius: The central patch radius resulting from un-delayed capture, Ruc(t), is plotted against the arctan curve and the data points. Their χ-squared values are given by the formula χ 2 = i = 1 258 ( R ( i ) μ ( i ) σ ( i ) ) 2 , where R ( i ) is the corresponding radius at time i .
Biophysica 02 00038 g007
Figure 8. Delayed capture radii: The central patch radii resulting from cross-correlation and convolution-type delay, R c c ( t ) and R c o ( t ) , respectively, are plotted against the arctan curve and the datapoints. Their χ-squared values are given by the formulae χ 2 ( t < 112 ) = i = 1 112 ( R ( i ) μ ( i ) σ ( i ) ) 2 and χ 2 ( t > 160 ) = i = 160 258 ( R ( i ) μ ( i ) σ ( i ) ) 2 .
Figure 8. Delayed capture radii: The central patch radii resulting from cross-correlation and convolution-type delay, R c c ( t ) and R c o ( t ) , respectively, are plotted against the arctan curve and the datapoints. Their χ-squared values are given by the formulae χ 2 ( t < 112 ) = i = 1 112 ( R ( i ) μ ( i ) σ ( i ) ) 2 and χ 2 ( t > 160 ) = i = 160 258 ( R ( i ) μ ( i ) σ ( i ) ) 2 .
Biophysica 02 00038 g008
Figure 9. Delayed cross correlation capture with time-dependent potential: The central patch radii resulting from cross-correlation-type delay capture and Fokker–Planck V-potential, Rccv(t), is plotted against the arctan curve and the data points. Their χ-squared values are given by the formula χ 2 = i = 1 258 ( R ( i ) μ ( i ) σ ( i ) ) 2 .
Figure 9. Delayed cross correlation capture with time-dependent potential: The central patch radii resulting from cross-correlation-type delay capture and Fokker–Planck V-potential, Rccv(t), is plotted against the arctan curve and the data points. Their χ-squared values are given by the formula χ 2 = i = 1 258 ( R ( i ) μ ( i ) σ ( i ) ) 2 .
Biophysica 02 00038 g009
Table 1. Time and strength constants for different types of capture delays.
Table 1. Time and strength constants for different types of capture delays.
c 2 (s−1) α (s−1) c 2 α (s−2) Γ (μm2 s−1)
Cross-Correlation0.300.0230.00690
Convolution10.0310.0310
Cross-Correlation and Potential0.260.0220.00570.19
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tsai, M.C.; Spendier, K. RBL-2H3 Mast Cell Receptor Dynamics in the Immunological Synapse. Biophysica 2022, 2, 428-439. https://doi.org/10.3390/biophysica2040038

AMA Style

Tsai MC, Spendier K. RBL-2H3 Mast Cell Receptor Dynamics in the Immunological Synapse. Biophysica. 2022; 2(4):428-439. https://doi.org/10.3390/biophysica2040038

Chicago/Turabian Style

Tsai, Ming Chih, and Kathrin Spendier. 2022. "RBL-2H3 Mast Cell Receptor Dynamics in the Immunological Synapse" Biophysica 2, no. 4: 428-439. https://doi.org/10.3390/biophysica2040038

APA Style

Tsai, M. C., & Spendier, K. (2022). RBL-2H3 Mast Cell Receptor Dynamics in the Immunological Synapse. Biophysica, 2(4), 428-439. https://doi.org/10.3390/biophysica2040038

Article Metrics

Back to TopTop