NIMROD: a program for inference via a normal approximation of the posterior in models with random effects based on ordinary differential equations

Pre-existing resistance in the latent reservoir can compromise VRC01 therapy during chronic HIV-1 infection

Passive immunization with broadly neutralizing antibodies (bNAbs) of HIV-1 appears a promising strategy for eliciting long-term HIV-1 remission. When administered concomitantly with the cessation of antiretroviral therapy (ART) to patients with established viremic control, bNAb therapy is expected to prolong remission. Surprisingly, in clinical trials on chronic HIV-1 patients, the bNAb VRC01 failed to prolong remission substantially. Identifying the cause of this failure is important for improving VRC01-based therapies and unraveling potential vulnerabilities of other bNAbs. In the trials, viremia resurged rapidly in most patients despite suppressive VRC01 concentrations in circulation, suggesting that VRC01 resistance was the likely cause of failure. ART swiftly halts viral replication, precluding the development of resistance during ART. If resistance were to emerge post ART, virological breakthrough would have taken longer than without VRC01 therapy. We hypothesized therefore that VRC01-resistant strains must have been formed before ART initiation, survived ART in latently infected cells, and been activated during VRC01 therapy, causing treatment failure. Current assays preclude testing this hypothesis experimentally. We developed a mathematical model based on the hypothesis and challenged it with available clinical data. The model integrated within-host HIV-1 evolution, stochastic latency reactivation, and viral dynamics with multiple-dose VRC01 pharmacokinetics. The model predicted that single but not higher VRC01-resistant mutants would pre-exist in the latent reservoir. We constructed a virtual patient population that parsimoniously recapitulated inter-patient variations. Model predictions with this population quantitatively captured data of VRC01 failure from clinical trials, presenting strong evidence supporting the hypothesis. We attributed VRC01 failure to single-mutant VRC01-resistant proviruses in the latent reservoir triggering viral recrudescence, particularly when VRC01 was at trough levels. Pre-existing resistant proviruses in the latent reservoir may similarly compromise other bNAbs. Our study provides a framework for designing bNAb-based therapeutic protocols that would avert such failure and maximize HIV-1 remission.

Antiretroviral therapy (ART) can control but not eradicate HIV-1. Stopping ART leads to rapid viral resurgence and progressive disease. ART is therefore administered lifelong. Tremendous efforts are ongoing to devise strategies that will enable stopping ART and yet prevent viral resurgence. One such strategy involves the administration of broadly neutralizing antibodies (bNAbs) of HIV-1 at the time of stopping ART. This strategy is expected to delay if not prevent viral resurgence. Surprisingly, treatment with VRC01, a potent bNAb, resulted in hardly any improvement in viral remission. In this study, we elucidate the cause of this failure. We hypothesized that VRC01-resistant strains may pre-exist in latently infected cells, which are unaffected by ART. They can thus outlast ART and get reactivated, triggering VRC01 failure. We built a detailed mathematical model based on this hypothesis and showed that it quantitatively captured observations of VRC01 failure in clinical trials on chronic HIV-1 patients. Our study thus identifies a potential vulnerability of bNAbs, namely, bNAb-resistant strains pre-existing in latently infected cells. Our model offers a framework for predicting bNAb-based treatment protocols that would preclude failure due to pre-existing resistance and maximally prolong remission.

Comparison of empirical and dynamic models for HIV viral load rebound after treatment interruption

OBJECTIVE

To compare empirical and mechanistic modeling approaches for describing HIV-1 RNA viral load trajectories after antiretroviral treatment interruption and for identifying factors that predict features of viral rebound process.

METHODS

We apply and compare two modeling approaches in analysis of data from 346 participants in six AIDS Clinical Trial Group studies. From each separate analysis, we identify predictors for viral set points and delay in rebound. Our empirical model postulates a parametric functional form whose parameters represent different features of the viral rebound process, such as rate of rise and viral load set point. The viral dynamics model augments standard HIV dynamics models-a class of mathematical models based on differential equations describing biological mechanisms-by including reactivation of latently infected cells and adaptive immune response. We use Monolix, which makes use of a Stochastic Approximation of the Expectation-Maximization algorithm, to fit non-linear mixed effects models incorporating observations that were below the assay limit of quantification.

