Erythrocytes are proposed to be involved in blood flow regulation through both shear- and oxygen-dependent mechanisms for the release of adenosine triphosphate (ATP), a potent vasodilator. In a recent study, the dynamics of shear-dependent ATP release from erythrocytes was measured using a microfluidic device with a constriction in the channel to increase shear stress. The brief period of increased shear stress resulted in ATP release within 25 to 75 milliseconds downstream of the constriction. The long-term goal of our research is to apply a similar approach to determine the dynamics of oxygen-dependent ATP release. In the place of the constriction, an oxygen permeable membrane would be used to decrease the hemoglobin oxygen saturation of erythrocytes flowing through the channel. This paper describes the first stage in achieving that goal, the development of a computational model of the proposed experimental system to determine the feasibility of altering oxygen saturation rapidly enough to measure ATP release dynamics. The computational model was constructed based on hemodynamics, molecular transport of oxygen and ATP, kinetics of luciferin/luciferase reaction for reporting ATP concentrations, light absorption by hemoglobin, and sensor characteristics. A linear model of oxygen saturation-dependent ATP release with variable time delay was used in this study. The computational results demonstrate that a microfluidic device with a 100 µm deep channel will cause a rapid decrease in oxygen saturation over the oxygen permeable membrane that yields a measurable light intensity profile for a change in rate of ATP release from erythrocytes on a timescale as short as 25 milliseconds. The simulation also demonstrates that the complex dynamics of ATP release from erythrocytes combined with the consumption by luciferin/luciferase in a flowing system results in light intensity values that do not simply correlate with ATP concentrations. A computational model is required for proper interpretation of experimental data.

The success of insects in terrestrial environments is due in large part to their ability to resist desiccation stress. Since the majority of water is lost across the cuticle, a relatively water-impermeable cuticle is a major component of insect desiccation resistance. Cuticular permeability is affected by the properties and mixing effects of component hydrocarbons, and changes in cuticular hydrocarbons can affect desiccation tolerance. A pre-exposure to a mild desiccation stress increases duration of desiccation survival in adult female Drosophila melanogaster, via a decrease in cuticular permeability. To test whether this acute response to desiccation stress is due to a change in cuticular hydrocarbons, we treated male and female D. melanogaster to a rapid desiccation hardening (RDH) treatment and used gas chromatography to examine the effects on cuticular hydrocarbon composition. RDH led to reduced proportions of unsaturated and methylated hydrocarbons compared to controls in females, but although RDH modified the cuticular hydrocarbon profile in males, there was no coordinated pattern. These data suggest that the phenomenon of RDH leading to reduced cuticular water loss occurs via an acute change in cuticular hydrocarbons that enhances desiccation tolerance in female, but not male, D. melanogaster.

BackgroundImmune checkpoint blockade therapy has clearly shown clinical activity in patients with triple-negative breast cancer, but less than half of the patients benefit from the treatments. While a number of ongoing clinical trials are investigating different combinations of checkpoint inhibitors and chemotherapeutic agents, predictive biomarkers that identify patients most likely to benefit remains one of the major challenges. Here we present a modular quantitative systems pharmacology (QSP) platform for immuno-oncology that incorporates detailed mechanisms of immune–cancer cell interactions to make efficacy predictions and identify predictive biomarkers for treatments using atezolizumab and nab-paclitaxel.MethodsA QSP model was developed based on published data of triple-negative breast cancer. With the model, we generated a virtual patient cohort to conduct in silico virtual clinical trials and make retrospective analyses of the pivotal IMpassion130 trial that led to the accelerated approval of atezolizumab and nab-paclitaxel for patients with programmed death-ligand 1 (PD-L1) positive triple-negative breast cancer. Available data from clinical trials were used for model calibration and validation.ResultsWith the calibrated virtual patient cohort based on clinical data from the placebo comparator arm of the IMpassion130 trial, we made efficacy predictions and identified potential predictive biomarkers for the experimental arm of the trial using the proposed QSP model. The model predictions are consistent with clinically reported efficacy endpoints and correlated immune biomarkers. We further performed a series of virtual clinical trials to compare different doses and schedules of the two drugs for simulated therapeutic optimization.ConclusionsThis study provides a QSP platform, which can be used to generate virtual patient cohorts and conduct virtual clinical trials. Our findings demonstrate its potential for making efficacy predictions for immunotherapies and chemotherapies, identifying predictive biomarkers, and guiding future clinical trial designs.

The survival rate of patients with breast cancer has been improved by immune checkpoint blockade therapies, and the efficacy of their combinations with epigenetic modulators has shown promising results in preclinical studies. In this prospective study, we propose an ordinary differential equation (ODE)-based quantitative systems pharmacology (QSP) model to conduct an in silico virtual clinical trial and analyze potential predictive biomarkers to improve the anti-tumor response in HER2-negative breast cancer. The model is comprised of four compartments: central, peripheral, tumor, and tumor-draining lymph node, and describes immune activation, suppression, T cell trafficking, and pharmacokinetics and pharmacodynamics (PK/PD) of the therapeutic agents. We implement theoretical mechanisms of action for checkpoint inhibitors and the epigenetic modulator based on preclinical studies to investigate their effects on antitumor response. According to model-based simulations, we confirm the synergistic effect of the epigenetic modulator and that pre-treatment tumor mutational burden, tumor-infiltrating effector T cell (Teff) density, and Teff to regulatory T cell (Treg) ratio are significantly higher in responders, which can be potential biomarkers to be considered in clinical trials. Overall, we present a readily reproducible modular model to conduct in silico virtual clinical trials on patient cohorts of interest, which is a step toward personalized medicine in cancer immunotherapy.

