Search

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 anEnd TB Strategythat 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 andin vivomouse studies suggest that regimens containing moxifloxacin have better bactericidal activity. However, testing every possible combination regimen with moxifloxacin eitherin vivoor 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 usedGranSim, 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 usingGranSimto 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.Author summaryTuberculosis (TB) is a top global health concern and treatment for TB requires multiple antibiotics taken for long periods of time, which is challenging for TB patients. Therefore, identifying regimens that are more effective and more patient-friendly than the standard treatment is urgently needed. It is also known that non-compliance leads to the development of drug resistant TB. In this work, we pair computational and experimental models to predict new regimens for the treatment of TB that optimize how fast bacteria are cleared using minimal dosage. We apply novel approaches to this goal and validate our predictions using a non-human primate model. Our findings suggest that systems pharmacological modeling should be employed as a method to narrow the design space for drug regimens for tuberculosis and other diseases as well.

INTRODUCTION: Mathematical and computational modeling have a long history of uncovering mechanisms and making predictions for biological systems. However, to create a model that can provide relevant quantitative predictions, models must first be calibrated by recapitulating existing biological datasets from that system. Current calibration approaches may not be appropriate for complex biological models because: 1) many attempt to recapitulate only a single aspect of the experimental data (such as a median trend) or 2) Bayesian techniques require specification of parameter priors and likelihoods to experimental data that cannot always be confidently assigned. A new calibration protocol is needed to calibrate complex models when current approaches fall short. METHODS: Herein, we develop CaliPro, an iterative, model-agnostic calibration protocol that utilizes parameter density estimation to refine parameter space and calibrate to temporal biological datasets. An important aspect of CaliPro is the user-defined pass set definition, which specifies how the model might successfully recapitulate experimental data. We define the appropriate settings to use CaliPro. RESULTS: We illustrate the usefulness of CaliPro through four examples including predator-prey, infectious disease transmission, and immune response models. We show that CaliPro works well for both deterministic, continuous model structures as well as stochastic, discrete models and illustrate that CaliPro can work across diverse calibration goals. CONCLUSIONS: We present CaliPro, a new method for calibrating complex biological models to a range of experimental outcomes. In addition to expediting calibration, CaliPro may be useful in already calibrated parameter spaces to target and isolate specific model behavior for further analysis. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s12195-020-00650-z) contains supplementary material, which is available to authorized users.

Tuberculosis (TB), caused by Mycobacterium tuberculosis, continues to be a major global health problem. Lung granulomas are organized structures of host immune cells that function to contain the bacteria. Cytokine expression is a critical component of the protective immune response, but inappropriate cytokine expression can exacerbate TB. Although the importance of proinflammatory cytokines in controlling M. tuberculosis infection has been established, the effects of anti-inflammatory cytokines, such as IL-10, in TB are less well understood. To investigate the role of IL-10, we used an Ab to neutralize IL-10 in cynomolgus macaques during M. tuberculosis infection. Anti–IL-10–treated nonhuman primates had similar overall disease outcomes compared with untreated control nonhuman primates, but there were immunological changes in granulomas and lymph nodes from anti–IL-10–treated animals. There was less thoracic inflammation and increased cytokine production in lung granulomas and lymph nodes from IL-10–neutralized animals at 3–4 wk postinfection compared with control animals. At 8 wk postinfection, lung granulomas from IL-10–neutralized animals had reduced cytokine production but increased fibrosis relative to control animals. Although these immunological changes did not affect the overall disease burden during the first 8 wk of infection, we paired computational modeling to explore late infection dynamics. Our findings support that early changes occurring in the absence of IL-10 may lead to better bacterial control later during infection. These unique datasets provide insight into the contribution of IL-10 to the immunological balance necessary for granulomas to control bacterial burden and disease pathology in M. tuberculosis infection.

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.

SUMMARYMycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), is a global health concern, yearly resulting in 10 million new cases of active TB. Immunologic investigation of lung granulomas is essential for understanding host control of bacterial replication. We identified and compared the pathological, cellular, and functional differences in granulomas at 4, 12, and 20 weeks post-infection in Chinese cynomolgus macaques. Original granulomas differed in transcription factor expression within adaptive lymphocytes with those at 12 weeks showing higher frequencies of CD8+T-bet+ T cells, while increases in CD4+T-bet+ T cells were observed at 20 weeks post-infection. The appearance of T-bet+ adaptive T cells at 12 and 20 weeks was coincident with a reduction in bacterial burden, suggesting their critical role in Mtb control. This study highlights the evolution of T cell responses within lung granulomas, suggesting that vaccines promoting the development and migration of T-bet+ T cells would enhance mycobacterial control.

According to the World Health Organization, tuberculosis (TB) is the leading cause of death from infectious disease worldwide (WHO, 2017). While there is no effective vaccine against adult pulmonary TB, more than a dozen vaccine candidates are in the clinical trial pipeline. These include both pre-exposure vaccines to prevent initial infections and post-exposure vaccines to prevent reactivation of latent disease. Many epidemiological models have been used to study TB, but most have not included a continuous age structure and the possibility of both pre- and post-exposure vaccination. Incorporating age-dependent death rates, disease properties, and social contact data allows for more realistic modeling of disease spread. We propose a continuous age-structured model for the epidemiology of tuberculosis with pre- and post-exposure vaccination. We use uncertainty and sensitivity analysis to make predictions about the efficacy of different vaccination strategies in a non-endemic setting (United States) and an endemic setting (Cambodia). In particular, we determine optimal age groups to target for pre-exposure and post-exposure vaccination in both settings. We find that the optimal age groups tend to be younger for Cambodia than for the US, and that post-exposure vaccination has a significantly larger effect than pre-exposure vaccination in the US.

Neutrophil infiltration into tuberculous granulomas is often associated with higher bacteria loads and severe disease but the basis for this relationship is not well understood. To better elucidate the connection between neutrophils and pathology in primate systems, we paired data from experimental studies with our next generation computational model GranSim to identify neutrophil-related factors, including neutrophil recruitment, lifespan, and intracellular bacteria numbers, that drive granuloma-level outcomes. We predict mechanisms underlying spatial organization of neutrophils within granulomas and identify how neutrophils contribute to granuloma dissemination. We also performed virtual deletion and depletion of neutrophils within granulomas and found that neutrophils play a nuanced role in determining granuloma outcome, promoting uncontrolled bacterial growth in some and working to contain bacterial growth in others. Here, we present three key results: We show that neutrophils can facilitate local dissemination of granulomas and thereby enable the spread of infection. We suggest that neutrophils influence CFU burden during both innate and adaptive immune responses, implying that they may be targets for therapeutic interventions during later stages of infection. Further, through the use of uncertainty and sensitivity analyses, we predict which neutrophil processes drive granuloma severity and structure.

Tuberculosis (TB) is the deadliest infectious disease worldwide. The design of new treatments for TB is hindered by the large number of candidate drugs, drug combinations, dosing choices, and complex pharmaco-kinetics/dynamics (PK/PD). Here we study the interplay of these factors in designing combination therapies by linking a machine-learning model, INDIGO-MTB, which predicts in vitro drug interactions using drug transcriptomics, with a multi-scale model of drug PK/PD and pathogen-immune interactions called GranSim. We calculate an in vivo drug interaction score (iDIS) from dynamics of drug diffusion, spatial distribution, and activity within lesions against various pathogen sub-populations. The iDIS of drug regimens evaluated against non-replicating bacteria significantly correlates with efficacy metrics from clinical trials. Our approach identifies mechanisms that can amplify synergistic or mitigate antagonistic drug interactions in vivo by modulating the relative distribution of drugs. Our mechanistic framework enables efficient evaluation of in vivo drug interactions and optimization of combination therapies.

Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis (TB), kills more individuals worldwide per year than any other infectious agent. As the hallmark of TB, lung granulomas are complex structures composed of immune cells that interact with and surround bacteria, infected cells, and a necrotic core. This interaction leads to diverse granuloma outcomes across time, ranging from bacterial sterilization to uncontrolled bacterial growth, as well as diverse spatial structures. At this time, there are no systematic quantitative methods to classify the formation, function, and spatial characteristics of granulomas. This type of analysis would enable better understanding and prediction of granuloma behaviors that have known associations with poor clinical outcomes for TB patients. Herein, we develop a temporal and spatial analysis framework for TB granulomas using a systems biology approach combining in silico granuloma modeling, geographic information systems, topological data analysis, and machine learning. We apply this framework to simulated granulomas to understand temporal granuloma dynamics, quantify granuloma spatial structure, and predict the relationship between granuloma structure and bacterial growth. As a proof-of-concept, we apply our in silico predictions to in vivo derived data to test our framework for future applications and as a personalized medicine intervention.

