Optimizing tuberculosis treatment efficacy: Comparing the standard regimen with Moxifloxacin-containing regimens.

Tuberculosis (TB) continues to be one of the deadliest infectious diseases in the world, causing ~1.5 million deaths every year. The World Health Organization initiated an End TB Strategy that aims to reduce TB-related deaths in 2035 by 95%. Recent research goals have focused on discovering more effective and more patient-friendly antibiotic drug regimens to increase patient compliance and decrease emergence of resistant TB. Moxifloxacin is one promising antibiotic that may improve the current standard regimen by shortening treatment time. Clinical trials and in vivo mouse studies suggest that regimens containing moxifloxacin have better bactericidal activity. However, testing every possible combination regimen with moxifloxacin either in vivo or clinically is not feasible due to experimental and clinical limitations. To identify better regimens more systematically, we simulated pharmacokinetics/pharmacodynamics of various regimens (with and without moxifloxacin) to evaluate efficacies, and then compared our predictions to both clinical trials and nonhuman primate studies performed herein. We used GranSim, our well-established hybrid agent-based model that simulates granuloma formation and antibiotic treatment, for this task. In addition, we established a multiple-objective optimization pipeline using GranSim to discover optimized regimens based on treatment objectives of interest, i.e., minimizing total drug dosage and lowering time needed to sterilize granulomas. Our approach can efficiently test many regimens and successfully identify optimal regimens to inform pre-clinical studies or clinical trials and ultimately accelerate the TB regimen discovery process.

Optimizing Solid Tumor Treatment with Antibody-drug Conjugates Using Agent-Based Modeling: Considering the Role of a Carrier Dose and Payload Class.

INTRODUCTION: Antibody-drug conjugates (ADCs) show significant clinical efficacy in the treatment of solid tumors, but a major limitation to their success is poor intratumoral distribution. Adding a carrier dose improves both distribution and overall drug efficacy of ADCs, but the optimal carrier dose has not been outlined for different payload classes. OBJECTIVE: In this work, we study two carrier dose regimens: 1) matching payload potency to cellular delivery but potentially not reaching cells farther away from blood vessels, or 2) dosing to tumor saturation but risking a reduction in cell killing from a lower amount of payload delivered per cell. METHODS: We use a validated computational model to test four different payloads conjugated to trastuzumab to determine the optimal carrier dose as a function of target expression, ADC dose, and payload potency. RESULTS: We find that dosing to tumor saturation is more efficacious than matching payload potency to cellular delivery for all payloads because the increase in the number of cells targeted by the ADC outweighs the loss in cell killing on targeted cells. An important exception exists if the carrier dose reduces the payload uptake per cell to the point where all cell killing is lost. Likewise, receptor downregulation can mitigate the benefits of a carrier dose. CONCLUSIONS: Because tumor saturation and in vitro potency can be measured early in ADC design, these results provide insight into maximizing ADC efficacy and demonstrate the benefits of using simulation to guide ADC design.

Simulating the Selection of Resistant Cells with Bystander Killing and Antibody Coadministration in Heterogeneous Human Epidermal Growth Factor Receptor 2-Positive Tumors.

Intratumoral heterogeneity is a leading cause of treatment failure resulting in tumor recurrence. For the antibody-drug conjugate (ADC) ado-trastuzumab emtansine (T-DM1), two major types of resistance include changes in human epidermal growth factor receptor 2 (HER2) expression and reduced payload sensitivity, which is often exacerbated by heterogenous HER2 expression and ADC distribution during treatment. ADCs with bystander payloads, such as trastuzumab-monomethyl auristatin E (T-MMAE), can reach and kill adjacent cells with lower receptor expression that cannot be targeted directly with the ADC. Additionally, coadministration of T-DM1 with its unconjugated antibody, trastuzumab, can improve distribution and minimize heterogeneous delivery. However, the effectiveness of trastuzumab coadministration and ADC bystander killing in heterogenous tumors in reducing the selection of resistant cells is not well understood. Here, we use an agent-based model to predict outcomes with these different regimens. The simulations demonstrate that both T-DM1 and T-MMAE benefit from trastuzumab coadministration for tumors with high average receptor expression (up to 70% and 40% decrease in average tumor volume, respectively), with greater benefit for nonbystander payloads. However, the benefit decreases as receptor expression is reduced, reversing at low concentrations (up to 360% and 430% increase in average tumor volume for T-DM1 and T-MMAE, respectively) for this mechanism that impacts both ADC distribution and efficacy. For tumors with intrinsic payload resistance, coadministration uniformly exhibits better efficacy than ADC monotherapy (50%-70% and 19%-36% decrease in average tumor volume for T-DM1 and T-MMAE, respectively). Finally, we demonstrate that several regimens select for resistant cells at clinical tolerable doses, which highlights the need to pursue other mechanisms of action for durable treatment responses. SIGNIFICANCE STATEMENT: Experimental evidence demonstrates heterogeneity in the distribution of both the antibody-drug conjugate and the target receptor in the tumor microenvironment, which can promote the selection of resistant cells and lead to recurrence. This study quantifies the impact of increasing the antibody dose and utilizing bystander payloads in heterogeneous tumors. Alternative cell-killing mechanisms are needed to avoid enriching resistant cell populations.

