The toxins associated with infectious diseases are potential targets for inhibitors which have the potential for prophylactic or therapeutic use. Many antibodies have been generated for this purpose, and the objective of this study was to develop a simple mathematical model that may be used to evaluate the potential protective effect of antibodies. This model was used to evaluate the contributions of antibody affinity and concentration to reducing antibody-receptor complex formation and internalization. The model also enables prediction of the antibody kinetic constants and concentration required to provide a specified degree of protection. We hope that this model, once validated experimentally, will be a useful tool for in vitro selection of potentially protective antibodies for progression to in vivo evaluation.
1. Introduction
Passive immunization using antibodies has been used successfully for treatment and prophylaxis of infectious disease in humans, and there is increasing interest in the use of antibodies for treatment of infectious diseases that may be used as terrorist weapons, but for which the risk is not sufficiently high to justify preventive vaccination of a large civilian population (see [1–4] and references therein). Toxins are an important potential target for designing therapies against these threats and a broad range of approaches have been taken to develop inhibitors that may be of prophylactic or therapeutic use [1, 5].
Antibody engineering techniques allow affinity maturation of antibodies, and these techniques are being exploited to produce inhibitors for a number of toxins [6, 7]. The emphasis of this approach is on producing reagents with high affinity, based on the proposition that higher affinity will provide better protection.
However affinity, by itself, is a poor predictor of protective or therapeutic potential. Antibodies with high in vitro affinity for toxins do not automatically confer protection in vivo [8, 9] and may exacerbate the toxicity [10, 11]. The effects of using multiple antibodies with high affinities may be additive [12] or synergistic [8] or without effect [9]. In addition, epitope specificity [13], antibody titre [14–18], and dissociation rate [19] have been correlated with protection.
Toxins are produced by a number of plants, animals and microorganisms. Toxins may act at the cell surface and either damage the cytoplasmic membrane or bind to a receptor and act via transmembrane signalling subsequent to that binding [20]. Alternatively, toxins may cross the cell membrane and act on intracellular targets [20]. For example, anthrax lethal toxin, ricin and cholera toxin bind to a cell surface receptor and make use of cellular membrane trafficking to enter the cell [21, 22].
The objective of this study is to develop a simple mathematical model that may be used to predict the optimum antibody parameters (kinetic constants and concentration) needed to inhibit the binding of the toxin to its receptor. These predictions may be used to select candidate antibodies for progression to in vivo evaluation and to assess the potential value of affinity enhancement.
This paper is an extension to our previous work [23]. In the model presented in the following we explicitly take into account the process of toxin internalization and diffusive fluxes around the cell.
2. Model
The kinetic model describing the interactions of toxins with cell receptors can be formulated based on the well-known analytical framework for ligand-receptor binding. The models of this process have been studied for many years and a vast amount of literature has accumulated on this subject (see [24–28] and references therein).
When a toxin diffuses in the extracellular environment and binds to the cell surface receptors, the toxin concentration will vary in both space and time. Any rigorous description of this process would entail a system of Partial Differential Equations (PDE), which couples extracellular diffusion with reaction kinetics of the cell surface. The resulting system of PDE is nonlinear and too complex to be treated analytically. This complexity makes any comprehensive study of parameter optimization unfeasible. From another perspective, it is well known that under some rather broad conditions (see [24–28] and references therein) the reaction-diffusion system of the ligand-receptor binding can be well approximated by a system of Ordinary Differential Equations in which the spatial variability of the process is simulated by different concentrations of species in initially predefined spatial domains (called compartments). Although this compartment model is significantly simpler than the initial reaction-diffusion system, it still allows a consistent description of reaction-diffusion transport in underlying system [25, 26, 28]. In the current paper we use the compartment-model approach for our analytical study and numerical simulations.
To begin, we consider the following simple model. The toxin, T, binds reversibly to cell surface receptors, R, with a forward rate k1 and a reverse rate k-1 to form the toxin-receptor complex CR which is then slowly internalized at a rate k3. The neutralizing antibody binds competitively to the toxin with on and off rates of k2 and k-2, respectively. The antibody-toxin complex, CA, remains in the extracellular space (see Figure 1).
Schematic representation of the model for receptor-toxin-antibody interaction.
We can easily write an equation for the toxin-receptor binding (namely, without antibody being present). For a spherical cell of radius a with the toxin binding to its surface [24–28],
(1)dCRdt=kfeRT+kreCR,
where CR is the concentration of the bound receptors (toxin-receptor complexes), R is the concentration of receptors, and T is the bulk toxin concentration (i.e., far from the cell surface) and is assumed to be spatially uniform. The effective forward and reverse rate coefficients are defined by [24–28]
(2)kfe=γk1,kre=γk-1,
where k1, k-1 are intrinsic reaction rates, kD=4πaD is the diffusion rate, D is the diffusivity of toxin in the extracellular space, and γ=1/(1+Rk1/kD)≤1 [25–27].
The bulk concentration of toxin T is mainly driven by the binding to antibody. Therefore, in this case we can write an equation system similar to (1) but without any “diffusive” modification of the intrinsic rate constants:
(3)dCAdt=k2AT+k-2CA,
where CA is the concentration of toxin-antibody complexes and A is the concentration of antibody.
The process of toxin internalization is phenomenologically introduced into our model by the following equation:
(4)dTidt=k3CR,
where Ti is the concentration of internalized toxin. The corresponding term should be included in (1), so we arrive at modified expression for kre(5)kre=γk-1-k3.
The systems (1), (3), and (4) should be supplemented with three conservation laws for concentrations of R, T, and A:
(6)R0=R+CR,(7)A0=A+CA,(8)T0=T+CT+CA+Ti,
where R0, T0, and A0 are the initial concentrations.
Equations (1), (3), (4), and (6)–(8) form a framework for our analysis. This is a system of nonlinear ODE (because of conservation laws (6)–(8) and because effective rates kfe, kre are functions of the receptor concentration). It can be easily solved numerically and also allows some analytical progress (see the following). If parameter γ≪1 (and this is the case in many practical situations), then this model can be reduced to the “well-mixed” kinetic model with constant kinetic rates [23].
It is worth emphasizing that the aim of our analytical framework is to develop a simple but scientifically rigorous model that may be used to predict the optimum antibody kinetic properties and concentration required to achieve a desired protective effect rather than develop a detailed, biologically accurate model that captures all the details of the toxin internalization process. Therefore, the model does not take into account the pharmacokinetics of the toxin-antibody complex [11] or receptor internalization and recycling [29, 30]. k3 is a lumped constant that describes all the processes that result in the appearance of the free toxin in the intracellular space [31]. Wiley and Cunningham [32] and Shankaran et al. [33] have also developed mathematical models of this type of process.
We are particularly interested in the behaviour of the model under conditions most likely to reflect the real biological situation, that is, toxin concentration much lower than the concentration of receptors (T0/R0≪1).
Testing of the model was carried out using COPASI (software application for simulation and analysis of biochemical networks and their dynamics [34]) and the kinetic parameters for the binding of ricin to its receptor and its internalization [35] and competition by the monoclonal antibody 2B11 [8]. The kinetic parameters used are shown in Table 1. The value of k3 used is that determined by Sandvig et al. [35] to be the rate of irreversible binding of ricin to HeLa cells. For simplicity, the simulation was carried out using all reactions taking place in the same compartment.
Kinetic constants used in numerical simulations (the binding of ricin to its receptor and the monoclonal antibody 2B11).
Reaction
Value
k1
1.3·105 M−1s−1
k-1
1.4·10-2 s−1
k2
1.25·105 M−1s−1
k-2
5.2·10-4 s−1
k3
3.3·10-5 s−1
To illustrate the model, we used toxin and receptor concentrations based on cell culture studies carried out in our laboratory. These typically use a cell concentration of 1·104 cells per 100 μL experiment and a ricin concentration of 10 pM. Assuming 3·107 receptors/cell [35], the receptor concentration is approximately 5 nM.
3. Analytical Results3.1. Cell Surface Binding
Initially we derive some analytical results for toxins that act at the cell surface and are not internalized; that is, we set k3=0 in (4). At equilibrium d/dt=0 and from (1) and (3) we can write
(9)CR=RTK1,CA=ATK2,
where K1=k1/k-1, K2=k2/k-2 are the association constants for the toxin binding to the receptor and antibody, respectively. It is worth noting that the parameter γ (diffusive correction of the intrinsic reaction rates) disappears from (9), so in this case the analytical results are identical to ones derived using the “well-mixed” approximation [23].
In order to simplify notations, we denote by z and y the equilibrium concentrations of the toxin-receptor and toxin-antibody complexes; that is,
(10)z=[CR]eq,y=[CA]eq.
From (9) and conservation laws (6)–(8) the following closed equation for z can be derived:
(11)(R0-z)(T0-z-y)-K1z=0,(12)y=A0ϵzR0-z(1-ϵ),
where ϵ=K1/K2.
Equation (11) can be written in a more conventional form of a cubic equation as follows:
(13)a3z3+a2z2+a1z+a0=0,
where
(14)a3=ϵ-1,a2=(1-ϵ)C0+ϵA0+R0,a1=-R0(C0+A0+(1-ϵ)T0),a0=T0R02,
and C0=R0+K1.
It is well known that (13) has a closed-form analytical solution (Cardano's formula [36]), which in our case provides a consistent way to derive exact solutions for the proposed model. Unfortunately these solutions still involve rather cumbersome expressions, which require further simplifications in order to be used in practical situations. In the following we present another approach that explicitly employs the smallness of ratio T0/R0≪1 and leads to a simple analytical expression for the protective properties of the antibody.
We observe that in the absence of antibody (i.e., A0=0), (11) is an elementary quadratic equation that has two roots. If we impose the obvious constraint z→0 as T0→0, then there is only one solution, which we designate as z0:
(15)z0=C02[1-(1-4R0T0C02)1/2].
Under the condition T0/R0≪1, this can be simplified to
(16)z0≈R0T0C0,C0=R0+K1.
Let us now evaluate the effect of adding an antibody. From a mathematical point of view this effect (i.e., change of z under condition A0>0) is captured entirely by the term y in (11), so our aim is to provide a reasonable analytical estimation of this term.
From (12) and based on our initial assumption of low toxin concentration (T0/R0≪1), we can deduce the following simple estimate y≈ϵzA0/R0. This then leads to a modified form of (11) as follows:
(17)(R0-z)(T0-z)-K*z=0,
where
(18)K*=K1+ϵA0.
We can see that this is the same form as the equation for z when A0=0, but now with K1 replaced with K*. This also implies that the analytical solution (16) is still valid but only with the substitution K1=K*.
In order to characterize the effect of an antibody on the binding of a toxin to its receptor, we introduce the nondimensional parameter Ψ, the relative reduction in CR due to the introduction of an antibody as follows:
(19)Ψ≡z(A0>0)z(A0=0).
The analytical results presented previously enable us easily to derive a simple formula for the antibody efficiency parameter Ψ. By using (10), (16), (18), and (19), we can readily deduce the following:
(20)Ψ=11+ϵλ,ϵ=K1K2,λ=A0C0.
This expression is the main result of the current paper and will be validated with numerical simulations.
To conclude this section let us briefly discuss some additional constraints for the parameters of our model in order for the expression (20) to be valid. As mentioned above the condition of low toxin concentration is always assumed in our study. Another simple condition can be derived from the constraint CR+CA≤T0 and by using (16):
(21)R0C0(1+ϵA0C0)≈ϵR0A0C02≤1,
since R0/C0≤1. This condition could always be checked retrospectively and always hold in our numerical simulations.
3.2. Toxin Internalization
For toxins that are internalized, the effect of antibodies that prevent receptor binding is to reduce the effective rate of internalization. To examine and evaluate this effect, we need to analyze the full systems (1), (3), and (4).
In order to characterize the effect of antibody concentration on the rate of toxin internalization, we introduce a new parameter as follows:
(22)G=Ti(A0>0)Ti(A0=0),
which is a function of time (i.e., G≡G(t)).
Our aim is to deduce function G based on the kinetic models (1), (3), and (4). It is evident that G≤1 for t>0 and G→1 as t→∞ (since in that case all toxin will be internalized).
For the toxins of interest, while the receptor binding is rapid (time sale ~1/(k1C0)) [24, 26], the subsequent internalization is much slower (time scale ~1/k3≫1/(k1C0)). This coupling of slow and fast processes in our system allows us to develop a simplified model of toxin internalization using the well-known framework of Quasi-Steady-State Approximation (QSSA); see [24–28, 37] and refs therein.
When applied to our system, QSSA elucidates the toxin internalization as a two-stage process. After the initial rapid binding of the toxin to the receptor we can simply set dCR/dt=0 in (1). The further slow evolution of T(t) (namely, quasi-steady state) is completely determined by the conservation laws (8) and (4) and spans a time scale of the order of the internalization time (~1/k3). In addition, for solving (4) at the initial stage of internalization, we can assume that Ti≪T0 and write
(23)Ti(t)=k3z0t,t≪1k3,
where z0 is given by expressions (15) and (16). The evolution of Ti(t) for the late stage of internalization can be readily derived from (4) and (6)–(8) by assuming [T0-Ti(t)]≪T0:
(24)Ti(t)=T0[1-exp(-k3t)],t≥1k3,
so Ti(t) exponentially approaches its saturation limit. A simulation of this process is shown in Figure 5, and the slow linear increase of Ti at the initial stage is clearly visible.
Now, consider the case where A0>0. According to (23) the main effect of the introduction of an antibody is to reduce the value of z0, as described in the previous section. Then, based on (22), (23), and (19) we can conclude that, during the quasi-equilibrium stage, the following approximation holds:
(25)G=Ti(A0>0)Ti(A0=0)≈Ψ,
where Ψ is given by expression (20).
The overall effect of introducing an antibody can be best described in terms of the internalization half-time, τi. Without antibody the latter can be estimated from (24) and condition Ti(τi)=T0/2. Thus from (23) we yield
(26)τi≈T02k3z0=C02k3R0.
For the internalization time with the presence of antibody we can apply reduced value of z0 and write the following simple formula:
(27)τiτi0≈1Ψ,
where τi0 is the internalization time in the absence of antibody (A0=0).
Equations (26) and (27) have a clear interpretation. As described in the previous section, the introduction of an antibody results in a decrease, at t≪τi, in the equilibrium value of CR (i.e., in z0). This can be related, in accordance with (23) and (26), to a corresponding decrease in the concentration of internalized toxin Ti and a consequent increase in the toxin internalization time (since it takes longer to achieve a give level of Ti). Since changes in z0 can be described comprehensively by means of the parameter Ψ, it still remains the only parameter needed to characterize the influence of an antibody on the concentration of internalized toxin (25), (27).
It is evident that the two main effects described above (reduction of the concentration of internalized toxin at a given time and increase in the time required for the internalized toxin to reach a given concentration) are not independent of each other. The linear relationships (25), (27) allow us to establish a general identity that relates these two effects for any time t.
Let us assume that for A0=0, τ0 is the time taken for the internalized toxin to reach a concentration Ti0 (i.e., τ0=Ti0/(k3z0); see (23)). The effect of introducing an antibody is to reduce the internalized toxin concentration to a value Ti≤Ti0. Then from (25), (27) we can derive the following identity:
(28)Tiτi=Ti0τi0,
where τi is the time required for the internalized toxin to reach Ti0 when A0>0. The identity (28) has no explicit dependency on antibody kinetic parameters or concentration and provides an easy way to calculate any of the parameters (Ti,Ti0,τi,τi0) if the other three are known.
4. Numerical Results and Discussion
We have derived an analytical expression for the parameter Ψ, the relative ability of an antibody to reduce the binding of a toxin to its receptor (20). Our derivation is based on the following assumptions:
toxin concentration is much lower than the receptor concentration,
for toxins that are internalized, the internalization rate is much slower than establishment of the receptor-toxin binding equilibrium.
Applying these assumptions, we found that parameter Ψ is independent of the toxin concentration (see (20)); that is, it is determined by the ratio of antibody to receptor concentration and not by the ratio of antibody to toxin concentration as commonly used. For the low toxin/receptor ratios likely to occur in biological situations, the condition (21) can be met by large range of antibody kinetic parameters. From this point of view (20) should be valid for most practical applications.
The implications of our analytical findings are illustrated by simulation of the complete kinetic models ((1),(3), (4), and (6)–(8)) using the kinetic constants for ricin and the anti-ricin antibody 2B11 (Table 1). Figure 2 is a simulation of the effect of the presence of an antibody on the binding of the toxin to its receptor (formation of CR). The antibody concentration is expressed as the dimensionless parameter λ=A0/C0. In this case, since R0 and T0≪K1, the parameter C0=R0+K1 is dominated by K1 (1.08·10-7).
Simulated effect of antibody concentration on formation of toxin-receptor complexes CR. Parameter λ=A0/C0, C0=R0+K1. The binding curves were created using the simulation package COPASI and the kinetic constants in Table 1. R0=5 nM, T0=10 pM, C0=1.15·10-7.
Figure 3 shows the effect of increasing antibody concentration on Ψ. There is a good agreement between the values of Ψ determined from (20) and from (19) using the equilibrium values of CR determined from simulation of the complete kinetic model (Figure 3). For instance, the results predict that, for this toxin and antibody combination, the additional protection provided by increasing the antibody concentration diminishes rapidly when λ exceeds 0.1.
Effect of antibody concentration on protection factor. Parameter Ψ (19) was determined from (20) (solid lines) and by using simulated values of CR from Figure 2 at 2500 sec (△), ϵ=25.9.
Figure 4 shows the relationship (20) between Ψ, antibody concentration and the toxin/antibody and the ratio of toxin/receptor dissociation constants (ϵ). This plot is valid for all combinations of toxin, receptor, and antibody consistent with the assumptions used to derive (20), principally T0≪R0. The antibody kinetic parameters and concentration required to provide a specified degree of protection may be determined from this plot. For example, any combination of ϵ and λ falling below the dashed line will reduce either CR or Ti by 80%.
Protection factor Ψ (19) as a function of parameter ϵ=K1/K2 and λ=A0/C0 ((20)): λ=0.01 (⚪); 0.025 (△); 0.05 (□); 0.1 (▽); 0.25 (◊). The range of values for λ and ϵ below dashed line corresponds to 80% protection.
Different time scales for formation of receptor-toxin complex CR (□) and associated toxin internalization Ti (solid lines). Results of COPASI simulation with kinetic constants from Table 1. λ=A0/C0, R0=5 nM, T0=10 pM, C0=1.15·10-7, ϵ=25.9.
This, in turn, enables important judgements to be made about antibody selection. For example, if an antibody concentration of 0.25C0 (λ=0.25) is achievable, then an antibody with an ϵ value of 50 will provide good protection (93% reduction in CR or Ti). If an antibody concentration of only 0.05C0 (λ=0.05) is achievable, then an ϵ value of 250 is required to achieve the same level of protection. The structure of (20) is such that a given increase in protection (Ψ or Γ) may be achieved by either an x-fold increase in ϵ or an x-fold increase in λ.
The effect of antibody on toxin internalization is simulated in Figure 5. Rapid equilibration of receptor and toxin is followed by slow accumulation of toxin within the cell. Equation (25) predicts that Ψ is the only parameter needed to characterize the influence of an antibody on toxin internalization. Figure 6 compares Γ calculated using (25), (20) with Γ determined using values of Ti and Ti0 at t=104 sec from this simulated data and shows good agreement between the two values under the condition T0≪R0, although the value of Γ is slightly greater than Ψ. The plot predicts the degree of protection provided by a given concentration of antibody and enables assessment of the value of increasing antibody concentration beyond a certain value. For example, to enhance the reduction of Ti from 90% to 95% requires doubling of A0.
Comparison of parameters Ψ and Γ. Γ (△) was determined using values of Ti and Ti0 at t=104 sec from toxin internalization time courses simulated using COPASI and the kinetic constants in Table 1. Parameter Ψ (solid line) was determined from (20). R0=5 nM, T0=10 pM, C0=1.15·10-7, ϵ=25.9.
The expression for Ψ, (20), assumes a quasi-equilibrium state in the system. In practice, this state may take significant time to be achieved. Figure 7 shows a simulation of the time taken by the ricin/receptor/2B11 system to reach the quasi-equilibrium state for λ=0.05. The value of Γ determined from the toxin internalization profiles (Figure 7) parallels this process; that is, experimental validation of Γ must allow sufficient time to elapse for the quasi-equilibrium state to be established.
Establishment of the quasi-equilibrium state in the presence of antibody. CR formation (△) was simulated using COPASI and the kinetic constants in Table 1. Γ (□) was determined using (25) and values Ti and Ti0 at t=104 sec using simulated toxin internalization time courses. R0=5 nM, T0=10 pM, C0=1.15·10-7, λ=0.05.
The relationship between the internalization time τi and Ψ described in (27) is shown in Figure 8. Ψ was determined from simulated toxin internalization time courses (Figure 5) as the time to internalize 5·10-14 M ricin. The slope of the fitted line is 1.07, close to the predicted value of 1.0.
Relationship between toxin internalization time τi and protection factor Ψ (19). Solid line is formula (27) and (□) is simulation with COPASI. τi was determined as the time to internalize 5·10-14 M of ricin. All other parameters are the same as in Figure 7.
In summary, the protection provided by an antibody against toxins that act either at the cell surface or after binding to the cell surface followed by internalization may be predicted from a simple kinetic model. Protection parameter Ψ is a simple function of antibody, receptor, and toxin concentrations and the kinetic parameters governing the binding of the toxin to the receptor and antibody:
(29)Ψ=11+(K1/K2)(A0/C0).
The calculated value of Ψ matches closely the degree of protection determined from numerical simulation of the binding and internalization reactions and provides a convenient method for predicting the optimum antibody parameters (concentration and dissociation constants) needed to provide effective treatment or prophylaxis for toxins.
Acknowledgments
The authors acknowledge helpful discussions with Dr. Chris Woodruff and Dr. Ralph Leslie.
ReislerR. B.SmithL. A.The need for continued development of ricin countermeasures20122012414973710.1155/2012/149737PrigentJ.PanigaiL.LamouretteP.SauvaireD.DevilliersK.PlaisanceM.VollandH.CréminonC.SimonS.Neutralising antibodies against ricin toxin2011652-s2.0-7995747632110.1371/journal.pone.0020166e20166OralH. B.OzakinC.AkdisC. A.Back to the future: antibody-based strategies for the treatment of infectious diseases200221225239KaufmannA. F.MeltzerM. I.SchmidG. P.The economic impact of a bioterrorist attack: are prevention and postattack intervention programs justifiable?19973283942-s2.0-0031112046RaineyG. J. A.YoungJ. A. T.Antitoxins: novel strategies to target agents of bioterrorism2004297217262-s2.0-444431779610.1038/nrmicro977NehalM.MichelleC.JuanL.StevenJ.StevenP.GiovanniD. A.LindaN.MichaelB.GeorgeS.LeslieC.A high-affinity monoclonal antibody to anthrax protective antigen passively protects rabbits before and after aerosolized200573795802RazaiA.Garcia-RodriguezC.LouJ.GerenI. N.ForsythC. M.RoblesY.TsaiR.SmithT. J.SmithL. A.SiegelR. W.FeldhausM.MarksJ. D.Molecular evolution of antibody affinity for sensitive detection of botulinum neurotoxin type A200535111581692-s2.0-2214447373110.1016/j.jmb.2005.06.003McGuinnessC. R.MantisN. J.Characterization of a novel high-affinity monoclonal immunoglobulin G antibody against the ricin B subunit2006746346334702-s2.0-3364690506010.1128/IAI.00324-06StaatsH. F.AlamS. M.ScearceR. M.KirwanS. M.ZhangJ. X.GwinnW. M.HaynesB. F.In vitro and in vivo characterization of anthrax anti-protective antigen and anti-lethal factor monoclonal antibodies after passive transfer in a mouse lethal toxin challenge model to define correlates of immunity20077511544354522-s2.0-3564894219410.1128/IAI.00529-07MaddaloniM.CookeC.WilkinsonR.StoutA. V.EngL.PincusS. H.Immunological characteristics associated with the protective efficacy of antibodies to ricin200417210622162282-s2.0-2442626550LoboE. D.HansenR. J.BalthasarJ. P.Antibody pharmacokinetics and pharmacodynamics20049311264526682-s2.0-864423491010.1002/jps.20178AlbrechtM. T.LiH.WilliamsonE. D.LeButtC. S.Flick-SmithH. C.QuinnC. P.WestraH.GallowayD.MateczunA.GoldmanS.GroenH.BaillieL. W. J.Human monoclonal antibodies against anthrax lethal factor and protective antigen act independently to protect against Bacillus anthracis infection and enhance endogenous immunity to anthrax20077511542554332-s2.0-3564901198410.1128/IAI.00261-07SmithT. J.LouJ.GerenI. N.ForsythC. M.TsaiR.LaPorteS. L.TeppW. H.BradshawM.JohnsonE. A.SmithL. A.MarksJ. D.Sequence variation within botulinum neurotoxin serotypes impacts antibody binding and neutralization2005739545054572-s2.0-2394447223510.1128/IAI.73.9.5450-5457.2005ReuvenyS.WhiteM. D.AdarY. Y.KafriY.AltboumZ.GozesY.KobilerD.ShaffermanA.VelanB.Search for correlates of protective immunity conferred by anthrax vaccine2001695288828932-s2.0-003505529510.1128/IAI.69.5.2888-2893.2001LittleS. F.IvinsB. E.FellowsP. F.PittM. L. M.NorrisS. L. W.AndrewsG. P.Defining a serological correlate of protection in rabbits for a recombinant anthrax vaccine2004223-44224302-s2.0-034509657510.1016/j.vaccine.2003.07.004PittM. L. M.LittleS. F.IvinsB. E.FellowsP.BarthJ.HewetsonJ.GibbsP.DertzbaughM.FriedlanderA. M.In vitro correlate of immunity in a rabbit model of inhalational anthrax20011932476847732-s2.0-003586057110.1016/S0264-410X(01)00234-1PeachmanK. K.RaoM.AlvingC. R.BurgeR.LepplaS. H.RaoV. B.MatyasG. R.Correlation between lethal toxin-neutralizing antibody titers and protection from intranasal challenge with Bacillus anthracis ames strain spores in mice after transcutaneous immunization with recombinant anthrax protective antigen20067417947972-s2.0-2964443835610.1128/IAI.74.1.794-797.2006MohamedN.ClagettM.LiJ.JonesS.PincusS.D'AliaG.NardoneL.BabinM.SpitalnyG.CaseyL.A high-affinity monoclonal antibody to anthrax protective antigen passively protects rabbits before and after aerosolized Bacillus anthracis spore challenge20057327958022-s2.0-1994443297210.1128/IAI.73.2.795-802.2005PapaliaG. A.BaerM.LuehrsenK.NordinH.FlynnP.MyszkaD. G.High-resolution characterization of antibody fragment/antigen interactions using Biacore T100200635911121192-s2.0-3375060607710.1016/j.ab.2006.08.032AloufJ. E.PopoffM. R.20063rdAmsterdam, The NetherlandsElsevierSandvigK.van DeursB.Membrane traffic exploited by protein toxins2002181242-s2.0-003644113510.1146/annurev.cellbio.18.011502.142107SandvigK.van DeursB.Transport of protein toxins into cells: pathways used by ricin, cholera toxin and Shiga toxin2002529149532-s2.0-003701014110.1016/S0014-5793(02)03182-4SkvortsovA.GrayP.Modeling and simulation of receptor-toxin-antibody interactionProceedings of the 18th World IMACS Congress and International Congress on Modelling and Simulation (MODSIM '09)July 2009Cairns, Australia1851912-s2.0-80052978606TruskeyG. A.YuanF.20092ndNew York, NY, USAPrentice HallLauffenburgerD. A.19932ndOxford, UKOxford Univercity PressGoldsteinB.DemboM.Approximating the effects of diffusion on reversible reactions at the cell surface: ligand-receptor kinetics1995684122212302-s2.0-0028960497CoppeyM.BerezhkovskiiA. M.SealfonS. C.ShvartsmanS. Y.Time and length scales of autocrine signals in three dimensions2007936191719222-s2.0-3454873965910.1529/biophysj.107.109736GhoshS.GopalakrishnanM.Forsten-WilliamsK.Self-consistent theory of reversible ligand binding to a spherical cell2007443443522-s2.0-38049103229MaxfieldF. R.McGrawT. E.Endocytic recycling2004521211322-s2.0-074230588010.1038/nrm1315SheffD. R.DaroE. A.HullM.MellmanI.The receptors recycling pathway contains two distinct populations of early endosomes with different sorting functions199914511231392-s2.0-003352593010.1083/jcb.145.1.123LamazeC.JohannesL.AloufJ. E.PopoffN. R.Intracellular trafficking of bacterial and plant protein toxins20063rdAmsterdam, The NetherlandsElsevier135153WileyH. S.CunninghamD. D.A steady state model for analyzing the cellular binding, internalization and degradation of polypeptide ligands19812524334402-s2.0-0019404817ShankaranH.ResatH.WileyH. S.Cell surface receptors for signal transduction and ligand transport: a design principles study200736, article e1012-s2.0-3434735557610.1371/journal.pcbi.0030101HoopsS.GaugesR.LeeC.PahleJ.SimusN.SinghalM.XuL.MendesP.KummerU.COPASI—a complex pathway simulator20062224306730742-s2.0-3384536851310.1093/bioinformatics/btl485SandvigK.OlsnesS.PihlA.Kinetics of binding of the toxic lectins abrin and ricin to surface receptors of human cells197625113397739842-s2.0-0017074398PolyaninA. D.ManzhirovA. V.2007New York, NY, USACRC PressTzafririA. R.EdelmanE. R.On the validity of the quasi-steady state approximation of bimolecular reactions in solution200523333433502-s2.0-1184427381810.1016/j.jtbi.2004.10.013