Macrophages respond to signals in the microenvironment by changing their functional phenotypes, a process known as polarization. Depending on the context, they acquire different patterns of transcriptional activation, cytokine expression and cellular metabolism which collectively constitute a continuous spectrum of phenotypes, of which the two extremes are denoted as classical (M1) and alternative (M2) activation. To quantitatively decode the underlying principles governing macrophage phenotypic polarization and thereby harness its therapeutic potential in human diseases, a systems-level approach is needed given the multitude of signaling pathways and intracellular regulation involved. Here we develop the first mechanism-based, multi-pathway computational model that describes the integrated signal transduction and macrophage programming under M1 (IFN-γ), M2 (IL-4) and cell stress (hypoxia) stimulation. Our model was calibrated extensively against experimental data, and we mechanistically elucidated several signature feedbacks behind the M1-M2 antagonism and investigated the dynamical shaping of macrophage phenotypes within the M1-M2 spectrum. Model sensitivity analysis also revealed key molecular nodes and interactions as targets with potential therapeutic values for the pathophysiology of peripheral arterial disease and cancer. Through simulations that dynamically capture the signal integration and phenotypic marker expression in the differential macrophage polarization responses, our model provides an important computational basis toward a more quantitative and network-centric understanding of the complex physiology and versatile functions of macrophages in human diseases.

Cancer immunotherapy has recently drawn remarkable attention as promising results in the clinic have shown its ability to improve the overall survival, and T cells are considered to be one of the primary effectors for cancer immunotherapy. Enhanced and restored T cell tumoricidal activity has shown great potential for killing cancer cells. Bispecific T cell engagers (TCEs) are a growing class of molecules that are designed to bind two different antigens on the surface of T cells and cancer cells to bring them in close proximity and selectively activate effector T cells to kill target cancer cells. New T cell engagers are being investigated for the treatment of solid tumors. The activity of newly developed T cell engagers showed a strong correlation with tumor target antigen expression. However, the correlation between tumor-associated antigen expression and overall response of cancer patients is poorly understood. In this study, we used a well-calibrated quantitative systems pharmacology (QSP) model extended to bispecific T cell engagers to explore their efficacy and identify potential biomarkers. In principle, patient-specific response can be predicted through this model according to each patient’s individual characteristics. This extended QSP model has been calibrated with available experimental data and provides predictions of patients’ response to TCE treatment.

Immunotherapy has shown great potential in the treatment of cancer; however, only a fraction of patients respond to treatment, and many experience autoimmune‐related side effects. The pharmaceutical industry has relied on mathematical models to study the behavior of candidate drugs and more recently, complex, whole‐body, quantitative systems pharmacology (QSP) models have become increasingly popular for discovery and development. QSP modeling has the potential to discover novel predictive biomarkers as well as test the efficacy of treatment plans and combination therapies through virtual clinical trials. In this work, we present a QSP modeling platform for immuno‐oncology (IO) that incorporates detailed mechanisms for important immune interactions. This modular platform allows for the construction of QSP models of IO with varying degrees of complexity based on the research questions. Finally, we demonstrate the use of the platform through two example applications of immune checkpoint therapy.

BackgroundT cells have been recognized as core effectors for cancer immunotherapy. How to restore the anti-tumor ability of suppressed T cells or improve the lethality of cytotoxic T cells has become the main focus in immunotherapy. Bispecific antibodies, especially bispecific T cell engagers (TCEs), have shown their unique ability to enhance the patient’s immune response to tumors by stimulating T cell activation and cytokine production in an MHC-independent manner. Antibodies targeting the checkpoint inhibitory molecules such as programmed cell death protein 1 (PD-1), PD-ligand 1 (PD-L1) and cytotoxic lymphocyte activated antigen 4 are able to restore the cytotoxic effect of immune suppressed T cells and have also shown durable responses in patients with malignancies. However, both types have their own limitations in treating certain cancers. Preclinical and clinical results have emphasized the potential of combining these two antibodies to improve tumor response and patients’ survival. However, the selection and evaluation of combination partners clinically is a costly endeavor. In addition, despite advances made in immunotherapy, there are subsets of patients who are non-responders, and reliable biomarkers for different immunotherapies are urgently needed to improve the ability to prospectively predict patients’ response and improve clinical study design. Therefore, mathematical and computational models are essential to optimize patient benefit, and guide combination approaches with lower cost and in a faster manner.MethodIn this study, we continued to extend the quantitative systems pharmacology (QSP) model we developed for a bispecific TCE to explore efficacy of combination therapy with an anti-PD-L1 monoclonal antibody in patients with colorectal cancer.ResultsPatient-specific response to TCE monotherapy, anti-PD-L1 monotherapy and the combination therapy were predicted using this model according to each patient’s individual characteristics.ConclusionsIndividual biomarkers for TCE monotherapy, anti-PD-L1 monotherapy and their combination have been determined based on the QSP model. Best treatment options for specific patients could be suggested based on their own characteristics to improve clinical trial efficiency. The model can be further used to assess plausible combination strategies for different TCEs and immune checkpoint inhibitors in different types of cancer.