Computational modeling implicates protein scaffolding in p38 regulation of Akt.

Cells process environmental cues by activating intracellular signaling pathways with numerous interconnections and opportunities for cross-regulation. We employed a systems biology approach to investigate intersections of kinase p38, a context-dependent tumor suppressor or promoter, with Akt and ERK, two kinases known to promote cell survival, proliferation, and drug resistance in cancer. Using live, single cell microscopy, multiplexed fluorescent reporters of p38, Akt, and ERK activities, and a custom automated image-processing pipeline, we detected marked heterogeneity of signaling outputs in breast cancer cells stimulated with chemokine CXCL12 or epidermal growth factor (EGF). Basal activity of p38 correlated inversely with amplitude of Akt and ERK activation in response to either ligand. Remarkably, small molecule inhibitors of p38 immediately decreased basal activities of Akt and ERK but increased the proportion of cells with high amplitude ligand-induced activation of Akt signaling. To identify mechanisms underlying cross-talk of p38 with Akt signaling, we developed a computational model incorporating subcellular compartmentalization of signaling molecules by scaffold proteins. Dynamics of this model revealed that subcellular scaffolding of Akt accounted for observed regulation by p38. The model also predicted that differences in the amount of scaffold protein in a subcellular compartment captured the observed single cell heterogeneity in signaling. Finally, our model predicted that reduction in kinase signaling can be accomplished by both scaffolding and direct kinase inhibition. However, scaffolding inhibition can potentiate future kinase activity by redistribution of pathway components, potentially amplifying oncogenic signaling. These studies reveal how computational modeling can decipher mechanisms of cross-talk between the p38 and Akt signaling pathways and point to scaffold proteins as central regulators of signaling dynamics and amplitude.

A virtual host model of Mycobacterium tuberculosis infection identifies early immune events as predictive of infection outcomes.

Tuberculosis (TB), caused by infection with Mycobacterium tuberculosis (Mtb), is one of the world's deadliest infectious diseases and remains a significant global health burden. TB disease and pathology can present clinically across a spectrum of outcomes, ranging from total sterilization of infection to active disease. Much remains unknown about the biology that drives an individual towards various clinical outcomes as it is challenging to experimentally address specific mechanisms driving clinical outcomes. Furthermore, it is unknown whether numbers of immune cells in the blood accurately reflect ongoing events during infection within human lungs. Herein, we utilize a systems biology approach by developing a whole-host model of the immune response to Mtb across multiple physiologic and time scales. This model, called HostSim, tracks events at the cellular, granuloma, organ, and host scale and represents the first whole-host, multi-scale model of the immune response following Mtb infection. We show that this model can capture various aspects of human and non-human primate TB disease and predict that biomarkers in the blood may only faithfully represent events in the lung at early time points after infection. We posit that HostSim, as a first step toward personalized digital twins in TB research, offers a powerful computational tool that can be used in concert with experimental approaches to understand and predict events about various aspects of TB disease and therapeutics.

Dynamic balance of pro- and anti-inflammatory signals controls disease and limits pathology.

