Search

The risk of active tuberculosis disease is 15-21 times higher in those coinfected with human immunodeficiency virus-1 (HIV) compared to tuberculosis alone, and tuberculosis is the leading cause of death in HIV+ individuals. Mechanisms driving synergy between Mycobacterium tuberculosis (Mtb) and HIV during coinfection include: disruption of cytokine balances, impairment of innate and adaptive immune cell functionality, and Mtb-induced increase in HIV viral loads. Tuberculosis granulomas are the interface of host-pathogen interactions. Thus, granuloma-based research elucidating the role and relative impact of coinfection mechanisms within Mtb granulomas could inform cohesive treatments that target both pathogens simultaneously. We review known interactions between Mtb and HIV, and discuss how the structure, function and development of the granuloma microenvironment create a positive feedback loop favoring pathogen expansion and interaction. We also identify key outstanding questions and highlight how coupling computational modeling with in vitro and in vivo efforts could accelerate Mtb-HIV coinfection discoveries.

Immunomodulatory drugs could contribute to a functional cure for Human Immunodeficiency Virus (HIV). Interleukin-15 (IL-15) promotes expansion and activation of CD8+ T cell and natural killer (NK) cell populations. In one study, an IL-15 superagonist, N-803, suppressed Simian Immunodeficiency Virus (SIV) in non-human primates (NHPs) who had received prior SIV vaccination. However, viral suppression attenuated with continued N-803 treatment, partially returning after long treatment interruption. While there is evidence of concurrent drug tolerance, immune regulation, and viral escape, the relative contributions of these mechanisms to the observed viral dynamics have not been quantified. Here, we utilize mathematical models of N-803 treatment in SIV-infected macaques to estimate contributions of these three key mechanisms to treatment outcomes: 1) drug tolerance, 2) immune regulation, and 3) viral escape. We calibrated our model to viral and lymphocyte responses from the above-mentioned NHP study. Our models track CD8+ T cell and NK cell populations with N-803-dependent proliferation and activation, as well as viral dynamics in response to these immune cell populations. We compared mathematical models with different combinations of the three key mechanisms based on Akaike Information Criterion and important qualitative features of the NHP data. Two minimal models were capable of reproducing the observed SIV response to N-803. In both models, immune regulation strongly reduced cytotoxic cell activation to enable viral rebound. Either long-term drug tolerance or viral escape (or some combination thereof) could account for changes to viral dynamics across long breaks in N-803 treatment. Theoretical explorations with the models showed that less-frequent N-803 dosing and concurrent immune regulation blockade (e.g. PD-L1 inhibition) may improve N-803 efficacy. However, N-803 may need to be combined with other immune therapies to countermand viral escape from the CD8+ T cell response. Our mechanistic model will inform such therapy design and guide future studies., Author summary Immune therapy may be a critical component in the functional cure for Human Immunodeficiency Virus (HIV). N-803 is an immunotherapeutic drug that activates antigen-specific CD8+ T cells of the immune system. These CD8+ T cells eliminate HIV-infected cells in order to limit the spread of infection in the body. In one study, N-803 reduced plasma viremia in macaques that were infected with Simian Immunodeficiency Virus, an analog of HIV. Here, we used mathematical models to analyze the data from this study to better understand the effects of N-803 therapy on the immune system. Our models indicated that inhibitory signals may be reversing the stimulatory effect of N-803. Results also suggested the possibilities that tolerance to N-803 could build up within the CD8+ T cells themselves and that the treatment may be selecting for virus strains that are not targeted by CD8+ T cells. Our models predict that N-803 therapy may be made more effective if the time between doses is increased or if inhibitory signals are blocked by an additional drug. Also, N-803 may need to be combined with other immune therapies to target virus that would otherwise evade CD8+ T cells.

Despite intensive research efforts, several fundamental disease processes for tuberculosis (TB) remain poorly understood. A central enigma is that host immunity is necessary to control disease yet promotes transmission by causing lung immunopathology. Our inability to distinguish these processes makes it challenging to design rational novel interventions. Elucidating basic immune mechanisms likely requires both in vivo and in vitro analyses, since Mycobacterium tuberculosis is a highly specialized human pathogen. The classic immune response is the TB granuloma organized in three dimensions within extracellular matrix. Several groups are developing cell culture granuloma models. In January 2018, NIAID convened a workshop, entitled "3-D Human in vitro TB Granuloma Model" to advance the field. Here, we summarize the arguments for developing advanced TB cell culture models and critically review those currently available. We discuss how integrating complementary approaches, specifically organoids and mathematical modeling, can maximize progress, and conclude by discussing future challenges and opportunities. Funding Agencies|Division of AIDS, National Institute of Allergy and Infectious Diseases, National Institutes of Health (NIH), Department of Health and Human Services [HHSN272201600001G]; Medical Research Council [MR/P023754/1]; NIH [1DP2HL127720-01, R33AI102239, UL1TR001070]; National Science Foundation [CMMI: 15-48571]; Swedish Research Council [2014-02289]

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.