RESULTS

Among the 346 participants, the median age at treatment interruption was 42. Ninety-three percent of participants were male and sixty-five percent, white non-Hispanic. Both models provided a reasonable fit to the data and can accommodate atypical viral load trajectories. The median set points obtained from two approaches were similar: 4.44 log10 copies/mL from the empirical model and 4.59 log10 copies/mL from the viral dynamics model. Both models revealed that higher nadir CD4 cell counts and ART initiation during acute/recent phase were associated with lower viral set points and identified receiving a non-nucleoside reverse transcriptase inhibitor (NNRTI)-based pre-ATI regimen as a predictor for a delay in rebound.

CONCLUSION

Although based on different sets of assumptions, both models lead to similar conclusions regarding features of viral rebound process.

A model for establishment, maintenance and reactivation of the immune response after vaccination against Ebola virus

The 2014-2016 Ebola outbreak in West Africa has triggered accelerated development of several preventive vaccines against Ebola virus. Under the EBOVAC1 consortium, three phase I studies were carried out to assess safety and immunogenicity of a two-dose heterologous vaccination regimen developed by Janssen Vaccines and Prevention in collaboration with Bavarian Nordic. To describe the immune response induced by the two-dose heterologous vaccine regimen, we propose a mechanistic ODE based model, which takes into account the role of immunological memory. We perform identifiability and sensitivity analysis of the proposed model to establish which kind of biological data are ideally needed in order to accurately estimate parameters, and additionally, which of those are non-identifiable based on the available data. Antibody concentrations data from phase I studies have been used to calibrate the model and show its ability in reproducing the observed antibody dynamics. Together with other factors, the establishment of an effective and reactive immunological memory is of pivotal importance for several prophylactic vaccines. We show that introducing a memory compartment in our calibrated model allows to evaluate the magnitude of the immune response induced by a booster dose and its long-term persistence afterwards.

Dynamics of the Humoral Immune Response to a Prime-Boost Ebola Vaccine: Quantification and Sources of Variation

The Ebola vaccine based on Ad26.ZEBOV/MVA-BN-Filo prime-boost regimens is being evaluated in multiple clinical trials. The long-term immune response to the vaccine is unknown, including factors associated with the response and variability around the response. We analyzed data from three phase 1 trials performed by the EBOVAC1 Consortium in four countries: the United Kingdom, Kenya, Tanzania, and Uganda. Participants were randomized into four groups based on the interval between prime and boost immunizations (28 or 56 days) and the sequence in which Ad26.ZEBOV and MVA-BN-Filo were administered. Consecutive enzyme-linked immunosorbent assay (ELISA) measurements of the IgG binding antibody concentrations against the Kikwit glycoprotein (GP) were available for 177 participants to assess the humoral immune response up to 1 year postprime. Using a mathematical model for the dynamics of the humoral response, from 7 days after the boost immunization up to 1 year after the prime immunization, we estimated the durability of the antibody response and the influence of different factors on the dynamics of the humoral response. Ordinary differential equations (ODEs) described the dynamics of antibody response and two populations of antibody-secreting cells (ASCs), short-lived (SL) and long-lived (LL). Parameters of the ODEs were estimated using a population approach. We estimated that half of the LL ASCs could persist for at least 5 years. The vaccine regimen significantly affected the SL ASCs and the antibody peak but not the long-term response. The LL ASC compartment dynamics differed significantly by geographic regions analyzed, with a higher long-term antibody persistence in European subjects. These differences could not be explained by the observed differences in cellular immune response.IMPORTANCE With no available licensed vaccines or therapies, the West African Ebola virus disease epidemic of 2014 to 2016 caused 11,310 deaths. Following this outbreak, the development of vaccines has been accelerated. Combining different vector-based vaccines as heterologous regimens could induce a durable immune response, assessed through antibody concentrations. Based on data from phase 1 trials in East Africa and Europe, the dynamics of the humoral immune response from 7 days after the boost immunization onwards were modeled to estimate the durability of the response and understand its variability. Antibody production is maintained by a population of long-lived cells. Estimation suggests that half of these cells can persist for at least 5 years in humans. Differences in prime-boost vaccine regimens affect only the short-term immune response. Geographical differences in long-lived cell dynamics were inferred, with higher long-term antibody concentrations induced in European participants.