Immune responses to pathogens are complex and not well understood in many diseases, and this is especially true for infections by persistent pathogens. One mechanism that allows for long-term control of infection while also preventing an over-zealous inflammatory response from causing extensive tissue damage is for the immune system to balance pro- and anti-inflammatory cells and signals. This balance is dynamic and the immune system responds to cues from both host and pathogen, maintaining a steady state across multiple scales through continuous feedback. Identifying the signals, cells, cytokines, and other immune response factors that mediate this balance over time has been difficult using traditional research strategies. Computational modeling studies based on data from traditional systems can identify how this balance contributes to immunity. Here we provide evidence from both experimental and mathematical/computational studies to support the concept of a dynamic balance operating during persistent and other infection scenarios. We focus mainly on tuberculosis, currently the leading cause of death due to infectious disease in the world, and also provide evidence for other infections. A better understanding of the dynamically balanced immune response can help shape treatment strategies that utilize both drugs and host-directed therapies.

A computational model tracks whole-lung Mycobacterium tuberculosis infection and predicts factors that inhibit dissemination.

Mycobacterium tuberculosis (Mtb), the causative infectious agent of tuberculosis (TB), kills more individuals per year than any other infectious agent. Granulomas, the hallmark of Mtb infection, are complex structures that form in lungs, composed of immune cells surrounding bacteria, infected cells, and a caseous necrotic core. While granulomas serve to physically contain and immunologically restrain bacteria growth, some granulomas are unable to control Mtb growth, leading to bacteria and infected cells leaving the granuloma and disseminating, either resulting in additional granuloma formation (local or non-local) or spread to airways or lymph nodes. Dissemination is associated with development of active TB. It is challenging to experimentally address specific mechanisms driving dissemination from TB lung granulomas. Herein, we develop a novel hybrid multi-scale computational model, MultiGran, that tracks Mtb infection within multiple granulomas in an entire lung. MultiGran follows cells, cytokines, and bacterial populations within each lung granuloma throughout the course of infection and is calibrated to multiple non-human primate (NHP) cellular, granuloma, and whole-lung datasets. We show that MultiGran can recapitulate patterns of in vivo local and non-local dissemination, predict likelihood of dissemination, and predict a crucial role for multifunctional CD8+ T cells and macrophage dynamics for preventing dissemination.

A Two-Pulse Cellular Stimulation Test Elucidates Variability and Mechanisms in Signaling Pathways.

Mammalian cells respond in a variable manner when provided with physiological pulses of ligand, such as low concentrations of acetylcholine present for just tens of seconds or TNFα for just tens of minutes. For a two-pulse stimulation, some cells respond to both pulses, some do not respond, and yet others respond to only one or the other pulse. Are these different response patterns the result of the small number of ligands being able to only stochastically activate the pathway at random times or an output pattern from a deterministic algorithm responding differently to different stimulation intervals? If the response is deterministic in nature, what parameters determine whether a response is generated or skipped? To answer these questions, we developed a two-pulse test that utilizes different rest periods between stimulation pulses. This "rest-period test" revealed that cells skip responses predictably as the rest period is shortened. By combining these experimental results with a mathematical model of the pathway, we further obtained mechanistic insight into potential sources of response variability. Our analysis indicates that in both intracellular calcium and NFκB signaling, response variability is consistent with extrinsic noise (cell-to-cell variability in protein levels), a short-term memory of stimulation, and high Hill coefficient processes. Furthermore, these results support recent works that have emphasized the role of deterministic processes for explaining apparently stochastic cellular response variability and indicate that even weak stimulations likely guide mammalian cells to appropriate fates rather than leaving outcomes to chance. We envision that the rest-period test can be applied to other signaling pathways to extract mechanistic insight.

Comparing efficacies of moxifloxacin, levofloxacin and gatifloxacin in tuberculosis granulomas using a multi-scale systems pharmacology approach.