The model herein aims to explore the dynamics of the spread of tuberculosis (TB) in an informal settlement or township. The population is divided into households of various sizes and also based on commuting status. The model dynamics distinguishes between three distinct social patterns: the exposure of commuters during travel, random diurnal interaction and familial exposure at night. Following the general SLIR models, the population is further segmented into susceptible (S), exposed/latently infected (L), active/infectious (I), and recovered (R) individuals. During the daytime, commuters travel on public transport, while non-commuters randomly interact in the community to mimic chance encounters with infectious persons. At night, each family interacts and sleeps together in the home. The risk of exposure to TB is based on the proximity, duration, and frequency of encounters with infectious persons. The model is applied to a hypothetical population to explore the effects of different intervention strategies including vaccination, wearing of masks or scarves during the commute, prophylactic treatment of latent infections and more effective case-finding and treatment. The most important findings of the model are: (1) members of larger families are responsible for more disease transmissions than those from smaller families, (2) daily commutes on public transport provide ideal conditions for transmission of the disease, (3) improved diagnosis and treatment has the greatest impact on the spread of the disease, and (4) detecting TB at the first clinic visit, when patients are still smear negative, is key.

Approximately one third of the world's population is infected with Mycobacterium tuberculosis. Limited information about how the immune system fights M. tuberculosis and what constitutes protection from the bacteria impact our ability to develop effective therapies for tuberculosis. We present an in vivo systems biology approach that integrates data from multiple model systems and over multiple length and time scales into a comprehensive multi-scale and multi-compartment view of the in vivo immune response to M. tuberculosis. We describe computational models that can be used to study (a) immunomodulation with the cytokines tumor necrosis factor and interleukin 10, (b) oral and inhaled antibiotics, and (c) the effect of vaccination.

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., Author summary Tuberculosis (TB) is caused by infection with the bacterium Mycobacterium tuberculosis (Mtb) and kills 1.5 million people each year. TB requires at least 6 months of treatment with up to four drugs, and is characterized by formation of granulomas in patient lungs. Granulomas are spherical collections of host cells and bacteria. Fluoroquinolones (FQs) are a class of drug that could help shorten TB treatment. Three FQs that are used to treat TB are: moxifloxacin (MXF), levofloxacin (LVX) or gatifloxacin (GFX). To date, it is unclear if one FQ is better than the others at treating TB, in part because little is known about how these drugs distribute and work inside the lung granulomas. We use computer simulations of Mtb infection and FQ treatment within granulomas to predict which FQ is better and why. Our computer model is calibrated to multiple experimental data sets. We compare the three FQs by multiple metrics, and predict that MXF is better than LVX and GFX because it kills bacteria more quickly, and it works better when patients miss doses. However, all three FQs are unable to kill a part of the bacterial population living in the center of granulomas. Our results can now inform future experimental studies.

Ethnopharmacological relevance The emergence of multidrug-resistant strains of Mycobacterium tuberculosis underscores the need for continuous development of new and efficient methods to determine the susceptibility of isolates of Mycobacterium tuberculosis in the search for novel antimycobacterial agents. Natural products constitute an important source of new drugs, and design and implementation of antimycobacterial susceptibility testing methods are necessary to evaluate the different extracts and compounds. In this study we have explored the antimycobacterial properties of 50 ethanolic extracts from different parts of 46 selected medicinal plants traditionally used in Sudan to treat infectious diseases. Materials and methods Plants were harvested and ethanolic extracts were prepared. For selected extracts, fractionation with hydrophilic and hydrophobic solvents was undertaken. A luminometry-based assay was used for determination of mycobacterial growth in broth cultures and inside primary human macrophages in the presence or absence of plant extracts and fractions of extracts. Cytotoxicity was also assessed for active fractions of plant extracts. Results Of the tested extracts, three exhibited a significant inhibitory effect on an avirulent strain of Mycobacterium tubercluosis (H37Ra) at the initial screening doses (125 and 6.25 µg/ml). These were bark and leaf extracts of Khaya senegalensis and the leaf extract of Rosmarinus officinalis L. Further fractions of these plant extracts were prepared with n-hexane, chloroform, ethyl acetate, n-butanol, ethanol and water, and the activity of these extracts was retained in hydrophobic fractions. Cytotoxicity assays revealed that the chloroform fraction of Khaya senegalensis bark was non-toxic to human monocyte-derived macrophages and other cell types at the concentrations used and hence, further analysis, including assessment of IC50 and intracellular activity was done with this fraction. Conclusion These results encourage further investigations to identify the active compound(s) within the chloroform fraction of Khaya senegalensis bark.

