High rifampicin peak plasma concentrations accelerate the slow phase of bacterial decline in tuberculosis patients: Evidence for heteroresistance

BACKGROUND

Antibiotic treatments are often associated with a late slowdown in bacterial killing. This separates the killing of bacteria into at least two distinct phases: a quick phase followed by a slower phase, the latter of which is linked to treatment success. Current mechanistic explanations for the in vitro slowdown are either antibiotic persistence or heteroresistance. Persistence is defined as the switching back and forth between susceptible and non-susceptible states, while heteroresistance is defined as the coexistence of bacteria with heterogeneous susceptibilities. Both are also thought to cause a slowdown in the decline of bacterial populations in patients and therefore complicate and prolong antibiotic treatments. Reduced bacterial death rates over time are also observed within tuberculosis patients, yet the mechanistic reasons for this are unknown and therefore the strategies to mitigate them are also unknown.

METHODS AND FINDINGS

We analyse a dose ranging trial for rifampicin in tuberculosis patients and show that there is a slowdown in the decline of bacteria. We show that the late phase of bacterial killing depends more on the peak drug concentrations than the total drug exposure. We compare these to pharmacokinetic-pharmacodynamic models of rifampicin heteroresistance and persistence. We find that the observation on the slow phase's dependence on pharmacokinetic measures, specifically peak concentrations are only compatible with models of heteroresistance and incompatible with models of persistence. The quantitative agreement between heteroresistance models and observations is very good ([Formula: see text]). To corroborate the importance of the slowdown, we validate our results by estimating the time to sputum culture conversion and compare the results to a different dose ranging trial.

CONCLUSIONS

Our findings indicate that higher doses, specifically higher peak concentrations may be used to optimize rifampicin treatments by accelerating bacterial killing in the slow phase. It adds to the growing body of literature supporting higher rifampicin doses for shortening tuberculosis treatments.

Current Approaches of Building Mechanistic Pharmacodynamic Drug-Target Binding Models

Mechanistic pharmacodynamic models that incorporate the binding kinetics of drug-target interactions have several advantages in understanding target engagement and the efficacy of a drug dose. However, guidelines on how to build and interpret mechanistic pharmacodynamic drug-target binding models considering both biological and computational factors are still missing in the literature. In this chapter, current approaches of building mechanistic PD models and their advantages are discussed. We also present a methodology on how to select a suitable model considering both biological and computational perspectives, as well as summarize the challenges of current mechanistic PD models.

Perspectives for systems biology in the management of tuberculosis

Standardised management of tuberculosis may soon be replaced by individualised, precision medicine-guided therapies informed with knowledge provided by the field of systems biology. Systems biology is a rapidly expanding field of computational and mathematical analysis and modelling of complex biological systems that can provide insights into mechanisms underlying tuberculosis, identify novel biomarkers, and help to optimise prevention, diagnosis and treatment of disease. These advances are critically important in the context of the evolving epidemic of drug-resistant tuberculosis. Here, we review the available evidence on the role of systems biology approaches - human and mycobacterial genomics and transcriptomics, proteomics, lipidomics/metabolomics, immunophenotyping, systems pharmacology and gut microbiomes - in the management of tuberculosis including prediction of risk for disease progression, severity of mycobacterial virulence and drug resistance, adverse events, comorbidities, response to therapy and treatment outcomes. Application of the Grading of Recommendations, Assessment, Development and Evaluation (GRADE) approach demonstrated that at present most of the studies provide "very low" certainty of evidence for answering clinically relevant questions. Further studies in large prospective cohorts of patients, including randomised clinical trials, are necessary to assess the applicability of the findings in tuberculosis prevention and more efficient clinical management of patients.

RESTAMP - Rate estimates by sequence-tag analysis of microbial populations