Highlights • Casual and workplace contacts are main drivers of disease spread. • The second wave of infections could be significantly larger than the first wave. • Reopening speed affects timing but not magnitude of second peak. • Reduced casual contacts decrease and delay the second wave of infections., The COVID-19 pandemic has highlighted the patchwork nature of disease epidemics, with infection spread dynamics varying wildly across countries and across states within the US. To explore this issue, we study and predict the spread of COVID-19 in Washtenaw County, MI, which is home to University of Michigan and Eastern Michigan University, and in close proximity to Detroit, MI, a major epicenter of the epidemic in Michigan. We apply a discrete and stochastic network-based modeling framework allowing us to track every individual in the county. In this framework, we construct contact networks based on synthetic population datasets specific for Washtenaw County that are derived from US Census datasets. We assign individuals to households, workplaces, schools, and group quarters (such as prisons or long term care facilities). In addition, we assign casual contacts to each individual at random. Using this framework, we explicitly simulate Michigan-specific government-mandated workplace and school closures as well as social distancing measures. We perform sensitivity analyses to identify key model parameters and mechanisms contributing to the observed disease burden in the three months following the first observed cases of COVID-19 in Michigan. We then consider several scenarios for relaxing restrictions and reopening workplaces to predict what actions would be most prudent. In particular, we consider the effects of 1) different timings for reopening, and 2) different levels of workplace vs. casual contact re-engagement. We find that delaying reopening does not reduce the magnitude of the second peak of cases, but only delays it. Reducing levels of casual contact, on the other hand, both delays and lowers the second peak. Through simulations and sensitivity analyses, we explore mechanisms driving the magnitude and timing of a second wave of infections upon re-opening. We find that the most significant factors are workplace and casual contacts and protective measures taken by infected individuals who have sought care. This model can be adapted to other US counties using synthetic population databases and data specific to those regions.

Current methods to optimize vaccine dose are purely empirically based, whereas in the drug development field, dosing determinations use far more advanced quantitative methodology to accelerate decision-making. Applying these established methods in the field of vaccine development may reduce the currently large clinical trial sample sizes, long time frames, high costs, and ultimately have a better potential to save lives. We propose the field of immunostimulation/immunodynamic (IS/ID) modelling, which aims to translate mathematical frameworks used for drug dosing towards optimizing vaccine dose decision-making. Analogous to PK/PD modelling, IS/ID modelling approaches apply mathematical models to describe the underlying mechanisms by which the immune response is stimulated by vaccination (IS) and the resulting measured immune response dynamics (ID). To move IS/ID modelling forward, existing datasets and further data on vaccine allometry and dose-dependent dynamics need to be generated and collate, requiring a collaborative environment with input from academia, industry, regulators, governmental and non-governmental agencies to share modelling expertise, and connect modellers to vaccine data., 11 pages, 1 Figure

