Translational modeling-based evidence for enhanced efficacy of standard-of-care drugs in combination with anti-microRNA-155 in non-small-cell lung cancer.

Elevated microRNA-155 (miR-155) expression in non-small-cell lung cancer (NSCLC) promotes cisplatin resistance and negatively impacts treatment outcomes. However, miR-155 can also boost anti-tumor immunity by suppressing PD-L1 expression. We developed a multiscale mechanistic model, calibrated with in vivo data and then extrapolated to humans, to investigate the therapeutic effects of nanoparticle-delivered anti-miR-155 in NSCLC, alone or in combination with standard-of-care drugs. Model simulations and analyses of the clinical scenario revealed that monotherapy with anti-miR-155 at a dose of 2.5 mg/kg administered once every three weeks has substantial anti-cancer activity. It led to a median progression-free survival (PFS) of 6.7 months, which compared favorably to cisplatin and immune checkpoint inhibitors. Further, we explored the combinations of anti-miR-155 with standard-of-care drugs, and found strongly synergistic two- and three-drug combinations. A three-drug combination of anti-miR-155, cisplatin, and pembrolizumab resulted in a median PFS of 13.1 months, while a two-drug combination of anti-miR-155 and cisplatin resulted in a median PFS of 11.3 months, which emerged as a more practical option due to its simple design and cost-effectiveness. Our analyses also provided valuable insights into unfavorable dose ratios for drug combinations, highlighting the need for optimizing dose regimen to prevent antagonistic effects. Thus, this work bridges the gap between preclinical development and clinical translation of anti-miR-155 and unravels the potential of anti-miR-155 combination therapies in NSCLC.

Mechanistic modeling for optimal design of dissolvable microneedle-based patches for transdermal drug delivery.

Polymeric microneedle (MN)-based patches are an efficient, non-invasive, and painless means of drug delivery through the skin to systemic circulation. The design of these MN-based patches can be customized for various drug delivery applications, particularly modified release of drugs. In this study, we developed a mathematical model of drug delivery via MN-based patches to study the effect of patch design properties on drug delivery kinetics and systemic pharmacokinetics (PK). We calibrated the model against two representative formulations: a rapid release patch of naloxone and a sustained-release patch of levonorgestrel. The model was then applied to assess the relative significance of model parameters in governing systemic PK of drugs and obtain optimal design parameters to achieve therapeutically meaningful drug levels in a clinical setting. We identified the importance of drug loading fraction, MN base radius, and MN height as the key control parameters responsible for drug PK.Clinical Relevance- Through the application of modeling and simulation, we can improve drug delivery from MN-based patches by identifying optimal design parameters to support the clinical translation of these novel drug delivery systems.

A modeling-based approach to optimize COVID-19 vaccine dosing schedules for improved protection.

While the development of different vaccines slowed the dissemination of SARS-CoV-2, the occurrence of breakthrough infections has continued to fuel the COVID-19 pandemic. To secure at least partial protection in the majority of the population through 1 dose of a COVID-19 vaccine, delayed administration of boosters has been implemented in many countries. However, waning immunity and emergence of new variants of SARS-CoV-2 suggest that such measures may induce breakthrough infections due to intermittent lapses in protection. Optimizing vaccine dosing schedules to ensure prolonged continuity in protection could thus help control the pandemic. We developed a mechanistic model of immune response to vaccines as an in silico tool for dosing schedule optimization. The model was calibrated with clinical data sets of acquired immunity to COVID-19 mRNA vaccines in healthy and immunocompromised participants and showed robust validation by accurately predicting neutralizing antibody kinetics in response to multiple doses of COVID-19 mRNA vaccines. Importantly, by estimating population vulnerability to breakthrough infections, we predicted tailored vaccination dosing schedules to minimize breakthrough infections, especially for immunocompromised individuals. We identified that the optimal vaccination schedules vary from CDC-recommended dosing, suggesting that the model is a valuable tool to optimize vaccine efficacy outcomes during future outbreaks.

A modeling-based approach to optimize COVID-19 vaccine dosing schedules for improved protection.