Microbial division rates determine the speed of mutation accumulation and thus the emergence of antimicrobial resistance. Microbial death rates are affected by antibiotic action and the immune system. Therefore, measuring these rates has advanced our understanding of host-pathogen interactions and antibiotic action. Several methods based on marker-loss or few inheritable neutral markers exist that allow estimating microbial division and death rates, each of which has advantages and limitations. Technical bottlenecks, i.e., experimental sampling events, during the experiment can distort the rate estimates and are typically unaccounted for or require additional calibration experiments. In this work, we introduce RESTAMP (Rate Estimates by Sequence Tag Analysis of Microbial Populations) as a method for determining bacterial division and death rates. This method uses hundreds of fitness neutral sequence barcodes to measure the rates and account for experimental bottlenecks at the same time. We experimentally validate RESTAMP and compare it to established plasmid loss methods. We find that RESTAMP has a number of advantages over plasmid loss or previous marker based techniques. (i) It enables to correct the distortion of rate estimates by technical bottlenecks. (ii) Rate estimates are independent of the sequence tag distribution in the starting culture allowing the use of an arbitrary number of tags. (iii) It introduces a bottleneck sensitivity measure that can be used to maximize the accuracy of the experiment. RESTAMP allows studying microbial population dynamics with great resolution over a wide dynamic range and can thus advance our understanding of host-pathogen interactions or the mechanisms of antibiotic action.

All nonadherence is equal but is some more equal than others? Tuberculosis in the digital era

Adherence to treatment for tuberculosis (TB) has been a concern for many decades, resulting in the World Health Organization's recommendation of the direct observation of treatment in the 1990s. Recent advances in digital adherence technologies (DATs) have renewed discussion on how to best address nonadherence, as well as offering important information on dose-by-dose adherence patterns and their variability between countries and settings. Previous studies have largely focussed on percentage thresholds to delineate sufficient adherence, but this is misleading and limited, given the complex and dynamic nature of adherence over the treatment course. Instead, we apply a standardised taxonomy - as adopted by the international adherence community - to dose-by-dose medication-taking data, which divides missed doses into 1) late/noninitiation (starting treatment later than expected/not starting), 2) discontinuation (ending treatment early), and 3) suboptimal implementation (intermittent missed doses). Using this taxonomy, we can consider the implications of different forms of nonadherence for intervention and regimen design. For example, can treatment regimens be adapted to increase the "forgiveness" of common patterns of suboptimal implementation to protect against treatment failure and the development of drug resistance? Is it reasonable to treat all missed doses of treatment as equally problematic and equally common when deploying DATs? Can DAT data be used to indicate the patients that need enhanced levels of support during their treatment course? Critically, we pinpoint key areas where knowledge regarding treatment adherence is sparse and impeding scientific progress.

Multi-scale modeling of drug binding kinetics to predict drug efficacy

Optimizing drug therapies for any disease requires a solid understanding of pharmacokinetics (the drug concentration at a given time point in different body compartments) and pharmacodynamics (the effect a drug has at a given concentration). Mathematical models are frequently used to infer drug concentrations over time based on infrequent sampling and/or in inaccessible body compartments. Models are also used to translate drug action from in vitro to in vivo conditions or from animal models to human patients. Recently, mathematical models that incorporate drug-target binding and subsequent downstream responses have been shown to advance our understanding and increase predictive power of drug efficacy predictions. We here discuss current approaches of modeling drug binding kinetics that aim at improving model-based drug development in the future. This in turn might aid in reducing the large number of failed clinical trials.

Quantifying the impact of treatment history on plasmid-mediated resistance evolution in human gut microbiota

To understand how antibiotic use affects the risk of a resistant infection, we present a computational model of the population dynamics of gut microbiota including antibiotic resistance-conferring plasmids. We then describe how this model is parameterized based on published microbiota data. Finally, we investigate how treatment history affects the prevalence of resistance among opportunistic enterobacterial pathogens. We simulate treatment histories and identify which properties of prior antibiotic exposure are most influential in determining the prevalence of resistance. We find that resistance prevalence can be predicted by 3 properties, namely the total days of drug exposure, the duration of the drug-free period after last treatment, and the center of mass of the treatment pattern. Overall this work provides a framework for capturing the role of the microbiome in the selection of antibiotic resistance and highlights the role of treatment history for the prevalence of resistance.

Reaction Kinetic Models of Antibiotic Heteroresistance

Bacterial heteroresistance (i.e., the co-existence of several subpopulations with different antibiotic susceptibilities) can delay the clearance of bacteria even with long antibiotic exposure. Some proposed mechanisms have been successfully described with mathematical models of drug-target binding where the mechanism's downstream of drug-target binding are not explicitly modeled and subsumed in an empirical function, connecting target occupancy to antibiotic action. However, with current approaches it is difficult to model mechanisms that involve multi-step reactions that lead to bacterial killing. Here, we have a dual aim: first, to establish pharmacodynamic models that include multi-step reaction pathways, and second, to model heteroresistance and investigate which molecular heterogeneities can lead to delayed bacterial killing. We show that simulations based on Gillespie algorithms, which have been employed to model reaction kinetics for decades, can be useful tools to model antibiotic action via multi-step reactions. We highlight the strengths and weaknesses of current models and Gillespie simulations. Finally, we show that in our models, slight normally distributed variances in the rates of any event leading to bacterial death can (depending on parameter choices) lead to delayed bacterial killing (i.e., heteroresistance). This means that a slowly declining residual bacterial population due to heteroresistance is most likely the default scenario and should be taken into account when planning treatment length.

