An integrated quantitative systems pharmacology virtual population approach for calibration with oncology efficacy endpoints.

In drug development, quantitative systems pharmacology (QSP) models are becoming an increasingly important mathematical tool for understanding response variability and for generating predictions to inform development decisions. Virtual populations are essential for sampling uncertainty and potential variability in QSP model predictions, but many clinical efficacy endpoints can be difficult to capture with QSP models that typically rely on mechanistic biomarkers. In oncology, challenges are particularly significant when connecting tumor size with time-to-event endpoints like progression-free survival while also accounting for censoring due to consent withdrawal, loss in follow-up, or safety criteria.

A multiscale model of the action of a capsid assembly modulator for the treatment of chronic hepatitis B.

Chronic hepatitis B virus (HBV) infection is strongly associated with increased risk of liver cancer and cirrhosis. While existing treatments effectively inhibit the HBV life cycle, viral rebound occurs rapidly following treatment interruption. Consequently, functional cure rates of chronic HBV infection remain low and there is increased interest in a novel treatment modality, capsid assembly modulators (CAMs).

Modeling resistance to the broadly neutralizing antibody PGT121 in people living with HIV-1.

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization.

A Continuation Technique for Maximum Likelihood Estimators in Biological Models.

Estimating model parameters is a crucial step in mathematical modelling and typically involves minimizing the disagreement between model predictions and experimental data. This calibration data can change throughout a study, particularly if modelling is performed simultaneously with the calibration experiments, or during an on-going public health crisis as in the case of the COVID-19 pandemic. Consequently, the optimal parameter set, or maximal likelihood estimator (MLE), is a function of the experimental data set.

Numerical methods and hypoexponential approximations for gamma distributed delay differential equations.

Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic gamma distributed DDEs have not previously been implemented. Modellers have therefore resorted to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations (ODEs).

A simple model of COVID-19 explains disease severity and the effect of treatments.

Considerable effort has been made to better understand why some people suffer from severe COVID-19 while others remain asymptomatic. This has led to important clinical findings; people with severe COVID-19 generally experience persistently high levels of inflammation, slower viral load decay, display a dysregulated type-I interferon response, have less active natural killer cells and increased levels of neutrophil extracellular traps. How these findings are connected to the pathogenesis of COVID-19 remains unclear.

Safety, pharmacokinetics and antiviral activity of PGT121, a broadly neutralizing monoclonal antibody against HIV-1: a randomized, placebo-controlled, phase 1 clinical trial.

Human immunodeficiency virus (HIV)-1-specific broadly neutralizing monoclonal antibodies are currently under development to treat and prevent HIV-1 infection. We performed a single-center, randomized, double-blind, dose-escalation, placebo-controlled trial of a single administration of the HIV-1 V3-glycan-specific antibody PGT121 at 3, 10 and 30 mg kg in HIV-uninfected adults and HIV-infected adults on antiretroviral therapy (ART), as well as a multicenter, open-label trial of one infusion of PGT121 at 30 mg kg in viremic HIV-infected adults not on ART (no. NCT02960581).

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance.

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy.

Immunological Correlates of the HIV-1 Replication-Competent Reservoir Size.

Understanding what shapes the latent human immunodeficiency virus type 1 (HIV-1) reservoir is critical for developing strategies for cure. We measured frequency of persistent HIV-1 infection after 5 years of suppressive antiretroviral therapy initiated during chronic infection. Pretreatment CD8+ T-cell activation, nadir CD4 count, and CD4:CD8 ratio predicted reservoir size.

In silico trials predict that combination strategies for enhancing vesicular stomatitis oncolytic virus are determined by tumor aggressivity.

Background: Immunotherapies, driven by immune-mediated antitumorigenicity, offer the potential for significant improvements to the treatment of multiple cancer types. Identifying therapeutic strategies that bolster antitumor immunity while limiting immune suppression is critical to selecting treatment combinations and schedules that offer durable therapeutic benefits. Combination oncolytic virus (OV) therapy, wherein complementary OVs are administered in succession, offer such promise, yet their translation from preclinical studies to clinical implementation is a major challenge.

Characterizing Chemotherapy-Induced Neutropenia and Monocytopenia Through Mathematical Modelling.

In spite of the recent focus on the development of novel targeted drugs to treat cancer, cytotoxic chemotherapy remains the standard treatment for the vast majority of patients. Unfortunately, chemotherapy is associated with high hematopoietic toxicity that may limit its efficacy. We have previously established potential strategies to mitigate chemotherapy-induced neutropenia (a lack of circulating neutrophils) using a mechanistic model of granulopoiesis to predict the interactions defining the neutrophil response to chemotherapy and to define optimal strategies for concurrent chemotherapy/prophylactic granulocyte colony-stimulating factor (G-CSF).

Determinants of combination GM-CSF immunotherapy and oncolytic virotherapy success identified through in silico treatment personalization.

Oncolytic virotherapies, including the modified herpes simplex virus talimogene laherparepvec (T-VEC), have shown great promise as potent instigators of anti-tumour immune effects. The OPTiM trial, in particular, demonstrated the superior anti-cancer effects of T-VEC as compared to systemic immunotherapy treatment using exogenous administration of granulocyte-macrophage colony-stimulating factor (GM-CSF). Theoretically, a combined approach leveraging exogenous cytokine immunotherapy and oncolytic virotherapy would elicit an even greater immune response and improve patient outcomes.

Equivalences between age structured models and state dependent distributed delay differential equations.

We use the McKendrick equation with variable ageing rate and randomly distributed mat-uration time to derive a state dependent distributed delay differential equation. We show that the resulting delay differential equation preserves non-negativity of initial conditions and we characterise local stability of equilibria. By specifying the distribution of maturation age, we recover state depen-dent discrete, uniform and gamma distributed delay differential equations.

A mathematical model of viral oncology as an immuno-oncology instigator.

We develop and analyse a mathematical model of tumour-immune interaction that explicitly incorporates heterogeneity in tumour cell cycle duration by using a distributed delay differential equation. We derive a necessary and sufficient condition for local stability of the cancer-free equilibrium in which the amount of tumour-immune interaction completely characterizes disease progression. Consistent with the immunoediting hypothesis, we show that decreasing tumour-immune interaction leads to tumour expansion.

Transit and lifespan in neutrophil production: implications for drug intervention.

A comparison of the transit compartment ordinary differential equation modelling approach to distributed and discrete delay differential equation models is studied by focusing on Quartino's extension to the Friberg transit compartment model of myelosuppression, widely relied upon in the pharmaceutical sciences to predict the neutrophil response after chemotherapy, and on a QSP delay differential equation model of granulopoiesis. An extension to the Quartino model is provided by considering a general number of transit compartments and introducing an extra parameter that allows for the decoupling of the maturation time from the production rate of cells. An overview of the well established linear chain technique, used to reformulate transit compartment models with constant transit rates as distributed delay differential equations (DDEs), is then given.