The architectural complexity of the lung is crucial to its ability to function as an organ of gas exchange; the branching tree structure of the airways transforms the tracheal cross-section of only a few square centimeters to a blood-gas barrier with a surface area of tens of square meters and a thickness on the order of a micron or less. Connective tissue comprised largely of collagen and elastic fibers provides structural integrity for this intricate and delicate system. Homeostatic maintenance of this connective tissue, via a balance between catabolic and anabolic enzyme-driven processes, is crucial to life. Accordingly, when homeostasis is disrupted by the excessive production of connective tissue, lung function deteriorates rapidly with grave consequences leading to chronic lung conditions such as pulmonary fibrosis. Understanding how pulmonary fibrosis develops and alters the link between lung structure and function is crucial for diagnosis, prognosis, and therapy. Further information gained could help elaborate how the healing process breaks down leading to chronic disease. Our understanding of fibrotic disease is greatly aided by the intersection of wet lab studies and mathematical and computational modeling. In the present review we will discuss how multi-scale modeling has facilitated our understanding of pulmonary fibrotic disease as well as identified opportunities that remain open and have produced techniques that can be incorporated into this field by borrowing approaches from multi-scale models of fibrosis beyond the lung.

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.

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.

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), continues to be a major global health problem. Lung granulomas are organized structures of host immune cells that function to contain the bacteria. Cytokine expression is a critical component of the protective immune response but inappropriate cytokine expression can exacerbate TB. While the importance of pro-inflammatory cytokines in controlling Mtb infection has been established, the effects of anti-inflammatory cytokines, such as IL-10, in TB are less well understood. To investigate the role of IL-10, we used an antibody to neutralize IL-10 in cynomolgus macaques during Mtb infection. Anti-IL-10 treated NHPs had similar overall disease outcomes compared to untreated control NHPs, but there were immunological changes in granulomas and lymph nodes from anti-IL-10 treated animals. There was less thoracic inflammation and increased cytokine production in lung granulomas and lymph nodes from IL-10 neutralized animals at 3-4 weeks post-infection compared to control animals. At 8 weeks post-infection, lung granulomas from IL-10 neutralized animals had reduced cytokine production but increased fibrosis relative to control animals. While these immunological changes did not affect the overall disease burden during the first 8 weeks of infection, we paired computational modeling to explore late infection dynamics. Our findings support that early changes occurring in the absence of IL-10 may lead to better bacterial control later during infection. These unique datasets provide insight into the contribution of IL-10 to the immunological balance necessary for granulomas to control bacterial burden and disease pathology in Mtb infection.

Fluoroquinolones represent the pillar of multidrug-resistant tuberculosis (MDR-TB) treatment, with moxifloxacin, levofloxacin, or gatifloxacin being prescribed to MDR-TB patients. Recently, several clinical trials of “universal” drug regimens, aiming to treat drug-susceptible and drug-resistant TB, have included a fluoroquinolone., Fluoroquinolones represent the pillar of multidrug-resistant tuberculosis (MDR-TB) treatment, with moxifloxacin, levofloxacin, or gatifloxacin being prescribed to MDR-TB patients. Recently, several clinical trials of “universal” drug regimens, aiming to treat drug-susceptible and drug-resistant TB, have included a fluoroquinolone. In the absence of clinical data comparing their side-by-side efficacies in controlled MDR-TB trials, a pharmacological rationale is needed to guide the selection of the most efficacious fluoroquinolone. The present studies were designed to test the hypothesis that fluoroquinolone concentrations (pharmacokinetics) and activity (pharmacodynamics) at the site of infection are better predictors of efficacy than the plasma concentrations and potency measured in standard growth inhibition assays and are better suited to determinations of whether one of the fluoroquinolones outperforms the others in rabbits with active TB. We first measured the penetration of these fluoroquinolones in lung lesion compartments, and their potency against bacterial populations that reside in each compartment, to compute lesion-centric pharmacokinetic-pharmacodynamic (PK/PD) parameters. PK modeling methods were used to quantify drug penetration from plasma to tissues at human-equivalent doses. On the basis of these metrics, moxifloxacin emerged with a clear advantage, whereas plasma-based PK/PD favored levofloxacin (the ranges of the plasma AUC/MIC ratio [i.e., the area under the concentration-time curve over 24 h in the steady state divided by the MIC] are 46 to 86 for moxifloxacin and 74 to 258 for levofloxacin). A comparative efficacy trial in the rabbit model of active TB demonstrated the superiority of moxifloxacin in reducing bacterial burden at the lesion level and in sterilizing cellular and necrotic lesions. Collectively, these results show that PK/PD data obtained at the site of infection represent an adequate predictor of drug efficacy against TB and constitute the baseline required to explore synergies, antagonism, and drug-drug interactions in fluoroquinolone-containing regimens.