Estimating treatment prolongation for persistent infections

Treatment of infectious diseases is often long and requires patients to take drugs even after they have seemingly recovered. This is because of a phenomenon called persistence, which allows small fractions of the bacterial population to survive treatment despite being genetically susceptible. The surviving subpopulation is often below detection limit and therefore is empirically inaccessible but can cause treatment failure when treatment is terminated prematurely. Mathematical models could aid in predicting bacterial survival and thereby determine sufficient treatment length. However, the mechanisms of persistence are hotly debated, necessitating the development of multiple mechanistic models. Here we develop a generalized mathematical framework that can accommodate various persistence mechanisms from measurable heterogeneities in pathogen populations. It allows the estimation of the relative increase in treatment length necessary to eradicate persisters compared to the majority population. To simplify and generalize, we separate the model into two parts: the distribution of the molecular mechanism of persistence in the bacterial population (e.g. number of efflux pumps or target molecules, growth rates) and the elimination rate of single bacteria as a function of that phenotype. Thereby, we obtain an estimate of the required treatment length for each phenotypic subpopulation depending on its size and susceptibility.

Smear positivity in paediatric and adult tuberculosis: systematic review and meta-analysis

BACKGROUND

Tuberculosis (TB) diagnosis continues to rely on sputum smear microscopy in many settings. We conducted a meta-analysis to estimate the percentage of children and adults with tuberculosis that are sputum smear positive.

METHODS

We searched PubMed, MEDLINE, Embase, and Global Health databases for studies that included both children and adults with all forms of active TB. The pooled percentages of children and adults with smear positive TB were estimated using the inverse variance heterogeneity model. This review was registered in the PROSPERO database under registration number CRD42015015331.

RESULTS

We identified 20 studies meeting our inclusion criteria that reported smear positivity for a total of 18,316 children and 162,574 adults from 14 countries. The pooled percentage of paediatric TB cases that were sputum smear positive was 6.8 % (95 % Confidence Interval (CI) 2.2-12.2 %), compared with 52.0 % (95 % CI 40.0-64.0 %) among adult cases. Eight studies reported data separately for children aged 0-4 and 5-14. The percentage of children aged 0-4 that were smear positive was 0.5 % (95 % CI 0.0-1.9 %), compared with 14.0 % (95 % CI 8.9-19.4 %) among children aged 5-14.

CONCLUSIONS

Children, especially those aged 0-4, are much less likely to be sputum smear positive than adults. National TB programs relying on sputum smear for diagnosis are at risk of under-diagnosing and underestimating the burden of TB in children.

Classic reaction kinetics can explain complex patterns of antibiotic action

Finding optimal dosing strategies for treating bacterial infections is extremely difficult, and improving therapy requires costly and time-intensive experiments. To date, an incomplete mechanistic understanding of drug effects has limited our ability to make accurate quantitative predictions of drug-mediated bacterial killing and impeded the rational design of antibiotic treatment strategies. Three poorly understood phenomena complicate predictions of antibiotic activity: post-antibiotic growth suppression, density-dependent antibiotic effects, and persister cell formation. We show that chemical binding kinetics alone are sufficient to explain these three phenomena, using single-cell data and time-kill curves of Escherichia coli and Vibrio cholerae exposed to a variety of antibiotics in combination with a theoretical model that links chemical reaction kinetics to bacterial population biology. Our model reproduces existing observations, has a high predictive power across different experimental setups (R(2) = 0.86), and makes several testable predictions, which we verified in new experiments and by analyzing published data from a clinical trial on tuberculosis therapy. Although a variety of biological mechanisms have previously been invoked to explain post-antibiotic growth suppression, density-dependent antibiotic effects, and especially persister cell formation, our findings reveal that a simple model that considers only binding kinetics provides a parsimonious and unifying explanation for these three complex, phenotypically distinct behaviours. Current antibiotic and other chemotherapeutic regimens are often based on trial and error or expert opinion. Our "chemical reaction kinetics"-based approach may inform new strategies, which are based on rational design.

