A reaction-diffusion model of the receptor-toxin-antibody interaction
- Vladas Skakauskas^{1},
- Pranas Katauskis^{1} and
- Alex Skvortsov^{2}Email author
https://doi.org/10.1186/1742-4682-8-32
© Skakauskas et al; licensee BioMed Central Ltd. 2011
Received: 15 April 2011
Accepted: 7 September 2011
Published: 7 September 2011
Abstract
Background
It was recently shown that the treatment effect of an antibody can be described by a consolidated parameter which includes the reaction rates of the receptor-toxin-antibody kinetics and the relative concentration of reacting species. As a result, any given value of this parameter determines an associated range of antibody kinetic properties and its relative concentration in order to achieve a desirable therapeutic effect. In the current study we generalize the existing kinetic model by explicitly taking into account the diffusion fluxes of the species.
Results
A refined model of receptor-toxin-antibody (RTA) interaction is studied numerically. The protective properties of an antibody against a given toxin are evaluated for a spherical cell placed into a toxin-antibody solution. The selection of parameters for numerical simulation approximately corresponds to the practically relevant values reported in the literature with the significant ranges in variation to allow demonstration of different regimes of intracellular transport.
Conclusions
The proposed refinement of the RTA model may become important for the consistent evaluation of protective potential of an antibody and for the estimation of the time period during which the application of this antibody becomes the most effective. It can be a useful tool for in vitro selection of potential protective antibodies for progression to in vivo evaluation.
Keywords
1. Background
The successful bio-medical application of antibodies is well-documented (see [1, 2] and references therein) and there is an ever-increasing interest in the application of antibodies for a mitigation of the effect of toxins associated with various biological threats (epidemic outbreaks or malicious releases) [3–5]. With the recent progress in bio-engineering, many antibodies with different affinity parameters have been generated. For a long time the main target of antibody design has been the antibody affinity. However, according to recent results [6], affinity, on its own, is a poor predictor of protective or therapeutic potential of an antibody. In fact, the treatment effect of an antibody can be described by a consolidated parameter which includes the reaction rates of the receptor-toxin-antibody kinetics and the relative concentration of reacting species [6]. As a result, any given value of this parameter determines an associated range of antibody kinetic properties and its relative concentration in order to achieve a desirable therapeutic effect. Analytical models, similar to those reported in [6], can be a useful tool for in vitro selection of potentially protective antibodies for progression to in vivo evaluation. They can significantly reduce the cost of research and development programs by optimizing associated experimental efforts. From this perspective, extension and validation of such models becomes an important goal for biomedical modelling which is partially addressed in the current study.
There are a number of ways of refining the simple kinetic model for the Receptor-Toxin-Antibody (RTA) system proposed in [6]. The possibilities include incorporating a mechanism of receptor recycling, complex pathways for toxin internalization or multiple receptor population [7]. The focus of our study is on incorporation of the diffusion effects in the theoretical framework of RTA, i.e. enhancement of the reaction RTA model [6] with the capability to account for the diffusion fluxes of reacting species [7]. Such enhancement not only enables the application of the RTA model in more realistic setting (i.e. instead of the simplified "well-mixed" approximation [6] the reaction-diffusive RTA model can describe propagation of toxin into a single cell or into a system of cells), but also provides a high fidelity estimation of the limiting uptake rate of toxin by a cell (especially when it is limited by diffusion). More importantly, the refined model allows consistent simulation of the so-called 'window of opportunity' (period of time after exposure to toxin when the application of an antibody is the most effective). We believe the two latter parameters (the limiting uptake rate and the 'window of opportunity') can become the key parameters in the optimization study for the future antibody design.
The incorporation of diffusion fluxes into the RTA model can be implemented based on a generalization of the well-known analytical framework for ligand-receptor binding [6–10]. From a mathematical point of view, the inclusion of diffusion terms into the RTA kinetic model leads to significant complications (system of nonlinear PDEs instead of system of ODEs), which usually prevent any analytical progress and implies numerical solutions. This was the main motivation for our approach to tackle the refined RTA model. The aim of this study is to numerically evaluate the protective properties of an antibody against a given toxin in the model of a spherical cell placed into a toxin-antibody solution. We consider the problem of the RTA interaction in the most general setting, when relative concentrations of species are arbitrary and all diffusive fluxes are taken into account (toxin, antibody and associated complexes). We calculate the antibody treatment efficiency parameter under various scenarios and identify the causes of time variation of this parameter.
We also study the RTA interaction in the 'Well-Mixed Solution' (WMS) model, i.e. when the solution of a toxin, antibody, and toxin-antibody complex is assumed to be uniformly mixed and homogeneously distributed in an extracellular space. In this case all diffusion fluxes disappear and the model can be described by Ordinary Differential Equations (ODE). It is worth noting that, since in such approach receptors are still confined to the single cell surface, our model is different from the "well-mixed" model proposed in [6] where all species are homogeneously distributed over the whole space. But in the case of a low internalization rate (i.e. low toxin inflow into a cell) the governing equations of these models are of the same type.
The paper is organized as follows. In Section 3 we introduce the reaction-diffusion model for RTA. The WMS model is presented in Section 4. The results are presented in Section 5. Conclusions and summarising remarks are presented in Section 6.
2. Notation
Ω - the extracellular domain, i.e. the problem domain where species diffuse and react (i.e. toxin, antibody, and toxin-antibody complex),
S_{ e } - the external surface of Ω,
S_{ c } - the cell surface (inner surface of Ω),
r_{0} - the concentration of receptors on the cell surface,
θ(t, x) - the the fraction of bounded receptors,
r_{0}θ - the concentration of the toxin-bound receptors (confined to S_{ c } ),
r_{0}(1 - θ) - the concentration of free receptors,
u_{ T } , u_{ A } , and u_{ C } - the concentrations of toxin, antibody, and toxin-antibody complex,
${u}_{T}^{0},{u}_{A}^{0},{u}_{C}^{0}$- the initial concentrations,
κ_{ T } , κ_{ A } , and κ_{ C } - the diffusivities of the toxin, antibody, and toxin-antibody complex,
k_{1}, k_{-1} - the forward and reverse constants of toxin-antibody reaction rate,
k_{2} and k_{-2} - the forward and reverse constants of toxin and receptor binding rate,
k_{3} - the rate constant of toxin internalization,
∂_{ n } - the outward normal derivative on S_{ e } or S_{ c } ,
∂_{ t } = ∂/∂t,
Δ - the Laplace operator,
ψ(t) - the antibody protection factor (a relative reduction of toxin inside a cell due to application of antibody).
3. Reaction-Diffusion Model for RTA Interaction
We disregard any excretion mechanism since we assume that it is nonsignificant over the time scales of interest (i.e. internalization time, time of toxin depletion etc).
The boundary conditions at the system above correspond to a case where initially the toxin and antibody are distributed homogeneously in the extracellular domain Ω. The boundary conditions on the outer boundary of the domain are assumed to be the constant concentrations of toxin and antibody and zero concentration of toxin-antibody complex. It is worth noting that in this case the gradients of u_{ T } , u_{ A } , u_{ C } are nonzero at the outer surface of the domain and they provide a time-dependent influx of species into Ω (with implication no conservation law for u_{ T } , u_{ A } , u_{ C } ). Indeed, in such an approach we disregard any depletion of toxin and antibody within Ω (the depletion will be taken into account in the WNS model, see below). In a practical experiment this setup can correspond to a single cell embedded into a large volume (compartment) of toxin-antibody solution, so toxin and antibody are in excess. In this context it is also worth noting that in the real biomedical scenarios the concentration of toxin is usually very low with respect to the concentration of receptor due to the high concentration of receptors on the surface of living cells and the high toxicological effect (lethal dose) of the most toxins of interest. This implies that the condition of the excess of antibody over toxin is practically relevant and are very easy to achieve (e.g. see experimental results of [11], where the concentration of ricin was about a thousand times less than the concentration of antibody), while the condition of the excess of toxin over receptor seems to be infeasible for any in vivo situation (but the latter condition still can be used in lab experiments for the model validation).
It is worth mentioning that models similar to (1)-(4) have been extensively studied in application to biouptake of pollutants by micro-organisms, cellular nutrition, heterogeneous catalysis and analytical instrumental measurements (for comprehensive review of these studies see [12–17], and references therein). Equations (1)-(4) can be presented in non-dimensional form by using scales of τ_{*} (time), l (length), and u_{*}(concentration). By substituting new variables, $x=l\stackrel{\u0304}{x}$, $t={\tau}_{*}\stackrel{\u0304}{t}$, ${r}_{0}=l{u}_{*}{\stackrel{\u0304}{r}}_{0}$, ${u}_{T}={u}_{*}{\u016b}_{T}$, ${u}_{A}={u}_{*}{\u016b}_{A}$, ${u}_{C}={u}_{*}{\u016b}_{C}$, ${u}_{T0}={u}_{*}{\u016b}_{T}^{0}$, ${u}_{A0}={u}_{*}{\u016b}_{A}^{0}$, ${\stackrel{\u0304}{k}}_{1}={\tau}_{*}{u}_{*}{k}_{1}$, ${\stackrel{\u0304}{k}}_{2}={\tau}_{*}{u}_{*}{k}_{2}$, ${\stackrel{\u0304}{k}}_{-1}={\tau}_{*}{k}_{-1}$, ${\stackrel{\u0304}{k}}_{-2}={\tau}_{*}{k}_{-2}$, ${\stackrel{\u0304}{k}}_{3}={\tau}_{*}{k}_{3}$, ${\stackrel{\u0304}{\kappa}}_{T}={\tau}_{*}{\kappa}_{T}{l}^{-2}$, ${\stackrel{\u0304}{\kappa}}_{A}={\tau}_{*}{\kappa}_{A}{l}^{-2}$, ${\stackrel{\u0304}{\kappa}}_{C}={\tau}_{*}{\kappa}_{C}{l}^{-2}$ into (1)-(4) we can deduce the same system, but only in non-dimensional variables. Therefore, for simplicity in what follows, we treat system (1)-(4) as non-dimensional.
By definition 0 ≤ ψ ≤ 1 with the lower values of ψ corresponding to the more profound therapeutic effect of antibody treatment.
So that ψ_{1} can be expressed in terms of only one 'bulk' variable ${u}_{T}^{\mathsf{\text{sat}}}\ge 0$. Indeed, the value of ψ^{sat} can be appreciably affected by the diffusivities of species, since κ_{ T } , κ_{ A } , κ_{ C } determine the saturation value ${u}_{T}^{\mathsf{\text{sat}}}$ by virtue of Eqs. (1)-(4).
4. WMS Model for RTA Interaction
Here a dot is placed over the variables to represent a time derivative; k_{4} = S_{ c }/V_{Ω}, where S_{ c } and V_{Ω} are the area of cell and the extracellular volume. For instance, for a spherical cell of radius ρ_{ c } , V_{Ω} is a domain between the cell and a concentric sphere of radius ρ_{ e } > ρ_{ c } , ${V}_{\Omega}=\frac{4}{3}\pi \left({\rho}_{e}^{3}-{\rho}_{c}^{3}\right)$, ${S}_{c}=4\pi {\rho}_{c}^{2}$, and ${k}_{4}=3{\rho}_{c}^{2}/\left({\rho}_{e}^{3}-{\rho}_{c}^{3}\right)$. For a simple model of cell culture (a uniformly distributed system of cells) the average density of cell distribution, n, is approximately equal to $3\u2215\left(4\pi {\rho}_{e}^{3}\right)$, so we can treat the 'external' scale ρ_{ e } as the size of a compartment occupied by an individual cell in the culture. From this perspective, the dependence of ψ(ρ_{ e } ) presented below can provide insight into the dependence of ψ on the cell packing density in the culture since ρ_{ e } ≈ [3/(4πn)]^{1/3} (see below).
The WMS model (8)-(11) is worth comparing with the model of the RTA interaction proposed in [6] (a kinetic model of uniformly distributed chemical species and cells). Despite these models being essentially different in their geometrical setting (in our case the receptors are still confined to a surface of a single cell), their governing equations become similar in the case when toxin inflow into a cell can be neglected (i.e. low internalization rate); the latter case seems to be very typical for many practical situations [7]. The WMS model (8)-(11) being a system of ODEs is much easier to analyze and solve numerically than the full RTA model (1)-(4) but indeed the WMS model cannot be used for estimating the effect of diffusivity of species on the protective properties of antibody (since it contains no diffusivity parameters).
where ε = K_{2}/K_{1}, K_{1} = k_{-1}/k_{1}, K_{2} = k_{-2}/k_{2}, R_{0} = r_{0}k_{4} (the same equation is given in [6] for the "well mixed" model). Then the solution of this equation enables the calculation of protection factor ψ_{2} = ψ^{sat} by means of Eq. (7).
where $F\left(x,y\right)=\left({q}_{1}-\sqrt{{q}_{1}^{2}-4{q}_{2}y}\right)\u2215\left(2{q}_{2}\right)$, q_{1} = K_{2} + εx - (ε - 2)y and q_{2} = q_{1} - (εK_{2} + y).
Comparison of saturation values of ψ for WMS model: ψ^{sat} = ψ_{2} (12) and (7), ψ^{sat} = ψ_{3} (13), and ψ^{sat} = ψ_{4}, where ψ_{4} is estimated from the solution of (8)-(11) and (7) at t = 10 000 s
k _{1} | k _{2} | ψ _{2} | ψ _{3} | ψ _{4} | |
---|---|---|---|---|---|
k_{3} = 0 | k_{3} = 0.000033 | ||||
0.013 | 0.0125 | 0.215524 | 0.215524 | 0.216026 | 0.206474 |
0.013 | 0.025 | 0.345686 | 0.345708 | 0.345903 | 0.332632 |
0.013 | 0.05 | 0.508760 | 0.508754 | 0.508767 | 0.493704 |
0.13 | 0.0125 | 0.027219 | 0.027220 | 0.027426 | 0.025913 |
5. Numerical Results
Values of parameters used in calculations
Parameter | Dimensional value | Non-dimensional value |
---|---|---|
k _{1} | 1.3 ·10^{5} M^{-1} s^{-1} | 1.3 · 10^{-2} |
k _{2} | 1.25 ·10^{5} M^{-1} s^{-1} | 1.25 · 10^{-2} |
k _{-1} | 1.4 ·10^{-4} s^{-1} | 1.4 · 10^{-4} |
k _{-2} | 5.2 ·10^{-4} s^{-1} | 5.2 · 10^{-4} |
k _{3} | 3.3 ·10^{-5} s^{-1} | 3.3 · 10^{-5} |
κ _{ T } | 10^{-6} cm^{2} s^{-1} | 10^{-2} |
κ _{ A } | 10^{-6} cm^{2} s^{-1} | 10^{-2} |
κ _{ C } | 10^{-6} cm^{2} s^{-1} | 10^{-2} |
ρ _{ c } | 10^{-3} cm | 10^{-1} |
ρ _{ e } | 2 · 10^{-3}, 5 · 10^{-3} cm | 2, 5 |
${u}_{A}^{0}$ | 6.02 · 10^{-13} cm^{-3} | 1 |
${u}_{T}^{0}$ | 3.01 · 10^{-13}, 6.02 · 10 ^{-14} cm^{-3} | 0.5, 0.1 |
Comparison of saturation values of ψ for model (1)-(4): ψ^{sat} = ψ_{1} and ψ^{sat} = ψ_{5}, where ψ_{1} is determined by (6) and (7) with ${u}_{T}^{\mathsf{\text{sat}}}$ estimated by (1)-(4) at t = 1000 s, while ψ_{5} is determined by (7) with θ^{sat} estimated by (1)-(4) at t = 10 000 s
κ_{ T }(κ_{ A }= 10^{-2}) | ψ _{1} | ψ _{5} | κ_{ A }(κ_{ T }= 10^{-3}) | ψ _{1} | ψ _{5} |
---|---|---|---|---|---|
10^{-2} | 0.9339 | 0.9339 | 10^{-1} | 0.1342 | 0.1345 |
5 · 10^{-3} | 0.8433 | 0.8433 | 10^{-2} | 0.1480 | 0.1483 |
10^{-3} | 0.1480 | 0.1483 | 10^{-3} | 0.3729 | 0.3726 |
10^{-4} | 0.0047 | 0.0034 | 10^{-4} | 0.9801 | 0.9801 |
where Λ is the conductance of the system (flux-concentration ratio), u_{ T } (t) is the concentration of toxin on the outer boundary of Ω, viz. ${u}_{T}\left(t\right)={u}_{T}^{0}$ for the boundary condition of constant concentration or ${u}_{T}\left(t\right)={u}_{T}^{0}\mathsf{\text{exp}}\left(-t\u2215{\tau}_{d}\right)$ for the no-flux boundary condition, κ_{*} is the effective diffusion of the toxin, τ_{ d } is the depletion time of toxin in the bulk, K_{*} = R_{0}/(R_{0} + K_{1}) [6]. It can be seen that the parameter κ_{*} and depletion time τ_{ d } (if the depletion of toxin is significant) become two 'aggregated' parameters that can be used to comprehensively characterize the influence of an antibody on toxin transport in the model of spherical cell.
which is called bioavalability number [15] and can be used to characterized the regime of toxin uptake by the cell [15, 16]. If L ≪ 1 the uptake flux is fully controlled by the internalization process, while in the opposite case L ≫1 it is controlled by diffusion. Note that for the case of ricin competitive binding to cell receptors and the mono-clonal antibody 2B11 the value of parameter L ≈ 10^{-2}, i.e. flux is mostly controlled by internalization process. Importantly, even in the case of diffusion dominated flux the transport of toxin can be characterized by a rich variety of regimes that are parameterized based on the so-called degree of lability, so these regimes correspond to the different asymptotical values of parameters κ_{*}, τ_{ d }[15–17].
where $\gamma =1\u2215{\tau}_{d}-1\u2215{\tau}_{d}^{0},{\tau}_{d}^{0}$ is the depletion time of toxin without antibody, L_{0} = K_{*}k_{3}ρ_{ c }/κ_{ T } .
Equation (16) is an analogue of expression (13) that accounts for diffusion effects and toxin depletion. We observe that, depending on a value of the parameter γ, the asymptotical behavior of the protection factor can be either zero (γ < 0), infinity (γ > 0) or non-zero constant (γ = 0). For the diffusion controlled flux, L_{0} ≫ 1 and ψ_{*} = p while for the kinetically controlled regime L_{0} ≪ 1 and ψ_{*} = 1. This implies that by changing the diffusivity of the reacting species (i.e. by introducing an antibody) it is possible (at least in theory) to control the behavior of the protection factor ψ(t). The case of the constant toxin influx (i.e. no depletion) simply corresponds to γ = 0. We observed most of these scenarios in our numerical simulations (see below).
that can be significantly different from the diffusivity of a toxin u_{ T } . It is worth noting that antibody diffusivity does not appear in this expression. The "limit of mean diffusivity" for κ_{*} given by (17) occurs only for the system with κ_{ T } ≠ κ_{ C }[16, 17] (which is usually the case because of a difference in molecular weights).
Some analytical models for the calculation of the toxin depletion times ${\tau}_{d}^{0}$, τ_{ d } have been proposed [16]. They are quite involved, and for details, we refer the reader to the original publications. The results [16] clearly demonstrate that the parameter $\gamma =1\u2215{\tau}_{d}-1\u2215{\tau}_{d}^{0}$ in (16) can depend on the 'external' scale ρ_{ e } (i.e. size of the cell 'compartment') in a quite convoluted way. As was mentioned above, the scale ρ_{ e } can be approximately related to the packing density of cells in a culture (ρ_{ e } ≈ [3/(4πn)]^{1/3}), so plots ψ(ρ_{ e } ) can be also interpreted as a simple qualitative illustration of the effect of variation in packing density n.
The plots ψ(ρ_{ e } ) in Figures 1, 2, 3 depict the dependence of the antibody protection factor ψ on the radius of external surface ρ_{ e } (i.e. a size of the cell compartment) and on toxin diffusivity κ_{ T } . We calculated ψ for two values of ρ_{ e } (ρ_{ e } = 2 and ρ_{ e } = 5), for two values of κ_{ T } (κ_{ T } = 10^{-2} and κ_{ T } = 10^{-3}), and three values ${u}_{T}^{0}$ (${u}_{T}^{0}=0.5$ Figure 1, ${u}_{T}^{0}=0.3$ Figure 2, ${u}_{T}^{0}=0.1$ Figure 3).
We believe that the analytical results (16) discussed above and the numerical examples similar to those presented in Figures 1, 2, 3 may be important for either the planning of experiments (especially in cell culture) or for the correct interpretation of experimental data, since they provide a simple estimation for the amplitude of the observable effect (protection factor) and for the timescale during which this effect can occur (~ 1/γ).
The results depicted in Figures 1, 2, 3 also provide an illustrative example of the main finding of our study: the time evolution of the protection factor ψ(t) may switch from monotonic to markedly non-monotonic behavior with a variation of diffusion parameters of the RTA model. This phenomenon is in line with the theoretical framework proposed in [15–17]) and was observed frequently in our simulations.
Following the well-established application of the concept of lability to the spherical cell model [15–17], an incorporation of diffusion effects into our model enabled the simulation of a new phenomenology, which may occur in the RTA system. For instance, with system (1)-(4) we were able to model competitive behavior of the reaction and diffusion fluxes. As was mentioned above, the latter often manifests itself in rather convoluted (non-monotonic) dependence of concentration of species and their diffusion coefficients, see Figures 4 and 5. The plots in Figures 4 and 5 depict a variety of scenarios for time evolution of ψ for the different diffusivity of toxin and antibody (other parameters were the same). We can clearly see a switch from monotonic to non-monotonic behavior as we decrease diffusivity of toxin κ_{ T } (Figure 5). The cases of non-monotonic behavior with a profound minimum of ψ(t) provide revealing examples of the practically important concept of a 'window of opportunity' discussed in the Background. Once the function ψ(t) moves far away from its minimal value, the 'blocking' effect of an antibody markedly decreases. We observe that the 'window of opportunity' is very scenario-dependent and it can be easily estimated from the plots similar to those presented in Figure 5.
The plots presented in Figures 1, 2, 3, 4, 5 illustrate detailed insights into the transport process associated with the different behavior of the protective function ψ(t) owing to the introduction of antibody into the RTA system. From these results we can see that the non-monotonic behavior is caused by the non-monotonic supply (transport) of toxin across the compartment to the cell surface; this transport can be affected by changing the diffusivity of species. The main conclusion from our numerical results is that the relative diffusivity of species can be used to control the effect of antibody treatment during a short time after the exposure to a toxin (usually a few minutes).
The plots presented in Figures 6 and 7 show the dependence of ψ on k_{1} and k_{2} for the WMS model and demonstrate a possible switch of ψ from monotonic to non-monotonic behavior as k_{2} grows. Calculations show that ψ(t) depends very weakly on ρ_{ e } and this an indication that the reactions in the surface layers around the cell (diffusive and reaction) provide the dominant contribution to the transport properties of the system. This result is depicted in Figures 6 and 7 where all curves are given for ρ_{ e } = 2 and ρ_{ e } = 5. As time increases, function ψ(t) tends to the diffusivity-dependent asymptote ψ^{sat} = ψ_{1} for system (1)-(4) and to the reaction-dependent asymptote ψ^{sat} = ψ_{3} for the WMS model irrespective of its short-time behavior (which is indeed controlled by (16)). The diffusive dependency of the saturation limit of the system (1)-(4) becomes evident if we recall that at the steady state, the flux of internalized toxin (i.e. flux across S_{ c } ) should be compensated by the diffusion influx across the outer surface S_{ e } . In order to validate this asymptotic behavior we computed a rich set of scenarios. These results are presented in Figures 4, 5, 6, 7 and in Tables 1 and 3. In Table 3 we compare ψ_{1} and ψ_{5} where ψ_{1} is determined by (6) and (7) with ${u}_{T}^{\mathsf{\text{sat}}}$ estimated by (1)-(4) at t = 1000 s, while ψ_{5} is determined by (7) with θ^{sat} estimated (1)-(4) at t = 10 000 s. We observe that function ψ(t) converges to an asymptotic value, but this convergence can be rather slow.
As was suggested by one of the anonymous referees, the observable strongly non-monotonic behavior of parameter ψ(t) in some of our modeling scenarios can possibly be explained by applying the concept of dynamic speciation to the formation of a toxin-antibody complex [15–17]). In the diffusion-controlled regime the dynamic speciation (i.e. the fast toxin-antibody kinetics over diffusion time) can lead to the significant contribution to the toxin flux towards the cell and (under condition κ_{ C } < κ_{ T } ) can even cause a 'retardation' effect [15]. After some estimations we found this hypothesis quite reasonable. For a cell size of ρ_{ c } ≈ 10^{-5} m the diffusion time is τ_{ κ } ≈ 0.2 s for κ ≈ 1 · 10^{-9} m^{2}s^{-1}. The estimation for equilibration time τ_{ e } was derived from the rigorous theoretical framework proposed in [23] for competitive binding system (application of this framework to the toxin-receptor and toxin-antibody binding can be found in [6]). Indeed the equilibration time τ_{ e } is a strong function of the toxin concentration; it rapidly decreases as the toxin concentration increases (reacting species can faster find each other to form a complex). If as a reference point we assume that the value of parameters correspond to the scenario of binding of ricin to receptor and to the antibody then for the toxin concentration T = 10 pM the reaction time is of order of 10 s. By further increasing the toxin concentration (five times in our simulations) it appears that we approach the transition threshold from the 'inert' to the 'dynamic' complex, so the toxin-antibody complex starts contributing to the diffusion flux. A more challenging task was to identify the scenarios where this additional contribution can be appreciable, since the total flux is mainly internalization-controlled. Nevertheless, if we recall that diffusivities and concentration of species are varied by an order of hundreds, then reaching the diffusion-controlled regime in some of our simulations looks quite feasible. A more detailed interpretation of our numerical results with the concept of dynamic speciation would require additional numerical calculations (careful estimations of equilibration time of complex for each scenario) and is outside of the scope of the current study.
6. Concluding Remarks
In summary, we have refined the RTA model developed in [6] by incorporating diffusion of reacting species in the extracellular space. By solving numerically the system of nonlinear PDEs of the model we managed to simulate a rich variety of reaction-diffusion processes that may occur in the RTA system. For various combinations of parameters (rates of reactions, diffusivity and initial concentrations) we estimated the effect of antibody on the toxin penetration into a cell and expressed the effect of the antibody treatment in terms of a non-dimensional protection factor (relative reduction of toxin concentration within a cell). We demonstrated that this factor can be a significantly non-monotonic function of time and its behavior depends on an interplay between diffusive and reaction processes in the RTA model. We also examined the time evolution of the protection factor and found that it eventually tends to a diffusivity-dependent asymptotic limit, but the convergence to this asymptote may take significant time. From this perspective, the refinement of the RTA model proposed in the present study becomes important for the consistent evaluation of protective potential of an antibody and for the estimation of the time period during which the application of this antibody becomes the most effective.
The selection of the rate constants for numerical simulations was motivated by data reported in the literature [11, 19–22], with the significant ranges of variability to provide a simple sensitivity analysis for the system under consideration. The values for other constants (i.e. diffusivity) were selected based on similarity with other models [8–10]. The chosen values of parameters enable an illustrative demonstration of a rich variety of regimes of the evolution that could occur in the RTA system. These regimes are similar to ones occurring in electrochemistry and ecological studies (performance of microelectrodes and toxin uptake by microorganisms). Further validation of the proposed model with a particular set of experimental data on toxin-neutralising antibodies (e.g. [11, 21]) would require a separate study. Such a study would include an application of a data fitting algorithm that accounts for the experimental data uncertainty as well as some additional assumptions about relationships of the model predictions (concentration of species, protection factor) with the observable quantities (i.e. cellular viability). The latter assumptions may significantly affect the experimental data fit and the evaluation of predictive skills of the proposed model. We will report on such study in a separate publication.
8. Authors contributions
AS suggested a simplified model of the RTA system. VS extended the model and proposed a consistent analytical framework for the modeling and simulation of the RTA interaction. PK implemented the implicit finite-difference scheme for the numerical solution of the refined PDE system. All authors contributed to the writing of the article, and read and approved the final manuscript.
Declarations
9. Acknowledgements
We thank Dr Peter Gray (DSTO) for valuable discussions and support and an anonymous referee who drew our attention to the problem of estimation of the biouptake of pollutants by micro-organisms which is treated by the similar analytical methods.
Authors’ Affiliations
References
- Oral HB, Ozakin C, Akdis CA: Back to the future: antibody-based strategies for the treatment of infectious diseases. Mol Biotechnol. 2002, 21: 225-239. 10.1385/MB:21:3:225.View ArticlePubMedGoogle Scholar
- Lobo ED, Hansen RJ, Balthasar JP: Antibody pharmacokinetics and pharmacodynamics. J Pharm Sci. 2004, 93: 2645-2668. 10.1002/jps.20178.View ArticlePubMedGoogle Scholar
- Kaufmann AF, Meltzer MI, Schmid GP: The economic impact of a bioterrorist attack: are prevention and postattack intervention programs justifiable?. Emerging Infectious Diseases. 1997, 3: 83-94. 10.3201/eid0302.970201.PubMed CentralView ArticlePubMedGoogle Scholar
- Rainey GJA, Young JAT: Antitoxins: Novel strategies to target agents of bioterrorism. Nat Rev Microbiol. 2004, 2: 721-726. 10.1038/nrmicro977.View ArticlePubMedGoogle Scholar
- Nehal M, Michelle C, Juan L, Steven J, Steven P, Giovanni DA, Linda N, Michael B, George S: Leslie, A high-affinity monoclonal antibody to anthrax protective antigen passively protects rabbits before and after aerosolized. Infect Immun. 2005, 73: 795-802. 10.1128/IAI.73.2.795-802.2005.View ArticleGoogle Scholar
- Skvortsov A, Gray P: Modeling and simulation of receptor-toxin-antibody interaction. Proc. 18th World IMACS/MODSIM Congress, Australia. 2009, 185-191.Google Scholar
- Truskey GA, Yuan F, Katz DF: Transport Phenomena in Biological Systems. 2009, Prentice Hall, SecondGoogle Scholar
- Goldstein B, Dembo M: Approximating the effects of diffusion on reversible reactions at the cell surface: Ligand-receptor kinetics. Biophys J. 1995, 68: 1222-1230. 10.1016/S0006-3495(95)80298-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Coppey M, Berezhkovskii AM, Sealfon SC, Shvartsman SY: Time and length scales of autocrine signals in three dimensions. Biophys J. 2007, 93: 1917-1922. 10.1529/biophysj.107.109736.PubMed CentralView ArticlePubMedGoogle Scholar
- Ghosh S, Gopalakrishnan M, Forsten-Williams K: Self-consistent theory of reversible ligand binding to a spherical cell. Phys Biol. 2007, 4: 344-View ArticleGoogle Scholar
- Prigent J, Panigai L, Lamourette P, Sauvaire D, Devilliers K: Neutralising antibodies against ricin toxin. PloS ONE. 2011, 6: e20166-10.1371/journal.pone.0020166.PubMed CentralView ArticlePubMedGoogle Scholar
- Berg HC: Random Walks in Biology. 1993, Princeton University PressGoogle Scholar
- Ruthven DM: Principles of Adsorption and Adsorption Processes. 1984, John Wiley & Sons, USAGoogle Scholar
- van Leeuwen HP: Revisited: The conception of lability of metal complexes. Electroanalysis. 2001, 13: 826-830. 10.1002/1521-4109(200106)13:10<826::AID-ELAN826>3.0.CO;2-J.View ArticleGoogle Scholar
- Pinheiro J, Galceran J, van Leeuwen HP: Metal speciation dynamics and bioavailability: Bulk depletion effects. Environ Sci Technol. 2004, 38: 2397-2405. 10.1021/es034579n.View ArticlePubMedGoogle Scholar
- van Leeuwen HP, Town RM, Buffle J, Cleven RFMJ, Davison W, Puy J, van Riemsdijk WH, Sigg L: Dynamic speciation analysis and bioavailability of metals in aquatic systems. Environ Sci Technol. 2005, 39: 8545-8556. 10.1021/es050404x.View ArticlePubMedGoogle Scholar
- Pinheiro JP, Minor M, van Leeuwen HP: Metal speciation dynamics in colloidal ligand dispersions. Part 2: Electrochemical lability. J Elec Chem. 2006, 587: 284-292. 10.1016/j.jelechem.2005.11.022.Google Scholar
- Samarskii AA: The Theory of Difference Schemes. 2001, Marcel Dekker, New YorkView ArticleGoogle Scholar
- Sandvig K, Olsnes S, Pihl A: Kinetics of binding of the toxic lectins abrin and ricin to surface receptors of human cells. J Biol Chem. 1976, 251: 3977-3984.PubMedGoogle Scholar
- McGuinness CR, Mantis NJ: Characterization of a novel high-affinity monoclonal immunoglobulin G antibody against the ricin B subunit. Infect Immun. 2006, 74: 3463-3470. 10.1128/IAI.00324-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu F, Fan S, Martiniuk F, Pincus S, Muller S, Kohler H, Tchou-Wong KM: Protective effects of anti-ricin A-chain antibodies delivered intracellularly against ricin-induced cytotoxicity. World J Biol Chem. 2010, 1: 188-195. 10.4331/wjbc.v1.i5.188.PubMed CentralView ArticlePubMedGoogle Scholar
- Lectures Notes in Immunology: Antigen-antibody interactions, University of Pavia. 2011, [http://nfs.unipv.it/nfs/minf/dispense/immunology/lectures/files/antigens_antibodies.html]
- McPherson RA, Zettner A: A mathematical analysis of the incubation time in competitive binding systems. Analyt Biochem. 1975, 64: 501-508. 10.1016/0003-2697(75)90460-1.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.