Predicting nitroimidazole antibiotic resistance mutations in Mycobacterium tuberculosis with protein engineering

Abstract

Our inability to predict which mutations could result in antibiotic resistance has made it difficult to rapidly identify the emergence of resistance, identify pre-existing resistant populations, and manage our use of antibiotics to effectively treat patients and prevent or slow the spread of resistance. Here we investigated the potential for resistance against the new antitubercular nitroimidazole prodrugs pretomanid and delamanid to emerge in Mycobacterium tuberculosis, the causative agent of tuberculosis (TB). Deazaflavin-dependent nitroreductase (Ddn) is the only identified enzyme within M. tuberculosis that activates these prodrugs, via an F420H2-dependent reaction. We show that the native menaquinone-reductase activity of Ddn is essential for emergence from hypoxia, which suggests that for resistance to spread and pose a threat to human health, the native activity of Ddn must be at least partially retained. We tested 75 unique mutations, including all known sequence polymorphisms identified among ~15,000 sequenced M. tuberculosis genomes. Several mutations abolished pretomanid and delamanid activation in vitro, without causing complete loss of the native activity. We confirmed that a transmissible M. tuberculosis isolate from the hypervirulent Beijing family already possesses one such mutation and is resistant to pretomanid, before being exposed to the drug. Notably, delamanid was still effective against this strain, which is consistent with structural analysis that indicates delamanid and pretomanid bind to Ddn differently. We suggest that the mutations identified in this work be monitored for informed use of delamanid and pretomanid treatment and to slow the emergence of resistance.

About the Authors:

Brendon M. Lee

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia

ORCID logo http://orcid.org/0000-0001-6476-5246

Liam K. Harold

Roles Data curation, Formal analysis, Investigation, Methodology, Visualization

Affiliation: Department of Microbiology and Immunology, University of Otago, Dunedin, New Zealand

ORCID logo http://orcid.org/0000-0002-7570-7348

Deepak V. Almeida

Roles Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Visualization

Affiliation: Center for Tuberculosis Research, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America

Livnat Afriat-Jurnou

Roles Conceptualization, Investigation, Methodology, Writing – review & editing

Affiliations Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia, MIGAL, Galilee Research Institute, Kiryat Shmona, Israel, Faculty of Sciences and Technology, Tel-Hai Academic College, Upper Galilee, Israel

ORCID logo http://orcid.org/0000-0002-4821-0328

Htin Lin Aung

Roles Data curation, Investigation, Methodology

Affiliation: Department of Microbiology and Immunology, University of Otago, Dunedin, New Zealand

ORCID logo http://orcid.org/0000-0002-0465-6256

Brian M. Forde

Roles Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – review & editing

Affiliation: School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

ORCID logo http://orcid.org/0000-0002-2264-4785

Kiel Hards

Roles Formal analysis, Investigation, Methodology, Supervision, Visualization, Writing – review & editing

Affiliation: Department of Microbiology and Immunology, University of Otago, Dunedin, New Zealand

ORCID logo http://orcid.org/0000-0002-7887-4977

Sacha J. Pidot

Roles Data curation, Formal analysis, Investigation, Methodology

Affiliation: Department of Microbiology and Immunology, University of Melbourne, Melbourne, Victoria, Australia

ORCID logo http://orcid.org/0000-0003-1202-6614

F. Hafna Ahmed

Roles Investigation, Methodology, Supervision, Writing – review & editing

Affiliation: Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia

ORCID logo http://orcid.org/0000-0002-9153-2699

A. Elaaf Mohamed

Roles Investigation, Methodology

Affiliation: Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia

ORCID logo http://orcid.org/0000-0002-0873-6990

Matthew C. Taylor

Roles Supervision, Writing – review & editing

Affiliation: Land and Water Flagship, The Commonwealth Scientific and Industrial Organisation, Canberra, Australian Capital Territory, Australia

ORCID logo http://orcid.org/0000-0001-9576-1555

Nicholas P. West

Roles Supervision, Writing – review & editing

Affiliation: School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

ORCID logo http://orcid.org/0000-0003-0955-9890

Timothy P. Stinear

Roles Data curation, Formal analysis, Methodology, Resources, Supervision, Writing – review & editing

Affiliation: Department of Microbiology and Immunology, University of Melbourne, Melbourne, Victoria, Australia

ORCID logo http://orcid.org/0000-0003-0150-123X

Chris Greening

Roles Formal analysis, Funding acquisition, Supervision, Visualization, Writing – review & editing

Affiliations Land and Water Flagship, The Commonwealth Scientific and Industrial Organisation, Canberra, Australian Capital Territory, Australia, School of Biological Sciences, Monash University, Clayton, Victoria, Australia

ORCID logo http://orcid.org/0000-0001-7616-0594

Scott A. Beatson

Roles Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Visualization, Writing – review & editing

Affiliation: School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia

ORCID logo http://orcid.org/0000-0002-1806-3283

Eric L. Nuermberger

Roles Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Visualization, Writing – review & editing

Affiliation: Center for Tuberculosis Research, Johns Hopkins University School of Medicine, Baltimore, Maryland, United States of America

Gregory M. Cook

Roles Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Visualization, Writing – review & editing

Affiliation: Department of Microbiology and Immunology, University of Otago, Dunedin, New Zealand

Colin J. Jackson

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Resources, Supervision, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliation: Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, Australia

ORCID logo http://orcid.org/0000-0001-6150-3822

Introduction