Computational and mathematical models in biology rely heavily on the parameters that characterize them. However, robust estimates for their values are typically elusive and thus a large parameter space becomes necessary for model study, particularly to make translationally impactful predictions. Sampling schemes exploring parameter spaces for models are used for a variety of purposes in systems biology, including model calibration and sensitivity analysis. Typically, random sampling is used; however, when models have a high number of unknown parameters or the models are highly complex, computational cost becomes an important factor. This issue can be reduced through the use of efficient sampling schemes such as Latin hypercube sampling (LHS) and Sobol sequences. In this work, we compare and contrast three sampling schemes – random sampling, LHS, and Sobol sequences – for the purposes of performing both parameter sensitivity analysis and model calibration. In addition, we apply these analyses to different types of computational and mathematical models of varying complexity: a simple ODE model, a complex ODE model, and an agent-based model. In general, the sampling scheme had little effect when used for calibration efforts, but when applied to sensitivity analyses, Sobol sequences exhibited faster convergence. While the observed benefit to convergence is relatively small, Sobol sequences are computationally less expensive to compute than LHS samples and also have the benefit of being deterministic, which allows for better reproducibility of results.

Within the first 2&ndash, 3 months of a Mycobacterium tuberculosis (Mtb) infection, 2&ndash, 4 mm spherical structures called granulomas develop in the lungs of the infected hosts. These are the hallmark of tuberculosis (TB) infection in humans and non-human primates. A cascade of immunological events occurs in the first 3 months of granuloma formation that likely shapes the outcome of the infection. Understanding the main mechanisms driving granuloma development and function is key to generating treatments and vaccines. In vitro, in vivo, and in silico studies have been performed in the past decades to address the complexity of granuloma dynamics. This study builds on our previous 2D spatio-temporal hybrid computational model of granuloma formation in TB (GranSim) and presents for the first time a more realistic 3D implementation. We use uncertainty and sensitivity analysis techniques to calibrate the new 3D resolution to non-human primate (NHP) experimental data on bacterial levels per granuloma during the first 100 days post infection. Due to the large computational cost associated with running a 3D agent-based model, our major goal is to assess to what extent 2D and 3D simulations differ in predictions for TB granulomas and what can be learned in the context of 3D that is missed in 2D. Our findings suggest that in terms of major mechanisms driving bacterial burden, 2D and 3D models return very similar results. For example, Mtb growth rates and molecular regulation mechanisms are very important both in 2D and 3D, as are cellular movement and modulation of cell recruitment. The main difference we found was that the 3D model is less affected by crowding when cellular recruitment and movement of cells are increased. Overall, we conclude that the use of a 2D resolution in GranSim is warranted when large scale pilot runs are to be performed and if the goal is to determine major mechanisms driving infection outcome (e.g., bacterial load). To comprehensively compare the roles of model dimensionality, further tests and experimental data will be needed to expand our conclusions to molecular scale dynamics and multi-scale resolutions.

Comprehensive assessment of the effectiveness of contact investigations for tuberculosis (TB) control is still lacking. In this study, we use a computational model, calibrated against notification data from Arkansas during the period 2001-2011, that reproduces independent data on key features of TB transmission and epidemiology. The model estimates that the Arkansas contact investigations program has avoided 18.6% (12.1-25.9%) of TB cases and 23.7% (16.4-30.6%) of TB deaths that would have occurred during 2001-2014 if passive diagnosis alone were implemented. If contacts of sputum smear-negative cases had not been included in the program, the percentage reduction would have been remarkably lower. In addition, we predict that achieving national targets for performance indicators of contact investigation programs has strong potential to further reduce TB transmission and burden. However, contact investigations are expected to have limited effectiveness on avoiding reactivation cases of latent infections over the next 60 years.

