Predicting rifampicin resistance in M. tuberculosis using machine learning informed by protein structural and chemical features.

Abstract

BackgroundRifampicin remains a key antibiotic in the treatment of tuberculosis. Despite advances in cataloguing resistance-associated variants (RAVs), novel and rare mutations in the relevent gene, rpoB, will be encountered in clinical samples, complicating the task of using genetics to predict whether a sample is resistant or not to rifampicin. We have trained a series of machine learning models with the aim of complementing genetics-based drug susceptibility testing. MethodsWe built a Test+Train dataset comprising 219 susceptible mutations and 46 RAVs. Features derived from the structure of the RNA polymerase or the change in chemistry introduced by the mutation were considered, however, only a few, notably the distance from the rifampicin binding site, were found to be predictive on their own. Due to the paucity of RAVs we used Monte Carlo cross-validation with 50 repeats to train four different machine learning models. ResultsAll four models behaved similarly with sensitivities and specificities in the range 0.84-0.88 and 0.94-0.97 although we preferred the ensemble of Decision Tree models as they are easy to inspect and understand. We showed that measuring distances from molecular dynamics simulations did not improve performance. ConclusionsIt is possible to predict whether a mutation in rpoB confers resistance to rifampicin using a machine learning model trained on a combination of structural, chemical and evolutionary features, however performance is moderate and training is complicated by the lack of data.