Granulomas are complex lung lesions that are the hallmark of tuberculosis (TB). Understanding antibiotic dynamics within lung granulomas will be vital to improving and shortening the long course of TB treatment. Three fluoroquinolones (FQs) are commonly prescribed as part of multi-drug resistant TB therapy: moxifloxacin (MXF), levofloxacin (LVX) or gatifloxacin (GFX). To date, insufficient data are available to support selection of one FQ over another, or to show that these drugs are clinically equivalent. To predict the efficacy of MXF, LVX and GFX at a single granuloma level, we integrate computational modeling with experimental datasets into a single mechanistic framework, GranSim. GranSim is a hybrid agent-based computational model that simulates granuloma formation and function, FQ plasma and tissue pharmacokinetics and pharmacodynamics and is based on extensive in vitro and in vivo data. We treat in silico granulomas with recommended daily doses of each FQ and compare efficacy by multiple metrics: bacterial load, sterilization rates, early bactericidal activity and efficacy under non-compliance and treatment interruption. GranSim reproduces in vivo plasma pharmacokinetics, spatial and temporal tissue pharmacokinetics and in vitro pharmacodynamics of these FQs. We predict that MXF kills intracellular bacteria more quickly than LVX and GFX due in part to a higher cellular accumulation ratio. We also show that all three FQs struggle to sterilize non-replicating bacteria residing in caseum. This is due to modest drug concentrations inside caseum and high inhibitory concentrations for this bacterial subpopulation. MXF and LVX have higher granuloma sterilization rates compared to GFX; and MXF performs better in a simulated non-compliance or treatment interruption scenario. We conclude that MXF has a small but potentially clinically significant advantage over LVX, as well as LVX over GFX. We illustrate how a systems pharmacology approach combining experimental and computational methods can guide antibiotic selection for TB.

Identifying mechanisms driving formation of granuloma-associated fibrosis during Mycobacterium tuberculosis infection.

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), is a pulmonary pathogen of major global concern. A key feature of Mtb infection in primates is the formation of granulomas, dense cellular structures surrounding infected lung tissue. These structures serve as the main site of host-pathogen interaction in TB, and thus to effectively treat TB we must clarify mechanisms of granuloma formation and their function in disease. Fibrotic granulomas are associated with both good and bad disease outcomes. Fibrosis can serve to isolate infected tissue from healthy tissue, but it can also cause difficulty breathing as it leaves scars. Little is known about fibrosis in TB, and data from non-human primates is just beginning to clarify the picture. This work focuses on constructing a hybrid multi-scale model of fibrotic granuloma formation, in order to identify mechanisms driving development of fibrosis in Mtb infected lungs. We combine dynamics of molecular, cellular, and tissue scale models from previously published studies to characterize the formation of two common sub-types of fibrotic granulomas: peripherally fibrotic, with a cuff of collagen surrounding granulomas, and centrally fibrotic, with collagen throughout granulomas. Uncertainty and sensitivity analysis, along with large simulation sets, enable us to identify mechanisms differentiating centrally versus peripherally fibrotic granulomas. These findings suggest that heterogeneous cytokine environments exist within granulomas and may be responsible for driving tissue scale morphologies. Using this model we are primed to better understand the complex structure of granulomas, a necessity for developing successful treatments for TB.

Computational and Empirical Studies Predict Mycobacterium tuberculosis-Specific T Cells as a Biomarker for Infection Outcome.

Identifying biomarkers for tuberculosis (TB) is an ongoing challenge in developing immunological correlates of infection outcome and protection. Biomarker discovery is also necessary for aiding design and testing of new treatments and vaccines. To effectively predict biomarkers for infection progression in any disease, including TB, large amounts of experimental data are required to reach statistical power and make accurate predictions. We took a two-pronged approach using both experimental and computational modeling to address this problem. We first collected 200 blood samples over a 2- year period from 28 non-human primates (NHP) infected with a low dose of Mycobacterium tuberculosis. We identified T cells and the cytokines that they were producing (single and multiple) from each sample along with monkey status and infection progression data. Machine learning techniques were used to interrogate the experimental NHP datasets without identifying any potential TB biomarker. In parallel, we used our extensive novel NHP datasets to build and calibrate a multi-organ computational model that combines what is occurring at the site of infection (e.g., lung) at a single granuloma scale with blood level readouts that can be tracked in monkeys and humans. We then generated a large in silico repository of in silico granulomas coupled to lymph node and blood dynamics and developed an in silico tool to scale granuloma level results to a full host scale to identify what best predicts Mycobacterium tuberculosis (Mtb) infection outcomes. The analysis of in silico blood measures identifies Mtb-specific frequencies of effector T cell phenotypes at various time points post infection as promising indicators of infection outcome. We emphasize that pairing wetlab and computational approaches holds great promise to accelerate TB biomarker discovery.

