Search

Trastuzumab CRD Calculation Details

Description of common pharmacokinetic parameters used in the CRD calculations.

Dasatinib CRD Calculation Details

Vismodegib CRD Calculation Details

Purpose: The translation of nonclinical oncology studies is a subject of continuous debate. We propose that translational oncology studies need to optimize both pharmacokinetic (drug exposure) and pharmacodynamic (xenograft model) aspects. While improvements in pharmacodynamic translatability can be obtained by choosing cell lines or patient-derived xenograft models closer to the clinical indication, significant ambiguity and variability exists when optimizing the pharmacokinetic translation of small molecule and biotherapeutic agents.Experimental Design and Results: In this work, we propose a pharmacokinetic-based strategy to select nonclinical doses for approved drug molecules. We define a clinically relevant dose (CRD) as the dosing regimen in mice that most closely approximates the relevant pharmacokinetic metric in humans. Such metrics include area under the time–concentration curve and maximal or minimal concentrations within the dosing interval. The methodology is applied to six drugs, including targeted agents and chemotherapeutics, small and large molecules (erlotinib, dasatinib, vismodegib, trastuzumab, irinotecan, and capecitabine). The resulting efficacy response at the CRD is compared with clinical responses.Conclusions: We conclude that nonclinical studies designed with the appropriate CRDs of approved drug molecules will maximize the translatability of efficacy results, which is critical when testing approved and investigational agents in combination. Clin Cancer Res; 23(4); 1080–90. ©2016 AACR.

Full Description with references for the data included in Figure 4.