While the development of different vaccines has slowed the dissemination of SARS-CoV-2, the occurrence of breakthrough infections continues to fuel the pandemic. As a strategy to secure at least partial protection, with a single dose of a given COVID-19 vaccine to maximum possible fraction of the population, delayed administration of subsequent doses (or boosters) has been implemented in many countries. However, waning immunity and emergence of new variants of SARS-CoV-2 suggest that such measures may jeopardize the attainment of herd immunity due to intermittent lapses in protection. Optimizing vaccine dosing schedules could thus make the difference between periodic occurrence of breakthrough infections or effective control of the pandemic. To this end, we have developed a mechanistic mathematical model of adaptive immune response to vaccines and demonstrated its applicability to COVID-19 mRNA vaccines as a proof-of-concept for future outbreaks. The model was thoroughly calibrated against multiple clinical datasets involving immune response to SARS-CoV-2 infection and mRNA vaccines in healthy and immunocompromised subjects (cancer patients undergoing therapy); the model showed robust clinical validation by accurately predicting neutralizing antibody kinetics, a correlate of vaccine-induced protection, in response to multiple doses of mRNA vaccines. Importantly, we estimated population vulnerability to breakthrough infections and predicted tailored vaccination dosing schedules to maximize protection and thus minimize breakthrough infections, based on the immune status of a sub-population. We have identified a critical waiting window for cancer patients (or, immunocompromised subjects) to allow recovery of the immune system (particularly CD4+ T-cells) for effective differentiation of B-cells to produce neutralizing antibodies and thus achieve optimal vaccine efficacy against variants of concern, especially between the first and second doses. Also, we have obtained optimized dosing schedules for subsequent doses in healthy and immunocompromised subjects, which vary from the CDC-recommended schedules, to minimize breakthrough infections. The developed modeling tool is based on generalized adaptive immune response to antigens and can thus be leveraged to guide vaccine dosing schedules during future outbreaks.

Dedifferentiation-mediated stem cell niche maintenance in early-stage ductal carcinoma in situ progression: insights from a multiscale modeling study.

We present a multiscale agent-based model of ductal carcinoma in situ (DCIS) to study how key phenotypic and signaling pathways are involved in the early stages of disease progression. The model includes a phenotypic hierarchy, and key endocrine and paracrine signaling pathways, and simulates cancer ductal growth in a 3D lattice-free domain. In particular, by considering stochastic cell dedifferentiation plasticity, the model allows for study of how dedifferentiation to a more stem-like phenotype plays key roles in the maintenance of cancer stem cell populations and disease progression. Through extensive parameter perturbation studies, we have quantified and ranked how DCIS is sensitive to perturbations in several key mechanisms that are instrumental to early disease development. Our studies reveal that long-term maintenance of multipotent stem-like cell niches within the tumor are dependent on cell dedifferentiation plasticity, and that disease progression will become arrested due to dilution of the multipotent stem-like population in the absence of dedifferentiation. We have identified dedifferentiation rates necessary to maintain biologically relevant multipotent cell populations, and also explored quantitative relationships between dedifferentiation rates and disease progression rates, which may potentially help to optimize the efficacy of emerging anti-cancer stem cell therapeutics.

Targeting a cell surface vitamin D receptor on tumor-associated macrophages in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive tumor with limited treatment options and poor prognosis. We applied the in vivo phage display technology to isolate peptides homing to the immunosuppressive cellular microenvironment of TNBC as a strategy for non-malignant target discovery. We identified a cyclic peptide (CSSTRESAC) that specifically binds to a vitamin D receptor, protein disulfide-isomerase A3 (PDIA3) expressed on the cell surface of tumor-associated macrophages (TAM), and targets breast cancer in syngeneic TNBC, non-TNBC xenograft, and transgenic mouse models. Systemic administration of CSSTRESAC to TNBC-bearing mice shifted the cytokine profile toward an antitumor immune response and delayed tumor growth. Moreover, CSSTRESAC enabled ligand-directed theranostic delivery to tumors and a mathematical model confirmed our experimental findings. Finally, in silico analysis showed PDIA3-expressing TAM in TNBC patients. This work uncovers a functional interplay between a cell surface vitamin D receptor in TAM and antitumor immune response that could be therapeutically exploited.