Computational modeling predicts IL-10 control of lesion sterilization by balancing early host immunity-mediated antimicrobial responses with caseation during mycobacterium tuberculosis infection.

Although almost a third of the world's population is infected with the bacterial pathogen Mycobacterium tuberculosis, our understanding of the functions of many immune factors involved in fighting infection is limited. Determining the role of the immunosuppressive cytokine IL-10 at the level of the granuloma has proven difficult because of lesional heterogeneity and the limitations of animal models. In this study, we take an in silico approach and, through a series of virtual experiments, we predict several novel roles for IL-10 in tuberculosis granulomas: 1) decreased levels of IL-10 lead to increased numbers of sterile lesions, but at the cost of early increased caseation; 2) small increases in early antimicrobial activity cause this increased lesion sterility; 3) IL-10 produced by activated macrophages is a major mediator of early antimicrobial activity and early host-induced caseation; and 4) increasing levels of infected macrophage derived IL-10 promotes bacterial persistence by limiting the early antimicrobial response and preventing lesion sterilization. Our findings, currently only accessible using an in silico approach, suggest that IL-10 at the individual granuloma scale is a critical regulator of lesion outcome. These predictions suggest IL-10-related mechanisms that could be used as adjunctive therapies during tuberculosis.

Multiscale Model of Mycobacterium tuberculosis Infection Maps Metabolite and Gene Perturbations to Granuloma Sterilization Predictions.

Granulomas are a hallmark of tuberculosis. Inside granulomas, the pathogen Mycobacterium tuberculosis may enter a metabolically inactive state that is less susceptible to antibiotics. Understanding M. tuberculosis metabolism within granulomas could contribute to reducing the lengthy treatment required for tuberculosis and provide additional targets for new drugs. Two key adaptations of M. tuberculosis are a nonreplicating phenotype and accumulation of lipid inclusions in response to hypoxic conditions. To explore how these adaptations influence granuloma-scale outcomes in vivo, we present a multiscale in silico model of granuloma formation in tuberculosis. The model comprises host immunity, M. tuberculosis metabolism, M. tuberculosis growth adaptation to hypoxia, and nutrient diffusion. We calibrated our model to in vivo data from nonhuman primates and rabbits and apply the model to predict M. tuberculosis population dynamics and heterogeneity within granulomas. We found that bacterial populations are highly dynamic throughout infection in response to changing oxygen levels and host immunity pressures. Our results indicate that a nonreplicating phenotype, but not lipid inclusion formation, is important for long-term M. tuberculosis survival in granulomas. We used virtual M. tuberculosis knockouts to predict the impact of both metabolic enzyme inhibitors and metabolic pathways exploited to overcome inhibition. Results indicate that knockouts whose growth rates are below ∼66% of the wild-type growth rate in a culture medium featuring lipid as the only carbon source are unable to sustain infections in granulomas. By mapping metabolite- and gene-scale perturbations to granuloma-scale outcomes and predicting mechanisms of sterilization, our method provides a powerful tool for hypothesis testing and guiding experimental searches for novel antituberculosis interventions.

Macrophage polarization drives granuloma outcome during Mycobacterium tuberculosis infection.

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), induces formation of granulomas, structures in which immune cells and bacteria colocalize. Macrophages are among the most abundant cell types in granulomas and have been shown to serve as both critical bactericidal cells and targets for M. tuberculosis infection and proliferation throughout the course of infection. Very little is known about how these processes are regulated, what controls macrophage microenvironment-specific polarization and plasticity, or why some granulomas control bacteria and others permit bacterial dissemination. We take a computational-biology approach to investigate mechanisms that drive macrophage polarization, function, and bacterial control in granulomas. We define a "macrophage polarization ratio" as a metric to understand how cytokine signaling translates into polarization of single macrophages in a granuloma, which in turn modulates cellular functions, including antimicrobial activity and cytokine production. Ultimately, we extend this macrophage ratio to the tissue scale and define a "granuloma polarization ratio" describing mean polarization measures for entire granulomas. Here we coupled experimental data from nonhuman primate TB granulomas to our computational model, and we predict two novel and testable hypotheses regarding macrophage profiles in TB outcomes. First, the temporal dynamics of granuloma polarization ratios are predictive of granuloma outcome. Second, stable necrotic granulomas with low CFU counts and limited inflammation are characterized by short NF-κB signal activation intervals. These results suggest that the dynamics of NF-κB signaling is a viable therapeutic target to promote M1 polarization early during infection and to improve outcome.

