Model parameters and fitted values.

Pyrazinamide (PZA) is a key drug in the treatment of Mycobacterium tuberculosis. Although not completely understood yet, the bactericidal mechanism of PZA starts with its diffusion into the cell and subsequent conversion into pyrazinoic acid (POA) after the hydrolysis of ammonia group. This leads to the acidification cycle, which involves: (1) POA extrusion into the extracellular environment, (2) reentry of protonated POA, and (3) release of a proton into the cytoplasm, resulting in acidification of the cytoplasm and accumulation of intracellular POA. To better understand this process, we developed a system of coupled non-linear differential equations, which successfully recapitulates the kinetics of PZA/POA observed in M. tuberculosis. The parametric space was explored, assessing the impact of different PZA and pH concentrations and variations in the kinetic parameters, finding scenarios of PZA susceptibility and resistance. Furthermore, our predictions show that the acidification cycle alone is not enough to result in significant intracellular accumulation of POA in experimental time scales when compared to other neutral pH scenarios. Thus, revealing the need of novel hypotheses and experimental evidence to determine the missing mechanisms that may explain the pH-dependent intracellular accumulation of POA and their subsequent effects.