Innate Immunity Plays a Key Role in Controlling Viral Load in COVID-19: Mechanistic Insights from a Whole-Body Infection Dynamics Model.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a pathogen of immense public health concern. Efforts to control the disease have only proven mildly successful, and the disease will likely continue to cause excessive fatalities until effective preventative measures (such as a vaccine) are developed. To develop disease management strategies, a better understanding of SARS-CoV-2 pathogenesis and population susceptibility to infection are needed. To this end, mathematical modeling can provide a robust in silico tool to understand COVID-19 pathophysiology and the in vivo dynamics of SARS-CoV-2. Guided by ACE2-tropism (ACE2 receptor dependency for infection) of the virus and by incorporating cellular-scale viral dynamics and innate and adaptive immune responses, we have developed a multiscale mechanistic model for simulating the time-dependent evolution of viral load distribution in susceptible organs of the body (respiratory tract, gut, liver, spleen, heart, kidneys, and brain). Following parameter quantification with in vivo and clinical data, we used the model to simulate viral load progression in a virtual patient with varying degrees of compromised immune status. Further, we ranked model parameters through sensitivity analysis for their significance in governing clearance of viral load to understand the effects of physiological factors and underlying conditions on viral load dynamics. Antiviral drug therapy, interferon therapy, and their combination were simulated to study the effects on viral load kinetics of SARS-CoV-2. The model revealed the dominant role of innate immunity (specifically interferons and resident macrophages) in controlling viral load, and the importance of timing when initiating therapy after infection.

A mathematical model for the quantification of a patient's sensitivity to checkpoint inhibitors and long-term tumour burden.

A large proportion of patients with cancer are unresponsive to treatment with immune checkpoint blockade and other immunotherapies. Here, we report a mathematical model of the time course of tumour responses to immune checkpoint inhibitors. The model takes into account intrinsic tumour growth rates, the rates of immune activation and of tumour-immune cell interactions, and the efficacy of immune-mediated tumour killing. For 124 patients, four cancer types and two immunotherapy agents, the model reliably described the immune responses and final tumour burden across all different cancers and drug combinations examined. In validation cohorts from four clinical trials of checkpoint inhibitors (with a total of 177 patients), the model accurately stratified the patients according to reduced or increased long-term tumour burden. We also provide model-derived quantitative measures of treatment sensitivity for specific drug-cancer combinations. The model can be used to predict responses to therapy and to quantify specific drug-cancer sensitivities in individual patients.

Intratumoral injection of hydrogel-embedded nanoparticles enhances retention in glioblastoma.

Intratumoral drug delivery is a promising approach for the treatment of glioblastoma multiforme (GBM). However, drug washout remains a major challenge in GBM therapy. Our strategy, aimed at reducing drug clearance and enhancing site-specific residence time, involves the local administration of a multi-component system comprised of nanoparticles (NPs) embedded within a thermosensitive hydrogel (HG). Herein, our objective was to examine the distribution of NPs and their cargo following intratumoral administration of this system in GBM. We hypothesized that the HG matrix, which undergoes rapid gelation upon increases in temperature, would contribute towards heightened site-specific retention and permanence of NPs in tumors. BODIPY-containing, infrared dye-labeled polymeric NPs embedded in a thermosensitive HG (HG-NPs) were fabricated and characterized. Retention and distribution dynamics were subsequently examined over time in orthotopic GBM-bearing mice. Results demonstrate that the HG-NPs system significantly improved site-specific, long-term retention of both NPs and BODIPY, with co-localization analyses showing that HG-NPs covered larger areas of the tumor and the peri-tumor region at later time points. Moreover, NPs released from the HG were shown to undergo uptake by surrounding GBM cells. Findings suggest that intratumoral delivery with HG-NPs has immense potential for GBM treatment, as well as other strategies where site-specific, long-term retention of therapeutic agents is warranted.