A computational tool integrating host immunity with antibiotic dynamics to study tuberculosis treatment.

While active tuberculosis (TB) is a treatable disease, many complex factors prevent its global elimination. Part of the difficulty in developing optimal therapies is the large design space of antibiotic doses, regimens and combinations. Computational models that capture the spatial and temporal dynamics of antibiotics at the site of infection can aid in reducing the design space of costly and time-consuming animal pre-clinical and human clinical trials. The site of infection in TB is the granuloma, a collection of immune cells and bacteria that form in the lung, and new data suggest that penetration of drugs throughout granulomas is problematic. Here we integrate our computational model of granuloma formation and function with models for plasma pharmacokinetics, lung tissue pharmacokinetics and pharmacodynamics for two first line anti-TB antibiotics. The integrated model is calibrated to animal data. We make four predictions. First, antibiotics are frequently below effective concentrations inside granulomas, leading to bacterial growth between doses and contributing to the long treatment periods required for TB. Second, antibiotic concentration gradients form within granulomas, with lower concentrations toward their centers. Third, during antibiotic treatment, bacterial subpopulations are similar for INH and RIF treatment: mostly intracellular with extracellular bacteria located in areas non-permissive for replication (hypoxic areas), presenting a slowly increasing target population over time. Finally, we find that on an individual granuloma basis, pre-treatment infection severity (including bacterial burden, host cell activation and host cell death) is predictive of treatment outcome.

Identifying Mechanisms of Homeostatic Signaling in Fibroblast Differentiation.

Fibroblasts play an important role in the wound-healing process by generating extracellular matrix (ECM) and undergoing differentiation into myofibroblasts, but these cells can also be involved in pathologic remodeling of tissue. Nascent ECM provides a substrate for re-epithelialization to occur, restoring damaged tissue to a functional state. Dysregulation of this process can result in fibrosis--stiffening and scarring of the tissue. Current treatments cannot halt or reverse this process. The molecular mechanisms underlying fibrotic dysregulation are poorly understood, providing an untapped pool of potential therapeutic targets. Transforming growth factor-ß (TGF-ß) and adhesion signaling are involved in inducing fibroblast differentiation into α-smooth muscle actin (αSMA) expressing myofibroblasts, while prostaglandin E2 (PGE2) has been shown to antagonize TGF-ß signaling; however, the temporal and mechanistic details of this relationship have not yet been fully characterized. We measured αSMA, a marker of fibroblast to myofibroblast differentiation, as a function of: TGF-ß1 receptor-ligand complex internalization, PGE2 binding, and adhesion signaling and developed a mathematical model capturing the molecular mechanisms of fibroblast differentiation. Using our model, we predict the following: Periodic dosing with PGE2 temporarily renders fibroblasts incapable of differentiation and refractory to additional TGF-ß1 stimulation; conversely, periodic dosing with TGF-ß1 in the presence of PGE2 induces a reduced signal response that can be further inhibited by the addition of more PGE2. Controlled fibroblast differentiation is necessary for effective wound healing; however, excessive accumulation of αSMA-expressing myofibroblasts can result in fibrosis. Homeostasis of αSMA in our model requires a balance of positive and negative regulatory signals. Sensitivity analysis predicts that PGE2 availability, TGF-ß1 availability, and the rate of TGF-ß1 receptor recycling each highly influence the rates of αSMA production. With this model, we are able to demonstrate that regulation of both TGF-ß1 and PGE2 signaling levels is essential for preventing fibroblast dysregulation.

In silico models of M. tuberculosis infection provide a route to new therapies.

Tuberculosis (TB) is a global health problem responsible for ~2 million deaths per year. Current antibiotic treatments are lengthy and fraught with compliance and resistance issues. There is a crucial need for additional approaches to provide a cost-effective means of exploring the design space for potential therapies. We discuss the use of mathematical and computational models in virtual experiments and virtual clinical trials both to develop new hypotheses regarding the disease and to provide a cost-effective means of discovering new treatment strategies.