Conventional oral therapies for the treatment of tuberculosis are limited by poor antibiotic distribution in granulomas, which contributes to lengthy treatment regimens and inadequate bacterial sterilization. Inhaled formulations are a promising strategy to increase antibiotic efficacy and reduce dose frequency. We develop a multiscale computational approach that accounts for simultaneous dynamics of a lung granuloma, carrier release kinetics, pharmacokinetics, and pharmacodynamics. Using this computational platform, we predict that a rationally designed inhaled formulation of isoniazid given at a significantly reduced dose frequency has better sterilizing capabilities and reduced toxicity than the current oral regimen. Furthermore, we predict that inhaled formulations of rifampicin require unrealistic carrier antibiotic loadings that lead to early toxicity concerns. Lastly, we predict that targeting carriers to macrophages has limited effects on treatment efficacy. Our platform can be extended to account for additional antibiotics and provides a new tool for rapidly prototyping the efficacy of inhaled formulations.

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.

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 M 1 polarization early during infection and to improve outcome.

Tuberculosis (TB), one of the most common infectious diseases, requires treatment with multiple antibiotics taken over at least 6 months. This long treatment often results in poor patient-adherence, which can lead to the emergence of multi-drug resistant TB. New antibiotic treatment strategies are sorely needed. New antibiotics are being developed or repurposed to treat TB, but as there are numerous potential antibiotics, dosing sizes and potential schedules, the regimen design space for new treatments is too large to search exhaustively. Here we propose a method that combines an agent-based multi-scale model capturing TB granuloma formation with algorithms for mathematical optimization to identify optimal TB treatment regimens. We define two different single-antibiotic treatments to compare the efficiency and accuracy in predicting optimal treatment regimens of two optimization algorithms: genetic algorithms (GA) and surrogate-assisted optimization through radial basis function (RBF) networks. We also illustrate the use of RBF networks to optimize double-antibiotic treatments. We found that while GAs can locate optimal treatment regimens more accurately, RBF networks provide a more practical strategy to TB treatment optimization with fewer simulations, and successfully estimated optimal double-antibiotic treatment regimens. Our results indicate surrogate-assisted optimization can locate optimal TB treatment regimens from a larger set of antibiotics, doses and schedules, and could be applied to solve optimization problems in other areas of research using systems biology approaches. Our findings have important implications for the treatment of diseases like TB that have lengthy protocols or for any disease that requires multiple drugs.

Tuberculosis (TB) is a world-wide health problem with approximately 2 billion people infected with Mycobacterium tuberculosis (Mtb, the causative bacterium of TB). The pathologic hallmark of Mtb infection in humans and Non-Human Primates (NHPs) is the formation of spherical structures, primarily in lungs, called granulomas. Infection occurs after inhalation of bacteria into lungs, where resident antigen-presenting cells (APCs), take up bacteria and initiate the immune response to Mtb infection. APCs traffic from the site of infection (lung) to lung-draining lymph nodes (LNs) where they prime T cells to recognize Mtb. These T cells, circulating back through blood, migrate back to lungs to perform their immune effector functions. We have previously developed a hybrid agent-based model (ABM, labeled GranSim) describing in silico immune cell, bacterial (Mtb) and molecular behaviors during tuberculosis infection and recently linked that model to operate across three physiological compartments: lung (infection site where granulomas form), lung draining lymph node (LN, site of generation of adaptive immunity) and blood (a measurable compartment). Granuloma formation and function is captured by a spatio-temporal model (i.e., ABM), while LN and blood compartments represent temporal dynamics of the whole body in response to infection and are captured with ordinary differential equations (ODEs). In order to have a more mechanistic representation of APC trafficking from the lung to the lymph node, and to better capture antigen presentation in a draining LN, this current study incorporates the role of dendritic cells (DCs) in a computational fashion into GranSim. Results: The model was calibrated using experimental data from the lungs and blood of NHPs. The addition of DCs allowed us to investigate in greater detail mechanisms of recruitment, trafficking and antigen presentation and their role in tuberculosis infection. Conclusion: The main conclusion of this study is that early events after Mtb infection are critical to establishing a timely and effective response. Manipulating CD8+ and CD4+ T cell proliferation rates, as well as DC migration early on during infection can determine the difference between bacterial clearance vs. uncontrolled bacterial growth and dissemination.