Tuberculosis (TB) is currently the leading cause of death from a single infectious agent [1]. The limitations of current treatment regimens, combined with the rapid emergence of multidrug-resistant tuberculosis (MDR-TB) strains, necessitate the development of new drugs. Three new antitubercular agents are now in advanced clinical development: bedaquiline and delamanid [2], which have been conditionally approved for MDR-TB treatment [3,4], and pretomanid [5], which is approved as part of a promising regimen containing bedaquiline and linezolid (http://www.newtbdrugs.org/pipeline/clinical). The nitroimidazoles, delamanid and pretomanid, are prodrugs that are reductively activated in an F420H2-dependent reaction in Mycobacterium tuberculosis by deazaflavin-dependent nitroreductase (Ddn) [6,7]. An initial hydride transfer step from the F420H2 cofactor leads to their decomposition into des-nitro products and releases reactive nitrogen species that elicit a bactericidal mode-of-action linked to respiratory poisoning and inhibition of mycolic acid synthesis [6–8].

Because pretomanid and delamanid are prodrugs that require activation, mutations that knock-out the activity of Ddn or the biosynthesis or reduction of the enzyme’s cofactor (F420), could confer resistance. However, the fitness cost of such knockouts may be considerable given that (i) F420 has been shown to be conditionally essential to the survival of M. tuberculosis, (ii) F420 is required by at least 28 different enzymes [9], and (iii) it plays important roles in hypoxic survival, protection against oxidative and nitrosative damage, and evasion of the host immune system [10–12]. Ddn is highly conserved across almost all species of mycobacteria (except Mycobacterium leprae), suggesting its physiological role is under strong evolutionary selection [13]. It has been hypothesised that Ddn serves as an F420H2-dependent menaquinone reductase given its membrane localisation and catalytic activity with the synthetic quinone analogue menadione [12,14], although further work is required to fully define its physiological role.

Despite the recent introduction of nitroimidazoles, cases of acquired clinical resistance, i.e. resistance that occurs during the long treatment of TB infection but is not necessarily transmissible, have already been reported [15,16]. Acquired resistance to pretomanid and delamanid can occur through genetic changes that cause loss of function within the biosynthetic pathway for F420 production or in the F420-dependent glucose 6-phosphate dehydrogenase (FGD) that catalyzes F420 reduction to F420H2 [17,18]. Laboratory studies have also shown that genetic changes that abolish Ddn activity can also confer resistance [18–20]. While such mutations could compromise treatments for individuals that are already infected, transmission of M. tuberculosis to healthy individuals after these genetic changes has never been documented. We speculate that in order to for tuberculosis to spread effectively and endanger health, these resistant strains would need to retain sufficient fitness to survive all stages of the lifecycle of M. tuberculosis, including recovery from hypoxia. Interestingly, delamanid-resistant isolates with mutations in ddn have been recovered from MDR-TB patients who never received delamanid or pretomanid [21,22], raising important questions about the fitness costs associated with such mutations and their potential impact on transmission.

In this study, we analysed Ddn orthologs from related mycobacteria to identify natural sequence variations that make mycobacteria resistant to pretomanid, as well as analysing the genomes of ~15,000 M. tuberculosis isolates to identify the spectrum of naturally occurring non-synonymous Ddn polymorphisms. Mutations were then made to Ddn at positions identified in orthologs and through analysis of non-synonymous polymorphisms to analyse their effect on the native activity (quinone reduction) and pretomanid activation. Altogether, 75 mutants, at 47 unique positions within the 151 amino acid Ddn protein were analysed. This identified a number of mutations that prevent pretomanid/delamanid activation, without full loss of the native menaquinone reductase activity. Analysis of complete and partial Ddn knock-outs demonstrated that it is essential for resuscitation from hypoxia, i.e. genetic variants that have lost their nitroimidazole-activation function, but retained their native activity, appear to be sufficiently fit to spread and cause disease in new patients. We examined a hypervirulent strain of M. tuberculosis from Vietnam that contains one such mutation (despite never being exposed to pretomanid), which is resistant to pretomanid. Curiously, sensitivity to delamanid activation was not affected in this strain, nor was delamanid activation by the Ddn variant in vitro, which is consistent with our structural analysis that suggests it binds to the active site of Ddn in an alternative orientation.

Results

The physiological role of Ddn

A role for Ddn as a quinone reductase was previously suggested based on activity with the synthetic menaquinone analogue menadione [12]. Here, we show that purified Ddn catalyzes the F420H2-dependent reduction of menaquinone, which is known to accept electrons from a variety of electron donors and transfer them to terminal oxidases or reductases in mycobacterial respiration [23]. Ddn catalyzed menaquinone-1 reduction in vitro with a kcat/KM value of ~103 M-1 s-1 and a KM of 22.4 ± 3.8 μM (Fig 1a, Table 1). It should be noted that the catalytic rate is likely to be higher in the crowded native environment of the M. tuberculosis cell where Ddn and menaquinone are co-localised at the cell membrane [24]. The enzyme kinetics of Ddn are comparable to another mycobacterial menaquinone oxidoreductase, MenJ (Rv0561c), which is found in the menaquinone biosynthetic pathway (kcat/KM of ~104 M-1 s-1 and KM of 30 μM [25]). The KM value is significantly lower than the median KM of enzymes with their native substrates (~130 μM), and the kcat/KM value is within the broad distribution of known enzymes (103−107 M-1 s-1), albeit at the lower end of the range [26,27]. Ddn orthologs encoded by Mycobacterium smegmatis (MSMEG_2027, MSMEG_5998) were also able to reduce menaquinone, suggesting that Ddn orthologs have similar physiological roles across the genus (Table 1). Indeed, Ddn and its orthologs are highly conserved and abundant throughout mycobacteria [13].

[Figure omitted. See PDF.]

Fig 1. Ddn is a respiratory primary dehydrogenase that couples F420H2 oxidation to menaquinone reduction.

A. Michaelis-Menten curve of the activity of purified Ddn with menaquinone-1 in the presence of F420H2. Rates were determined by measuring change of absorbance at 420 nm. Error bars indicate SEM (n = 3). B. Relative oxygen consumption of M. smegmatis membranes. F420, NADH, and glucose 6-phosphate were present in all assays, while the availability of the F420-dependent glucose 6-phosphate dehydrogenase (FGD) was varied. Membranes were purified from four different genetic backgrounds: wild-type, a transformant with the empty complementation vector pMV261, a transformant with pMV261 expressing the ddn gene, and a strain containing a chromosomal cytochrome bd oxidase deletion (ΔcydA). Error bars indicate SEM (n = 3). Stars above each column indicate a statistically significant difference compared to the first column. Other significant differences are shown by starred lines between columns. * = p ≤ 0.05, ** = p ≤ 0.01, **** = p ≤ 0.0001 n.s = p > 0.05 (one-way ANOVA followed by Tukey’s multiple comparison test). C. Rate of NADH oxidation in the presence and absence of F420H2. Rates were determined by following change of absorbance at 340 nm. n.s indicates p > 0.05, two tailed t-test, error bars indicate SEM (n = 3).

https://doi.org/10.1371/journal.ppat.1008287.g001

[Figure omitted. See PDF.]

Table 1. Kinetic parameters of native F420H2-dependent menaquinone reductase activities of Ddn and M. smegmatis orthologs in in vitro assays.

https://doi.org/10.1371/journal.ppat.1008287.t001

Having established that Ddn can reduce menaquinone with efficiency that could be physiologically relevant, along with the fact that menaquinone is the only native quinone in M. tuberculosis [28–30], we then investigated whether the menaquinone reductase activity of Ddn was coupled to the mycobacterial respiratory chain by comparing the rates of respiratory oxygen consumption of M. smegmatis in the presence and absence of a complementation vector expressing ddn. We observed that the addition of glucose 6-phosphate (G6P), F420, and F420 dependent glucose-6-phosphate dehydrogenase (FGD), which catalyzes the reduction of F420 to F420H2 for use by F420H2-dependent enzymes such as Ddn [31], resulted in a statistically significant increase in oxygen consumption by mycobacterial membranes when the membranes are activated by NADH (Fig 1b). This was dependent on the presence of cytochrome bd oxidase, which utilizes reduced menaquinone (i.e. menaquinol) as an electron source (Fig 1b). The rate of NADH oxidation was not affected by the addition of F420H2 in the absence of ddn, which implies endogenous NADH-dependent oxidases do not contribute to the increased oxygen consumption (Fig 1c). This supports the idea that mycobacteria can couple F420H2 oxidation to O2 reduction through the respiratory chain via cytochrome bd oxidase and demonstrates that bacteria can use F420H2 as a respiratory electron donor. We also observed that the extent of oxygen consumption was significantly higher in membranes purified from the strain that expressed ddn, relative to empty vector controls (Fig 1b), which suggest Ddn serves as a menaquinone reductase in the respiratory chain, with the remaining stimulation attributable to F420H2 oxidation by native Ddn homologs of M. smegmatis. Altogether, these results suggest that Ddn likely plays a role in bacterial respiration through reduction of menaquinone.

Natural sequence variation in Ddn can result in loss of pretomanid activation with retention of the native activity

Having established that Ddn is a menaquinone reductase, we then investigated the effects of natural sequence variation within Ddn orthologs from mycobacteria on the native, and drug-activating, activities. We expressed and purified Ddn orthologs from M. tuberculosis, M. marinum, M. smegmatis, M. vanbaalenii, M. avium and M. ulcerans. All Ddn orthologs catalyzed the reduction of menadione, suggesting that this is a shared physiological function that has been under consistent selective pressure. However, only Ddn from M. tuberculosis and M. marinum could activate pretomanid (Fig 2a, Table 2). This demonstrates that the native Ddn activity can exist in in the absence of nitroimidazole reductase activity, i.e. the drug-activating activity is not coupled to the native activity. The binding site of Ddn has been defined over several studies, identifying a number of polar amino acids within the substrate binding site that contribute to activity, particularly Y65, S78, Y130, Y133, Y136 (Fig 2b). Three of these binding site residues are fully conserved among the Ddn orthologs tested here, whereas Y65 and Y133 are more variable (Fig 2c). These results show that sequence variation among Ddn orthologs in mycobacteria can result in the loss of the promiscuous, nitroimidazole-activating activity.

[Figure omitted. See PDF.]

Fig 2. The activity of Ddn and orthologs from other mycobacterial species and key substrate binding pocket residues.

(a) The activity of Ddn and orthologs with menadione, pretomanid, and delamanid. Error bars show standard deviations from three independent replicates. (b) Multiple sequence alignment of Ddn and orthologs. Highlighted residues indicate active site residues and numbers indicate their residue position in Ddn. (c) The binding pocket of Ddn consists of F420/F420H2, Tyr65, Ser78, Tyr130 and Tyr136.

https://doi.org/10.1371/journal.ppat.1008287.g002

[Figure omitted. See PDF.]

Table 2. Reduction rates of F420H2-dependent menadione and nitroimidazole reductase activities of Ddn mutants.

https://doi.org/10.1371/journal.ppat.1008287.t002

Having demonstrated that mutations can cause the loss of nitroimidazole-reducing activity in vitro, we investigated whether this corresponded to differences in nitroimidazole susceptibility in vivo. M. tuberculosis (H37Rv) and M. marinum, which were the only two species that encoded Ddn orthologs with in vitro pretomanid activation activity, were found to be susceptible to pretomanid and delamanid treatment (Table 3). In contrast, species in which the Ddn orthologs did not exhibit pretomanid activation activity (M. smegmatis, M. ulcerans, M. avium) have been shown to be naturally resistant to pretomanid [5,32,33]. This analysis confims that the prodrug-activating activity of Ddn from M. tuberculosis H37Rv and M. marinum must result from sequence differences to the other Ddn orthologs tested.

[Figure omitted. See PDF.]

Table 3. MICs of pretomanid and delamanid against different mycobacterial strains.

https://doi.org/10.1371/journal.ppat.1008287.t003

Genome sequences of M. tuberculosis were then searched to identify nonsynonymous sequence polymorphisms within the ddn gene. We found that around 1.5% (219/14,876) of presumptive M. tuberculosis genomes screened encoded a non-synonymous mutation in ddn, including several that have arisen independently in unrelated strains (S1 Table). Altogether, we identified 46 non-synonymous substitutions and 2 deletions in ddn, distributed throughout the M. tuberculosis phylogeny (Table 4).

[Figure omitted. See PDF.]

Table 4. Distribution of Ddn alleles in the major M. tuberculosis lineages.

Distribution of Ddn alleles in five major lineages of M. tuberculosis.

https://doi.org/10.1371/journal.ppat.1008287.t004

All 46 mutants found in our genomic survey were expressed, purified, and assayed with menadione, pretomanid, and delamanid (Table 2, Figs 3 and 4). Every mutant tested was able to reduce menadione, which is consistent with selective pressure for maintenance of the physiological function of Ddn. The majority of the mutants could also reduce/activate pretomanid and delamanid. This is not unexpected, as these M. tuberculosis strains have not (to our knowledge) been exposed to either drug and have therefore not been under selective pressure to develop resistance. However, of the 46 mutants, several mutants did not activate pretomanid (L49P, S78Y, K79Q, W88R, and Y133C) indicating that any strain of M. tuberculosis with these mutations would be unable to activate the drug. Indeed, an M. tuberculosis strain haboring the W88R mutation has been shown to be resistant to pretomanid in vitro [20], and the L49P mutant was obtained from an in vivo resistance selection experiment [34]. It is notable that the S78Y and Y133C mutants retained the ability to activate delamanid.

[Figure omitted. See PDF.]

Fig 3. Mutations to Ddn that were made in this study.

The holo-enzyme structure of Ddn with a reconstructed N-terminus and F420 and pretomanid bound is shown, with mutations highlighted from mutants selected for pretomanid resistance in laboratory (orange spheres), non-synonymous mutations found in sequenced M. tuberculosis strains (green spheres), and potential spontaneous mutations within the substrate binding site of Ddn (blue spheres).

https://doi.org/10.1371/journal.ppat.1008287.g003

[Figure omitted. See PDF.]

Fig 4. Natural sequence polymorphisms found in M. tuberculosis stains.

The activity of Ddn mutants with menadione, pretomanid, and delamanid. Bar color represents level of activity related to the rest of the graph with the highest activity colored red and the lowest activity colored blue. Error bars show standard deviations from three independent replicates. Data from this Figure is also present in Tables 2 and 3.

https://doi.org/10.1371/journal.ppat.1008287.g004

The S78Y polymorphism is found in the genome of N0008 [35], a clinical isolate of the hypervirulent Beijing family [15], and in two other genomes (SRA accessions:ERR718320 and ERR751847) that are phylogenetically closely related, indicating a shared evolutionary history and suggesting that there was no substantial selective pressure to eliminate this mutation, i.e. its fitness cost must have been relatively small (Fig 5). We obtained the Beijing strain N0008 to investigate whether the S78Y mutation in the ddn gene, which results in loss of pretomanid reduction in vitro (Table 3), corresponded to resistance to pretomanid in vivo. No other genetic changes previously observed to cause nitroimidazole resistance (e.g. loss of F420 biosynthetic genes) were apparent in the genome of N0008. We observed a 64-fold increase in minimum inhibitory concentration (MIC) of pretomanid against N0008 (256 μg mL-1) compared to H37Rv (4 μg mL-1). This suggests SNPs, such as L49P, S78Y and W88R, in the ddn gene can confer resistance to pretomanid and confirms that transmissible pretomanid-resistant populations of M. tuberculosis already exist.

[Figure omitted. See PDF.]

Fig 5. Phylogeny of M. tuberculosis strains with sequence polymorphisms in Ddn.

Maximum-likelihood phylogenetic tree of 332 M. tuberculosis strains built using 2,099 non-recombinant core genome SNPs (relative to H37Rv). The tree shows all M. tuberculosis strains for which public genome data is available in Genbank or the short read archive (SRA) that have either synonymous and non-synonymous mutations in Ddn (relative to H37Rv). The amino acid changes of the non-synonymous mutations are indicated on the branch tips. Blank tips represent strains with synonymous mutations in Ddn. Branches are coloured by Lineage: Lineage 1 (Indo-Oceanic), yellow; Lineage 2 (East Asian-Beijing), red; Lineage 3 (East African-Indian), purple; Lineage 4 (Euro-American), green; Lineage 5 (West African), teal. The scale bar indicates branch length in number of SNPs. Genome alignments, recombination filtering and phylogenetic reconstruction were done using Parsnp, Gubbins and RaxML, respectively. The phylogenetic tree was visualised using Figtree version 1.4.3 (https://tree.bio.ed.ac.uk/software/figtree).

https://doi.org/10.1371/journal.ppat.1008287.g005

The potential for spontaneous pretomanid resistance mutations to arise in M. tuberculosis

Previous studies have used laboratory evolution or engineering to investigate the potential for pathogens to evolve resistance to antibiotics [36,37]. We took a similar approach here, using structure-guided mutagenesis to investigate the robustness of the nitroreductase activity of Ddn to mutations. The binding site of Ddn has been defined over several studies, identifying a number of polar amino acids within the substrate binding site that contribute to activity, particularly Y65, S78, Y130, Y133, Y136. Three of these binding site residues are fully conserved among the Ddn orthologs tested here, whereas Y65 and Y133 are more variable (Fig 2c). We made, expressed and purified 26 mutants (including the Y65S, S78Y, and Y133C sequence differences observed some of the other Ddn orthologs) of these five key substrate binding residues (Y65, S78, Y130, Y133, Y136) (Fig 3), to test how spontaneous mutations at these positions affect nitroimidazole activation. Across the 26 mutants tested, all retained significant levels of the native activity, while 16 did not display detectable pretomanid activation (Table 2, Fig 6). In other words, while the native activity was not greatly affected by sequence variation, the promiscuous nitroreductase activity was extremely sensitive to mutation.

[Figure omitted. See PDF.]

Fig 6. Kinetic activities of mutants with potential spontaneous mutations within the substrate binding site of Ddn.

The activity of Ddn mutants with menadione, pretomanid, and delamanid. Bar color represents level of activity related to the rest of the graph with the highest activity colored red and the lowest activity colored blue. Error bars show standard deviations from three independent replicates. Data from this Figure is also present in Tables 2 and 3.

https://doi.org/10.1371/journal.ppat.1008287.g006

To estimate the rate at which mutations could spontaneously arise at one of the positions identified in this work, we used our genomic survey of ddn in M. tuberculosis to estimate the level of sequence variation in ddn. This revealed that virtually every distinct lineage of M. tuberculosis contains SNPs in ddn; given the burden of TB in the world and the potential for widespread pretomanid administration, alongside the apparent sensitivity of the promiscuous pretomanid activation activity to such SNPs in contrast to the native activity, it appears that spontaneous resistance enabling mutations of ddn could readily occur and spread with sufficient selective pressure.

The molecular basis of resistance

The effects of the mutations shown in Table 2 are generally consistent with our mechanistic understanding of Ddn from mutagenesis and computational simulation [7,38]. For example, S78 is thought to interact with the nitro-moiety of pretomanid and to stabilize the transition state; none of the S78 mutants retained pretomanid nitroreductase activity, suggesting this interaction is particularly important. Y65, Y130, Y133 and Y136 are known to form a hydrophobic wall in the binding site, which can move during the catalytic cycle to shield the active site from solvent and thereby facilitate pretomanid reduction [38]. This is in keeping with the observation that tyrosine to phenylalanine mutations at positions 65, 133 and 136 were essentially neutral, whereas substitution with other residues (Met, Leu, Cys, Trp, Thr, Ser, Glu) led to loss of pretomanid activation. Menadione reduction by Ddn is less susceptible to loss of activity through mutation, which is consistent with work showing native activities are substantially more robust to mutation than promiscuous activities (such as pretomanid activation) [39], as well as the observation that menadione is more chemically labile.

Delamanid activation is less susceptible to resistance mutations

Of the 75 mutants we made and tested, 25 did not reduce pretomanid at detectable levels but only 10 lost the ability to reduce delamanid. Thus, although pretomanid and delamanid are superficially similar, they must interact with Ddn differently. These enzymatic data extend to whole cell activity, as we observed that the N0008 strain harbouring the S78Y mutation was resistant to pretomanid but remained susceptible to delamanid (Table 3). To better understand the molecular basis for the differential effects of mutations on pretomanid and delamanid activation, we used molecular docking (we were unable to capture ternary crystal structures with ligand and cofactor) to obtain low-energy poses of these substrates in the crystal structure of Ddn. The binding of delamanid is predicted to be different than that of pretomanid (Fig 7B and 7C), with the dual methyl and phenoxy-methyl substituents on the oxazole ring preventing delamanid binding above the deazaflavin ring of F420 in a ring-stacked orientation, as is seen in pretomanid, which has an oxazine ring with a single substituent in the analogous position. This results in a change in the angle at which the nitroimidazole group interacts with the cofactor that results in an increase in the distance to S78Y. Thus, there appears to be a plausible molecular basis for the differing effects of the S78Y mutation on pretomanid and delamanid activity in vitro and in vivo.

[Figure omitted. See PDF.]

Fig 7. The substrate-binding pocket of Ddn (PDB: 3R5R).

(a) Menadione (docking score = -6.584 kJ/mol) binds above the deazaflavin ring of F420H2 in a complementary pocket formed by the aromatic rings of the tyrosine residues and the hydroxyl group of Ser78. (b) Pretomanid (docking score = -5.995 kJ/mol) docked into wild-type Ddn with F420H2 bound. The nitroimidazo-oxazine moiety can bind parallel to the deazaflavin group of F420, extending the nitro-group deep into the active site, towards Ser78. (c) Delamanid (docking score = -4.890 kJ/mol) docked into wild-type Ddn with F420H2 bound. The dual methyl/phenoxy-methyl substitution at the C6 position of the oxazole ring creates steric hindrance with the deazaflavin ring of F420, causing it to bind above F420 in a perpendicular orientation. This results in additional distance to Ser78.

https://doi.org/10.1371/journal.ppat.1008287.g007

Ddn mutants are virulent in mice but show defective recovery from hypoxic stress in vitro

To understand the fitness costs of loss of Ddn activity through mutation, we investigated the in vitro enzymatic properties of Ddn variants that have been identified in M. tuberculosis. Several isogenic pretomanid-resistant M. tuberculosis mutants have recently been selected by pretomanid monotherapy in infected mice and shown to harbour mutations in Ddn (M1T, L49P, L64P, R112W, C149Y and from insertion of the IS6110 transposable element at D108) [34]. We also investigated two mutants from a previous study in which resistant strains were identified from an in vitro selection experiment (S22L, W88R). We tested the ability of these eight mutants to reduce native (menadione) and drug (pretomanid and delamanid) substrates in vitro, finding that all mutations resulted in loss of detectable pretomanid activation, consistent with their selection in resistant strains. We observed that the M1T mutant (loss of start codon) and the mutants harbouring the IS6110 insertions did not produce any functional protein and therefore had no detectable native, nor prodrug-activating activity. In contrast, many of the point mutants, such as the L64P mutant, retained a small amount of activity with menadione, suggesting that it might retain a fraction of its native function (Table 2). Because many of these mutants were found in M. tuberculosis living in mouse lungs, the fitness cost of these mutations is clearly low enough to allow them to survive in that environment.

To understand the fitness costs of loss of Ddn activity in more detail, we investigated the ability of M. tuberculosis strains with mutations in Ddn to survive in different conditions. We therefore tested a selection of variants: wild-type H37Rv (full native activity), L64P (partial loss of native activity), as well as M1T and IS6110 (complete loss of native activity) to examine whether these mutants were attenuated for multiplication and survival in mice. These mutants showed no difference in growth or survival compared to wild type after low-dose aerosol infection of mice (Fig 8a), suggesting that these Ddn mutants would be transmissible.

[Figure omitted. See PDF.]

Fig 8. Fitness analysis of Ddn mutants.

(a) M. tuberculosis Ddn mutants are not attenuated in mice. Mean lung colony-forming unit (CFU) counts of the H37Rv strain and isogenic Ddn mutants are shown following low-dose aerosol infection. Error bars represent SEM biological triplicate. Competition of the Δddn mutant compared to isogenic wild type mc2 6230. (b) Growth and survival of the Δddn mutant compared to isogenic wild type in aerobic conditions. (c) Growth and survival of the Δddn mutant compared to isogenic wild type in hypoxic conditions. (A and B) Arrows indicate samples taken for competition assay. (d) Percentage of the Δddn mutant when grown in competition with the isogenic wild type. Two-way ANOVA with Tukey multiple comparisons test where *, p<0.05; **, p<0.01. Red * is comparison between aerobic (Day 7) and Hypoxic (Day 37) at given timepoint, Black * is comparison between aerobic (Day 7) and Hypoxic (Day 70) at given timepoint. Error bars for b, c and d representative of standard deviation biological triplicate.

https://doi.org/10.1371/journal.ppat.1008287.g008

A previous study showed that M. tuberculosis mutants that were unable to biosynthesise F420 have a survival defect when recovering from hypoxia [12]. To better understand the competitive cost that loss of Ddn function causes to M. tuberculosis during recovery from hypoxia, we created a hygromycin marked deletion mutant in the mc26230 strain of Mycobacterium tuberculosis using phage transduction. The suspected Δddn mutants with hygromycin resistance displayed MIC values against delamanid of greater than 512 μM. The deletion mutants were then confirmed using whole genome sequencing (S1 Fig). The Δddn mutant displayed no growth or survival defects compared to the isogenic wild type when grown in aerobic or hypoxic conditions (Fig 8b and 8c) consistent with the previous results on ddn point mutants (Fig 8a). We then tested the marked Δddn strain in competition with the wild type isogenic strain when reaerated from hypoxia. As a control we took the day 7 aerobic growth sample of each strain and grew them in competition. We found that the Δddn mutant had lost no competitiveness when taken from aerobic conditions with the proportional percentage of Δddn remaining at about ~50% (Fig 8d). We then took equal proportions of each strain from hypoxia at 37 days and 70 days, which we grew in competition (Fig 8c). The proportion of Δddn was then measured at days 0, 7 and 17. We found that with the Day 37 hypoxic sample the Δddn mutant had an initial drop in proportion to 28% at 7 days but by day 17 the proportion of the Δddn mutant had returned to 57% (Fig 8d). The day 70 hypoxic sample had an initial drop off Δddn proportion to 33% which was maintained through until day 17 (Fig 8d). These data demonstrate that the Δddn mutant cannot compete with wild type after exit from prolonged hypoxia and that Ddn contributes to fitness under these conditions, perhaps because they are more susceptible to oxidative stress or have some defect in respiration.

Discussion

The fitness trade-off between resistance and native function

Given the conservation of Ddn-like genes throughout mycobacteria, its activity with menaquinone, and the inability of Ddn knock-outs to recover from hypoxia, we suggest that mutations that completely knock-out Ddn activity (including the loss of F420 biosynthesis through mutations to F420 biosynthetic genes, loss of F420 reductase activity through knockout of FGD, and introduction of stop codons or large genetic insertions/deletions in ddn) could result in substantial loss of fitness given that the ability to recover from hypoxia is an important aspect of M. tuberculosis pathogenesis [3,23,40–42]. Thus, for nitroimidazole resistance to spread and endanger health, the activity of Ddn must be either retained or otherwise compensated for. This means that mutations that are neutral or near-neutral for the wild-type activity, but knock out nitroimidazole-activation, will be of great clinical concern.

The binding of the prodrugs pretomanid and delamanid by Ddn, which leads to the activation of these prodrugs, is a promiscuous activity that is not fully coupled to its ability to bind and reduce its native substrate; i.e. loss of pretomanid activation is possible without total loss of the native activity. It has long been known that promiscuous activities of enzymes are less robust to the effects of mutations than native enzyme functions [39], meaning that it is possible that mutations could cause loss of prodrug activation without significant loss of native function. In this study we tested this assumption for Ddn, mutating ~1/3 of the amino acid positions in M. tuberculosis Ddn, including all known naturally occurring polymorphisms, showing that many such mutations (including several already present in clinical isolates) can cause loss of delamanid and/or pretomanid activation without loss of the native activity. Extending our enzymatic measurements to measurement of MICs revealed that related species with sequence differences in ddn, and clinical isolates of M. tuberculosis that harbor mutations within the active site, are resistant to pretomanid. These results have important implications for the clinical use of nitroimidazoles in TB treatment in order to optimize treatment outcomes and to prevent or slow the development of resistance. For instance, pretomanid will not be effective in patients infected with naturally occurring M. tuberculosis variants harbouring an S78Y mutation and indiscriminant use of pretomanid against such variants (such as in regions in which the N0008 strain is endemic) could drive selective amplification and spread of pretomanid resistance. Moreover, the chances of spontaneous mutation of Ddn in currently sensitive strains is significant, given our observation of the number of single nucleotide mutations that can knock out nitroreductase activity. Given that the clinical isolate N0008 is highly transmissible, it is likely that other Ddn variants in which nitroreductase activity can be abolished without substantial loss of the native activity will also be infectious owing to the minimal loss of fitness these mutations cause with the native activity. Thus, fitness-neutral Ddn mutations are likely to be a prominent route through which clinical pretomanid and delamanid resistance spreads.

Our findings have broad implications for the continued clinical development and usage of nitroimidazole antitubercular agents. Through multiple ongoing phase II/III clinical trials, the TB Alliance and others are evaluating novel combination therapies including pretomanid. The ongoing ZeNix and TB-PRACTECAL trials are extending the study of the BPaL (bedaquiline, pretomanid, linezolid) regimen that showed highly promising results against MDR- and XDR-TB in the Nix-TB trial, while the SimpliciTB trial studies the BPaMZ regimen (bedaquiline, pretomanid, moxifloxacin, pyrazinamide) against MDR- and drug-susceptible TB. Meanwhile, comparable delamanid-containing regimens are being studied in the endTB and MDR-END trials. Now that both delamanid and pretomanid are approved for clinical use, intensive monitoring of ddn polymorphisms is recommended to ensure informed regimen selection and allow interventions that will reduce spreading of transmissible resistance. Our findings indicate that delamanid binds to Ddn in a different conformation than pretomanid does, suggesting that it could be effective against some pretomanid-resistant isolates. It is also possible that combination therapy with both nitroimidazoles could help prevent the evolution and spread of resistance (since the simultaneous loss of both activities will be less likely to result from single amino acid substitution). Further studies on how delamanid and pretomanid are activated in the M. tuberculosis cell will inform the development of improved nitroimidazole therapies and testing of a broad range of nitroimidazole analogs against a panel of Ddn variants could help to identify compounds for which resistance is less likely to evolve.

Materials and methods

Plasmid construction and point mutation

The E. coli codon optimised sequences for Ddn from M. tuberculosis, MSMEG_5998 from M. smegmatis, MMAR_5035 from M. marinum, MVAN_5261 from M. vanbaalenii, MAV_0613 from M. avium, and MUL_4109 from M. ulcerans were purchased as gene strings from ThermoFisher Scientific (Massachusetts, USA) and cloned into the expression vector pMAL-c2X using Gibson assembly [43]. All mutations to Ddn were made by site-directed mutagenesis using Gibson assembly [43]. Construction of MSMEG_2027 and FGD has been described previously [44,45].

Protein expression and purification

MSMEG_2027 was expressed and purified as previously described [13]. Ddn, Ddn mutants, and Ddn orthologs were transformed into E. coli BL21 (DE3) cells and grown on LB agar containing 100 μg/ml ampicillin. Single colonies were picked and inoculated in LB media with 100 μg/ml ampicillin. Starter cultures were grown overnight and diluted 1/100 and grown at 37 °C until OD = 0.4. Cultures were induced with IPTG to a final concentration of 0.3 mM, and grown for 3 h at 25 °C. Cells were harvested by centrifugation at 8,500 × g for 20 minutes at 4 °C and resuspended in lysis buffer (20 mM Tris-Cl, pH 7.5, 200 mM NaCl) and lysed by sonication using an Omni Sonicator Ruptor 400 (2 x 6 min. at 50% power). The soluble extract was obtained by centrifugation at 13,500 × g for 1 h at 4 °C. The protein was purified using amylose resin (NEB) using the provided protocol. Briefly the lysate was passed over the amylose resin and washed with 12 column volumes of lysis buffer. The protein was eluted using elution buffer (same as lysis buffer but with 10 mM maltose). Samples were frozen at -80 °C in 20 mM tris pH 7.5, 200 mM NaCl, 10 mM maltose, and 10% glycerol.

FGD was expressed and purified as described by Bashiri et al. [45], with minor modifications. FGD was transformed into E. coli BL21 (DE3) cells and grown on LB agar containing 100 μg/ml ampicillin. Single colonies were picked and inoculated in Terrific Broth (TB) [46] with 100 μg/ml. ampicillin. Starter cultures were grown overnight and diluted 1/100 into auto-induction media (20 g/l tryptone, 5 g/l yeast extract, 5 g/l NaCl, 6 g/l Na2HPO4, 3 g/l KH2PO4, 6 ml/l glycerol, 2 g/l lactose, 0.5 g/l glucose, 100 mg/ml ampicillin) and grown at 30 °C for 24 h. The cells were harvested by centrifugation at 8500 g for 20 min at 4 °C and resuspended in lysis buffer (20 mM NaPO4 pH 8, 300 mM NaCl, 25 mM imidazole) and lysed by sonication using an Omni Sonicator Ruptor 400 (2 x 6 min at 50% power). The soluble extract was obtained by centrifugation at 13500 g for 1 h at 4 °C. The soluble fraction was filtered and loaded 5-ml HisTrap HP column (GE Healthcare) and washed with lysis buffer. The protein was eluted with elution buffer (lysis buffer with 250 mM imidazole). The purified protein was dialyzed in 50 mM Tris-Cl, pH 7, 200 mM ammonium sulphate. Samples were frozen at -80 °C in 50 mM Tris-Cl, pH 7, 200 mM ammonium sulphate, and 10% glycerol.

Enzymatic assays

F420 was purified from M. smegmatis mc24517 as described by Ahmed et al.[13]. F420 was reduced overnight with 10 μM FGD and 10 mM glucose-6-phosphate in 20 mM Tris-CL, pH 7.5 under anaerobic conditions. FGD was removed by spin filtration in a 1 mL 10 K MWCO spin filter (Millipore) and used F420H2 within 8 hours. Enzyme assays were performed according to Ahmed et al. [13,47] in 200 mM Tris-HCl, pH 7.5, 0.1% Triton X-100, 25 μM F420H2, and 25 μM of substrate at room temperature. Enzyme concentrations used were between 0.1 μM and 1 μM. Activity was monitored following the oxidation of F420H2 which was measured spectrophotometrically at 420 nm (ε = 41,400 M−1 cm−1) [48]. Specific enzymatic activity with delamanid was measured using fluorescence (excitation/emission: 400 nm/470 nm).

Computational analysis

Sequences of Ddn and orthologs were obtained from the NCBI sequence database. Alignment of the sequences were performed using MUSCLE [49] via the EMBL-EBI web services [50]. Autodock Vina [51] was used to dock menadione, pretomanid and delamanid into Ddn (PDB ID: 3R5R) [7]). The protein and ligand were prepared using in Autodock tools with default settings [52] and visualized using Pymol [53]. Substrate structures were obtained from the ZINC database [54].

Drug susceptibility testing

Minimum inhibitory concentration (MIC) testing of M. marinum was performed according to the method of Wiegand et al. [55] with some minor modifications. A culture was grown on Brown and Buckle media and colonies were scraped and diluted to an OD600 of 0.2, and was used as the inoculum for the MIC assay. Plates were incubated at 30 °C in a humid environment and were read after 5 days of incubation. MICs of M. tuberculosis isolates were determined as described previously [56] with slight modifications. Briefly, U-bottomed 96-well microtitre plates containing Middlebrook 7H9 broth supplemented with oleic acid, albumin, dextrose and catalase were inoculated with M. tuberculosis H37Rv or N0008 (OD600 0.02) in the presence of delamanid and pretomanid (0.015–512 μg/ml), and incubated at 37 °C for 5 days. The presence/absence of a cell pellet was checked visually. Resazurin (0.03% w/v) was then added, and plates were further incubated for 7 days. Viable cells reduce resazurin (blue) to resorufin (pink). Bedaquiline (0.015–32 μg/ml) was used as a drug susceptibility control.

Dataset

To determine the frequency of non-synonymous sequence polymorphisms in the ddn gene we assembled a collection of 5,184 publicly available complete and draft M. tuberculosis genomes from Genbank [57] on the 7th of April 2017. Additionally, using the structured query https://www.ncbi.nlm.nih.gov/sra/(Illumina[Platform]) AND “Mycobacterium tuberculosis”[orgn:__txid1773]) to search the SRA (https://www.ncbi.nlm.nih.gov/sra/) opened on the 28th of April 2017, we identified 9,692 unassembled M. tuberculosis genome datasets.

Lineage-typing of M. tuberculosis genomes

Classification of M. tuberculosis phylogenetic lineages was performed using KvarQ version 0.12.3a1 and default parameters [58]. Using KvarQ scan, unassembled sequencing data were screened against a database of previously defined sequence polymorphisms known to delineate the 7 major M. tuberculosis phylogenetic lineages (S1 Table).

Characterisation of ddn sequence polymorphisms in M. tuberculosis

Using BLASTn (with default parameters) complete and draft M. tuberculosis genome assemblies from Genbank were queried with the ddn gene from M. tuberculosis H37Rv (Genbank accession:AL123456; [59]), identifying 145 M. tuberculosis genomes carrying alternative ddn alleles (relative to H37Rv). For these 145 genomes the Ddn sequence was extracted and aligned using Clustal Omega version 1.2.4 [60] to identify non-synonymous mutation and other sequence polymorphism which could impact the function of Ddn (S1 Table).

For the unassembled genomes datasets, raw sequencing reads were aligned to the reference genome M. tuberculosis H37Rv using Snippy version 2.9 (BWA-mem version 0.7.12 [61])(https://github.com/tseeman/snippy). Variant calling was performed using Snippy (Freebayes version 0.9.21 [62] with default parameters (minimum read coverage of 10x and 90% read concordance at the variant locus). 244 unassembled genome datasets were identified as having a mutation in ddn, only three of which were represented amongst the afore-mentioned 145 assembled genomes available in GenBank.

Taxonomic profiling of unassembled sequencing data

Kraken [63], an ultrafast sequence classification tool, was used to assign taxonomic labels to the unassembled M. tuberculosis genomes from the SRA. Sequence reads for each genome from the SRA were screened against a NCBI refseq database (https://www.ncbi.nlm.nih.gov/refseq) using Kraken version 0.10.5 to identify contaminants. Of the 244 M. tuberculosis genomes with sequence polymorphisms in ddn, 67 (27.5%) were found to be contaminated with DNA from unidentified organisms or from Mycobacterium species other than M. tuberculosis (S2 Table). These genomes were excluded from all further analysis.

Phylogenetic reconstruction of M. tuberculosis strains with sequencing polymorphisms in Ddn

To determine the phylogenetic distribution of ddn alleles in M. tuberculosis we carried out a phylogenomic analysis of M. tuberculosis strains with sequence polymorphism in ddn. Firstly, unassembled genomes from the SRA were assembled de novo using SPAdes version 3.9.0 [64]. Next, all 322 complete and draft genomes (S3 Table) were aligned with Parsnp version 1.2 [65] using the genome of H37Rv as a reference to produce a core genome alignment of 848,476 bp. Core genome SNPs were identified and recombinant regions removed using Gubbins version 2.2.0 [66]. Finally, a maximum-likelihood phylogenetic tree was estimated using RAxML version 8.2.9 [67] under the GTRGAMMA nucleotide substitution model.

Evaluation of ddn mutant fitness in low-dose aerosol mouse infection model

Isogenic pretomanid-resistant mutants were selected in the M. tuberculosis H37Rv strain background and subjected to whole genome sequencing, as previously described [34]. In the present study, female BALB/c mice were aerosol-infected with the H37Rv parent strain or one of the following isogenic ddn mutants, implanting approximately 2 log10 CFU/lung in each case: L64P, M1T, or IS6110 insertion at D108. Three mice from each group were sacrificed on the day after infection to confirm the number of CFU implanted. The multiplication and persistence of each strain was evaluated by sacrificing 3 mice/group and quantifying the CFU counts by plating serial dilutions of lung homogenate on Middlebrook selective 7H11 agar at 2, 4, 8, 13 and 23 weeks post-infection.

Evaluation of ddn mutant fitness during and after hypoxic stress

The H37Rv parent and isogenic ddn mutants used for mouse infection were inoculated into tubes containing Dubos Tween Albumin Broth (BD Difco) supplemented with the hypoxia indicator dye, methylene blue (500 μg/ml). Tubes were sealed with rubber stoppers and incubated at 37 °C with slow magnetic stirring. The reduction and decolorization of methylene blue served as a visual indicator of oxygen levels corresponding to NRP stage 2 [68,69], which was attained in all cultures by day 24. Aliquots were removed for quantitative culture on 7H11 agar after 15, 24, 29, 43 and 57 days of incubation. Samples were collected by piercing the rubber stoppers with a syringe and sampling did not introduce atmospheric oxygen into the tubes, as shown by maintenance of the decolorized dye state. After day 57, each strain was returned to normoxic conditions by transferring the contents of the tubes to sterile Erlenmeyer flasks and incubated at 37 °C in a shaker for another 14 days [70]. Aliquots were removed for quantitative culture after 7 and 14 days of normoxia.

DDN marked deletion construction

To generate the DDN Hygromycin marked deletion mutant we used a high throughput phage-based method based on Jain et al (2014) [71]. The phasmid construct for this deletion was a gift from William R. Jacobs. The DDN deletion phasmid DNA was electroporated into electrocompetent M. smegmatis. The M. smegmatis was recovered for two hours in 7H9, then split and mixed into two 7H9 agar overlay mixtures (7H9 powder + 0.3% agar) and spread on 7H11 plates. The plates were incubated for two to three days at 30°C to allow phage growth and lysis. Resulting plaques were picked using a pasteur pipette and resuspended in MP buffer (50 mM Tris, 150 mM NaCl, 10 mM MgSO4, 2 mM CaCl2, pH 7.6–7.8). The MP buffer containing the phage was split and mixed with two M. smegmatis cultures, followed by mixing into two 7H9 agar overlay mixtures. The overlay mixes were spread on 7H11 plates and examined for plaques after two to three days at 30°C. This resulted in an increased number of plaques, which were harvested by pipetting MP buffer on the of the plate and shaking at 100 rpm for 3 hours. The MP buffer containing the phage was harvested from the plate with a syringe before being passed through a 20 μm filter to remove the bacteria. The resulting sterile MP buffer with phage titre was quantified by serial dilution followed by spot plating onto dry M. smegmatis overlay plates. The phage amplification was repeated until a titre of 109, 1010 phage was achieved. mc2 6230 Mycobacterium tuberculosis was grown up in ink wells to an optical density (600 nm) of 0.6 to 0.8, before harvest at 3,126 x g. Cells were then transduced by resuspension in the MP buffer containing the phage titre at 37 °C (to prevent phage lysing the cell) overnight followed by plating on 7H11 with 50 μg/mL hygromycin. Colonies with hygromycin resistance were indicative of successful double homologous recombination with the gene being replaced with hygR-sacB DNA fragment. The Δddn mutant was screened for resistance to delamanid. The Δddn mutant was confirmed through whole genome sequencing (Illumina 150 bp paired end reads) with reads mapped using the Bowtie2 plugin for Geneious to the wild type H37RV genome (S1a Fig) and a H37RV genome with the expected hygromycin marker (S1b Fig).

Competition assays

For the competition assay aerobically grown cultures the Δddn mutant and wild type mc2 6230 were grown in 7H9 in separate inkwells and monitored by spot platting 5 μL of a serial dilution in triplicate on 7H11 agar plates. At 7 days a new combined inkwell was inoculated with 100μL of each strain. For the competition assay from hypoxically grown cultures the Δddn mutant and wild type mc2 6230 were grown in separate 30 mL serum vials in 7H9 and monitored by spot platting 5 μL of a serial dilution in triplicate on 7H11 agar plates. Oxygen depletion was estimated by de-colourisation of 1.5 μg/mL methylene blue which occurred between 17–27 days. At 37 days 100 μL of the Δddn mutant and wild type mc2 6230 was taken and used to inoculate a single inkwell with 10 mL 7H9. At 70 days 200 μL of the Δddn mutant and wild type mc2 6230 was taken and used to inoculate a single inkwell with 10 mL 7H9. The combined inkwells with both strains were assayed for the proportion of the Δddn mutant by spot platting 5 μL of a serial dilution in triplicate on 7H11 agar plates with and without 50 μg/mL hygromycin. The Δddn mutant proportion of total cells was established from the number of colonies on the hygromycin agar plates and expressed as a percentage of total colonies on the plates with no hygromycin.

Supporting information

S1 Fig.

(A) Map of the Δddn mutant reads using Bowtie 2 to H37Rv genome. (B) Map of the Δddn mutant reads using Bowtie 2 to H37Rv genome with expected hygromycin marker.

https://doi.org/10.1371/journal.ppat.1008287.s001

(TIF)

S1 Table. Distribution of Ddn alleles in the major TB lineages.

https://doi.org/10.1371/journal.ppat.1008287.s002

(XLSX)

S2 Table. Taxonomic classification of sequence reads from bacterial genome data sets described as M. tuberculosis on the sequence read archive (SRA) using Kraken: Red shading indicates genome datasets contaminated with DNA from unknown organisms (Reads classified (%)) or datasets where the majority of reads, identified as genus Mycobacterium, match to a Mycobacterium species other than M. tuberculosis.

https://doi.org/10.1371/journal.ppat.1008287.s003

(XLSX)

S3 Table. M. tuberculosis genome datasets used in this study.

https://doi.org/10.1371/journal.ppat.1008287.s004

(XLSX)

Acknowledgments

We acknowledge Thomas Cuddihy (QFAB Bioinformatics) for assistance with genome data retrieval. We thank William R. Jacobs Jr. for the gift of the phasmid used to make the Δddn mutant. We thank Sebastien Gagneux for the gift of the M. tuberculosis N0008 strain.

Citation: Lee BM, Harold LK, Almeida DV, Afriat-Jurnou L, Aung HL, Forde BM, et al. (2020) Predicting nitroimidazole antibiotic resistance mutations in Mycobacterium tuberculosis with protein engineering. PLoS Pathog 16(2): e1008287. https://doi.org/10.1371/journal.ppat.1008287

References

1. WHO. Global Tuberculosis Report 2016. Cdc 2016. 2016; 214. ISBN 978 92 4 156539 4

2. Matsumoto M, Hashizume H, Tomishige T, Kawasaki M, Tsubouchi H, Sasaki H, et al. OPC-67683, a nitro-dihydro-imidazooxazole derivative with promising action against tuberculosis in vitro and in mice. PLoS Med. 2006;3: 2131–2144. pmid:17132069

3. Diacon AH, Pym A, Grobusch M, Patientia R, Rustomjee R, Page-Shipp L, et al. The Diarylquinoline TMC207 for Multidrug-Resistant Tuberculosis. N Engl J Med. 2009;360: 2397–2405. pmid:19494215

4. Gupta R, Gao M, Cirule A, Xiao H, Geiter LJ, Wells CD. Delamanid for Extensively Drug-Resistant Tuberculosis. N Engl J Med. 2015;373: 291–292. pmid:26176402

5. Stover CK, Warrener P, VanDevanter DR, Sherman DR, Arain TM, Langhorne MH, et al. A small-molecule nitroimidazopyran drug candidate for the treatment of tuberculosis. Nature. 2000;405: 962–966. pmid:10879539

6. Singh R, Manjunatha U, Boshoff HIM, Ha YH, Niyomrattanakit P, Ledwidge R, et al. PA-824 kills nonreplicating Mycobacterium tuberculosis by intracellular NO release. Science. 2008;322: 1392–5. pmid:19039139

7. Cellitti SE, Shaffer J, Jones DH, Mukherjee T, Gurumurthy M, Bursulaya B, et al. Structure of Ddn, the deazaflavin-dependent nitroreductase from Mycobacterium tuberculosis involved in bioreductive activation of PA-824. Structure. 2012;20: 101–112. pmid:22244759

8. Manjunatha U, Boshoff HIM, Barry CE. The mechanism of action of PA-824: Novel insights from transcriptional profiling. Commun Integr Biol. 2009;2: 215–218. pmid:19641733

9. Greening C, Ahmed FH, Mohamed AE, Lee BM, Pandey G, Warden AC, et al. Physiology, Biochemistry, and Applications of F420—and Fo -Dependent Redox Reactions. Microbiol Mol Biol Rev. 2016;80: 451–493. pmid:27122598

10. Purwantini E, Mukhopadhyay B. Conversion of NO2 to NO by reduced coenzyme F420 protects mycobacteria from nitrosative damage. Proc Natl Acad Sci U S A. 2009;106: 6333–6338. pmid:19325122

11. Hasan MR, Rahman M, Jaques S, Purwantini E, Daniels L. Glucose 6-phosphate accumulation in mycobacteria: implications for a novel F420-dependent anti-oxidant defense system. J Biol Chem. 2010;285: 19135–19144. pmid:20075070

12. Gurumurthy M, Rao M, Mukherjee T, Rao SPS, Boshoff HI, Dick T, et al. A novel F420-dependent anti-oxidant mechanism protects Mycobacterium tuberculosis against oxidative stress and bactericidal agents. Mol Microbiol. 2013;87: 744–755. pmid:23240649

13. Ahmed FH, Carr PD, Lee BM, Afriat-Jurnou L, Mohamed AE, Hong NS, et al. Sequence-Structure-Function Classification of a Catalytically Diverse Oxidoreductase Superfamily in Mycobacteria. J Mol Biol. 2015;427: 3554–3571. pmid:26434506

14. de Souza GA, Leversen NA, Målen H, Wiker HG. Bacterial proteins with cleaved or uncleaved signal peptides of the general secretory pathway. J Proteomics. 2011;75: 502–510. pmid:21920479

15. Hoffmann H, Kohl TA, Hofmann-Thiel S, Merker M, Beckert P, Jaton K, et al. Delamanid and bedaquiline resistance in MTB ancestral Beijing genotype causing extensive drug-resistant TB in a tibetan refugee. Am J Respir Crit Care Med. 2016;193: 337–340. pmid:26829425

16. Bloemberg GV., Keller PM, Stucki D, Trauner A, Borrell S, Latshang T, et al. Acquired Resistance to Bedaquiline and Delamanid in Therapy for Tuberculosis. N Engl J Med. 2015;373: 1986–1988. pmid:26559594

17. Choi KP, Bair TB, Bae YM, Daniels L. Use of transposon Tn5367 mutagenesis and a nitroimidazopyran-based selection system to demonstrate a requirement for fbiA and fbiB in coenzyme F420 biosynthesis by Mycobacterium bovis BCG. J Bacteriol. 2001;183: 7058–7066. pmid:11717263

18. Manjunatha UH, Boshoff H, Dowd CS, Zhang L, Albert TJ, Norton JE, et al. Identification of a nitroimidazo-oxazine-specific protein involved in PA-824 resistance in Mycobacterium tuberculosis. Proc Natl Acad Sci U S A. 2006;103: 431–436. pmid:16387854

19. Hoffmann H, Borroni E, Schena E, Nedialkova L, Hofmann-Thiel S, Cirillo D. Delamanid susceptibility testing of Mycobacterium tuberculosis using the resazurin microtitre assay and the BACTEC™ MGIT™ 960 system—authors’ response. J Antimicrob Chemother. 2016;71: 3625–3625. pmid:27621178

20. Haver H, Chua A, Ghode P. Mutations in genes for the F420 biosynthetic pathway and a nitroreductase enzyme are the primary resistance determinants in spontaneous in vitro-selected PA-824-. Antimicrob agents. 2015;59: 5316–5323.

21. Fujiwara M, Kawasaki M, Hariguchi N, Liu Y, Matsumoto M. Mechanisms of resistance to delamanid, a drug for Mycobacterium tuberculosis. Tuberculosis. 2018;108: 186–194. https://doi.org/10.1016/j.tube.2017.12.006 pmid:29523322

22. Schena E, Nedialkova L, Borroni E, Battaglia S, Cabibbe AM, Niemann S, et al. Delamanid susceptibility testing of Mycobacterium tuberculosis using the resazurin microtitre assay and the BACTEC™ MGIT™ 960 system. J Antimicrob Chemother. 2016;71: 1532–1539. pmid:27076101

23. Cook GM, Greening C, Hards K, Berney M. Chapter One-Energetics of Pathogenic Bacteria and Opportunities for Drug Development. Adv Microb Physiol. 2014;65: 1–62. pmid:25476763

24. Sinha S, Kosalai K, Arora S, Namane A, Sharma P, Gaikwad AN, et al. Immunogenic membrane-associated proteins of Mycobacterium tuberculosis revealed by proteomics. Microbiology. 2005;151: 2411–2419. pmid:16000731

25. Upadhyay A, Kumar S, Rooker SA, Koehn JT, Crans DC, McNeil MR, et al. Mycobacterial MenJ: an oxidoreductase involved in menaquinone biosynthesis. ACS Chem Biol. 2018;13: 2498–2507. pmid:30091899

26. Bar-Even A, Milo R, Noor E, Tawfik DS. The moderately efficient enzyme: futile encounters and enzyme floppiness. Biochemistry. 2015;54: 4969–4977. pmid:26219075

27. Bar-Even A, Noor E, Savir Y, Liebermeister W, Davidi D, Tawfik DS, et al. The moderately efficient enzyme: evolutionary and physicochemical trends shaping enzyme parameters. Biochemistry. 2011;50: 4402–4410. pmid:21506553

28. Hards K, Rodriguez SM, Cairns C, Cook GM. Alternate quinone coupling in a new class of succinate dehydrogenase may potentiate mycobacterial respiratory control. FEBS Lett. 2019;593: 475–486. pmid:30675730

29. Holsclaw CM, Sogi KM, Gilmore SA, Schelle MW, Leavell MD, Bertozzi CR, et al. Structural characterization of a novel sulfated menaquinone produced by stf3 from Mycobacterium tuberculosis. ACS Chem Biol. 2008;3: 619–624. pmid:18928249

30. Hatfull GF, Jacobs WR Jr. Molecular genetics of mycobacteria. American Society for Microbiology Press; 2014.

31. Oyugi MA, Bashiri G, Baker EN, Johnson-Winters KL. Investigating the reaction mechanism of F420-dependent glucose-6-phosphate dehydrogenase from Mycobacterium tuberculosis: kinetic analysis of the wild-type and mutant enzymes. Biochemistry. 2016.

32. Ji B, Lefrançois S, Robert J, Chauffour A, Truffot C, Jarlier V. In vitro and in vivo activities of rifampin, streptomycin, amikacin, moxifloxacin, R207910, linezolid, and PA-824 against Mycobacterium ulcerans. Antimicrob Agents Chemother. 2006;50: 1921–1926. pmid:16723546

33. Upton AM, Cho S, Yang TJ, Kim Y, Wang Y, Lu Y, et al. In vitro and in vivo activities of the nitroimidazole TBA-354 against Mycobacterium tuberculosis. Antimicrob Agents Chemother. 2015;59: 136–144. pmid:25331696

34. Rifat D, Li S-Y, Ioerger T, Lanoix J-P, Lee J, Bashiri G, et al. Mutations in Rv2983 as a novel determinant of resistance to nitroimidazole drugs in Mycobacterium tuberculosis. bioRxiv. 2018; 457754.

35. Comas I, Coscolla M, Luo T, Borrell S, Holt KE, Kato-Maeda M, et al. Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans. Nat Genet. 2013;45: 1176–1182. pmid:23995134

36. Hart KM, Ho CMW, Dutta S, Gross ML, Bowman GR. Modelling proteins’ hidden conformations to predict antibiotic resistance. Nat Commun. 2016;7: 12965. pmid:27708258

37. Orencia MC, Yoon JS, Ness JE, Stemmer WPC, Stevens RC. Predicting the emergence of antibiotic resistance by directed evolution and structural analysis. Nat Struct Mol Biol. 2001;8: 238–242.

38. Mohamed AE, Condic-Jurkic K, Ahmed FH, Yuan P, O’Mara ML, Jackson CJ, et al. Hydrophobic shielding drives catalysis of hydride transfer in a family of F420H2-dependent enzymes. Biochemistry. 2016. pmid:27951661

39. Aharoni A, Gaidukov L, Khersonsky O, Gould SMQ, Roodveldt C, Tawfik DS. The’evolvability’of promiscuous protein functions. Nat Genet. 2005;37: 73–76. pmid:15568024

40. Hards K, Robson JR, Berney M, Shaw L, Bald D, Koul A, et al. Bactericidal mode of action of bedaquiline. J Antimicrob Chemother. 2015;70: 2028–2037. pmid:25754998

41. Lamprecht DA, Finin PM, Rahman MA, Cumming BM, Russell SL, Jonnala SR, et al. Turning the respiratory flexibility of Mycobacterium tuberculosis against itself. Nat Commun. 2016;7: 12393. pmid:27506290

42. Ai J-W, Ruan Q-L, Liu Q-H, Zhang W-H. Updates on the risk factors for latent tuberculosis reactivation and their managements. Emerg Microbes Infect. 2016;5: 1–8.

43. Gibson DG, Young L, Chuang R-Y, Venter JC, Hutchison Ca, Smith HO, et al. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat Methods. 2009;6: 343–5. pmid:19363495

44. Lapalikar GV., Taylor MC, Warden AC, Scott C, Russell RJ, Oakeshott JG. F 420H 2-dependent degradation of aflatoxin and other furanocoumarins is widespread throughout the Actinomycetales. PLoS One. 2012;7: e30114. pmid:22383957

45. Bashiri G, Squire CJ, Moreland NJ, Baker EN. Crystal structures of F420-dependent glucose-6-phosphate dehydrogenase FGD1 involved in the activation of the anti-tuberculosis drug candidate PA-824 reveal the basis of coenzyme and substrate binding. J Biol Chem. 2008;283: 17531–17541. pmid:18434308

46. Sambrook J, Fritsch EF, Maniatis T. Molecular cloning: a laboratory manual. Cold spring harbor laboratory press; 1989.

47. Jirapanjawat T, Ney B, Taylor MC, Warden AC, Afroze S, Russell RJ, et al. The redox cofactor F420 protects mycobacteria from diverse antimicrobial compounds and mediates a reductive detoxification system. Appl Environ Microbiol. 2016;82: 6810–6818. pmid:27637879

48. Purwantini E, Mukhopadhyay B, Spencer RW, Daniels L. Effect of temperature on the spectral properties of coenzyme F420 and related compounds. Anal Biochem. 1992;205: 342–350. pmid:1443583

49. Edgar RC. MUSCLE : multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2008;32: 1792–1797. pmid:15034147

50. Li W, Cowley A, Uludag M, Gur T, McWilliam H, Squizzato S, et al. The EMBL-EBI bioinformatics web and programmatic tools framework. Nucleic Acids Res. 2015;43: W580–W584. pmid:25845596

51. Trott O, Olson AJ. AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, EfficientOptimization, and Multithreading. J Comput Chem. 2010;31: 455–61. pmid:19499576

52. Morris G, Huey R. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J …. 2009;30: 2785–2791.

53. DeLano WL. The PyMOL molecular graphics system. http://pymol.org. 2002.

54. Irwin JJ, Shoichet BK. ZINC − A Free Database of Commercially Available Compounds for Virtual Screening ZINC—A Free Database of Commercially Available Compounds for Virtual Screening. J Chem Inf Model. 2005;45: 177–182. pmid:15667143

55. Wiegand I, Hilpert K, Hancock REW. Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat Protoc. 2008;3: 163–75. pmid:18274517

56. Heikal A, Hards K, Cheung CY, Menorca A, Timmer MSM, Stocker BL, et al. Activation of type II NADH dehydrogenase by Quinolinequinones mediates antitubercular cell death. J Antimicrob Chemother. 2016;71: 2840–2847. pmid:27365187

57. Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, et al. GenBank. Nucleic Acids Res. 2013;41: D36–D42. pmid:23193287

58. Steiner A, Stucki D, Coscolla M, Borrell S, Gagneux S. KvarQ: targeted and direct variant calling from fastq reads of bacterial genomes. BMC Genomics. 2014;15: 881. pmid:25297886

59. Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, et al. Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998;393: 537–544. pmid:9634230

60. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high‐quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7.

61. Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv Prepr arXiv13033997. 2013.

62. Garrison E, Marth G. Haplotype-based variant detection from short-read sequencing. arXiv Prepr arXiv12073907. 2012.

63. Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15: R46–R46. pmid:24580807

64. Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. J Comput Biol. 2012;19: 455–477. pmid:22506599

65. Treangen TJ, Ondov BD, Koren S, Phillippy AM. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 2014;15: 524. pmid:25410596

66. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015;43: e15–e15. pmid:25414349

67. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30: 1312–1313. pmid:24451623

68. Klinkenberg LG, Sutherland LA, Bishai WR, Karakousis PC. Metronidazole lacks activity against Mycobacterium tuberculosis in an in vivo hypoxic granuloma model of latency. J Infect Dis. 2008;198: 275–283. pmid:18491971

69. Wayne LG. In vitro model of hypoxically induced nonreplicating persistence of Mycobacterium tuberculosis. Mycobacterium tuberculosis protocols. Springer; 2001. pp. 247–269.

70. Sherrid AM, Rustad TR, Cangelosi GA, Sherman DR. Characterization of a Clp protease gene regulator and the reaeration response in Mycobacterium tuberculosis. PLoS One. 2010;5: e11622. pmid:20661284

71. Jain P, Hsu T, Arai M, Biermann K, Thaler DS, Nguyen A, et al. Specialized Transduction Designed for Precise High-Throughput Unmarked Deletions in Mycobacterium tuberculosis. Rubin EJ, editor. MBio. 2014;5: e01245–14. pmid:24895308