Differential risk of tuberculosis reactivation among anti-TNF therapies is due to drug binding kinetics and permeability.

Increased rates of tuberculosis (TB) reactivation have been reported in humans treated with TNF-α (TNF)-neutralizing drugs, and higher rates are observed with anti-TNF Abs (e.g., infliximab) as compared with TNF receptor fusion protein (etanercept). Mechanisms driving differential reactivation rates and differences in drug action are not known. We use a computational model of a TB granuloma formation that includes TNF/TNF receptor dynamics to elucidate these mechanisms. Our analyses yield three important insights. First, drug binding to membrane-bound TNF critically impairs granuloma function. Second, a higher risk of reactivation induced from Ab-type treatments is primarily due to differences in TNF/drug binding kinetics and permeability. Apoptotic and cytolytic activities of Abs and pharmacokinetic fluctuations in blood concentration of drug are not essential to inducing TB reactivation. Third, we predict specific host factors that, if augmented, would improve granuloma function during anti-TNF therapy. Our findings have implications for the development of safer anti-TNF drugs to treat inflammatory diseases.

Characterizing heterogeneous single-cell dose responses computationally and experimentally using threshold inhibition surfaces and dose-titration assays.

Single cancer cells within a tumor exhibit variable levels of resistance to drugs, ultimately leading to treatment failures. While tumor heterogeneity is recognized as a major obstacle to cancer therapy, standard dose-response measurements for the potency of targeted kinase inhibitors aggregate populations of cells, obscuring intercellular variations in responses. In this work, we develop an analytical and experimental framework to quantify and model dose responses of individual cancer cells to drugs. We first explore the connection between population and single-cell dose responses using a computational model, revealing that multiple heterogeneous populations can yield nearly identical population dose responses. We demonstrate that a single-cell analysis method, which we term a threshold inhibition surface, can differentiate among these populations. To demonstrate the applicability of this method, we develop a dose-titration assay to measure dose responses in single cells. We apply this assay to breast cancer cells responding to phosphatidylinositol-3-kinase inhibition (PI3Ki), using clinically relevant PI3Kis on breast cancer cell lines expressing fluorescent biosensors for kinase activity. We demonstrate that MCF-7 breast cancer cells exhibit heterogeneous dose responses with some cells requiring over ten-fold higher concentrations than the population average to achieve inhibition. Our work reimagines dose-response relationships for cancer drugs in an emerging paradigm of single-cell tumor heterogeneity.

In vivo Auto-tuning of Antibody-Drug Conjugate Delivery for Effective Immunotherapy using High-Avidity, Low-Affinity Antibodies.

Antibody-drug conjugates (ADCs) have experienced a surge in clinical approvals in the past five years. Despite this success, a major limitation to ADC efficacy in solid tumors is poor tumor penetration, which leaves many cancer cells untargeted. Increasing antibody doses or co-administering ADC with an unconjugated antibody can improve tumor penetration and increase efficacy when target receptor expression is high. However, it can also reduce efficacy in low-expression tumors where ADC delivery is limited by cellular uptake. This creates an intrinsic problem because many patients express different levels of target between tumors and even within the same tumor. Here, we generated High-Avidity, Low-Affinity (HALA) antibodies that can automatically tune the cellular ADC delivery to match the local expression level. Using HER2 ADCs as a model, HALA antibodies were identified with the desired HER2 expression-dependent competitive binding with ADCs in vitro. Multi-scale distribution of trastuzumab emtansine and trastuzumab deruxtecan co-administered with the HALA antibody were analyzed in vivo, revealing that the HALA antibody increased ADC tumor penetration in high-expression systems with minimal reduction in ADC uptake in low-expression tumors. This translated to greater ADC efficacy in immunodeficient mouse models across a range of HER2 expression levels. Furthermore, Fc-enhanced HALA antibodies showed improved Fc-effector function at both high and low expression levels and elicited a strong response in an immunocompetent mouse model. These results demonstrate that HALA antibodies can expand treatment ranges beyond high expression targets and leverage strong immune responses.