The efficacy of antiretroviral therapy is significantly compromised by medication non-adherence. Long-acting enteral systems that can ease the burden of daily adherence have not yet been developed. Here we describe an oral dosage form composed of distinct drug–polymer matrices which achieved week-long systemic drug levels of the antiretrovirals dolutegravir, rilpivirine and cabotegravir in a pig. Simulations of viral dynamics and patient adherence patterns indicate that such systems would significantly reduce therapeutic failures and epidemiological modelling suggests that using such an intervention prophylactically could avert hundreds of thousands of new HIV cases. In sum, weekly administration of long-acting antiretrovirals via a novel oral dosage form is a promising intervention to help control the HIV epidemic worldwide.

No abstract

The 2019 novel coronavirus, SARS-CoV-2, is an emerging pathogen of critical significance to international public health. Knowledge of the interplay between molecular-scale virus-receptor interactions, single-cell viral replication, intracellular-scale viral transport, and emergent tissue-scale viral propagation is limited. Moreover, little is known about immune system-virus-tissue interactions and how these can result in low-level (asymptomatic) infections in some cases and acute respiratory distress syndrome (ARDS) in others, particularly with respect to presentation in different age groups or pre-existing inflammatory risk factors like diabetes. A critical question for treatment and protection is why it appears that the severity of infection may correlate with the initial level of virus exposure. Given the nonlinear interactions within and among each of these processes, multiscale simulation models can shed light on the emergent dynamics that lead to divergent outcomes, identify actionable "choke points" for pharmacologic interactions, screen potential therapies, and identify potential biomarkers that differentiate response dynamics. Given the complexity of the problem and the acute need for an actionable model to guide therapy discovery and optimization, we introduce a prototype of a multiscale model of SARS-CoV-2 dynamics in lung and intestinal tissue that will be iteratively refined. The first prototype model was built and shared internationally as open source code and interactive, cloud-hosted executables in under 12 hours. In a sustained community effort, this model will integrate data and expertise across virology, immunology, mathematical biology, quantitative systems physiology, cloud and high performance computing, and other domains to accelerate our response to this critical threat to international health.

We develop a physiological model of granulopoiesis which includes explicit modelling of the kinetics of the cytokine granulocyte colony-stimulating factor (G-CSF) incorporating both the freely circulating concentration and the concentration of the cytokine bound to mature neutrophils. G-CSF concentrations are used to directly regulate neutrophil production, with the rate of differentiation of stem cells to neutrophil precursors, the effective proliferation rate in mitosis, the maturation time, and the release rate from the mature marrow reservoir into circulation all dependent on the level of G-CSF in the system. The dependence of the maturation time on the cytokine concentration introduces a state-dependent delay into our differential equation model, and we show how this is derived from an age-structured partial differential equation model of the mitosis and maturation, and also detail the derivation of the rest of our model. The model and its estimated parameters are shown to successfully predict the neutrophil and G-CSF responses to a variety of treatment scenarios, including the combined administration of chemotherapy and exogenous G-CSF. This concomitant treatment was reproduced without any additional fitting to characterise drug-drug interactions.

No abstract

To understand the diversity of immune responses to SARS-CoV-2 and distinguish features that predispose individuals to severe COVID-19, we developed a mechanistic, within-host mathematical model and virtual patient cohort. Our results suggest that virtual patients with low production rates of infected cell derived IFN subsequently experienced highly inflammatory disease phenotypes, compared to those with early and robust IFN responses. In these in silico patients, the maximum concentration of IL-6 was also a major predictor of CD8+ T cell depletion. Our analyses predicted that individuals with severe COVID-19 also have accelerated monocyte-to-macrophage differentiation mediated by increased IL-6 and reduced type I IFN signalling. Together, these findings suggest biomarkers driving the development of severe COVID-19 and support early interventions aimed at reducing inflammation.

We develop a mathematical model of platelet, megakaryocyte, and thrombopoietin dynamics in humans. We show that there is a single stationary solution that can undergo a Hopf bifurcation, and use this information to investigate both normal and pathological platelet production, specifically cyclic thrombocytopenia. Carefully estimating model parameters from laboratory and clinical data, we then argue that a subset of parameters are involved in the genesis of cyclic thrombocytopenia based on clinical information. We provide model fits to the existing data for both platelet counts and thrombopoietin levels by changing four parameters that have physiological correlates. Our results indicate that the primary change in cyclic thrombocytopenia is an interference with, or destruction of, the thrombopoietin receptor with secondary changes in other processes, including immune-mediated destruction of platelets and megakaryocyte deficiency and failure in platelet production. This study contributes to the understanding of the origin of cyclic thrombocytopenia as well as extending the modeling of thrombopoiesis.

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. However, regimen scheduling of combination immunostimulation and T-VEC therapy has yet to be established. Here, we calibrate a computational biology model of sensitive and resistant tumour cells and immune interactions for implementation into an in silico clinical trial to test and individualize combination immuno- and virotherapy. By personalizing and optimizing combination oncolytic virotherapy and immunostimulatory therapy, we show improved simulated patient outcomes for individuals with late-stage melanoma. More crucially, through evaluation of individualized regimens, we identified determinants of combination GM-CSF and T-VEC therapy that can be translated into clinically-actionable dosing strategies without further personalization. Our results serve as a proof-of-concept for interdisciplinary approaches to determining combination therapy, and suggest promising avenues of investigation towards tailored combination immunotherapy/oncolytic virotherapy.