article: A few broadly neutralizing antibodies, isolated from HIV-1 infected individuals, recognize epitopes in the membrane proximal external region (MPER) of gp41 that are transiently exposed during viral entry. The best characterized, 4E10 and 2F5, are polyreactive, binding to the viral membrane and their epitopes in the MPER. We present a model to calculate, for any antibody concentration, the probability that during the pre-hairpin intermediate, the transient period when the epitopes are first exposed, a bound antibody will disable a trivalent gp41 before fusion is complete. When 4E10 or 2F5 bind to the MPER, a conformational change is induced that results in a stably bound complex. The model predicts that for these antibodies to be effective at neutralization, the time to disable an epitope must be shorter than the time the antibody remains bound in this conformation, about five minutes or less for 4E10 and 2F5. We investigate the role of avidity in neutralization and show that 2F5 IgG, but not 4E10, is much more effective at neutralization than its Fab fragment. We attribute this to 2F5 interacting more stably than 4E10 with the viral membrane. We use the model to elucidate the parameters that determine the ability of these antibodies to disable epitopes and propose an extension of the model to analyze neutralization data. The extended model predicts the dependencies of for neutralization on the rate constants that characterize antibody binding, the rate of fusion of gp41, and the number of gp41 bridging the virus and target cell at the start of the pre-hairpin intermediate. Analysis of neutralization experiments indicate that only a small number of gp41 bridges must be disabled to prevent fusion. However, the model cannot determine the exact number from neutralization experiments alone.
  A small but growing number of HIV-1 broadly neutralizing monoclonal antibodies (mAbs) have been isolated from HIV-1 infected individuals (reviewed in [1]). A subset of these (4E10,2F5, m66. 6 and Z13e1) bind to epitopes at the base of gp41, along a segment adjacent to the viral membrane known as the membrane proximal external region (MPER) that links the gp41 transmembrane domain to its ectodomain 2,3. Two, 4E10 and 2F5, have been extensively characterized. Each has been shown to be polyreactive [4]–[6], with an elongated heavy chain complementarity determining region 3 (CDRH3) that has a hydrophobic surface that does not contact the peptide antigen [7], [8], but that contributes to the ability of 4E10 and 2F5 to bind to the virion membrane and to broadly neutralise [9]–[13]. Recently, self-antigens have been identified that these antibodies recognize [14]. Of the two mAbs, 2F5 interacts more strongly with membrane as demonstrated by its ability to bind to virus-like particles bearing no envelop glycoprotein (Env spikes) [6]. Detailed structural studies and surface plasmon resonance measurements (SPR) have shown that the epitopes for 4E10 and 2F5 are within the MPER, lying along the viral surface with parts immersed in the viral membrane [15]–[17]. The binding of these antibodies to their epitopes on the viral membrane appears to be step-wise [4], [17], [18] with the membrane associated antibody and the MPER undergoing conformaional changes [15], [16], [18] that result in a more stably bound complex [4], [8], [17], [18]. For both mAbs and their Fabs, the rate constants that describe the kinetics of the two step process have been determined for binding to liposome-gp41 peptide conjugates [4]. The epitopes recognized by 4E10 and 2F5 on gp41 are displayed on the viral membrane after the binding of HIV-1 to the target cell has occurred and the processes that result in membrane fusion have begun [2]. This stage, between pre- and postfusion, called the pre-hairpin intermediate [19], has a lifetime that has been estimated to range from a few seconds to many minutes [20]–[23]. The tenet that 4E10 and 2F5 only recognize epitopes on the MPER during the pre-hairpin intermediate appears to strictly hold only for HIV-1 isolates that are highly resistant to neutralization [24]. For HIV-1 isolates that are less resistant to neutralization, the MPER epitopes may be accessible to antibody binding before the pre-hairpin intermediate stage has been initiated. Alternatively, enhanced susceptibility to MPER neutralizing antibodies could be a result of prolonged exposure of the intermediate state during viral entry [25]. Although numerous studies have determined rate constants for the binding of 4E10 and 2F5 to peptide epitopes in solution and peptides conjugated to liposomes [2], [4], [10], [12], [13], [15], [17], [18], and estimates of the lifetime of the exposure of the MPER on HIV-1 have been obtained [20]–[22], [26], a quantitative description of how these parameters determine the mAb concentrations that are effective in inhibiting fusion has been lacking. Here we present a model that allows us to calculate, as a function of the serum antibody concentration, the probability of a trivalent MPER being disabled by polyreactive mAbs. A question we also consider is how the probability of disabling a single epitope is related to the probability of blocking membrane fusion of HIV-1 and target cell when, at the start of the pre-hairpin intermediate stage, epitopes bridge the two membranes. We present a model that predicts the form of a neutralization curve, i. e. , the percent neutralization versus the antibody or Fab concentration, and how the for neutralization depends on, as well as on the rate constants that characterize antibody binding and the rate of fusion of gp41. We use the model to fit neutralization experiments where neutralization is achieved with either 4E10 or 2F5 IgG or Fab. Our analysis of these experiments indicates that at the start of the pre-hairpin intermediate only a few gp41 trimers are involved in bridging the virion and target cell. As the bonds break sequentially those that remain come under additional strain and are disabled more easily.
  We first consider a polyreactive Fab at low concentration that can bind to both a site in the MPER and to lipids in the membrane. The binding reactions that the Fab can undergo in the model and their rate constants are illustrated in Figure 1. The parameters of the model are summarized in Table 1. By the low concentration limit we mean that at most one Fab can be bound to a spike, even though the spike contains three gp41 binding sites for the Fab. After we explore this limit we treat the general case where up to three Fab can be bound to one spike. We assume it takes only one Fab bound for a sufficient time to an MPER to disable a spike, i. e. we take the stoichiometry of trimer neutralization to be one [32]. We use the term spike and epitope interchangeably in what follows, so in the model an epitope is composed of three gp41 with each containing a single binding site for the Fab. A free epitope has all three sites free. In the model free epitope with surface concentration is lost through fusion with rate. Free epitope is also lost when a Fab, with solution concentration, binds to the epitope. Free epitope is gained when a Fab bound to an epitope with surface concentration dissociates from the epitope. A bound Fab can induce a conformational change leading to the formation of long-lived complexes with surface concentration [4], [15], [17]. A Fab bound in this way for a sufficient time can render the epitope nonfunctional [30]. This conversion to a nonfunctional state takes a mean time equal to. Since the mean lifetime of an epitope bound to a 4E10 or 2F5 Fab or IgG before it undergoes a conformational change is short 5–20 s ([4], Table 2), we expect there to be negligible disabling of epitopes before the conformational change occurs. We therefore take disabling of the epitope to occur only after a stable complex has formed. In the low Fab concentration limit the following set of ordinary differential equations (ODEs) describes the kinetics of the epitope in the presence of Fab. The initial conditions for these equations are that at t = 0, when the MPER first becomes accessible, and. (4) (5) (6) We only need the solution of these equations at after all epitopes have either fused or been disabled (see Eq. (3) ). The solution of Eqs (4) – (6) for is given in SI2 of Text S1, the Supporting information. Substitution of the expression we obtain for into Eq. (3) gives us the desired expression for (7) The parameter is an effective three-dimensional forward rate constant given in terms of the rate constants for the reversible binding of the Fab to the epitope, and, and the rate constants for the reversible transition induced by the Fab binding to the epitope, and. In addition, depends on the mean time that an Fab is required to remain bound to an epitope to disable it, . (8) Note that when, , and therefore. This is because in the model the Fab or mAb can only disable the epitope after it makes a conformational change to a stably bound configuration, and the rate at which this conformational change occurs is. The faster the epitope is disabled the higher will be the probability of preventing fusion. From Eq. (8), a sufficient requirement to achieve maximal disabling of the epitope by the antibody is that or equivalently that (9) For all the mAbs and Fabs in Table 1, , so the condition for maximal disabling, Eq. (9), reduces to requiring that. For a Fab to be effective in neutralizing, as Fabs of 4E10 and 2F5 are [8], [9], the model requires that, i. e. , the time it takes to disable an epitope must be much shorter than the time a Fab remains bound to the epitope. From Table 2 we see this requires that, i. e. , disabling must occur in five minutes or less. In Figure 2 we plot as a function of, for two different Fab concentrations and a series of different values of. Figure 2 demonstrates that when, becomes independent of. When, Eq. (8) becomes (10) Eq. (10) is the limit of Eq. (8) as, i. e. , the instant disabling limit, and therefore, as long as, there is no gain achieved in the Fab' s ability to disable an epitope by decreasing the time required to disrupt fusion () or by increasing the time the Fab remains bound to the MPER in the stable state (). In this limit, the rate-limiting step in the disabling of an epitope is the formation of the stably bound antibody-epitope complex, i. e. , the bound complex after the epitope has undergone a conformational change. Eq. (10) shows that equals the product of the rate of binding to the epitope, , and the probability that the bound complex will undergo a conformational change rather than breakup, . As can be seen from Eq. (7), is determined by the competition between the rate of formation of the stably bound complex, , and the rate of fusion, . Increasing the affinity of the Fab for the MPER will increase. Lowering will slow the rate of breakup of the stably bound complex, while increasing will enhance its formation. In SI3 of Text S1 we write down and solve the equations for the model when we allow up to three Fab to bind to a single epitope. We show that the results for in the low concentration limit, Eqs. (7) and (10), hold when the Fab concentration satisfies the inequality: (11) For the rate constants in Table 2 we have that this inequality holds when M for 4E10 Fab and M for 2F5 Fab. Since 50% neutralization () is typically in the range of 1–50 nM [8]–[10], we expect the simple expression we previously derived for, Eqs. (7) and (10), to give an excellent description up to these concentrations. Figure 3, for 4E10 Fab, shows a comparison between the predictions for the low concentration limit and the full model. Agreement is excellent for Fab concentrations less than M. In the model we assume that at most one IgG can bind an epitope even though it is composed of three gp41 molecules [28]. In addition, we assume that an IgG can' t cross-link two different epitopes. On average the distance between epitopes on HIV-1 is too far for an antibody to span [30]. It has been argued that this allows HIV-1 to evade the effects of antibody avidity and prevent antibodies from forming long-lived attachments by cross-linking different epitopes on the viral surface [31]. However, recent studies show that Env spikes are clustered on mature HIV-1 particles [33]. In the contact region between simian immunodeficiency virus (SIV) and a CD4+ T cell in the prefusion stage, about seven Env spikes are clustered in the contact region with a center to center distance between 140–170 Å [34]. This was closer than the average spacing of Env spikes in unbound SIV (Å) and comparable to the span between two Fab sites on an IgG (Å [35], [36]). The binding reactions that the IgG can undergo in the model, and the rate constants that determine the kinetics of the model, are shown in Figure 4. In Figure 4A the IgG cannot bind membrane alone. In Figure 4B the IgG first binds to a lipid site in the membrane with equilibrium constant and then, through a surface reaction, binds to the MPER. The IgG can have only one of its Fab bound to lipid or it can attach to lipid with both its Fab binding sites (in Figure 4B). It forms the second bound site with an equilibrium cross-linking constant. In addition, if it is bound to the MPER with one site free, that site can also bind to the lipid with equilibrium cross-linking constant (and in Figure 4B). The equations that describe the model in Figure 4B are given in SI4 of Text S1. There we show that these equations can be reduced to the same set of equations that describe the model for Fab, Eqs. (7) and (10), with the one change, that for IgG the forward rate for binding (12) where is the forward rate constant for an Fab to bind to an MPER and is the concentration of lipid sites in the membrane. For models where the Fab cannot bind to lipid in the absence of the MPER, , and therefore. The factor of two arises because for the same molar concentration an IgG has twice as many binding sites as a Fab. In the model, because cross-linking between epitopes is prohibited, the only avidity effects come from the formation of doubly bound antibodies with either both sites bound to lipid or one site bound to a MPER and one to lipid. The avidity effect is characterized by the term. For the values in Table 1, for 4E10, indicating that there are essentially no 4E10 IgG that have both sites bound to lipid on the liposome-peptide conjugates used in the experiments of Table 2 [4], nor are there any IgG with one site bound to a MPER and the other to lipid. For 2F5, or about 40% of the 2F5 IgG that are bound to lipid and not MPER peptide are doubly bound. (This follows because in Figure 4.) In addition we would expect in the absence of fusion that there would be 1. 3 times more 2F5 IgG bridging an MPER and a lipid site than bound to MPER alone. We note that the value of may be an over estimate since these experiments were done on peptide-liposome conjugates and some cross-linking of the MPER peptides may have occurred. The picture that 2F5 interacts more stably than 4E10 with membrane lipid is in agreement with previous experiments. SPR studies of 4E10 and 2F5 binding to liposomes without peptide (Figure 3A and 3B in [4]) indicate that bound 2F5 dissociates more slowly from the liposomes than bound 4E10. Also 2F5 was found to bind to virus-like particles bearing no envelop glycoprotein while such binding could not be detected with 4E10 [6]. This is consistent with 4E10 being described by the wedging model and 2F5 by the encounter model [13]. From the values we obtained for, we predict that in a neutralization experiment 4E10 IgG and Fab should be similarly effective when the concentration per mole of Fab sites is the same, while we predict the opposite for 2F5. Shown in Figure 5 are neutralization experiments for 4E10 and 2F5 Fab and IgG, plotted as a function of the molar concentration per Fab site. Figure 5A and 5B show that unlike 2F5,4E10 IgG is only marginally more effective than its Fab fragment. The inserts in A and B of Figure 5 are the predicted values of for the parameters given in Table 1. These parameters were determined from SPR binding studies to peptide-liposome conjugates, while the neutralization studies used pseudotyped virus, and therefore neither the membrane composition nor the manner in which the epitope was displayed on the membrane was the same. None-the-less, the qualitative prediction is born out. For the same concentration of Fab sites, in a neutralizing experiment 2F5 IgG is much more effective than its Fab compared with 4E10 IgG and its Fab as has previously been observed [36]. From Table 1 we see the value of for 2F5 is greater than for 4E10. Therefore, the model predicts that 2F5 will be more effective at neutralization than 4E10. Neutralization experiments that use either 2F5 or 4E10 mAb, shown in Figure 5C, agree with this prediction. For mAbs and Fab that bind epitopes within the MPER we have derived an expression for the probability, , of disabling an epitope on a virion that is bound to a target cell, Eqs (7) and (10), as a function of the Fab or IgG concentration. In deriving we assumed that there was a uniform distribution of gp41 trimers bridging the virus and target cell and that was sufficiently large that was independent of. An obvious question is, how is this probability related to the fraction of target cells prevented from being infected in a neutralization assay? For a virion that initially has epitopes bound to a target cell, and to prevent fusion all must be disabled, we propose that the probability of preventing fusion is (13) where is the probability a bond in the presence of other bonds is disabled before it fuses. Writing as a product of different probabilities assumes that the individual probabilities are independent, i. e. , each individual probability depends only on the number of bonds present that bridge the virus and target cell and not on what has happened previously. Writing in this way also recognizes the possibility that the number of gp41 trimers bridging the region between virion and target cell influences the probability of a bond rupturing before fusion. In general cellular membranes repel each other as a result of nonspecific forces that fall off with the distance between the membranes [37]. For the membranes to attach, specific bonding is required to overcome the repulsive forces. The bonds that hold the two membranes together are under stress with the bonds at the edge of the contact region between the two membranes experiencing larger stresses than those in the interior of the contact region [38]. For this reason we expect that when there is one or only a few bonds present there will be increased strain on a bond and as a result it will be more easily disabled by a mAb then when a bond is surrounded by many other bonds as illustrated in Figure 6. Models have been presented to determine the stoichiometry of entry, defined as the number of attached gp41 trimers required to infect the target cell [32], [39], [40]. Although this parameter is not defined in our model, our model encompasses the stoichiometry of entry. For example, if in Eq. (13) we have that and, then. If the virion and target cell are attached by three trimers, if one of those trimers is disable, the two remaining trimers will spontaneously disable. The stoichiometry of entry has yet to be rigorously determined, but it is most probably greater than one [32]. To analyze neutralization experiments we treat the population of virion-target cell complexes as if they were homogeneous, each held together with the same number of bonds at the start of the experiment. We also assume that each target cell has at most one HIV-1 bound to it. Since is the probability of a target cell not becoming infected when bound to a virion in the presence of an antibody concentration, it also equals the fraction of target cells in a neutralization assay that are prevented from being infected at antibody concentration. To proceed we need to obtain an expression for as function of and of. We do this by assuming a form for that has the following limiting behavior: when the number of bonds is large, is independent of and equal to; When the number of bonds is small, is elevated (more easily disabled) compared to. (14) where. When, and is independent of the number of bonds present, while when, and the th bond breaks spontaneously. To help us interpret our results in the next section when we analyze neutralization experiments, it will be useful to have a concrete expression for. A reasonable choice because of its flexibility, few parameters, and its range between zero and one, is a Hill function which we write as follows: (15) In Figure 6, and are plotted for, and. For any set of parameters, when, and. When, approaches as increases. It can be seen from Figure 6 and Table 3 that for these parameters and are both greater than 0. 9. Therefore, in this example, the last two trimers between virion and target cell have a high probability of spontaneously disabling indicating the stoichiometry of entry. From Eqs. (13), (14) and (15) we can write in the following form: (16) where (17) and. It is useful to rewrite in terms of the, the antibody concentration at which 50% of the virion-target cell complexes do not go on to fuse. In terms of the, Eq. (16) becomes (18) where (19) and (20) The parameter was chosen so that when. It is a decreasing function of that is equal to one when. Eq. (18) predicts the following dependencies of on the parameters of the model: The faster fusion occurs (larger) the higher; the more bonds required to be disabled (larger and therefore smaller) the higher; and the more rapidly the antibody binds to the epitope and induces a stable complex (larger) the lower. Note that the is inversely proportional to. Since we have an expression for how depends on the binding parameters, Eq. (10), using the values obtained from SPR data (Table 2), we can predict how the will shift as a result of modifying the rate constants. For example, increasing the mAb affinity by decreases by a factor of ten is predicted to lower the for mAbs 4E10 and 2F5 by factors of 4. 86 and 3. 1 respectively. Since is directly proportional to, increasing the mAb affinity by increasing by a factor of ten is predicted to lower the for both mAbs by a factor of ten. Eq. (18) predicts that if a series of neutralizing experiments are performed under the same conditions with different mAbs or Fabs, and the fraction of target cells neutralized is plotted against, all the experiments should lie on the same curve. To test this prediction we use the four neutralization experiments shown in Figure 5. We simultaneously fit the four neutralization experiments with five free parameters, four and, that we take to be the same in all experiments. As can be seen in Figure 7A, when plotted as a function of the four experiments lie on a universal curve. Figure 7A gives a best fit value for, where the error bounds are the 95% confidence interval obtained from bootstrapping. Using other parameter estimation methods to simultaneously fit the four neutralization experiments, as discussed in SI5 of Text S1, give similar values for. When we fit each experiment separately ranged from 0. 45 to 1. 15 with an average. Figures 7B () and 7C () show that the fits get progressively worse with increasing. To interpret the value of that we have obtained, it is instructive to consider the example in Figure 6. As can be seen from Table 3, a in this example would mean that initially there were between three () and four () bonds attaching the virus and target cell with the last bond, , breaking spontaneously (). Since for our experimental system we don' t know the values of and, we can' t calculate the values of and determine the closest to our best fit value of 0. 64. However, we do know that when, so for our value, the values for must be small. This suggests that as in our example, corresponds to the total number of gp41 attachments at the start of the experiment. The bond marks a transition. For the bonds are highly strained and as a result more easily disabled than when.
  The probability of disabling a trivalent gp41 in the pre-hairpin intermediate depends on a number of factors including the concentration of mAb available to bind to gp41 and the number of gp41 bridging the virus and target cell. We expect that for a fixed antibody concentration, as the number of gp41 attachments increases, the probability decreases (it becomes harder to disable a bond) reaching a plateau where the probability (in our model) is independent of the number of bonds present. In this limit, increasing the number of bonds increases the contact area between virus and target cell while the bond density remains uniform [37]. For this limit we presented a kinetic model and used it to calculate, for a given concentration of polyreactive mAbs or Fabs, the probability that during the pre-hairpin intermediate, a trivalent gp41 will be rendered incapable of completing the fusion process when a single Fab or mAb binds to one of the three exposed MPERs. The model predicted that for the antibody to be effective in neutralizing HIV-1, the time it takes to disable an epitope must be shorter than the time the epitope remains in the long-lived conformation induced by antibody binding. For 4E10 and 2F5, whose dissociation rate constants for the breakup of the long-lived complex s−1 (Table 1 and [4], [17]), the time to disable an epitope must be about five minutes or less. Increasing the time the antibodies remain bound in the long-lived conformation (lowering) would not improve their neutralization capabilities. Rather, it is the rate of formation of the long-lived conformation that is the limiting step in determining the ability of these antibodies to disable epitopes. The effective forward rate constant for this step, the parameter in the model, given by Eq. (10), equals the product of the forward rate constant for binding of the antibody in solution to the epitope, , and the probability that the bound complex undergoes a conformational change rather than dissociating, . In the encounter model, which we believe well characterizes the kinetics of binding of 2F5, is the overall rate constant for the binding of the antibody from solution to the membrane, diffusion of the membrane bound antibody to the epitope, and binding of the antibody to the epitope through a surface reaction. Increasing the antibody' s affinity for the membrane or its affinity for the MPER is predicted to result in higher values of, higher values for the probability of disabling an epitope and lower values for the. Increasing the affinity of the antibody by lowering its dissociation rate constant for the MPER, , will increase the probability of forming a long-lived complex, while increasing the affinity by increasing the forward rate constants will increase the rate of formation of the initial antibody-MPER complex. From previous SPR studies of 2F5 and 4E10 Fab and IgG binding to peptide-lipisome conjugates [4] we determined the values of for these four ligands. These values led us to predict that 2F5 IgG would exhibit a much stronger avidity effect than 4E10 IgG. The prediction was tested by comparing neutralization experiments using Fab and IgG and plotting the ligand concentration in Fab sites per mole. If there were no avidity effect, such a plot should yield two identical curves. In the plots in Figure 5A and 5B the prediction was born out. The mAb 2F5 was much more effective than its Fab, while the mAb 4E10 was only marginally better than its Fab. In the model the strong avidity of 2F5 IgG arises as result of its being able to form cross-links between the MPER and the viral membrane. This result is consistent with experiments showing that 2F5 interacts more stably with the viral membrane than does 4E10 [6]. We cannot rule out that an additional contribution to the avidity comes from the mAb cross-linking different MPERs on the same virion, something not in the model. Indeed, in light of recent studies showing that Env spikes are clustered on mature HIV-1 [33], and because 2F5 and 4E10 interact with gp41 after it has attached to the target cell membrane and clustered in the contact region between the two membranes, it is to be expected that some cross-linking of MPERs occurs. Still, 4E10 shows little improvement over its Fab in neutralization, indicating that for 4E10 there is little or no cross-linking of MPERs. In contrast, 2F5 shows a high avidity, but what contribution, if any, comes from 2F5 cross-linking two MPERs rather than MPER and viral membrane, is unknown. Although, the model we presented to calculate delineated the roles of the rate constants in determining the ability of 2F5 and 4E10 to disable an epitope, it couldn' t be directly applied to analyze neutralization experiments. To do this we considered an HIV and target cell bridged by N gp41 at the start of a neutaliation experiment and asked the question, what is the probability of sequentially disabling all N gp41s before one of them fused? To answer this question required a model that quantified the influence of the number of bridging bonds on bond breakage. We proposed an expression for, the probability of disabling a bond when bonds are present, given by Eq. (14). The expression was chosen so that it had the correct limiting behavior, that when was large was independent of and when was small the probability of disabling increased as the bonds came under additional strain. Using this expression for we were able to obtain an expression for the probability of disabling all bonds before fusion occurred, Eq. (18). One prediction of the model, Eq. (20) is how the depends on the model parameters. For a mean time of fusion of 10–30 min () [20]–[22], , the value for 4E10 in Table 2, and, nM. For, the value for 2F5 from Table 1, nM. These predicted values, obtained for values of from SPR experiments with peptide-vesicle conjugates, are consistent with the values determined from neutralization experiments [8]–[11], [13], [24], [28], [36], [41] and lend some support to the model. The model also predicted that for a series of neutralization experiments with different mAbs and Fabs, a plot of the fraction of virus neutralized in each experiment versus, should lie on a single curve. After we fit Eq. (18) to four experiments that used 2F5 and 4E10 mAb and Fab, and replotted the data as a function of, the data did lie on a universal curve (Figure 7A). From the analysis of these data we could not determine the number of gp41 bridging the virus and target cell at the start of the experiment, but the analysis indicated that these bonds are few in number and under considerably more strain then if there were a large number of attachments in the contact region. There are relatively few Env spikes on HIV-1. The average number has been reported to be between seven and 14 spikes per HIV-1 [30], [33], [42], [43]. Zhu et al. found spikes per HIV-1 with a range of four to 35 [30]. Using electron tomography about five to seven spikes have been observed in the contact region between virion and CD+ T cell in the prefusion stage [34]. This was the case when the virion was either SIV or HIV. However, it is not known how many gp41 membrane spanning trimers are established after the loss of gp120 and the entrance into the pre-hairpin intermediate. Attempts have been made to determine the stoichiometry of HIV entry, but as yet reliable estimates of the minimum number of attachments required for fusion to occur have not been established [32], [39], [40], [44]. It therefore remains to be seen whether this picture, where just a few weak attachments are mediating fusion and neutralization can be achieved by disabling these, is a reasonable description of what is occurring in the neutralization experiments we have analyzed.
  Neutralization was determined by monitoring reductions in luciferase reporter-gene expression after a single round of infection by pseudotyped HIV-1 virus in TZM/bl cells, as described previously [45] The Tzm/bl cell neutralization assay was done using 2F5 and 4E10 mAbs and their Fab fragments. Antibodies were preincubated with the pseudovirus for 1 h using a dose titration starting at and a series (seven) of 3-fold dilutions. Best fit values for the parameters where obtained using the Common Los Alamos Software Library fitting routine DNLSI which is based on a finite difference Levenberg-Marquardt algorithm for solving nonlinear least-squares problems. Estimates of the 95% confidence intervals were obtained by using the bootstrap method described by Efron and Tibshirani [46]. For each estimate 500 simulations were performed. Increasing the number of simulations caused no significant change in the error estimate. The program that runs under Linux and Mac OS, and is called FUI, can be downloaded from Professor Rob J. de Boer' s website, http: //www-binf. bio. uu. nl/rdb/software. html. Additional methods, discussed in SI5 of Text S1, were used to check the parameter estimations.
summary: Most people who become infected with HIV generate a strong antibody response to the infecting virus population. Unfortunately, the protection offered by the antibody is short lived as the virus rapidly mutates and renders the antibodies impotent in preventing further infection. There are a few antibodies, however, that have been isolated from infected individuals that can block infection by many different viral strains. Among these are several that target sites on the HIV that are exposed only after the virus has attached to a cell. These antibodies have a brief window of time to prevent fusion of the virus and cell. They are special in that they bind both to the viral membrane and to sequences on the gp41 protein that lie along the viral surface. Here, we present a model that predicts the concentrations at which these antibodies effectively neutralize the virus. The model tells us what properties of antibody binding are key in determining efficient neutralization and what properties have little influence. A prediction of the model is that in a standard neutralization assay there are only a small number of attachments between virus and cell and disabling these is sufficient to prevent infection.
keywords: medicine
  infectious diseases
  hiv
  immunology
  biology
  viral diseases
  immunoglobulins
title: Estimating the Probability of Polyreactive Antibodies 4E10 and 2F5 Disabling a gp41 Trimer after T Cell-HIV Adhesion