On being the right size: the impact of population size and stochastic effects on the evolution of drug resistance in hospitals and the community

The evolution of drug resistant bacteria is a severe public health problem, both in hospitals and in the community. Currently, some countries aim at concentrating highly specialized services in large hospitals in order to improve patient outcomes. Emergent resistant strains often originate in health care facilities, but it is unknown to what extent hospital size affects resistance evolution and the resulting spillover of hospital-associated pathogens to the community. We used two published datasets from the US and Ireland to investigate the effects of hospital size and controlled for several confounders such as antimicrobial usage, sampling frequency, mortality, disinfection and length of stay. The proportion of patients acquiring both sensitive and resistant infections in a hospital strongly correlated with hospital size. Moreover, we observe the same pattern for both the percentage of resistant infections and the increase of hospital-acquired infections over time. One interpretation of this pattern is that chance effects in small hospitals impede the spread of drug-resistance. To investigate to what extent the size distribution of hospitals can directly affect the prevalence of antibiotic resistance, we use a stochastic epidemiological model describing the spread of drug resistance in a hospital setting as well as the interaction between one or several hospitals and the community. We show that the level of drug resistance typically increases with population size: In small hospitals chance effects cause large fluctuations in pathogen population size or even extinctions, both of which impede the acquisition and spread of drug resistance. Finally, we show that indirect transmission via environmental reservoirs can reduce the effect of hospital size because the slow turnover in the environment can prevent extinction of resistant strains. This implies that reducing environmental transmission is especially important in small hospitals, because such a reduction not only reduces overall transmission but might also facilitate the extinction of resistant strains. Overall, our study shows that the distribution of hospital sizes is a crucial factor for the spread of drug resistance.

The increasing spread of bacteria, which are resistant to antibiotics, is a serious threat to clinical care. Currently, several countries aim at concentrating highly specialized services in large hospitals in order to improve patient outcomes. However, empirical studies have shown that resistance levels correlate with hospital size. To illustrate this correlation, we analyze two published datasets from the US and Ireland and controlled for antimicrobial usage, disinfection and length of stay. The proportion of patients acquiring both sensitive and resistant infections in hospitals strongly correlated with hospital size. Moreover, we observe the same pattern for both the percentage of resistant infections and the temporal increase of hospital-acquired infections. To investigate to what extent hospital size can directly affect the prevalence of antibiotic resistance, we use mathematical models describing the epidemic spread of resistance in hospitals and the community. We find that small hospitals typically lead to considerably lower resistance levels than large hospitals. However, this beneficial effect of small hospital size may be reduced if bacteria are transmitted indirectly via the environment. Therefore, reducing environmental transmission might be particularly important in small hospitals. Overall, our findings suggest that the short-term benefits of larger hospitals may come at the price of increasing resistance in the long term.

Informed switching strongly decreases the prevalence of antibiotic resistance in hospital wards

Antibiotic resistant nosocomial infections are an important cause of mortality and morbidity in hospitals. Antibiotic cycling has been proposed to contain this spread by a coordinated use of different antibiotics. Theoretical work, however, suggests that often the random deployment of drugs ("mixing") might be the better strategy. We use an epidemiological model for a single hospital ward in order to assess the performance of cycling strategies which take into account the frequency of antibiotic resistance in the hospital ward. We assume that information on resistance frequencies stems from microbiological tests, which are performed in order to optimize individual therapy. Thus the strategy proposed here represents an optimization at population-level, which comes as a free byproduct of optimizing treatment at the individual level. We find that in most cases such an informed switching strategy outperforms both periodic cycling and mixing, despite the fact that information on the frequency of resistance is derived only from a small sub-population of patients. Furthermore we show that the success of this strategy is essentially a stochastic phenomenon taking advantage of the small population sizes in hospital wards. We find that the performance of an informed switching strategy can be improved substantially if information on resistance tests is integrated over a period of one to two weeks. Finally we argue that our findings are robust against a (moderate) preexistence of doubly resistant strains and against transmission via environmental reservoirs. Overall, our results suggest that switching between different antibiotics might be a valuable strategy in small patient populations, if the switching strategies take the frequencies of resistance alleles into account.