Innate immunity plays a key role in controlling viral load in COVID-19: mechanistic insights from a whole-body infection dynamics model.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a pathogen of immense public health concern. Efforts to control the disease have only proven mildly successful, and the disease will likely continue to cause excessive fatalities until effective preventative measures (such as a vaccine) are developed. To develop disease management strategies, a better understanding of SARS-CoV-2 pathogenesis and population susceptibility to infection are needed. To this end, physiologically-relevant mathematical modeling can provide a robust in silico tool to understand COVID-19 pathophysiology and the in vivo dynamics of SARS-CoV-2. Guided by ACE2-tropism (ACE2 receptor dependency for infection) of the virus, and by incorporating cellular-scale viral dynamics and innate and adaptive immune responses, we have developed a multiscale mechanistic model for simulating the time-dependent evolution of viral load distribution in susceptible organs of the body (respiratory tract, gut, liver, spleen, heart, kidneys, and brain). Following calibration with in vivo and clinical data, we used the model to simulate viral load progression in a virtual patient with varying degrees of compromised immune status. Further, we conducted global sensitivity analysis of model parameters and ranked them for their significance in governing clearance of viral load to understand the effects of physiological factors and underlying conditions on viral load dynamics. Antiviral drug therapy, interferon therapy, and their combination was simulated to study the effects on viral load kinetics of SARS-CoV-2. The model revealed the dominant role of innate immunity (specifically interferons and resident macrophages) in controlling viral load, and the importance of timing when initiating therapy following infection.

Image-guided mathematical modeling for pharmacological evaluation of nanomaterials and monoclonal antibodies.

While plasma concentration kinetics has traditionally been the predictor of drug pharmacological effects, it can occasionally fail to represent kinetics at the site of action, particularly for solid tumors. This is especially true in the case of delivery of therapeutic macromolecules (drug-loaded nanomaterials or monoclonal antibodies), which can experience challenges to effective delivery due to particle size-dependent diffusion barriers at the target site. As a result, disparity between therapeutic plasma kinetics and kinetics at the site of action may exist, highlighting the importance of target site concentration kinetics in determining the pharmacodynamic effects of macromolecular therapeutic agents. Assessment of concentration kinetics at the target site has been facilitated by non-invasive in vivo imaging modalities. This allows for visualization and quantification of the whole-body disposition behavior of therapeutics that is essential for a comprehensive understanding of their pharmacokinetics and pharmacodynamics. Quantitative non-invasive imaging can also help guide the development and parameterization of mathematical models for descriptive and predictive purposes. Here, we present a review of the application of state-of-the-art imaging modalities for quantitative pharmacological evaluation of therapeutic nanoparticles and monoclonal antibodies, with a focus on their integration with mathematical models, and identify challenges and opportunities. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease Diagnostic Tools > in vivo Nanodiagnostics and Imaging Nanotechnology Approaches to Biology > Nanoscale Systems in Biology.

A mathematical model to predict nanomedicine pharmacokinetics and tumor delivery.

Towards clinical translation of cancer nanomedicine, it is important to systematically investigate the various parameters related to nanoparticle (NP) physicochemical properties, tumor characteristics, and inter-individual variability that affect the tumor delivery efficiency of therapeutic nanomaterials. Comprehensive investigation of these parameters using traditional experimental approaches is impractical due to the vast parameter space; mathematical models provide a more tractable approach to navigate through such a multidimensional space. To this end, we have developed a predictive mathematical model of whole-body NP pharmacokinetics and their tumor delivery in vivo, and have conducted local and global sensitivity analyses to identify the factors that result in low tumor delivery efficiency and high off-target accumulation of NPs. Our analyses reveal that NP degradation rate, tumor blood viscosity, NP size, tumor vascular fraction, and tumor vascular porosity are the key parameters in governing NP kinetics in the tumor interstitium. The impact of these parameters on tumor delivery efficiency of NPs is discussed, and optimal values for maximizing NP delivery are presented.

A modeling platform for the lymphatic system.

We present a physiologically-based pharmacokinetic modeling platform capable of simulating the biodistribution of different therapeutic agents, including cells, their interactions within the immune system, redistribution across lymphoid compartments, and infiltration into tumor tissues. This transport-based platform comprises a distinctive implementation of a tumor compartment with spatial heterogeneity which enables the modeling of tumors of different size, necrotic state, and agent infiltration capacity. We provide three validating and three exploratory examples that illustrate the capabilities of the proposed approach. The results show that the model can recapitulate immune cell balance across different compartments, respond to antigen stimulation, simulate immune vaccine effects, and immune cell infiltration to tumors. Based on the results, the model can be used to study problems pertinent to current immunotherapies and has the potential to assist medical techniques that rely on the transport of biological species.