Adaptive protocols based on predictions from a mechanistic model of the effect of IL7 on CD4 counts

In human immunodeficiency virus-infected patients, antiretroviral therapy suppresses the viral replication, which is followed in most patients by a restoration of CD4+ T cells pool. For patients who fail to do so, repeated injections of exogenous interleukin 7 (IL7) are experimented. The IL7 is a cytokine that is involved in the T cell homeostasis and the INSPIRE study has shown that injections of IL7 induced a proliferation of CD4+ T cells. Phase I/II INSPIRE 2 and 3 studies have evaluated a protocol in which a first cycle of three IL7 injections is followed by a new cycle at each visit when the patient has less than 550 CD4 cells/μL. Restoration of the CD4 concentration has been demonstrated, but the long-term best adaptive protocol is yet to be determined. A mechanistic model of the evolution of CD4 after IL7 injections has been developed, which is based on a system of ordinary differential equations and includes random effects. Based on the estimation of this model, we use a Bayesian approach to forecast the dynamics of CD4 in new patients. We propose four prediction-based adaptive protocols of injections to minimize the time spent under 500 CD4 cells/μL for each patient, without increasing the number of injections received too much. We show that our protocols significantly reduce the time spent under 500 CD4 over a period of two years, without increasing the number of injections. These protocols have the potential to increase the efficiency of this therapy.

Controlling IL-7 Injections in HIV-Infected Patients

Immune interventions consisting in repeated injections are broadly used as they are thought to improve the quantity and the quality of the immune response. However, they also raise several questions that remain unanswered, in particular the number of injections to make or the delay to respect between different injections to achieve this goal. Practical and financial considerations add constraints to these questions, especially in the framework of human studies. We specifically focus here on the use of interleukin-7 (IL-7) injections in HIV-infected patients under antiretroviral treatment, but still unable to restore normal levels of [Formula: see text] T lymphocytes. Clinical trials have already shown that repeated cycles of injections of IL-7 could help maintaining [Formula: see text] T lymphocytes levels over the limit of 500 cells/[Formula: see text]L, by affecting proliferation and survival of [Formula: see text] T cells. We then aim at answering the question: how to maintain a patients level of [Formula: see text] T lymphocytes by using a minimum number of injections (i.e., optimizing the strategy of injections)? Based on mechanistic models that were previously developed for the dynamics of [Formula: see text] T lymphocytes in this context, we model the process by a piecewise deterministic Markov model. We then address the question by using some recently established theory on impulse control problem in order to develop a numerical tool determining the optimal strategy. Results are obtained on a reduced model, as a proof of concept: the method allows to define an optimal strategy for a given patient. This method could be applied to optimize injections schedules in clinical trials.

Dynamic models for estimating the effect of HAART on CD4 in observational studies: Application to the Aquitaine Cohort and the Swiss HIV Cohort Study

Highly active antiretroviral therapy (HAART) has proved efficient in increasing CD4 counts in many randomized clinical trials. Because randomized trials have some limitations (e.g., short duration, highly selected subjects), it is interesting to assess the effect of treatments using observational studies. This is challenging because treatment is started preferentially in subjects with severe conditions. This general problem had been treated using Marginal Structural Models (MSM) relying on the counterfactual formulation. Another approach to causality is based on dynamical models. We present three discrete-time dynamic models based on linear increments models (LIM): the first one based on one difference equation for CD4 counts, the second with an equilibrium point, and the third based on a system of two difference equations, which allows jointly modeling CD4 counts and viral load. We also consider continuous-time models based on ordinary differential equations with non-linear mixed effects (ODE-NLME). These mechanistic models allow incorporating biological knowledge when available, which leads to increased statistical evidence for detecting treatment effect. Because inference in ODE-NLME is numerically challenging and requires specific methods and softwares, LIM are a valuable intermediary option in terms of consistency, precision, and complexity. We compare the different approaches in simulation and in illustration on the ANRS CO3 Aquitaine Cohort and the Swiss HIV Cohort Study.

