Advances in Computational ImmunologyView this Special Issue
Research Article | Open Access
Liusong Yin, Xiaoying Chen, Abhinav Tiwari, Paolo Vicini, Timothy P. Hickling, "The Role of Aggregates of Therapeutic Protein Products in Immunogenicity: An Evaluation by Mathematical Modeling", Journal of Immunology Research, vol. 2015, Article ID 401956, 14 pages, 2015. https://doi.org/10.1155/2015/401956
The Role of Aggregates of Therapeutic Protein Products in Immunogenicity: An Evaluation by Mathematical Modeling
Therapeutic protein products (TPP) have been widely used to treat a variety of human diseases, including cancer, hemophilia, and autoimmune diseases. However, TPP can induce unwanted immune responses that can impact both drug efficacy and patient safety. The presence of aggregates is of particular concern as they have been implicated in inducing both T cell-independent and T cell-dependent immune responses. We used mathematical modeling to evaluate several mechanisms through which aggregates of TPP could contribute to the development of immunogenicity. Modeling interactions between aggregates and B cell receptors demonstrated that aggregates are unlikely to induce T cell-independent immune responses by cross-linking B cell receptors because the amount of signal transducing complex that can form under physiologically relevant conditions is limited. We systematically evaluate the role of aggregates in inducing T cell-dependent immune responses using a recently developed multiscale mechanistic mathematical model. Our analysis indicates that aggregates could contribute to T cell-dependent immune response by inducing high affinity epitopes which may not be present in the nonaggregated TPP and/or by enhancing danger signals to break tolerance. In summary, our computational analysis is suggestive of novel insights into the mechanisms underlying aggregate-induced immunogenicity, which could be used to develop mitigation strategies.
Therapeutic protein products (TPP) from nonhuman, humanized, and human origins include monoclonal antibodies (mAbs), Fc fusion proteins, blood factors, hormones, cytokines, chemokines, and engineered protein scaffolds . They have been widely used to treat a variety of human diseases, including cancer, anemia, hemophilia, rheumatoid arthritis, multiple sclerosis, and inflammatory bowel diseases [1, 2]. Their large success is mainly due to increased target specificity, decreased intrinsic toxicity, and longer half-lives compared with small molecule drugs . These advantages have led to the expansion of TPP in the drug market, with annual revenues of over 100 billion US dollars [1, 2]. However, unwanted immune responses against TPP, such as generation of anti-drug antibodies (ADA), have raised concerns on both drug efficacy and patient safety [4–8]. The effect of ADA on clinical outcomes ranges from no obvious impact to severe loss of efficacy and adverse effects such as infusion reactions . The mechanisms leading to the generation of immunogenicity are yet to be established, but several risk factors have been proposed [9–12], which can be classified as follows: (i) patient-related: genetic background, immunological status, and prior exposure , (ii) treatment-related: route, dose, and frequency of administration [7, 13], and (iii) product-related: drug origins, characteristics such as protein structures and aggregates, and formulations .
Among these risk factors, aggregates of TPP are of particular concern due to their potential role in inducing both T cell-independent and T cell-dependent immune responses [14–17] (Figure 1). It has been previously found that aggregated recombinant human interferon alpha2b generated by thermal stress, low pH, or oxidization stress is more immunogenic in mice compared with nonaggregated product [18–20]. High immunogenicity in mice has also been observed for aggregates of other TPP, such as human mAbs [21–23], human epoetin alfa , human factor VIII [25, 26], human interferon beta , and murine growth hormone . In the clinic, the different ADA incidence rates for several recombinant human interferon beta drugs have been attributed to the differences in aggregation levels . However, the detailed mechanism by which aggregates increase immunogenicity, especially in humans, is yet to be established. For example, it is unknown whether aggregates increase immunogenicity through a T cell-dependent or T cell-independent pathway; and which processes of ADA production could be altered by aggregates is also unknown. In the case of TPP, immunogenicity could be induced through both T cell-dependent and T cell-independent pathways [9, 12]. In the T cell-dependent pathway, antigenic peptides derived from TPP could be presented by major histocompatibility complex class II molecules (MHC II) on antigen-presenting cells (APC) that have been matured by danger signal to stimulate antigen-specific CD4+ T cells. Activated CD4+ T cells would then stimulate antigen-specific B cells that will be responsible for the production of ADA, which are usually affinity matured IgG. It has been found that, in comparison with the nonaggregated form, aggregated mAb results in an increase in the amount of total peptides and the number of epitopes eluted from MHC II . This suggests that aggregates may increase immunogenicity by enhancing antigen processing and presentation in the T cell-dependent pathway. Aggregates could also contribute to T cell-dependent immunogenicity by increasing the danger signal for dendritic cell maturation. Consistent with this, a recent study suggested that aggregated mAb induces significantly higher dendritic cell maturation compared with unstressed mAb . Lastly, aggregates could form repetitively arranged B cell epitopes in a paracrystalline manner to cross-link B cell receptors (BCRs), which in turn will activate antigen-specific B cells to generate ADA, mostly IgM, via the T cell-independent pathway . However, the scarcity of clinical data and the difficulty to isolate the impact of aggregates from other immunogenicity risk factors are major impediments to understand the mechanisms of aggregate-induced immunogenicity.
Mathematical modeling offers the advantage of fast and cost-effective assessment and so it can be used in complement with experimental analysis to study immune responses [31–34]. It also provides quantitative means to dissect each component of a complex response for a deeper understanding of the mechanisms underlying aggregate-induced immunogenicity. Multiple mechanistic mathematical models have been previously developed to study immune responses against various pathogens. For example, antigen processing and presentation by APC and the activation of T helper cells by interactions between T cell receptors and MHC II-peptide complexes have been modeled and the simulation results agree with a variety of experimental data . A mathematical model was also developed for predicting the clonal selection of B cells and antibody production by plasma cells . The role of activation threshold and infections in the dynamics of autoimmune diseases has been studied mathematically as well [37, 38]. Mathematical models have been proposed and experimentally validated for T cell-dependent antibody responses to a wide range of antigens, including Haemophilus influenzae type b, hepatitis B virus, cancer antigens, and influenza A virus [39–43]. The T cell-independent activation of B cells by multivalent hapten-polymer has been modeled, where fitting to experimental data revealed that a minimum number of BCRs, in the range of 7 to 15, need to be cross-linked by a single multivalent ligand to stimulate a B cell [44, 45]. With regards to TPP-induced immunogenicity, several pharmacokinetics (PK) models have been developed to study the impact of ADA on mAb therapy . For example, by incorporating ADA-drug interactions into empirical PK modeling, we developed a PK/ADA model to quantitatively assess the extent and timing of ADA generation, affinity maturation, and ADA-mediated TPP elimination . More recently, we built a mechanistic, multiscale mathematical model of TPP-induced immunogenicity, recapitulating the key processes underlying T cell-dependent generation of ADA, such as antigen presentation, activation of immune cells, and production of ADA as well as in vivo disposition of ADA and TPP [47, 48]. This system-level model consists of a subcellular module for antigen presentation, a cellular module for immune system activation and antibody production, and a whole-body module for drug disposition. The model is able to reproduce key immunological phenomena such as antibody affinity maturation and enhanced secondary response [47, 48]. More importantly, a case study on immune response against adalimumab (a fully human anti-TNF alpha IgG1 mAb) showed reasonable agreement between model simulations and experimental observations [47, 48]. Owing to its flexibility and comprehensiveness this system-level model provides us with an ideal platform to probe mechanisms through which aggregates could generate immunogenicity.
In this study, we evaluate whether aggregates could induce T cell-independent or T cell-dependent immune response. In the former case, we model the interactions between multivalent aggregates and BCRs and examine the formation of signal-transducing complex (STC) under physiologically relevant conditions. For the latter case, we use our previously developed system-level model to investigate the impact of antigen processing and presentation, number and affinity of epitopes, and danger signal on ADA production due to aggregates.
2. Materials and Methods
2.1. Aggregates in the T Cell-Independent Pathway: Interactions between Multivalent Aggregates and BCRs
An aggregate () is assumed to be a homogeneous product formed by the combination of monomers, which gives it a valency of . The binding of to BCR is assumed to be sequential (see Figure 2(a) for an example with ) and can be represented by the following second-order reactions:where and are the th reaction’s binding and dissociation rates, respectively, and is the complex formed by binding of to BCRs. It is assumed that a BCR could bind to any free site on and dissociate from any bound site on . The above reactions can be described by the following ordinary differential equations that govern the time evolution of , , and BCR:We selected three (low, medium, and high) physiologically relevant levels for input parameters association constant () and initial concentration ().  is concentration at , as an initial condition for ordinary differential equations, which is estimated using the following equation: where [Ag] is the total TPP concentration, is the aggregation percentage in TPP, and is the valency of aggregates. [Ag] ranges from 500 to 105 pM based on 30 μg dose of interferon beta 1b and 40 mg dose of anti-TNF mAb adalimumab, respectively [29, 47–49]; spans from 2 to 15% based on a previous report on the characterization and quantitation of aggregates in recombinant human interferon beta drug products ; and varies from 10 to 100 based on the sizes of nonaggregated and aggregated TPP [18, 23, 29, 50, 51]. Taken together, the low and high levels of are 0.1 and 1500 pM, respectively. The association constant has been previously reported to be 10−7 pM−1 for antibodies with low intrinsic affinities and 10−3 pM−1 for affinity matured antibodies, and hence these were selected as low and high levels . The middle levels for total (12 pM) and (10−5 pM−1) are the geometric means of corresponding low and high levels. The rate of binding of an antigen to its corresponding BCR, , is relatively constant [52, 53], so we fixed it to 8.64 × 10−3 pM−1 day−1. By contrast, the rate of dissociation () is expected to increase with because the resistance of against torsion and bending grows due to the steric hindrance from progressive binding of BCRs . For simplicity we assume that decreases exponentially with and the base for exponential decay is 0.5 as previously identified while modeling interactions between multivalent hapten-polymer and BCRs . The initial BCR concentration is the product of number of BCRs per cell, B cell concentration, and percentage of antigen-specific B cells. The number of BCRs per cell and B cell concentration have been previously reported as ~105 and ~108 L−1, respectively [41, 44, 45, 47, 48]. Studies on the percentage of antigen-specific B cells are limited, but it has been reported to be <0.002% for vaccinia virus  and <1% for individual antigens . The above estimates were used to define the input range of BCR concentration at as an initial condition for the ordinary differential equations in the simulation.
In the model, the STC is the number of that cross-links at least BCRs as defined in [44, 45]:The model was simulated using the ordinary differential equation solver ode15s in MATLAB (The MathWorks, Inc., Natick, MA).
2.2. Aggregates in the T Cell-Dependent Pathway: Impact on Antigen Processing and Presentation and Danger Signal
For this analysis, we use our previously developed mechanistic, multiscale mathematical model for T cell-dependent ADA production [47, 48]. In this system-level model aggregates could contribute to increased ADA production by enhancing either the antigen processing and presentation or the danger signal for dendritic cell maturation (denoted by red arrows in Figure 3). We simulate the impact of aggregates by increasing (i) the rate of internalization of TPP into the endosome, (ii) the rate of degradation/processing of TPP into antigenic peptides, (iii) the number of epitopes generated, (iv) the affinity of epitopes to MHC II, and (v) the level of danger signal. Subsequently, for each of these conditions, we examine the endosomal levels of aggregates and epitope, the number of MHC II-peptide complexes on APC, and the levels of ADA production. To simulate B cell clonal selection and antibody affinity maturation, B cells and ADA are divided into 17 subgroups based on the binding affinity to antigen [36, 47, 48]. In our analysis, we define ADA production as the sum of the 17 subgroups.
3.1. Aggregates Are Unlikely to Induce T Cell-Independent Immune Response because the Number of STC Formed Is Limited
To evaluate whether aggregates could induce T cell-independent antibody responses through BCR cross-linking, we examine the number of STC formed per B cell for different parameter combinations (see Section 2 for details). The model output for interactions between aggregates and BCR is the STC formed per B cell, which was previously defined as the number of which cross-links a minimum number of BCRs [44, 45]. It has been reported that a multivalent ligand stimulates B cell activation only if it cross-links a minimum number () of BCRs, which is usually between 7 and 15 [44, 45]. We calculated the number of STC for , 5, and 10 under different total , , and BCR levels. Surprisingly, our computer simulation analysis showed that if or 5, no more than one STC per cell could be observed under physiological levels of total , BCR, and (data not shown). Even if is lowered to 2, more than one STC per cell can form only under limited conditions, when the sensitive parameters are near the upper limits of the physiologically plausible ranges (Figure 2(b)). In the case of pM−1, no more than one STC could form (Figure 2(b), left panel). For pM−1, more than one STC could form at high levels of total (1.5 × 10−3 pM) but only near the upper limit of antigen-specific B cells percentage (1%) (Figure 2(b), middle panel). Finally, when pM−1, more than one STC could form at all total levels but only with antigen-specific B cell percentage >0.006% (Figure 2(b), right panel). These results from our computer simulation showed that STC per cell is very sensitive to and total concentrations of and BCRs (but not to binding rate , data not shown). Overall, this analysis suggests that aggregates are unlikely to induce T cell-independent activation of B cells and consequent ADA production under physiologically plausible conditions. Therefore, aggregates may only contribute to ADA production through a T cell-dependent pathway, which we explore next.
3.2. Aggregates Could Enhance ADA Production by Increasing the Danger Signal to Maturate Dendritic Cells
To evaluate the T cell-dependent effect of aggregates on ADA production, we modulated those parameters in our system-level immunogenicity model [47, 48] that may be impacted by aggregation. This model consists of a subcellular module for antigen presentation in APC, a cellular module for immune cell activation and ADA production, and a whole-body module for drug and ADA disposition (Figure 3). Aggregates have been previously shown to increase danger signal for dendritic cell maturation and T cell activation [12, 22, 30, 56]. Specifically, aggregated mAb upregulated the dendritic cell maturation marker CD83 and CD4+ T cell costimulatory molecules CD80 and CD86 as well as cytokines produced by CD4+ T cells, such as IL-2 and IL-10 [30, 56]. Due to the complexity of dendritic cell maturation by danger signal and the unavailability of many parameters associated with this process, it is simply modeled as being driven by endotoxin lipopolysaccharide (LPS) [47, 48]. LPS is widely used in immunological studies for dendritic cell maturation [57–61] and is present in many TPP . The cytokine profiles induced by LPS and aggregates of mAb are very similar [22, 63]. Using our system-level model, we previously simulated ADA production induced by adalimumab, a fully anti-TNF alpha IgG1 mAb used to treat various inflammatory and autoimmune diseases, with a danger signal of 350 ng LPS  (Figure 4(e)). If aggregates increase the danger signal by 5-fold, ADA production is increased by 20-fold (Figure 4(f)). We also simulated ADA production for low danger signal levels (Figures 4(a)–4(d)) as the actual amount induced by nonaggregated TPP is unknown. In essence, ADA production depends on the level of danger signal (Figures 4(a)–4(f)). Therefore, our simulations suggest aggregates could enhance ADA production by increasing danger signal to enhance maturation of dendritic cells and subsequently activate T cells.
3.3. Aggregates Could Not Enhance ADA Production by Increasing Antigen Processing and Presentation If High Affinity Epitopes Are Already Present in Nonaggregated TPP
Antigen processing and presentation are the key events in T cell-dependent immunogenicity of TPP . Previous studies have demonstrated that aggregation enhances antigen’s uptake, processing, and presentation by APC [12, 22, 30, 56, 64]. More recently, a study showed that aggregated mAb could directly increase the total number of different peptides and the number of epitopes presented by MHC II compared with nonaggregated mAb . To evaluate whether aggregation-enhanced antigen processing and presentation could increase ADA production, we simulated these effects of aggregates in our model by changing its internalization or degradation rate or the number and affinity of epitopes generated and assessing their impact on final ADA production.
We previously simulated ADA production induced by adalimumab with an internalization rate of 14.4 day−1 (IR0), a degradation rate of 17.28 day−1 (DR0), and two predicted adalimumab epitopes with high binding affinities of 123 and 85 nM to common MHC II allele DRB104:01 . To model the aggregates’ effect on antigen processing, we increased either internalization or degradation rate by 16.6-fold based on a previous study which reported that aggregated mAb resulted in a 16.6-fold increase in total peptides associated with MHC II  and then assessed the levels of endosomal aggregates and epitopes, MHC II-peptide complexes on cell surface, and ADA production. As expected, conditional on the parameters and structure of the model simulation, increasing internalization rate by 16.6-fold resulted in a similar fold increase in aggregates internalized into endosome and epitopes generated by its degradation (Figures 5(a)-5(b) and 5(e)-5(f)). Increasing degradation rate by 16.6-fold resulted in the same fold decrease in endosomal aggregates, but the levels of epitopes were unchanged, which suggested that epitope generation was limited by the amount of aggregates internalized and not by the degradation rate (Figures 5(a)-5(b) and 5(i)-5(j)). Moreover, increasing internalization or degradation rate by 16.6-fold did not significantly change the number of MHC II-peptide complexes presented on the surface of APC (Figures 5(c), 5(g), and 5(k)). Aggregates could also impact the FcR binding and potentially affect the antigen uptake . We therefore evaluated a larger range of internalization and degradation rate. Our conclusions were unaffected by larger increases (200-fold) in internalization or degradation rate (data not shown). Consistent with MHC II-peptide complex presentation levels, increasing internalization or degradation rate by 16.6-fold had little impact on final ADA production (Figures 5(d), 5(h), and 5(l)). We next modeled the effect of aggregates on the number of epitopes presented. As expected, including aggregate-induced generation of new epitopes led to the surface presentation of corresponding MHC II-peptide complexes whose levels depend on the binding affinity of epitope to MHC II (Figures 6(a)–6(c), 6(e)–6(g), and 6(i)–6(k)). Surprisingly, if two high affinity epitopes are already present, then the inclusion of new epitopes did not increase ADA production (Figures 6(d), 6(h), and 6(l)). Taken together, these analyses suggest that aggregate-induced high antigen processing and presentation cannot enhance ADA production if high affinity epitopes are already present.
3.4. Aggregates Could Enhance ADA Production by Inducing the Presentation of Epitopes with Higher Affinities than Those from Nonaggregated TPP
MHC II-restricted epitopes are generated with μM to nM affinity range . We next evaluated whether aggregate-induced high antigen processing and presentation could increase immunogenicity when nonaggregated TPP present low affinity (μM range) epitopes. We started with 40 mg dose of nonaggregated TPP administered biweekly and two epitopes with of 1230 and 850 nM representing low affinity epitopes of μM range [65, 66] and monitored the number of MHC II-peptide complexes on surface of APC and ADA production (Figures 7(a)–7(d)). We next increased the internalization rate by 16.6-fold to mimic the effect of aggregates and again saw no increase in antigen presentation and ADA production (Figures 7(e)–7(h)). Notably, when aggregates induced the presentation of a high affinity epitope ( nM), ADA production increased by >4-fold (Figure 7(l)) due to enhanced antigen presentation (Figures 7(i)–7(k)). We further evaluated the effect of aggregate-induced high affinity epitopes on ADA production under different dose levels, all of which demonstrated that induction of a high affinity epitope could significantly increase ADA production (compare top and bottom rows in Figure 8), whereas increase in internalization rate had no effect (compare top and middle rows in Figure 8). These computational modeling results indicate that aggregates could contribute to ADA production by inducing the presentation of high affinity epitopes that may not be present in nonaggregated TPP.
In this study, we used mathematical modeling to comprehensively evaluate mechanisms through which aggregates of TPP could contribute to immunogenicity. By modeling the interactions between aggregates and BCRs, we find that aggregates are unlikely to induce T cell-independent antibody responses through BCR cross-linking due to the limited number of STC that could form under physiologically plausible conditions. Thereafter, using our previously developed multiscale, mechanistic mathematical model for the T cell-dependent induction of ADA by TPP, we systematically evaluated the potential roles of aggregates in ADA generation by dissecting the individual steps leading to it. Our analyses indicate that aggregates could contribute to immunogenicity by increasing the danger signal to maturate dendritic cells and activate T cells and/or by inducing the presentation of high affinity epitopes that may not be present in nonaggregated TPP.
TPP could aggregate during manufacturing, storage, handling, or delivery to patients due to agitation, light exposure, temperature elevation, oxidation, pH change, and leaching [12, 17, 23, 24, 29, 30, 56, 67]. Aggregation has been proposed as a strong risk factor for TPP-induced immunogenicity due to its potential role in both T cell-independent and T cell-dependent antibody responses [10, 12, 14, 16, 17]. Several previous studies in mice have demonstrated that for different TPP aggregates induced a stronger ADA production compared with nonaggregated forms [18–21, 25, 27, 28]. However, the mechanisms underlying aggregate-induced ADA production are not clear. A recent study in mice transgenic for human IgG demonstrated that only light-induced oligomers of IgG induced an immune response, which was ablated by the depletion of CD4+ cells . The data from this mouse model are in agreement with the mathematical model in which aggregates induce immune responses in a T cell-dependent manner.
Repetitively arranged epitopes in a paracrystalline structure of viral particles could cross-link BCRs to induce T cell-independent IgM or in some cases IgG3 responses [68–72]. It is expected that aggregates of TPP, potentially resembling the structure of highly repetitive epitopes, could induce T cell-independent antibody responses in a similar way [12, 16, 17]. The model does not directly consider the nature of a polyclonal B cell response, but it is consistent with that. Specifically, multiple epitopes from the aggregates being bound by the BCR can be represented by the differential binding rate constants in the model, and the different number of B cell epitopes on aggregates can be captured by the complex forming between aggregates and various number of BCRs. Surprisingly, by modeling the interactions between aggregates and BCRs, we find that aggregates are unlikely to induce T cell-independent antibody responses because only a few STC can form under physiologically plausible conditions for antigen-specific B cells, antigen dose, and binding affinity (Figure 2(b), left and center panels). This is consistent with previous studies in mice that showed no significant T cell-independent IgG3 antibody response against aggregated recombinant murine growth hormone  or anti-TNFα murine mAb , although IgM production was not evaluated in either case. However, it should be noted that, under conditions of high binding affinity and BCR concentration and appropriate antigen concentrations, significant number of STC could form, with a potential to induce T cell-independent antibody response (Figure 2(b), right panel). High BCR concentration can be achieved through high percentage of antigen-specific B cells for particular TPP or through B cell proliferation due to lowering of activation threshold by cytokines , second messenger diacylglycerol , costimulatory signal , or Bruton’s tyrosine kinase . Appropriate antigen concentration can result from specific dosing strategies. Therefore, particular attention should be given while administering TPP to patients in those conditions. Future experiments directly investigating the downstream signaling of BCR cross-linking in the presence of aggregated and nonaggregated TPP and studies evaluating whether T cell-independent IgM is induced in response to aggregates can further elucidate the role of aggregates in T cell-independent ADA production.
T cell-dependent ADA production is thought to be the major pathway through which TPP induce immunogenicity as in the case of IgG1 and IgG4 generated against anti-TNFα mAb to treat a variety of inflammatory and autoimmune diseases . Antigen processing and presentation by professional APC, such as dendritic cells, macrophages, and B cells, play a key role in T cell-dependent antibody responses . It has been shown that aggregates could enhance antigen uptake by APC thereby increasing peptides associated with MHC II and could induce dendritic cell maturation and T cell activation [22, 30, 56]. However, human data directly ascribing ADA levels to aggregates are still lacking. In this study, we systematically evaluated whether aggregate-enhanced antigen processing and presentation could increase ADA production. Our computer simulation suggests that the amount of antigenic peptides in endosome is limited by antigen internalization rate, not degradation rate, and the number of MHC II-peptide complexes presented on cell surface is mainly restricted by the binding affinity of epitopes. Our modeling analyses indicate that induction of high affinity epitopes by aggregates that may not be present in nonaggregated TPP and increased danger signal by aggregates to maturate dendritic cells could result in increased ADA production (Figures 4–8). A specifically designed experimental study that examines the binding affinities of peptides to MHC II derived from dendritic cells treated with aggregated or nonaggregated TPP would verify whether aggregates can induce the presentation of high affinity epitopes not present in nonaggregated TPP. In this work, we modeled aggregates-induced danger signal as LPS. However, it should be noted that aggregates have the potential to bind to a variety of pattern recognition receptors as well as FcR. Therefore the kinetics, activation thresholds, and receptors engaged by aggregates are more diverse and complicated than those of LPS and need further investigation.
This work improves our understanding of aggregate-induced immunogenicity and could be utilized to develop prediction and mitigation strategies. Overall, this modeling study suggests that aggregates could enhance immunogenicity; therefore enough attention should be given to reduce aggregation during manufacturing, storage, handling, and administration. In particular, potential high affinity CD4+ T cell epitopes are of great concern because their presentation in nonaggregated TPP will result in high levels of immunogenicity regardless of aggregation. On the other hand, even if they are not presented in nonaggregated TPP, an aggregation-induced presentation will also result in enhanced immunogenicity. Thus, efforts should be made towards experimental identification or in silico prediction of high affinity epitopes during immunogenicity assessment, and potential high affinity epitopes should be avoided while designing novel TPP as they carry a strong risk for ADA generation.
Our recently developed mechanistic system-level mathematical model for ADA production is a useful tool to evaluate human immunogenicity against TPP as it incorporates protein-specific antigenic properties and host-specific immunological characteristics, although further experimental validation is needed to increase confidence in ADA predictions [47, 48]. Multiple product- and patient-related risk factors have been proposed to impact immunogenicity of TPP [7, 8, 10–14, 77, 78]. As confidence in its properties increases, this system-level model could potentially be used to design new hypotheses and study other risk factors besides aggregation. For example, though the model is developed for healthy subjects, it can be easily modified to account for the effect of different disease statuses. For example, the profile of ADA generation observed in autoimmune patients [79, 80] can be simulated by including either a lower activation threshold for immune cells [37, 38] or preexisting immunity against TPP [79, 80]. Also, peptide editor HLA-DM plays a key role in MHC II antigen presentation and CD4+ T cell epitope selection by favoring the presentation of peptides with higher kinetic stabilities [65, 81–84]. To evaluate the effect of HLA-DM-mediated epitope selection on ADA production, it could be included in the subcellular module of antigen processing and presentation to select the epitopes presented based on peptide susceptibility to HLA-DM-mediated peptide exchange . Other ADA production impact factors that could be evaluated by this system-level model include time delays between administration, immune cell activation and migration from tissue to lymphoid compartments , contraction of effector B cells and T cells [85, 86], effect of immunomodulators through comedication , and different antibody isotypes generated by short- and long-lived plasma cells [42, 88, 89]. Therefore, this model could generate new hypotheses about immunogenicity and could be used with experiments to decipher the mechanisms underlying immunogenicity of TPP and develop corresponding mitigation strategies.
In summary, our computational analyses suggest that aggregates are unlikely to induce T cell-independent antibody responses through BCR cross-linking due to limited formation of STC under physiologically plausible conditions. In contrast, aggregates could contribute to immunogenicity via the T cell-dependent pathway by inducing the presentation of high affinity epitopes that may not be present in nonaggregated TPP and/or by enhancing danger signal to maturate dendritic cells and activate T cells. This study provides novel insights into how aggregates could contribute to overall immunogenicity and suggests novel mechanistic hypotheses eventually suitable for experimental testing.
Paolo Vicini is currently working at Clinical Pharmacology and DMPK, MedImmune, Cambridge CB21 6GH, UK.
Conflict of Interests
All authors are current or former employees of Pfizer Inc.
Xiaoying Chen and Abhinav Tiwari contributed to the work equally.
This work was supported by a Pfizer Worldwide Research and Development Postdoctoral Fellowship.
- D. S. Dimitrov, “Therapeutic proteins,” in Therapeutic Proteins, vol. 899 of Methods in Molecular Biology, pp. 1–26, Humana Press, 2012.
- G. Walsh, “Biopharmaceutical benchmarks 2014,” Nature Biotechnology, vol. 32, no. 10, pp. 992–1000, 2014.
- V. Brinks, D. Weinbuch, M. Baker et al., “Preclinical models used for immunogenicity prediction of therapeutic proteins,” Pharmaceutical Research, vol. 30, no. 7, pp. 1719–1728, 2013.
- R. T. Purcell and R. F. Lockey, “Immunologic responses to therapeutic biologic agents,” Journal of Investigational Allergology & Clinical Immunology, vol. 18, no. 5, pp. 335–342, 2008.
- V. Jawa, L. P. Cousens, M. Awwad, E. Wakshull, H. Kropshofer, and A. S. De Groot, “T-cell dependent immunogenicity of protein therapeutics: preclinical assessment and mitigation,” Clinical Immunology, vol. 149, no. 3, pp. 534–555, 2013.
- J. R. Maneiro, E. Salgado, and J. J. Gomez-Reino, “Immunogenicity of monoclonal antibodies against tumor necrosis factor used in chronic immune-mediated inflammatory conditions: systematic review and meta-analysis,” JAMA Internal Medicine, vol. 173, no. 15, pp. 1416–1428, 2013.
- P. A. van Schouwenburg, T. Rispens, and G. J. Wolbink, “Immunogenicity of anti-TNF biologic therapies for rheumatoid arthritis,” Nature Reviews Rheumatology, vol. 9, no. 3, pp. 164–172, 2013.
- G. Shankar, S. Arkin, L. Cocea et al., “Assessment and reporting of the clinical immunogenicity of therapeutic proteins and peptides—harmonized terminology and tactical recommendations,” The AAPS Journal, vol. 16, no. 4, pp. 658–673, 2014.
- A. S. De Groot and D. W. Scott, “Immunogenicity of protein therapeutics,” Trends in Immunology, vol. 28, no. 11, pp. 482–490, 2007.
- S. K. Singh, “Impact of product-related factors on immunogenicity of biotherapeutics,” Journal of Pharmaceutical Sciences, vol. 100, no. 2, pp. 354–387, 2011.
- C. Krieckaert, T. Rispens, and G. Wolbink, “Immunogenicity of biological therapeutics: from assay to patient,” Current Opinion in Rheumatology, vol. 24, no. 3, pp. 306–311, 2012.
- S. Sethu, K. Govindappa, M. Alhaidari, M. Pirmohamed, K. Park, and J. Sathish, “Immunogenicity to biologics: mechanisms, prediction and reduction,” Archivum Immunologiae et Therapiae Experimentalis, vol. 60, no. 5, pp. 331–344, 2012.
- A. C. Moss, V. Brinks, and J. F. Carpenter, “Review article: immunogenicity of anti-TNF biologics in IBD—the role of patient, product and prescriber factors,” Alimentary Pharmacology & Therapeutics, vol. 38, no. 10, pp. 1188–1197, 2013.
- S. Kumar, S. K. Singh, X. Wang, B. Rup, and D. Gill, “Coupling of aggregation and immunogenicity in biotherapeutics: T- and B-cell immune epitopes may contain aggregation-prone regions,” Pharmaceutical Research, vol. 28, no. 5, pp. 949–961, 2011.
- S. Kumar, M. A. Mitchell, B. Rup, and S. K. Singh, “Relationship between potential aggregation-prone regions and HLA-DR-binding T-cell immune epitopes: implications for rational design of novel and follow-on therapeutic antibodies,” Journal of Pharmaceutical Sciences, vol. 101, no. 8, pp. 2686–2701, 2012.
- K. D. Ratanji, J. P. Derrick, R. J. Dearman, and I. Kimber, “Immunogenicity of therapeutic proteins: influence of aggregation,” Journal of Immunotoxicology, vol. 11, no. 2, pp. 99–109, 2014.
- M. Sauerborn, V. Brinks, W. Jiskoot, and H. Schellekens, “Immunological mechanism underlying the immune response to recombinant human protein therapeutics,” Trends in Pharmacological Sciences, vol. 31, no. 2, pp. 53–59, 2010.
- S. Hermeling, L. Aranha, J. M. A. Damen et al., “Structural characterization and immunogenicity in wild-type and immune tolerant mice of degraded recombinant human interferon alpha2b,” Pharmaceutical Research, vol. 22, no. 12, pp. 1997–2006, 2005.
- S. Hermeling, H. Schellekens, C. Maas, M. F. B. G. Gebbink, D. J. A. Crommelin, and W. Jiskoot, “Antibody response to aggregated human interferon alpha2b in wild-type and transgenic immune tolerant mice depends on type and level of aggregation,” Journal of Pharmaceutical Sciences, vol. 95, no. 5, pp. 1084–1096, 2006.
- P. Human, H. Ilsley, C. Roberson et al., “Assessment of the immunogenicity of mechanically induced interferon aggregates in a transgenic mouse model,” Journal of Pharmaceutical Sciences, vol. 104, no. 2, pp. 722–730, 2015.
- V. Filipe, W. Jiskoot, A. H. Basmeleh, A. Halim, H. Schellekens, and V. Brinks, “Immunogenicity of different stressed IgG monoclonal antibody formulations in immune tolerant transgenic mice,” mAbs, vol. 4, no. 6, pp. 740–752, 2012.
- M. K. Joubert, M. Hokom, C. Eakin et al., “Highly aggregated antibody therapeutics can enhance the in vitro innate and late-stage T-cell immune responses,” The Journal of Biological Chemistry, vol. 287, no. 30, pp. 25266–25279, 2012.
- A. J. Freitag, M. Shomali, S. Michalakis et al., “Investigation of the immunogenicity of different types of aggregates of a murine monoclonal antibody in mice,” Pharmaceutical Research, vol. 32, no. 2, pp. 430–444, 2015.
- A. Seidl, O. Hainzl, M. Richter et al., “Tungsten-induced denaturation and aggregation of epoetin alfa during primary packaging as a cause of immunogenicity,” Pharmaceutical Research, vol. 29, no. 6, pp. 1454–1467, 2012.
- V. S. Purohit, C. R. Middaugh, and S. V. Balasubramanian, “Influence of aggregation on immunogenicity of recombinant human factor VIII in hemophilia A mice,” Journal of Pharmaceutical Sciences, vol. 95, no. 2, pp. 358–371, 2006.
- D. S. Pisal, M. P. Kosloski, C. R. Middaugh, R. B. Bankert, and S. V. Balu-Iyer, “Native-like aggregates of factor VIII are immunogenic in von Willebrand factor deficient and hemophilia a mice,” Journal of Pharmaceutical Sciences, vol. 101, no. 6, pp. 2055–2065, 2012.
- M. M. C. van Beers, M. Sauerborn, F. Gilli, V. Brinks, H. Schellekens, and W. Jiskoot, “Aggregated recombinant human interferon beta induces antibodies but no memory in immune-tolerant transgenic mice,” Pharmaceutical Research, vol. 27, no. 9, pp. 1812–1824, 2010.
- A. H. Fradkin, J. F. Carpenter, and T. W. Randolph, “Glass particles as an adjuvant: a model for adverse immunogenicity of therapeutic proteins,” Journal of Pharmaceutical Sciences, vol. 100, no. 11, pp. 4953–4964, 2011.
- J. G. Barnard, K. Babcock, and J. F. Carpenter, “Characterization and quantitation of aggregates and particles in interferon-β products: potential links between product quality attributes and immunogenicity,” Journal of Pharmaceutical Sciences, vol. 102, no. 3, pp. 915–928, 2013.
- V. Rombach-Riegraf, A. C. Karle, B. Wolf et al., “Aggregation of human recombinant monoclonal antibodies influences the capacity of dendritic cells to stimulate adaptive T-cell responses in vitro,” PLoS ONE, vol. 9, no. 1, Article ID e86322, 2014.
- C. Lundegaard, O. Lund, C. Keşmir, S. Brunak, and M. Nielsen, “Modeling the adaptive immune system: predictions and simulations,” Bioinformatics, vol. 23, no. 24, pp. 3265–3275, 2007.
- J. D. Gómez-Mantilla, I. F. Trocóniz, Z. Parra-Guillén, and M. J. Garrido, “Review on modeling anti-antibody responses to monoclonal antibodies,” Journal of Pharmacokinetics and Pharmacodynamics, vol. 41, no. 5, pp. 523–536, 2014.
- T. P. Hickling, X. Chen, P. Vicini, and S. Nayak, “A review of quantitative modeling of B cell responses to antigenic challenge,” Journal of Pharmacokinetics and Pharmacodynamics, vol. 41, no. 5, pp. 445–459, 2014.
- S. Palsson, T. P. Hickling, E. L. Bradshaw-Pierce et al., “The development of a fully-integrated immune response model (FIRM) simulator of the immune response through integration of multiple subset models,” BMC Systems Biology, vol. 7, article 95, 2013.
- N. G. B. Agrawal and J. J. Linderman, “Mathematical modeling of helper T lymphocyte/antigen-presenting cell interactions: analysis of methods for modifying antigen processing and presentation,” Journal of Theoretical Biology, vol. 182, no. 4, pp. 487–504, 1996.
- G. I. Bell, “Mathematical model of clonal selection and antibody production,” Journal of Theoretical Biology, vol. 29, no. 2, pp. 191–232, 1970.
- K. B. Blyuss and L. B. Nicholson, “The role of tunable activation thresholds in the dynamics of autoimmunity,” Journal of Theoretical Biology, vol. 308, pp. 45–55, 2012.
- K. B. Blyuss and L. B. Nicholson, “Understanding the roles of activation threshold and infections in the dynamics of autoimmune disease,” Journal of Theoretical Biology, vol. 375, pp. 13–20, 2015.
- M. Oprea and A. S. Perelson, “Exploring the mechanisms of primary antibody responses to T cell-dependent antigens,” Journal of Theoretical Biology, vol. 181, no. 3, pp. 215–236, 1996.
- A. Rundell, R. DeCarlo, H. HogenEsch, and P. Doerschuk, “The humoral immune response to Haemophilus influenzae type b: a mathematical model based on T-zone and germinal center B-cell dynamics,” Journal of Theoretical Biology, vol. 194, no. 3, pp. 341–381, 1998.
- F. Castiglione, F. Toschi, M. Bernaschi et al., “Computational modeling of the immune response to tumor antigens,” Journal of Theoretical Biology, vol. 237, no. 4, pp. 390–400, 2005.
- H. Y. Lee, D. J. Topham, S. Y. Park et al., “Simulation and prediction of the adaptive immune response to influenza A virus infection,” Journal of Virology, vol. 83, no. 14, pp. 7151–7165, 2009.
- S. M. Ciupe, R. M. Ribeiro, and A. S. Perelson, “Antibody responses during hepatitis B viral infection,” PLoS Computational Biology, vol. 10, no. 7, Article ID e1003730, 2014.
- B. Sulzer and A. S. Perelson, “Equilibrium binding of multivalent ligands to cells: effects of cell and receptor density,” Mathematical Biosciences, vol. 135, no. 2, pp. 147–185, 1996.
- B. Sulzer and A. S. Perelson, “Immunons revisited: binding of multivalent antigens to b cells,” Molecular Immunology, vol. 34, no. 1, pp. 63–74, 1997.
- X. Chen, T. Hickling, E. Kraynov, B. Kuang, C. Parng, and P. Vicini, “A mathematical model of the effect of immunogenicity on therapeutic protein pharmacokinetics,” The AAPS Journal, vol. 15, no. 4, pp. 1141–1154, 2013.
- X. Chen, T. P. Hickling, and P. Vicini, “A mechanistic, multiscale mathematical model of immunogenicity for therapeutic proteins: part 2—model applications,” CPT: Pharmacometrics & Systems Pharmacology, vol. 3, no. 9, article e134, 10 pages, 2014.
- X. Chen, T. P. Hickling, and P. Vicini, “A mechanistic, multiscale mathematical model of immunogenicity for therapeutic proteins: part 1—theoretical model,” CPT: Pharmacometrics & Systems Pharmacology, vol. 3, no. 9, pp. 1–9, 2014.
- S. E. Grossberg, J. Oger, L. D. Grossberg, A. Gehchan, and J. P. Klein, “Frequency and magnitude of interferon beta neutralizing antibodies in the evaluation of interferon beta immunogenicity in patients with multiple sclerosis,” Journal of Interferon & Cytokine Research, vol. 31, no. 3, pp. 337–344, 2011.
- L. O. Narhi, J. Schmit, K. Bechtold-Peters, and D. Sharma, “Classification of protein aggregates,” Journal of Pharmaceutical Sciences, vol. 101, no. 2, pp. 493–498, 2012.
- M. Reth, “Matching cellular dimensions with molecular sizes,” Nature Immunology, vol. 14, no. 8, pp. 765–767, 2013.
- J. Foote and H. N. Eisen, “Kinetic and affinity limits on antibodies produced during immune responses,” Proceedings of the National Academy of Sciences of the United States of America, vol. 92, no. 5, pp. 1254–1256, 1995.
- J. Foote and C. Milstein, “Kinetic maturation of an immune response,” Nature, vol. 352, no. 6335, pp. 530–532, 1991.
- S. Crotty, P. Felgner, H. Davies, J. Glidewell, L. Villarreal, and R. Ahmed, “Cutting edge: long-term B cell memory in humans after smallpox vaccination,” The Journal of Immunology, vol. 171, no. 10, pp. 4969–4973, 2003.
- A. P. Kodituwakku, C. Jessup, H. Zola, and D. M. Roberton, “Isolation of antigen-specific B cells,” Immunology and Cell Biology, vol. 81, no. 3, pp. 163–170, 2003.
- M. Ahmadi, C. J. Bryson, E. A. Cloake et al., “Small amounts of sub-visible aggregates enhance the immunogenic potential of monoclonal antibody therapeutics,” Pharmaceutical Research, vol. 32, no. 4, pp. 1383–1394, 2015.
- C. Buelens, V. Verhasselt, D. De Groote, K. Thielemans, M. Goldman, and F. Willems, “Human dendritic cell responses to lipopolysaccharide and CD40 ligation are differentially regulated by interleukin-10,” European Journal of Immunology, vol. 27, no. 8, pp. 1848–1852, 1997.
- F. Granucci, E. Ferrero, M. Foti, D. Aggujaro, K. Vettoretto, and P. Ricciardi-Castagnoli, “Early events in dendritic cell maturation induced by LPS,” Microbes and Infection, vol. 1, no. 13, pp. 1079–1084, 1999.
- R. M. Cisco, Z. Abdel-Wahab, J. Dannull et al., “Induction of human dendritic cell maturation using transfection with RNA encoding a dominant positive toll-like receptor 4,” Journal of Immunology, vol. 172, no. 11, pp. 7162–7168, 2004.
- K. Kawamura, K. Iyonaga, H. Ichiyasu, J. Nagano, M. Suga, and Y. Sasaki, “Differentiation, maturation, and survival of dendritic cells by osteopontin regulation,” Clinical and Diagnostic Laboratory Immunology, vol. 12, no. 1, pp. 206–212, 2005.
- Y.-C. Wang, X.-B. Hu, F. He et al., “Lipopolysaccharide-induced maturation of bone marrow-derived dendritic cells is regulated by notch signaling through the up-regulation of CXCR4,” The Journal of Biological Chemistry, vol. 284, no. 23, pp. 15993–16003, 2009.
- D. Verthelyi and V. Wang, “Trace levels of innate immune response modulating impurities (IIRMIs) synergize to break tolerance to therapeutic proteins,” PLoS ONE, vol. 5, no. 12, Article ID e15252, 2010.
- M. Rossol, H. Heine, U. Meusch et al., “LPS-induced cytokine production in human monocytes and macrophages,” Critical Reviews in Immunology, vol. 31, no. 5, pp. 379–446, 2011.
- J. C. Jones, E. W. Settles, C. R. Brandt, and S. Schultz-Cherry, “Virus aggregating peptide enhances the cell-mediated response to influenza virus vaccine,” Vaccine, vol. 29, no. 44, pp. 7696–7703, 2011.
- L. Yin, J. M. Calvo-Calle, O. Dominguez-Amorocho, and L. J. Ster, “HLA-DM constrains epitope selection in the human CD4 T cell response to vaccinia virus by favoring the presentation of peptides with longer HLA-DM-mediated half-lives,” Journal of Immunology, vol. 189, no. 8, pp. 3983–3994, 2012.
- J. Bessa, S. Boeckle, H. Beck et al., “The immunogenicity of antibody aggregates in a novel transgenic mouse model,” Pharmaceutical Research, vol. 32, no. 7, pp. 2344–2359, 2015.
- W. Wang, “Protein aggregation and its inhibition in biopharmaceutics,” International Journal of Pharmaceutics, vol. 289, no. 1-2, pp. 1–30, 2005.
- M. F. Bachmann and R. M. Zinkernagel, “Neutralizing antiviral B cell responses,” Annual Review of Immunology, vol. 15, pp. 235–270, 1997.
- E. Szomolanyi-Tsuda and R. M. Welsh, “T-cell-independent antiviral antibody responses,” Current Opinion in Immunology, vol. 10, no. 4, pp. 431–435, 1998.
- C. Babin, N. Majeau, and D. Leclerc, “Engineering of papaya mosaic virus (PapMV) nanoparticles with a CTL epitope derived from influenza NP,” Journal of Nanobiotechnology, vol. 11, article 10, 2013.
- C. M. Snapper, T. M. McIntyre, R. Mandler et al., “Induction of IgG3 secretion by interferon γ: a model for T cell-independent class switching in response to T cell-independent type 2 antigens,” The Journal of Experimental Medicine, vol. 175, no. 5, pp. 1367–1371, 1992.
- T. Fehr, M. F. Bachmann, E. Bucher et al., “Role of repetitive antigen patterns for induction of antibodies against antibodies,” The Journal of Experimental Medicine, vol. 185, no. 10, pp. 1785–1792, 1997.
- P. K. A. Mongini, P. F. Highet, and J. K. Inman, “Human B cell activation: effect of T cell cytokines on the physicochemical binding requirements for achieving cell cycle progression via the membrane IgM signaling pathway,” Journal of Immunology, vol. 155, no. 7, pp. 3385–3400, 1995.
- M. L. Wheeler, M. B. Dong, R. Brink, X.-P. Zhong, and A. L. DeFranco, “Diacylglycerol kinase zeta limits B cell antigen receptor-dependent activation of ERK Signaling to inhibit early antibody responses,” Science Signaling, vol. 6, no. 297, article ra91, 2013.
- G. G. B. Klaus, M. Holman, C. Johnson-Léger, J. R. Christenson, and M. R. Kehry, “Interaction of B cells with activated T cells reduces the threshold for CD40-mediated B cell activation,” International Immunology, vol. 11, no. 1, pp. 71–79, 1999.
- L. P. Kil, M. J. W. de Bruijn, M. van Nimwegen et al., “Btk levels set the threshold for B-cell activation and negative selection of autoreactive B cells in mice,” Blood, vol. 119, no. 16, pp. 3744–3756, 2012.
- D. W. Scott and A. S. De Groot, “Can we prevent immunogenicity of human protein drugs?” Annals of the Rheumatic Diseases, vol. 69, supplement 1, pp. i72–i76, 2010.
- V. Brinks, W. Jiskoot, and H. Schellekens, “Immunogenicity of therapeutic proteins: the use of animal models,” Pharmaceutical Research, vol. 28, no. 10, pp. 2379–2385, 2011.
- L. Xue, M. Fiscella, M. Rajadhyaksha et al., “Pre-existing biotherapeutic-reactive antibodies: survey results within the american association of pharmaceutical scientists,” The AAPS Journal, vol. 15, no. 3, pp. 852–855, 2013.
- L. Xue and B. Rup, “Evaluation of pre-existing antibody presence as a risk factor for posttreatment anti-drug antibody induction: analysis of human clinical study data for multiple biotherapeutics,” The AAPS Journal, vol. 15, no. 3, pp. 893–896, 2013.
- V. S. Sloan, P. Cameron, G. Porter et al., “Mediation by HLA-DM of dissociation of peptides from HLA-DR,” Nature, vol. 375, no. 6534, pp. 802–806, 1995.
- N. K. Nanda and A. J. Sant, “DM determines the cryptic and immunodominant fate of T cell epitopes,” The Journal of Experimental Medicine, vol. 192, no. 6, pp. 781–788, 2000.
- A. J. Sant, F. A. Chaves, S. A. Jenks et al., “The relationship between immunodominance, DM editing, and the kinetic stability of MHC class II:peptide complexes,” Immunological Reviews, vol. 207, pp. 261–278, 2005.
- L. Yin, P. Trenh, A. Guce et al., “Susceptibility to HLA-DM protein is determined by a dynamic conformation of major histocompatibility complex class II molecule bound with peptide,” The Journal of Biological Chemistry, vol. 289, no. 34, pp. 23449–23464, 2014.
- K. R. Garrod, H. D. Moreau, Z. Garcia et al., “Dissecting T cell contraction in vivo using a genetically encoded reporter of apoptosis,” Cell Reports, vol. 2, no. 5, pp. 1438–1447, 2012.
- H. Lee, S. Haque, J. Nieto et al., “A p53 axis regulates B cell receptor-triggered, innate immune system-driven B cell clonal expansion,” Journal of Immunology, vol. 188, no. 12, pp. 6093–6108, 2012.
- C. L. M. Krieckaert, G. M. Bartelds, W. F. Lems, and G. J. Wolbink, “The effect of immunomodulators on the immunogenicity of TNF-blocking therapeutic monoclonal antibodies: a review,” Arthritis Research & Therapy, vol. 12, no. 5, article 217, 2010.
- A. Radbruch, G. Muehlinghaus, E. O. Luger et al., “Competence and competition: the challenge of becoming a long-lived plasma cell,” Nature Reviews Immunology, vol. 6, no. 10, pp. 741–750, 2006.
- C. H. Rozanski, R. Arens, L. M. Carlson et al., “Sustained antibody responses depend on CD28 function in bone marrow-resident plasma cells,” The Journal of Experimental Medicine, vol. 208, no. 7, pp. 1435–1446, 2011.
Copyright © 2015 Liusong Yin et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.