Tremelimumab, an anti‐cytotoxic T‐lymphocyte antigen‐4 monoclonal antibody that enhances T‐cell activation, was evaluated in a randomized, double‐blind, placebo‐controlled, phase IIb study (NCT01843374) in patients with unresectable malignant mesothelioma. The study demonstrated no clinically meaningful differences in overall survival (OS). The objective of this analysis was to evaluate the relationship of exposure with OS. A population pharmacokinetic (PK) model adequately described the PK data. Three factors (sex, C‐reactive protein, and baseline tumor size) were identified as statistically significant PK predictors (P

Linear models are generally reliable methods for analyzing tumor growth in vivo, with drug effectiveness being represented by the steepness of the regression slope. With immunotherapy, however, not all tumor growth follows a linear pattern, even after log transformation. Tumor kinetics models are mechanistic models that describe tumor proliferation and tumor killing macroscopically, through a set of differential equations. In drug combination studies, although an additional drug-drug interaction term can be added to such models, however, the drug interactions suggested by tumor kinetics models cannot be translated directly into synergistic effects. We have developed a novel statistical approach that simultaneously models tumor growth in control, monotherapy, and combination therapy groups. This approach makes it possible to test for synergistic effects directly and to compare such effects among different studies.

A sequential pharmacokinetic (PK) and pharmacodynamic (PD) model was built with Nonlinear Mixed Effects Modelling based on data from a first-in-human trial of a novel biologic, MEDI7836. MEDI7836 is a human immunoglobulin G1 lambda (IgG1λ-YTE) monoclonal antibody, with an Fc modification to reduce metabolic clearance. MEDI7836 specifically binds to, and functionally neutralizes interleukin-13. Thirty-two healthy male adults were enrolled into a dose-escalation clinical trial. Four active doses were tested (30, 105, 300, and 600 mg) with 6 volunteers enrolled per cohort. Eight volunteers received placebo as control. Following single subcutaneous administration (SC), individual time courses of serum MEDI7836 concentrations, and the resulting serum IL13 modulation in vivo, were quantified. A binding pharmacokinetic-pharmacodynamic (PK-PD) indirect response model was built to characterize the exposure-driven modulation of the target over time by MEDI7836. While the validated bioanalytical assay specification quantified the level of free target (i.e., a free IL13 assay), emerging clinical data suggested dose-dependent increase in systemic IL13 concentration over time, indicative of a total IL13 assay. The target time course was modelled as a linear combination of free target and a percentage of the drug-target complex to fit the clinical data. This novel PK-PD modelling approach integrates independent knowledge about the assay characteristics to successfully elucidate apparently complex observations.

A 2‐day meeting was held by members of the UK Quantitative Systems Pharmacology Network (http://www.qsp‐uk.net/) in November 2018 on the topic of Translational Challenges in Oncology. Participants from a wide range of backgrounds were invited to discuss current and emerging modeling applications in nonclinical and clinical drug development, and to identify areas for improvement. This resulting perspective explores opportunities for impactful quantitative pharmacology approaches. Four key themes arose from the presentations and discussions that were held, leading to the following recommendations: • Evaluate the predictivity and reproducibility of animal cancer models through precompetitive collaboration. • Apply mechanism of action (MoA) based mechanistic models derived from nonclinical data to clinical trial data. • Apply MoA reflective models across trial data sets to more robustly quantify the natural history of disease and response to differing interventions. • Quantify more robustly the dose and concentration dependence of adverse events through mathematical modelling techniques and modified trial design.

Drug development is a lengthy, costly process with low probability of success. Biopharmaceuticals are highly specific molecules, with efficacy and safety closely tied to target biology and pharmacology. The "learning-predicting-confirming" continuum by translational and clinical modeling and simulation (M&S) was implemented at every decision point for mavrilimumab, a human monoclonal antibody in development for rheumatoid arthritis (RA). This tutorial uses mavrilimumab as an example to demonstrate rational discovery, preclinical development, clinical study design, and dose selection of biotherapeutics by M&S.

The use of peptides as therapeutic agents is undergoing a renaissance with the expectation of new drugs with enhanced levels of efficacy and safety. Their clinical potential will be only fully realised once their physicochemical and pharmacokinetic properties have been precisely controlled. Here we demonstrate a reversible peptide self-assembly strategy to control and prolong the bioactivity of a native peptide hormone in vivo. We show that oxyntomodulin, a peptide with potential to treat obesity and diabetes, self-assembles into a stable nanofibril formulation which subsequently dissociates to release active peptide and produces a pharmacological effect in vivo. The subcutaneous administration of the nanofibrils in rats results in greatly prolonged exposure, with a constant oxyntomodulin bioactivity detectable in serum for at least 5 days as compared to free oxyntomodulin which is undetectable after only 4 h. Such an approach is simple, cost-efficient and generic in addressing the limitations of peptide therapeutics., The clinical potential of peptide therapeutic agents can only be fully realised once their physicochemical and pharmacokinetic properties are precisely controlled. Here the authors show a reversible peptide self-assembly strategy to control and prolong the bioactivity of a native peptide hormone in vivo.

Mutant epidermal growth factor receptor (EGFR) is a major driver of non-small-cell lung cancer (NSCLC). Marketed first generation inhibitors, such as erlotinib, effect a transient beneficial response in EGFR mutant NSCLC patients before resistance mechanisms render these inhibitors ineffective. Secondary oncogenic EGFR mutations account for approximately 50% of relapses, the most common being the gatekeeper T790M substitution that renders existing therapies ineffective. The discovery of PF-06459988 (1), an irreversible pyrrolopyrimidine inhibitor of EGFR T790M mutants, was recently disclosed.1 Herein, we describe our continued efforts to achieve potency across EGFR oncogenic mutations and improved kinome selectivity, resulting in the discovery of clinical candidate PF-06747775 (21), which provides potent EGFR activity against the four common mutants (exon 19 deletion (Del), L858R, and double mutants T790M/L858R and T790M/Del), selectivity over wild-type EGFR, and desirable ADME properties. Compound 21 is ...

Essentials Baseline coagulation activity can be detected in non-bleeding state by in vivo biomarker levels. A detailed mathematical model of coagulation was developed to describe the non-bleeding state. Optimized model described in vivo biomarkers with recombinant activated factor VII treatment. Sensitivity analysis predicted prothrombin fragment 1 + 2 and D-dimer are regulated differently. SummaryBackground Prothrombin fragment 1 + 2 (F1 + 2), thrombin–antithrombin III complex (TAT) and D-dimer can be detected in plasma from non-bleeding hemostatically normal subjects or hemophilic patients. They are often used as safety or pharmacodynamic biomarkers for hemostatis-modulating therapies in the clinic, and provide insights into in vivo coagulation activity. Objectives To develop a quantitative systems pharmacology (QSP) model of the blood coagulation network to describe in vivo biomarkers, including F1 + 2, TAT, and D-dimer, under non-bleeding conditions. Methods The QSP model included intrinsic and extrinsic coagulation pathways, platelet activation state-dependent kinetics, and a two-compartment pharmacokinetics model for recombinant activated factor VII (rFVIIa). Literature data on F1 + 2 and D-dimer at baseline and changes with rFVIIa treatment were used for parameter optimization. Multiparametric sensitivity analysis (MPSA) was used to understand key proteins that regulate F1 + 2, TAT and D-dimer levels. Results The model was able to describe tissue factor (TF)-dependent baseline levels of F1 + 2, TAT and D-dimer in a non-bleeding state, and their increases in hemostatically normal subjects and hemophilic patients treated with different doses of rFVIIa. The amount of TF required is predicted to be very low in a non-bleeding state. The model also predicts that these biomarker levels will be similar in hemostatically normal subjects and hemophilic patients. MPSA revealed that F1 + 2 and TAT levels are highly correlated, and that D-dimer is more sensitive to the perturbation of coagulation protein concentrations. Conclusions A QSP model for non-bleeding baseline coagulation activity was established with data from clinically relevant in vivo biomarkers at baseline and changes in response to rFVIIa treatment. This model will provide future mechanistic insights into this system.

Immunotherapy and immune checkpoint blocking antibodies such as anti-PD-1 are approved and significantly improve the survival of advanced non-small cell lung cancer (NSCLC) patients, but there has been little success in identifying biomarkers capable of separating the responders from non-responders before the onset of the therapy. In this study, we developed a quantitative system pharmacology (QSP) model to represent the anti-tumor immune response in human NSCLC that integrated our knowledge of tumor growth, antigen processing and presentation, T cell activation and distribution, antibody pharmacokinetics, and immune checkpoint dynamics. The model was calibrated with the available data and was used to identify potential biomarkers as well as patient-specific response based on the patient parameters. The model predicted that in addition to tumor mutational burden (TMB), a known biomarker for anti-PD-1 therapy in NSCLC, the number of effector T cells and regulatory T cells in the tumor and blood is a predictor of the responders. Furthermore, the model simulated a set of 12 patients with known TMB and MHC/antigen-binding affinity from a recent clinical trial (ClinicalTrials.gov number, NCT02259621) on neoadjuvant nivolumab therapy in resectable lung cancer and predicted an augmented durable response in patients with adjuvant nivolumab treatment in addition to the clinical trial protocol of neoadjuvant nivolumab treatment followed by resection. Overall, the model provides a valuable framework to model tumor immunity and response to immune checkpoint blockers to enhance biomarker discovery and performing virtual clinical trials to aid in design and interpretation of the current trials with fewer patients. Electronic supplementary material The online version of this article (10.1208/s12248-019-0350-x) contains supplementary material, which is available to authorized users.

Precision medicine aims to use patient genomic, epigenomic, specific drug dose, and other data to define disease patterns that may potentially lead to an improved treatment outcome. Personalized dosing regimens based on tumor drug penetration can play a critical role in this approach. State-of-the-art techniques to measure tumor drug penetration focus on systemic exposure, tissue penetration, cellular or molecular engagement, and expression of pharmacological activity. Using in silico methods, this information can be integrated to bridge the gap between the therapeutic regimen and the pharmacological link with clinical outcome. These methodologies are described, and challenges ahead are discussed. Supported by many examples, this review shows how the combination of these techniques provides enhanced patient-specific information on drug accessibility at the tumor tissue level, target binding, and downstream pharmacology. Our vision of how to apply tumor drug penetration measurements offers a roadmap for the clinical implementation of precision dosing., SCOPUS: re.j, info:eu-repo/semantics/published

The low response rate of immune checkpoint blockade in breast cancer has highlighted the need for predictive biomarkers to identify responders. While a number of clinical trials are ongoing, testing all possible combinations is not feasible. In this study, a quantitative systems pharmacology model is built to integrate immune–cancer cell interactions in patients with breast cancer, including central, peripheral, tumour-draining lymph node (TDLN) and tumour compartments. The model can describe the immune suppression and evasion in both TDLN and the tumour microenvironment due to checkpoint expression, and mimic the tumour response to checkpoint blockade therapy. We investigate the relationship between the tumour response to checkpoint blockade therapy and composite tumour burden, PD-L1 expression and antigen intensity, including their individual and combined effects on the immune system, using model-based simulations. The proposed model demonstrates the potential to make predictions of tumour response of individual patients given sufficient clinical measurements, and provides a platform that can be further adapted to other types of immunotherapy and their combination with molecular-targeted therapies. The patient predictions demonstrate how this systems pharmacology model can be used to individualize immunotherapy treatments. When appropriately validated, these approaches may contribute to optimization of breast cancer treatment.

Over the past decade, several immunotherapies have been approved for the treatment of melanoma. The most prominent of these are the immune checkpoint inhibitors, which are antibodies that block the inhibitory effects on the immune system by checkpoint receptors, such as CTLA-4, PD-1 and PD-L1. Preclinically, blocking these receptors has led to increased activation and proliferation of effector cells following stimulation and antigen recognition, and subsequently, more effective elimination of cancer cells. Translation from preclinical to clinical outcomes in solid tumors has shown the existence of a wide diversity of individual patient responses, linked to several patient-specific parameters. We developed a quantitative systems pharmacology (QSP) model that looks at the mentioned checkpoint blockade therapies administered as mono-, combo- and sequential therapies, to show how different combinations of specific patient parameters defined within physiological ranges distinguish different types of virtual patient responders to these therapies for melanoma. Further validation by fitting and subsequent simulations of virtual clinical trials mimicking actual patient trials demonstrated that the model can capture a wide variety of tumor dynamics that are observed in the clinic and can predict median clinical responses. Our aim here is to present a QSP model for combination immunotherapy specific to melanoma.

The lack of standardization in the way that quantitative and systems pharmacology (QSP) models are developed, tested, and documented hinders their reproducibility, reusability, and expansion or reduction to alternative contexts. This in turn undermines the potential impact of QSP in academic, industrial, and regulatory frameworks. This article presents a minimum set of recommendations from the UK Quantitative and Systems Pharmacology Network (UK QSP Network) to guide QSP practitioners seeking to maximize their impact, and stakeholders considering the use of QSP models in their environment.

International audience; Hierarchical processes spanning several orders of magnitude of both space and time underlie nearly all cancers. Multi-scale statistical, mathematical, and computational modeling methods are central to designing, implementing and assessing treatment strategies that account for these hierarchies. The basic science underlying these modeling efforts is maturing into a new discipline that is close to influencing and facilitating clinical successes. The purpose of this review is to capture the state-of-the-art as well as the key barriers to success for multi-scale modeling in clinical oncol- ogy. We begin with a summary of the long-envisioned promise of multi-scale modeling in clinical oncology, including the synthesis of disparate data types into models that reveal underlying mechanisms and allow for experimental testing of hypotheses. We then evaluate the mathematical techniques employed most widely and present several examples illustrat- ing their application as well as the current gap between pre- clinical and clinical applications. We conclude with a discus- sion of what we view to be the key challenges and opportunities for multi-scale modeling in clinical oncology.

Recent advances in molecular biology have enabled personalized cancer therapies with molecularly targeted agents (MTAs), which offer a promising future for cancer therapy. Dynamic modeling and simulation (MS) is a powerful mathematical approach linking drug exposures to pharmacological responses, providing a quantitative assessment of in vivo drug potency. Accordingly, a growing emphasis is being placed upon MS to quantitatively understand therapeutic exposure-response relationships of MTAs in nonclinical models.An overview of MS approaches for MTAs in nonclinical models is presented with discussion about mechanistic extrapolation of antitumor efficacy from bench to bedside. Emphasis is placed upon recent advances in MS approaches linking drug exposures, biomarker responses (e.g. target modulation) and pharmacological outcomes (e.g. antitumor efficacy).For successful personalized cancer therapies with MTAs, it is critical to mechanistically and quantitatively understand their exposure-response relationships in nonclinical models, and to logically and properly apply such knowledge to the clinic. Particularly, MS approaches to predict pharmacologically active concentrations of MTAs in patients based upon nonclinical data would be highly valuable in guiding the design and execution of clinical trials. Proactive approaches to understand their exposure-response relationships could substantially increase probability of achieving a positive proof-of-concept in the clinic.

High throughput molecular and functional profiling of patients is a key driver of precision medicine. DNA and RNA characterization has been enabled at unprecedented cost and scale through rapid, disruptive progress in sequencing technology, but challenges persist in data management and interpretation. We analyze the state-of-the-art of large-scale unbiased sequencing in drug discovery and development, including technology, application, ethical, regulatory, policy and commercial considerations, and discuss issues of LUS implementation in clinical and regulatory practice.

Background In drug development, in vivo assessment of target engagement provides confidence when testing the drug's mechanism of action and improves the likelihood of clinical success. For biologics, receptor occupancy (RO) determined from circulating cells can provide evidence of target engagement. Integrating this information with mathematical modeling can further enhance the understanding of drug-target interactions and the biological factors that are critical to the successful modulation of the target and ultimately the disease state. Methods This mini-review presents two specific types of mathematical models used to describe antibody-receptor systems and highlights how experimental data can inform the model parameters. Simulations are used to illustrate how various mechanisms influence RO, PK and total cellular receptor profiles. Results The simulations demonstrate the effect antibody-receptor internalization, affinity and receptor turnover have on commonly acquired data in drug development. Conclusions Integrating RO data with mathematical models such as the two presented here (target-mediated drug disposition and site-of-action models) can provide a more comprehensive view of the biological system, which can be used to test hypotheses, extrapolate preclinical findings to humans and impact clinical study designs and risk assessments for the successful development of biotherapeutics. © 2015 International Clinical Cytometry Society

A number of therapeutics have been developed or are under development aiming to modulate the coagulation network to treat various diseases. We used a systems model to better understand the effect of modulating various components on blood coagulation. A computational model of the coagulation network was built to match in-house in vitro thrombin generation and activated Partial Thromboplastin Time (aPTT) data with various concentrations of recombinant factor VIIa (FVIIa) or factor Xa added to normal human plasma or factor VIII-deficient plasma. Sensitivity analysis applied to the model revealed that lag time, peak thrombin concentration, area under the curve (AUC) of the thrombin generation profile, and aPTT show different sensitivity to changes in coagulation factors' concentrations and type of plasma used (normal or factor VIII-deficient). We also used the model to explore how variability in concentrations of the proteins in coagulation network can impact the response to FVIIa treatment.

The objectives of this analysis were to develop a population pharmacokinetics (PK) model of durvalumab, an anti-PD-L1 antibody, and quantify the impact of baseline and time-varying patient/disease characteristics on PK. Pooled data from two studies (1,409 patients providing 7,407 PK samples) were analyzed with nonlinear mixed effects modeling. Durvalumab PK was best described by a two-compartment model with both linear and nonlinear clearances. Three candidate models were evaluated: a time-invariant clearance (CL) model, an empirical time-varying CL model, and a semimechanistic time-varying CL model incorporating longitudinal covariates related to disease status (tumor shrinkage and albumin). The data supported a slight decrease in durvalumab clearance with time and suggested that it may be associated with a decrease in nonspecific protein catabolic rate among cancer patients who benefit from therapy. No covariates were clinically relevant, indicating no need for dose adjustment. Simulations indicated similar overall PK exposures following weight-based and flat-dosing regimens.

When the immune system responds to tumour development, patterns of immune infiltrates emerge, highlighted by the expression of immune checkpoint-related molecules such as PDL1 on the surface of cancer cells. Such spatial heterogeneity carries information on intrinsic characteristics of the tumour lesion for individual patients, and thus is a potential source for biomarkers for anti-tumour therapeutics. We developed a systems biology multiscale agent-based model to capture the interactions between immune cells and cancer cells, and analysed the emergent global behaviour during tumour development and immunotherapy. Using this model, we are able to reproduce temporal dynamics of cytotoxic T cells and cancer cells during tumour progression, as well as three-dimensional spatial distributions of these cells. By varying the characteristics of the neoantigen profile of individual patients, such as mutational burden and antigen strength, a spectrum of pretreatment spatial patterns of PDL1 expression is generated in our simulations, resembling immuno-architectures obtained via immunohistochemistry from patient biopsies. By correlating these spatial characteristics with in silico treatment results using immune checkpoint inhibitors, the model provides a framework for use to predict treatment/biomarker combinations in different cancer types based on cancer-specific experimental data.

Mutant epidermal growth factor receptor (EGFR) is a major driver of non-small-cell lung cancer (NSCLC). Marketed first generation inhibitors, such as erlotinib, effect a transient beneficial response in EGFR mutant NSCLC patients before resistance mechanisms render these inhibitors ineffective. Secondary oncogenic EGFR mutations account for approximately 50% of relapses, the most common being the gatekeeper T790M substitution that renders existing therapies ineffective. The discovery of PF-06459988 (1), an irreversible pyrrolopyrimidine inhibitor of EGFR T790M mutants, was recently disclosed.1 Herein, we describe our continued efforts to achieve potency across EGFR oncogenic mutations and improved kinome selectivity, resulting in the discovery of clinical candidate PF-06747775 (21), which provides potent EGFR activity against the four common mutants (exon 19 deletion (Del), L858R, and double mutants T790M/L858R and T790M/Del), selectivity over wild-type EGFR, and desirable ADME properties. Compound 21 is currently being evaluated in phase-I clinical trials of mutant EGFR driven NSCLC.

Interleukin (IL)-13 is involved in the pathogenesis of some types of asthma. Tralokinumab is a human immunoglobulin G4 monoclonal antibody that specifically binds to IL-13. Two placebo-controlled phase II studies (phase IIa, NCT00873860 and phase IIb, NCT01402986) have been conducted in which tralokinumab was administered subcutaneously. This investigation aimed to characterize tralokinumab's dose-exposure-response (forced expiratory volume in 1 s (FEV1 )) relationship in patients with asthma and to predict the most appropriate dose for phase III. An integrated population pharmacokinetic-pharmacodynamic (PK/PD) modeling analysis was required for phase III dose selection, due to differing phase II patient populations, designs, and regimens. Analysis of combined datasets enabled the identification of tralokinumab's dose-exposure-FEV1 response relationship in patients with asthma. Near-maximal FEV1 increase was predicted at a dose of 300 mg SC once every 2 weeks (Q2W). This dose was chosen for tralokinumab in the phase III clinical development program for treatment of severe, uncontrolled asthma.

Item does not contain fulltext Patient groups prone to polypharmacy and special subpopulations are susceptible to suboptimal treatment. Refined dosing in special populations is imperative to improve therapeutic response and/or lowering the risk of toxicity. Model-informed precision dosing (MIPD) may improve treatment outcomes by achieving the optimal dose for an individual patient. There is, however, relatively little published evidence of large-scale utility and impact of MIPD, where it is often implemented as local collaborative efforts between academia and healthcare. This article highlights some successful applications of bringing MIPD to clinical care and proposes strategies for wider integration in healthcare. Considerations are brought up herein that will need addressing to see MIPD become "widespread clinical practice," among those, wider interdisciplinary collaborations and the necessity for further evidence-based efficacy and cost-benefit analysis of MIPD in healthcare. The implications of MIPD on regulatory policies and pharmaceutical development are also discussed as part of the roadmap.

Quantitative relationships between 9-β-d-arabinofuranosyl-2-fluoroadenine (F-ara-A) concentrations and lymphosuppression have not been reported, but would be useful for regimen design. A population pharmacokinetic/pharmacodynamic model was constructed in this study using data from 41 hematopoietic cell transplant (HCT) recipients conditioned with busulfan in combination with fludarabine (total dose 120 mg/m2, Protocol 1519) or with fludarabine (total dose 250 mg/m2) with rabbit antithymocyte globulin (rATG, Protocol 2041). Individual pharmacokinetic parameters were fixed to post hoc Bayesian estimates, and circulating absolute lymphocyte counts (ALC) were obtained during the 3 weeks prior to graft infusion. A semi-physiological cell-kill model with three lymphocyte transit compartments was applied and aptly characterized the time course of suppression of circulating ALC by fludarabine administration. Drug- and system-specific parameters were estimated using a maximum likelihood expectation maximization algorithm, and the final model was qualified using an internal visual predictive check. The final model successfully characterized the time course and variability in ALC. Pharmacodynamic parameters exhibited considerable between subject variability (38.9–211 %). The HCT protocol was the only covariate associated with the pharmacodynamic parameters, specifically the lymphocyte kill rate, the transit rate between lymphocyte compartments, and the baseline ALC. This model can be used to simulate the degree of lymphosuppression for design of future fludarabine-based conditioning regimens.

A key role of B cells in the mammalian immune response is the generation of antibodies that serve to protect the organism against antigenic challenges. The same process may also be detrimental in the context of autoimmunity. Several modeling approaches have been applied to this aspect of the immune response, from predicting potential epitopes to describing B cells progress through developmental models and simulating antibody production. Here we review some of the modeling techniques, and summarize models that describe different activation mechanisms for B cells, including T cell dependent and independent models. We focus on viral infection as a prototype system, and briefly describe case studies in other disease areas such as bacterial infection and oncology. We single out aspects of the B cell response for which there are current knowledge gaps. We outline areas in need of further research in modeling applications to ultimately produce a "B cell module" for a complete immune response model.