Permeability of TB drugs through the mycolic acid monolayer: a tale of two force fields.

Tuberculosis (TB) treatment becomes challenging due to the unique cell wall structure of (M. tb). Among various components of the M.tb cell wall, mycolic acid (MA) is of particular interest because it is speculated to exhibit extremely low permeability for most of the drug molecules, thus helping M.tb to survive against medical treatment. However, no quantitative assessment of the thermodynamic barrier encountered by various well-known TB drugs in the mycolic acid monolayer has been performed so far using computational tools. On this premise, our present work aims to probe the permeability of some first and second line TB drugs, namely ethambutol, ethionamide, and isoniazid, through the modelled mycolic acid monolayer, using molecular dynamics (MD) simulation with two sets of force field (FF) parameters, namely GROMOS 54A7-ATB (GROMOS) and CHARMM36 (CHARMM) FFs. Our findings indicate that both FFs provide consistent results in terms of the mode of drug-monolayer interactions but significantly differ in the drug permeability through the monolayer. The mycolic acid monolayer generally exhibited a higher free energy barrier of crossing with CHARMM FF, while with GROMOS FF, better stability of drug molecules on the monolayer surface was observed, which can be attributed to the greater electrostatic potential at the monolayer-water interface, found for the later. Although both the FF parameters predicted the highest resistance against ethambutol (permeability values of 8.40 × 10 cm s and 9.61 × 10 cm s for the CHARMM FF and the GROMOS FF, respectively), results obtained using GROMOS were found to be consistent with the water solubility of drugs, suggesting it to be a slightly better FF for modelling drug-mycolic acid interactions. Therefore, this study enhances our understanding of TB drug permeability and highlights the potential of the GROMOS FF in simulating drug-mycolic acid interactions.