Mathematical modeling of pulmonary tuberculosis therapy: Insights from a prototype model with rifampin

There is a critical need for improved and shorter tuberculosis (TB) treatment. Current in vitro models of TB, while valuable, are poor predictors of the antibacterial effect of drugs in vivo. Mathematical models may be useful to overcome the limitations of traditional approaches in TB research. The objective of this study was to set up a prototype mathematical model of TB treatment by rifampin, based on pharmacokinetic, pharmacodynamic and disease submodels. The full mathematical model can simulate the time-course of tuberculous disease from the first day of infection to the last day of therapy. Therapeutic simulations were performed with the full model to study the antibacterial effect of various dosage regimens of rifampin in lungs. The model reproduced some qualitative and quantitative properties of the bactericidal activity of rifampin observed in clinical data. The kill curves simulated with the model showed a typical biphasic decline in the number of extracellular bacteria consistent with observations in TB patients. Simulations performed with more simple pharmacokinetic/pharmacodynamic models indicated a possible role of a protected intracellular bacterial compartment in such a biphasic decline. This modeling effort strongly suggests that current dosage regimens of RIF may be further optimized. In addition, it suggests a new hypothesis for bacterial persistence during TB treatment.