Pulmonary fibrosis is pathologic remodeling of lung tissue that can result in difficulty breathing, reduced quality of life, and a poor prognosis for patients. Fibrosis occurs as a result of insult to lung tissue, though mechanisms of this response are not well-characterized. The disease is driven in part by dysregulation of fibroblast proliferation and differentiation into myofibroblast cells, as well as pro-fibrotic mediator-driven epithelial cell apoptosis. The most well-characterized pro-fibrotic mediator associated with pulmonary fibrosis is TGF-β1. Excessive synthesis of, and sensitivity to, pro-fibrotic mediators as well as insufficient production of and sensitivity to anti-fibrotic mediators has been credited with enabling fibroblast accumulation. Available treatments neither halt nor reverse lung damage. In this study we have two aims: to identify molecular and cellular scale mechanisms driving fibroblast proliferation and differentiation as well as epithelial cell survival in the context of fibrosis, and to predict therapeutic targets and strategies. We combine in vitro studies with a multi-scale hybrid agent-based computational model that describes fibroblasts and epithelial cells in co-culture. Within this model TGF-β1 represents a pro-fibrotic mediator and we include detailed dynamics of TGF-β1 receptor ligand signaling in fibroblasts. PGE2 represents an anti-fibrotic mediator. Using uncertainty and sensitivity analysis we identify TGF-β1 synthesis, TGF-β1 activation, and PGE2 synthesis among the key mechanisms contributing to fibrotic outcomes. We further demonstrate that intervention strategies combining potential therapeutics targeting both fibroblast regulation and epithelial cell survival can promote healthy tissue repair better than individual strategies. Combinations of existing drugs and compounds may provide significant improvements to the current standard of care for pulmonary fibrosis. Thus, a two-hit therapeutic intervention strategy may prove necessary to halt and reverse disease dynamics.

Drug resistant tuberculosis is increasing world-wide. Resistance against isoniazid (INH), rifampicin (RIF), or both (multi-drug resistant TB, MDR-TB) is of particular concern, since INH and RIF form part of the standard regimen for TB disease. While it is known that suboptimal treatment can lead to resistance, it remains unclear how host immune responses and antibiotic dynamics within granulomas (sites of infection) affect emergence and selection of drug-resistant bacteria. We take a systems pharmacology approach to explore resistance dynamics within granulomas. We integrate spatio-temporal host immunity, INH and RIF dynamics, and bacterial dynamics (including fitness costs and compensatory mutations) in a computational framework. We simulate resistance emergence in the absence of treatment, as well as resistance selection during INH and/or RIF treatment. There are four main findings. First, in the absence of treatment, the percentage of granulomas containing resistant bacteria mirrors the non-monotonic bacterial dynamics within granulomas. Second, drug-resistant bacteria are less frequently found in non-replicating states in caseum, compared to drug-sensitive bacteria. Third, due to a steeper dose response curve and faster plasma clearance of INH compared to RIF, INH-resistant bacteria have a stronger influence on treatment outcomes than RIF-resistant bacteria. Finally, under combination therapy with INH and RIF, few MDR bacteria are able to significantly affect treatment outcomes. Overall, our approach allows drug-specific prediction of drug resistance emergence and selection in the complex granuloma context. Since our predictions are based on pre-clinical data, our approach can be implemented relatively early in the treatment development process, thereby enabling pro-active rather than reactive responses to emerging drug resistance for new drugs. Furthermore, this quantitative and drug-specific approach can help identify drug-specific properties that influence resistance and use this information to design treatment regimens that minimize resistance selection and expand the useful life-span of new antibiotics.

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.

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.