Using the SAEM algorithm for mechanistic joint models characterizing the relationship between nonlinear PSA kinetics and survival in prostate cancer patients

Joint modeling is increasingly popular for investigating the relationship between longitudinal and time-to-event data. However, numerical complexity often restricts this approach to linear models for the longitudinal part. Here, we use a novel development of the Stochastic-Approximation Expectation Maximization algorithm that allows joint models defined by nonlinear mixed-effect models. In the context of chemotherapy in metastatic prostate cancer, we show that a variety of patterns for the Prostate Specific Antigen (PSA) kinetics can be captured by using a mechanistic model defined by nonlinear ordinary differential equations. The use of a mechanistic model predicts that biological quantities that cannot be observed, such as treatment-sensitive and treatment-resistant cells, may have a larger impact than PSA value on survival. This suggests that mechanistic joint models could constitute a relevant approach to evaluate the efficacy of treatment and to improve the prediction of survival in patients.

The stochastic system approach for estimating dynamic treatments effect

The problem of assessing the effect of a treatment on a marker in observational studies raises the difficulty that attribution of the treatment may depend on the observed marker values. As an example, we focus on the analysis of the effect of a HAART on CD4 counts, where attribution of the treatment may depend on the observed marker values. This problem has been treated using marginal structural models relying on the counterfactual/potential response formalism. Another approach to causality is based on dynamical models, and causal influence has been formalized in the framework of the Doob-Meyer decomposition of stochastic processes. Causal inference however needs assumptions that we detail in this paper and we call this approach to causality the "stochastic system" approach. First we treat this problem in discrete time, then in continuous time. This approach allows incorporating biological knowledge naturally. When working in continuous time, the mechanistic approach involves distinguishing the model for the system and the model for the observations. Indeed, biological systems live in continuous time, and mechanisms can be expressed in the form of a system of differential equations, while observations are taken at discrete times. Inference in mechanistic models is challenging, particularly from a numerical point of view, but these models can yield much richer and reliable results.

Quantifying and predicting the effect of exogenous interleukin-7 on CD4+ T cells in HIV-1 infection

Exogenous Interleukin-7 (IL-7), in supplement to antiretroviral therapy, leads to a substantial increase of all CD4⁺ T cell subsets in HIV-1 infected patients. However, the quantitative contribution of the several potential mechanisms of action of IL-7 is unknown. We have performed a mathematical analysis of repeated measurements of total and naive CD4⁺ T cells and their Ki67 expression from HIV-1 infected patients involved in three phase I/II studies (N = 53 patients). We show that, besides a transient increase of peripheral proliferation, IL-7 exerts additional effects that play a significant role in CD4⁺ T cell dynamics up to 52 weeks. A decrease of the loss rate of the total CD4⁺ T cell is the most probable explanation. If this effect could be maintained during repeated administration of IL-7, our simulation study shows that such a strategy may allow maintaining CD4⁺ T cell counts above 500 cells/µL with 4 cycles or fewer over a period of two years. This in-depth analysis of clinical data revealed the potential for IL-7 to achieve sustained CD4⁺ T cell restoration with limited IL-7 exposure in HIV-1 infected patients with immune failure despite antiretroviral therapy.

HIV infection is characterized by a decrease of CD4⁺ T-lymphocytes in the blood. Whereas antiretroviral treatment succeeds to control viral replication, some patients fail to reconstitute their CD4⁺ T cell count to normal value. IL-7 is a promising cytokine under evaluation for its use in HIV infection, in supplement to antiretroviral therapy, as it increases cell proliferation and survival. Here, we use data from three clinical trials testing the effect of IL-7 on CD4⁺ T-cell recovery in treated HIV-infected individuals and use a simple mathematical model to quantify IL-7 effects by estimating the biological parameters of the model. We show that the increase of peripheral proliferation could not explain alone the long-term dynamics of T cells after IL-7 injections underlining other important effects such as the improvement of cell survival. We also investigate the feasibility and the efficiency of repetitions of IL-7 cycles and argue for further evaluation through clinical trials.