The standard treatment of tuberculosis (TB) takes six to nine months to complete and this lengthy therapy contributes to the emergence of drug-resistant TB. TB is caused by Mycobacterium tuberculosis (Mtb) and the ability of this bacterium to switch to a dormant phenotype has been suggested to be responsible for the slow clearance during treatment. A recent study showed that the replication rate of a non-virulent mycobacterium, Mycobacterium smegmatis, did not correlate with antibiotic susceptibility. However, the question whether this observation also holds true for Mtb remains unanswered. Here, in order to mimic physiological conditions of TB infection, we established a protocol based on long-term infection of primary human macrophages, featuring Mtb replicating at different rates inside the cells. During conditions that restricted Mtb replication, the bacterial phenotype was associated with reduced acid-fastness. However, these phenotypically altered bacteria were as sensitive to isoniazid, pyrazinamide and ethambutol as intracellularly replicating Mtb. In support of the recent findings with M. smegmatis, we conclude that replication rates of Mtb do not correlate with antibiotic tolerance. Funding Agencies|Bill and Melinda Gates Foundation; Swedish Research Council [2009-3821, 2012-3349]; Swedish International Development Cooperation Agency; Swedish Heart-Lung Foundation; King Oscar II Foundation; Carl Trygger Foundation; Clas Groschinsky Foundation

There is a large body of literature describing molecular level interactions between Mycobacterium tuberculosis (Mtb) and macrophages. Macrophages initiate a range of anti-bacterial mechanisms in response to infection, and Mtb is capable of surviving and circumventing many of these responses. We apply a computational approach to ask: what are the effects on the cellular level of these opposing interactions? The model considers the interplay between bacterial killing and the pathogen's interference with macrophage function. The results reveal an oscillating balance between host and pathogen, but the balance is transient and varies in length, indicating that stochasticity in the bacterial population or host response could contribute to the diverse incubation periods observed in exposed individuals. The model captures host and strain variation and gives new insight into host-pathogen compatibility and co-evolution.

An epidemiological model is presented that considers five possible states of a population: susceptible (S), exposed (W), infectious (Y), in treatment (Z) and recovered (R). In certain instances transition rates (from one state to another) depend on the time spent in the state; therefore the states W, Y and Z depend on time and length of stay in that state - similar to age-structured models. The model is particularly amenable to describe delays of exposed persons to become infectious and re-infection of exposed persons. Other transitions that depend on state time include the case finding and diagnosis, increased death rate and treatment interruption. The mathematical model comprises of a set of partial differential and ordinary differential equations. Non-steady state solutions are first presented, followed by a bifurcation study of the stationary states.

It is widely accepted that β-lactam antimicrobials cause cell death through a mechanism that interferes with cell wall synthesis. Later studies have also revealed that β-lactams modify the autolysis function (the natural process of self-exfoliation of the cell wall) of cells. The dynamic equilibrium between growth and autolysis is perturbed by the presence of the antimicrobial. Studies with Staphylococcus aureus to determine the minimum inhibitory concentration (MIC) have revealed complex responses to methicillin exposure. The organism exhibits four qualitatively different responses: homogeneous sensitivity, homogeneous resistance, heterogeneous resistance and the so-called ‘Eagle-effect’. A mathematical model is presented that links antimicrobial action on the molecular level with the overall response of the cell population to antimicrobial exposure. The cell population is modeled as a probability density function F ( x , t ) that depends on cell wall thickness x and time t . The function F ( x,t ) is the solution to a Fokker–Planck equation. The fixed point solutions are perturbed by the antimicrobial load and the advection of F ( x , t ) depends on the rates of cell wall synthesis, autolysis and the antimicrobial concentration. Solutions of the Fokker–Planck model are presented for all four qualitative responses of S. aureus to methicillin exposure.

The tri-frame model gives mathematical expression to the transcription and translation processes, and considers all three reading frames (RFs). RNA polymerases transcribe DNA in single nucleotide increments, but ribosomes translate mRNA in pairings of three (triplets or codons). The set of triplets in the mRNA, starting with the initiation codon (usually AUG) defines the open reading frame (ORF). Since ribosomes do not always translocate three nucleotide positions, two additional RFs are accessible. The −1 RF and the +1 RF are triplet pairings of the mRNA, which are accessed by shifting one nucleotide position in the 5′ and 3′ directions, respectively. Transcription is modeled as a linear operator that maps the initial codons in all three frames into other codon sets to account for possible transcriptional errors. Translational errors (missense errors) originate from misacylation of tRNAs and misreading of aa-tRNAs by the ribosome. Translation is modeled as a linear mapping from codons into aa-tRNA species, which includes misreading errors. A final transformation from aa-tRNA species into amino acids provides the probability distributions of possible amino acids into which the codons in all three frames could be translated. An important element of the tri-frame model is the ribosomal occupancy probability. It is a vector in R 3 that gives the probability to find the ribosome in the ORF, −1 or +1 RF at each codon position. The sequence of vectors, from the first to the final codon position, gives a history of ribosome frameshifting. The model is powerful: it provides explicit expressions for (1) yield of error-free protein, (2) fraction of prematurely terminated polypeptides, (3) number of transcription errors, (4) number of translation errors and (5) mutations due to frameshifting. The theory is demonstrated for the three genes rpsU , dnaG and rpoD of Escherichia coli , which lie on the same operon, as well as for the prfB gene.