- Open Access
Three different mutations in the DNA topoisomerase 1B in Leishmania infantum contribute to resistance to antitumor drug topotecan
Parasites & Vectors volume 14, Article number: 438 (2021)
The evolution of drug resistance is one of the biggest challenges in leishmaniasis and has prompted the need for new antileishmanial drugs. Repurposing of approved drugs is a faster and very attractive strategy that is gaining supporters worldwide. Different anticancer topoisomerase 1B (TOP1B) inhibitors have shown strong antileishmanial activity and promising selective indices, supporting the potential repurposing of these drugs. However, cancer cells and Leishmania share the ability to become rapidly resistant. The aim of this study was to complete a whole-genome exploration of the effects caused by exposure to topotecan in order to highlight the potential mechanisms deployed by Leishmania to favor its survival in the presence of a TOP1B inhibitor.
We used a combination of stepwise drug resistance selection, whole-genome sequencing, functional validation, and theoretical approaches to explore the propensity of and potential mechanisms deployed by three independent clones of L. infantum to resist the action of TOP1B inhibitor topotecan.
We demonstrated that L. infantum is capable of becoming resistant to high concentrations of topotecan without impaired growth ability. No gene deletions or amplifications were identified from the next-generation sequencing data in any of the three resistant lines, ruling out the overexpression of efflux pumps as the preferred mechanism of topotecan resistance. We identified three different mutations in the large subunit of the leishmanial TOP1B (Top1BF187Y, Top1BG191A, and Top1BW232R). Overexpression of these mutated alleles in the wild-type background led to high levels of resistance to topotecan. Computational molecular dynamics simulations, in both covalent and non-covalent complexes, showed that these mutations have an effect on the arrangement of the catalytic pentad and on the interaction of these residues with surrounding amino acids and DNA. This altered architecture of the binding pocket results in decreased persistence of topotecan in the ternary complex.
This work helps elucidate the previously unclear potential mechanisms of topotecan resistance in Leishmania by mutations in the large subunit of TOP1B and provides a valuable clue for the design of improved inhibitors to combat resistance in both leishmaniasis and cancer. Our data highlights the importance of including drug resistance evaluation in drug discovery cascades.
Despite 1 million new cases of leishmaniasis being declared every year, there is still no effective vaccine available for humans. Moreover, available treatments for this major neglected tropical disease (NTD) are limited and outdated . Clinical manifestations vary in severity and include cutaneous, mucocutaneous, and visceral leishmaniasis. The latter form, caused by Leishmania infantum, is fatal within 2 years if left untreated . Leishmaniasis is also responsible for a significant health, psychosocial, and economic burden around the world [3, 4]. Due to the limited pharmacopeia, most available antileishmanial drugs are used in both humans and dogs. Of note, dogs constitute the main reservoir for the zoonotic life cycle of Leishmania [5, 6].
The rapid emergence and spread of resistant parasite strains has prompted the need for new intervention pathways and antileishmanial drugs. In this way, repurposing of approved drugs has become a very attractive strategy to tackle leishmaniases and is gaining supporters worldwide [1, 7, 8]. DNA topoisomerases (TOP) have garnered attention since their discovery by James Wang in 1971 . TOP are key enzymes for many essential biological functions such as DNA replication, transcription, recombination, DNA repair, and DNA segregation . TOP enzymes are preferential targets against rapid-dividing, highly proliferative eukaryotic cells, such as those responsible for various types of cancer (e.g. ovarian cancer or metastatic carcinoma of the colon) [11, 12]. Type 1B TOP (TOP1B) are ATP-independent enzymes that relax topological tensions in supercoiled DNA during replication and transcription processes mediated by DNA and RNA polymerases. TOP1B enzymes produce topological changes (e.g. relaxation) in supercoiled DNA substrates through a five-step coordinated sequence: (i) DNA binding; (ii) DNA cleavage and formation of a transient protein–DNA covalent interaction; (iii) DNA controlled rotation to release torsional stress; (iv) DNA religation; and (v) DNA release. Eukaryotic TOP1B are monomeric enzymes, except in kinetoplastid organisms (e.g. Leishmania) that rely on heterodimeric TOP1B enzymes, with coding genes located in different chromosomes, 34 and 4, respectively [13, 14]. The major structural and functional differences between Homo sapiens TOP1B (hTOP1B) and L. infantum TOP1B (LiTOP1B) [15, 16], coupled with Leishmania’s rapid proliferation rate, offer great potential for selective chemotherapy . Different anticancer TOP1B inhibitors, including the water-soluble camptothecin derivative topotecan (TPT) currently used to treat ovarian carcinoma , have been successfully tested against Leishmania parasites in in vitro, ex vivo, and in vivo murine models, demonstrating strong antileishmanial activity and promising selective indices, supporting the potential repurposing of these drugs [8, 17, 18]. However, tumor cells and Leishmania parasites share an undesirable and important feature in being prone to becoming resistant to drug treatment. L. infantum relies on DNA copy number variations (CNVs) for regulating the expression of drug targets [19, 20]) and drug resistance genes [20, 21]. In addition to CNVs, single-nucleotide polymorphisms (SNPs) in drug target genes or in transporters can lead to drug resistance without the need for altering gene content [22,23,24]. Markedly, recent studies point to the need for experimental generation of drug resistance to promising compounds in order to clinically evaluate and eventually circumvent the phenomenon [25,26,27].
The goal of this study was to complete a whole-genome exploration of the effects caused by prolonged exposure to TPT to highlight the different potential drug resistance mechanisms (changes in ploidy, CNVs, SNPs, etc.) deployed by L. infantum to favor its survival in the presence of an FDA-approved TOP1B inhibitor. This information is critical to the fight against drug-resistant Leishmania parasites, furthering knowledge of the mechanisms these parasites use to persist in the presence of drugs, as well as serving as a foundation for improved drug repurposing strategies.
Leishmania infantum (MHOM/MA/67/ITMAP-263) wild-type (WT) promastigotes and mutants resistant to > 700 μM TPT (TPT700.1, TPT700.2, and TPT700.3) generated in vitro in a stepwise manner were grown in M199 medium at 25 °C supplemented with 10% fetal bovine serum, and 5 μg/ml of hemin at pH 7.0. In addition, 700 μM of TPT (topotecan hydrochloride hydrate, Sigma-Aldrich) was added to the media for the maintenance of the endpoint mutants. Growth curves were performed in 25 cm2 cell culture flasks by seeding 1 × 106 parasites/ml, and the number of parasites was determined daily—up to 7 days—by manual counting using the Neubauer hemocytometer. Growth assays were performed with at least three biological replicates from independent cultures (n = 3), each of which included three technical replicates.
Three L. infantum WT independent clones, which were obtained by plating the MHOM/MA/67/ITMAP-263 strain onto solid M199, were independently selected in 25 cm2 flasks containing 5 ml M199 medium supplemented with 10% heat-inactivated fetal bovine serum (FBS) and 5 μg/ml hemin in the presence of increasing TPT concentrations, as described previously . Briefly, the stepwise drug selection ranged from 1 × the EC50 of TPT (24 μM as determined in the present study) up to 16 × the EC50 of TPT (384 μM), with a twofold increase in drug concentration (24, 48, 96, 192, and 384 μM) every three sub-culturing passages. Topotecan hydrochloride hydrate (Sigma-Aldrich, St. Louis, MO, USA) was used as the source of TPT.
Drug susceptibility assays
Antileishmanial values in promastigotes were determined by monitoring the growth of parasites after 72 h of incubation at 25 °C in the presence of increasing concentrations of TPT, by measuring A600 using a Cytation 5 multimode reader (BioTek, USA). Drug efficacy assays were performed with at least three biological replicates from independent cultures (n = 3). EC50 values were calculated based on dose–response curves analyzed by nonlinear regression with GraphPad Prism 8.4.3 software (GraphPad Software, La Jolla, CA, USA). Statistical analyses were performed using unpaired two-tailed t-tests. A P-value < 0.05 was considered statistically significant.
Whole-genome sequencing of TPT-resistant mutants
Comparative whole-genome sequencing (WGS) was performed as described previously . Briefly, genomic DNA was prepared from a mid-log phase clonal culture of each TPT-resistant mutant. DNA was quantified fluorometrically, and 50 ng of material was used for library preparation using a Nextera™ DNA Sample Preparation Kit (Illumina) according to the manufacturer’s instructions. The size distribution of Nextera™ libraries was validated using an Agilent 2100 Bioanalyzer and High Sensitivity DNA chips (Agilent Technologies). Sequencing libraries were quantified with the QuantiFluor® dsDNA System and sequenced using an Illumina MiSeq platform with 250-nucleotide paired-ends reads. An average genome coverage of over 50-fold was obtained for the mutants. This approach allowed identification of SNPs when compared with the reference genome sequence of L. infantum JPCM5 (TriTrypDB v9.0)  and L. infantum 263 WT . Sequence reads were aligned to the L. infantum JPCM5 genome and L. infantum 263 WT using the software bwa-mem . The maximum number of mismatches was four, the seed length was 32, and two mismatches were allowed within the seed. Read duplicates were marked using Picard (http://broadinstitute.github.io/picard), and GATK was applied for InDel realignment and SNP and InDel discovery in the three TPT mutants. PCR amplification and conventional DNA sequencing verified SNPs of interest detected by WGS. CNVs were derived from read depth coverage by comparing the coverage of uniquely mapped reads between each of the three TPT mutants and the L. infantum 263 WT in 5 kb non-overlapping genomic windows for the 36 chromosomes (normalized to the total number of uniquely mapped reads for each strain) . The sequence data for the L. infantum TPT-resistant mutants is available at the NCBI BioProject (https://www.ncbi.nlm.nih.gov/bioproject/) under study accession PRJNA647847 and sample accessions SAMN15599759, SAMN15599760, and SAMN15599761, corresponding to TPT-resistant clones TPT700.1, TPT700.2, and TPT700.3, respectively.
Quantitative real-time RT-PCR
RNAs from the WT and the three mutants were extracted using the RNeasy Mini Kit (Qiagen) according to the manufacturer's recommendations. The cDNA was synthesized using the iScript™ Reverse Transcription Supermix (Bio-Rad), and amplified in the iTaq™ universal SYBR® Green Supermix Kit (Bio-Rad) using a CFX Opus Real-Time PCR System (Bio-Rad). The expression levels of LinJ.34.3220 (Fw: 5′-CGACTTCGAGCCCATTTATCA-3′; Rv: 5′-ACCTTCAGCTTGCCCATTAG-3′) and LinJ.04.0070 (Fw: 5′-CCCTCCGTCAAGAAAGTTGT-3′; Rv: 5′-CTCTTCCACATTGCCCAGT-3′) were derived from three technical and three biological replicates and were normalized to constitutively expressed mRNA encoding glyceraldehyde-3-phosphate dehydrogenase (LinJ.36.2480; Fw: 5′-GTACACGGTGGAGGCTGTG-3′; Rv: 5′-CCCTTGATGTGGCCCTCGG-3′).
DNA constructs and nucleofection
The WT allele of LinJ.34.3220 (DNA topoisomerase IB large subunit) and its three mutant variants were amplified from genomic DNA derived from either L. infantum or mutants TPT700.1, TPT700.2, and TPT700.3 (mutations T560A, G572C, and T724C, respectively) using compatible primer pairs (Fw: 5′-XbaI-ATGAAGGTGGAGAA-3′; Rv: 5′-HindIII-TACACCCTCAAAGC-3′). PCR fragments were ligated into a pGEM T-easy vector (Promega) to confirm the quality of the insert by Sanger sequencing. PCR fragments were then cloned in the Leishmania expression vector pSP72αhygα, which contains the gene hygromycin phosphotransferase (hyg), a selectable marker in Leishmania . A total of 20 μg of plasmid DNA for episomal expression, either the empty vector (mock) or carrying the genes of interest, were delivered into L. infantum WT promastigotes by nucleofection, as described previously . Selection was achieved in the presence of a final concentration of 300 μg/ml hygromycin.
Computational methods for the in silico study of Leishmania topoisomerase 1B
The three-dimensional structure of the L. donovani TOP1LS complex (PDB ID: 2B9S ) contains residues 27–456 and 221–262 of the large and small subunits, respectively. Residues missing in the crystal structure, 427–430 of the large subunit, have been modelled as reported in Roy et al. . Following this procedure, a non-covalent interaction of the LiTOP1 bi-subunit with a 22-base-pair DNA substrate was modelled using the program UCSF Chimera  by fitting the structure of the protein backbone atoms on the coordinates of the human enzyme trapped in complex with the DNA (PDB ID 1A36). The same procedure was carried out in order to model an LiTOP1-DNA covalent interaction, using the human TOP1-DNA-topotecan ternary complex, extrapolated from the PDB structure 1K4T .
Once the two systems (wild-type LiTOP1–DNA covalent and non-covalent complexes) were obtained, the program UCSF Chimera  was used to introduce the three point mutations identified in our whole-genome sequencing experiments (F187Y, G191A, and W242R). All systems were subjected to molecular dynamics (MD) simulations aimed at understanding the structural/dynamical effect of the mutations, thus explaining their impact on TPT sensitivity. Briefly, each covalent and non-covalent complex was placed in a triclinic simulative box filled with a water molecule TIP3P model , and the resulting systems were rendered electroneutral by the addition of sodium counterions; this step, i.e. the topology building, was performed using the Amber 14 all-atoms force field, highly suitable for simulating nucleic acids . The systems were first subjected to a round of minimization of 10,000 steps of steepest descent followed by 10,000 steps of conjugate gradient. Relaxation of water molecules and thermalization in NPT environment were carried out for 1.2 ns in time steps of 1 fs. Six runs of 200 ps each were carried out while increasing the temperature by 50 K at each step, ranging from 50 to 300 K.
The systems were then simulated with a 2 fs time step for 300 ns in periodic boundary conditions, using a cut-off of 8 Å for the evaluation of short-range non-bonded interactions and the particle mesh Ewald method for long-range electrostatic interactions . The temperature was kept constant at 300 K with Langevin dynamics , whereas pressure was fixed at 1 atmosphere through the Langevin piston method . The bond lengths of solute and water molecules were restrained with the SHAKE  and SETTLE  algorithms, respectively. Atomic positions were saved every 250 steps (i.e. 0.5 ps) for analysis with the Gromacs 4.6 package .
In vitro resistance selection and characterization of the TPT-resistant mutants
Selection for resistance to TPT was performed to evaluate the possible propensity for drug resistance to this anticancer TOP1B inhibitor if repurposed against visceral leishmaniasis. Experimental concentration–response assays with L. infantum WT promastigotes revealed a TPT-sensitive phenotype within the low micromolar range (Fig. 1a; EC50 = 25 µM). Selection of independent TPT-resistant mutants began at 1 × EC50 up to 16 × EC50 of the WT parental line. The selection procedure was fast (three subsequent passages per drug concentration), and cells rapidly adapted to growing concentrations of TPT. All three clones reached the final selection step at 16 × EC50 (passage 3). At this point, clones were re-evaluated in terms of growth ability and sensitivity to TPT (Fig. 1a). All three clones were able to survive concentrations higher than 700 µM (ca. 28 × EC50 of the WT), demonstrating very high EC50 values (612, 519, and 602 µM for TPT700.1, TPT700.2, and TPT700.3, respectively) compared to the parental WT strain (~ 25 µM). Once the final selection step was reached, the growth rate of each TPT mutant was determined in order to verify their fitness compared to the WT parent line. As depicted in Fig. 1b, the TPT-resistant phenotype experimentally induced in promastigotes did not lead to any significant difference in the growth rate of the parasite when compared to the WT cells (Fig. 1b). Resistant phenotypes were stable after growth in media lacking TPT for over 60 passages (Fig. 1c).
Whole-genome sequencing was performed to compare the three independent TPT-resistant strains to the WT line to pinpoint and further explore potential mechanisms underlying the resistant phenotype. WGS was conducted by Illumina next-generation sequencing on the three independent L. infantum TPT700-resistant lines selected at 16 × EC50, as well as the isogenic L. infantum WT counterpart line. For all strains, this produced genome assemblies of 31 Mb with a coverage depth of at least 50-fold.
No significant deletions or amplifications were identified; however, the TPT-resistant lines demonstrated aneuploidy compared to the WT parental line. Two cases of reduction in ploidy were observed in the TPT-resistant mutants (Fig. 2a; Additional file 1: Dataset S1). These losses affected chromosome 12 in clone TPT700.3 (Fig. 2a, b) and chromosome 32 in the three resistant mutants (Fig. 2c). The log2 TPT700.3/WT read ratio for chromosome 12 (Fig. 2b) was close to −0.5, which should represent a loss of one allele compared to WT parasites (going from four to three chromosome copies). In the same way, the log2 TPT/WT read ratio (Fig. 2c) for chromosome 32 was close to −1, pinpointing a trisomic-to-disomic shift in all three TPT-resistant mutants (Fig. 2c), which could correlate with reduced TPT sensitivity. No cases of supernumerary chromosomes were observed in the TPT-resistant mutants (Fig. 2a; Additional file 1: Dataset S1).
A search for point mutations revealed several homozygous and heterozygous SNPs in the different TPT-resistant mutants (Additional file 2: Dataset S2). Nineteen common genes were mutated in the three TPT-resistant mutants (Table 1). In order to find potential SNP candidates, the different clones were examined for SNPs occurring in the same ORF but at different positions (to avoid mutations that are most likely to occur as natural polymorphisms in the parental strain). SNPs fulfilling these characteristics occurred in eight genes, including two coding for hypothetical proteins (LinJ.19.1680 and LinJ.19.1690), three proteophosphoglycan-coding genes (LinJ.35.0490, LinJ.35.0510, and LinJ.35.0550) often found mutated in our various sequencing screens, LinJ.24.1120 (encoding a putative pre-mRNA splicing factor), LinJ.34.0710 (encoding a putative flagellar attachment zone protein), and the gene coding for the large subunit of the heterodimeric LiTOP1B (LinJ.34.3220). The three different mutations identified in the LiTOP1B ORF by next-generation sequencing (NGS) and confirmed by targeted PCR were homozygous mutations top1BF187Y and top1BW232R for TPT700.1 and TPT700.3 mutants, respectively, and heterozygous SNP top1BG191A for TPT700.2. The levels of top1B (LinJ.34.3220 and LinJ.04.0070) mRNA did not vary between the different clones and the WT parental cell line (Additional file 3: Table S1).
Impact of LiTOP1B mutant variants in L. infantum TPT drug resistance
Because top1B is an essential gene in Leishmania, our preferred method for studying the role of TOP1B mutations was episomal transfection of WT and mutated forms of the gene in the WT background (Fig. 3). To this end, the WT copy as well as the three mutated forms of top1B were recovered by PCR, confirmed by Sanger sequencing, and episomally overexpressed in the WT background. Markedly, while episomal overexpression of gene LinJ.34.3220 (WT or mutated) increases the potential amount of TOP1B large subunit available, the amount of functional TOP1B heterodimer in the cell remains stable (independent of the overexpression), as it is limited by the basal amount of TOP1B small subunit. An increase in the total amount of functional WT TOP1B heterodimer in the presence of TPT would lead to an increase in the sensitivity to TPT. Overexpression of the leishmanial top1BWT allele did not confer TPT sensitivity in the WT cell line when compared with the mock control, confirming that the amount of TOP1B heterodimer remains stable. On the other hand, if the SNPs found in the mutated variants of the TOPIB large subunit play a role in TPT resistance and are able to bind to the small subunit (competing with the WT form), these would lead to a reduced sensitivity towards TPT. As depicted in Fig. 3, overexpression of all three mutated variants of the top1B gene resulted in a significant decrease in TPT sensitivity in the WT background: top1BF187Y (2.62-fold; P < 0.0001), top1BG191A (1.96-fold; P = 0.0006), and top1BW232R (2.28-fold; P = 0.0001).
Structural effect of the SNPs and impact on TPT sensitivity
The mechanism of inhibition exerted by TPT against TOP1B is strictly dependent on the enzyme catalytic cycle. Once the DNA cleavage has occurred, TPT is able to intercalate the DNA bases and stabilize the protein–DNA covalent complex, thus blocking the cycle and stalling the enzyme. If one of the catalytic steps is compromised, this can affect TPT sensitivity.
In order to evaluate the impact of the point mutations on the structure of the protein and its interaction with the DNA substrate for a correct progression of the catalytic cycle, we performed molecular dynamics (MD) simulations of the WT and mutant enzymes (F187Y, W232R, and G191A). In particular, we carried out MD simulations both of the WT and mutant non-covalent and covalent TOP1B–DNA complexes (for a total of eight systems), allowing detection of the mutations’ effect in any step of the two phases of the catalytic cycle, i.e. cleavage and religation. This gives important insight into the mechanism(s) underlying resistance due to the fact that these mutations are located in close proximity to the catalytic pentad and to the TPT binding site (Fig. 4).
In both groups of simulations (covalent and non-covalent complexes), the mutation of each residue has an effect on the arrangement of the catalytic pentad (R314, K352, R410, and H453 of the large subunit and Y222 of the small subunit) (Figs. 5, 6) and on the interaction of these residues with the surrounding amino acids and the DNA substrate (Tables 2, 3). This likely alters the cleavage/religation equilibrium, which is at the basis of an altered drug sensitivity. Of note, the profile of interaction at the cleavage site (−1/+1 bases) of the protein–DNA in the non-covalent complex is hardly affected in all mutants, making us hypothesize that due to the improper arrangement of the catalytic pentad, the reaction rate is affected. Similarly, the catalytic pentad arrangement is impacted, in particular at the level of residue K352, which is crucial for the religation reaction to proceed and for the stabilization of TPT in the DNA binding pocket.
Topotecan is a semi-synthetic, water-soluble analogue of camptothecin (CPT), as well as the first FDA-approved oral TOP1B inhibitor for the treatment of several types of cancer. CPT derivatives have been repeatedly suggested as a good source of repurposed drugs for the treatment of a variety of infectious diseases caused by protozoan parasites [1, 13, 46]. With the alarming decrease in effectiveness of first-line drugs in areas where L. infantum is endemic, repurposed drugs could represent a faster solution at lower cost . However, the predisposition of Leishmania to develop drug resistance should be addressed when repurposing a drug [1, 25, 27]. In this study, we used a combination of stepwise drug resistance selection, whole-genome sequencing, and theoretical approaches to explore the propensity of and potential mechanisms deployed by three independent clones of L. infantum to resist the activity of the TOP1B inhibitor TPT. One of the major strengths of this approach is that both the parent cell line and three directly derived drug-resistant lines are studied together, and thus, any confounding factor derived from strain-related heterogeneity is excluded from the analysis [23, 48].
Firstly, we demonstrated that L. infantum is able to become resistant to high concentrations of TPT. While the mechanisms involved in TPT resistance have not been fully elucidated in tumor cells, several studies have shown the implication of different drug transporters, such as multidrug resistance-associated protein 1 (ABCC1/MRP1) or the ABCG2 transporter [49, 50]. Likewise, Leishmania parasites rely on the amplification of ABC transporter MRPA (ABCC) and inactivation of the aquaglyceroporin 1 (AQP1) gene to counter the action of antimonial drugs [21, 23, 51]. As gene expression in Leishmania is regulated predominantly by gene dosage [20, 21], we proceeded to search for large-scale copy number variations (deletions and duplications) in chromosomes of the three selected clones for comparison with the unselected parental line. In the past, different ABCG and ABCC efflux-pump gene clusters were identified in L. infantum as part of chromosomes 6, 23, and 31 . In addition, overexpression (with no changes in gene dosage) of the ABCG6 transporter is known to be involved in CPT resistance in Leishmania parasites . Of note, none of these regions was found amplified in any of the three TPT-resistant clones in our whole-genome comparative analysis.
The absence of significant amplifications, coupled with the very unusual fact that no prominent changes in ploidy were observed for any of the studied clones , led us to suspect the potential implication of SNPs and small nucleotide insertions or deletions (indels) in the TPT-resistant phenotypes. Several SNPs have been shown to contribute to drug resistance (e.g. miltefosine, antimonial drugs, etc.) in Leishmania parasites by altering the activity of specific transporters or modifying different detoxification pathways [22, 23, 29, 55, 56]. Here we focused on SNPs and indels present in the three TPT-resistant clones and, at the same time, occurring in the same ORF but at different positions. Among the eight genes fulfilling these criteria, we identified the ORF coding for the large subunit of the DNA topoisomerase IB, which is the main target of CPT derivatives once bound to the DNA during cell replication [10, 11]. Clones TPT700.1 and TPT700.3 displayed homozygous mutations in the top1B gene, while the SNP identified in TPT700.2 was heterozygous. Although previously observed in diploid Leishmania parasites, homozygous mutations are rare because of their “nonreversible” nature. The homozygous mutations in clones TPT700.1 and TPT700.3 may have originated from loss of heterozygosity, a well-described phenomenon in Leishmania [57, 58]. Importantly, these results reinforce previous works demonstrating the possibility that, although rare, Leishmania can generate SNPs associated with drug resistance without the need for alteration of its genomic architecture and gene expression .
Due to the impossibility of generating a null mutant, our preferred method for studying the role of mutated variants of top1B consisted of episomal transfection of the mutated forms into a WT strain . Since these transfected parasites still carry the top1BWT allele, the mutated, overexpressed forms of the protein are in competition with the WT large subunit in their binding to the small subunit (in order to make up the functional heterodimer). As such, we were able to only partly recreate the highly resistant phenotype observed in the TPT700 mutants. However, the relative strength of each mutation followed the same drug resistance trend in both the original mutants (TPT700.1 > TPT700.3 > TPT700.2) and the episomal transfectants (top1BF187Y > top1BW232R > top1BG191A).
To better understand the potential contribution of these three SNPs to TPT resistance in the mutants, we performed several MD simulations. All three residues identified in the TPT700 mutants (F187, G191, and W232) were conserved between the human and parasitic enzyme and can be structurally aligned (Fig. 4). They were located in proximity of the TPT binding site, in close proximity to residues found to be crucial in the human enzyme for the interaction and stabilization of TPT with residues R364 and N722 once intercalated between DNA bases , corresponding to leishmanial residues R190 and N221. Thus, we can hypothesize that a change in one of the residues of this cluster may influence the arrangement of the TPT binding site. Although these residues are in proximity of the catalytic pentad [11, 14, 61], the mutations identified in this study are likely affecting the ability of binding of TPT to the DNA-TOP1B complex without altering the global catalytic function of the enzyme. Since TPT hinders DNA rotation within the covalent complex, a reduced binding of this drug would result in a faster religation step. Indeed, it has been very well established that a malfunction of cleavage/religation reactions will be reflected through altered protein drug sensitivity [11, 60]. Of note, two of the three SNPs identified in the TPT700 mutants (F187Y and G191A) were located within the conserved region corresponding to amino acids 361–365 in the hTOP1B enzyme. These results confirm the findings of Rubin et al. showing that a substitution of residue F361 can induce high levels of resistance against a CPT derivate (e.g. 9-nitro-20(S)camptothecin) in human U-937 myeloid leukemia cells . Likewise, Li et al. showed that certain substitutions in the 361–364 region affect DNA cleavage/ligation by the enzyme, as well as contribute to resistance against CPT since they may be included in the CPT-binding domain . These results suggest that these mutations are able to modify the architecture of the binding site, decreasing the persistence of TPT in the binding pocket, as well indicate that CPT and TPT may share binding sites in the LiTOP1B–DNA complex.
Furthermore, in the covalent complexes, K352 and R410 demonstrate a changed profile of interaction in all three TPT700 mutants. In particular, the hydrogen bond between K352 and D353 is lost. This interaction is crucial for the correct position of K352 (known to be a key player in the religation reaction), and when incorrectly positioned affects the religation rate and thus TPT sensitivity [15, 64]. Moreover, D353 is itself involved in the network of residues and TPT interaction. As such, the lack of the K352-D353 hydrogen bond and side chains orientation may be a main cause for rearrangement of the TPT binding site and lowered stabilization of the drug in the binding pocket, further explaining the observed resistance.
Importantly, the three TPT700 mutants became resistant to TPT without impairing their ability to proliferate in vitro. The process of becoming drug-resistant can lead to different evolutionary disadvantages (“fitness cost”), such as reduced survival . However, this concept remains controversial in Leishmania and is highly dependent on the parasite’s genetic and environmental context . Likewise, drug-resistant cancer cell lines exhibit different fitness–cost profiles, including subpopulations with increased fitness when compared to their sensitive counterparts . The absence of fitness cost in vitro in the TPT700 mutants could be due in part to the fact that these cells do not use ATP-dependent drug efflux pumps to resist treatment (e.g. MRPA in antimony resistance), for which they would have to divert energy away from proliferation towards running of the pumps. Moreover, the absence of a cost in terms of growth would also explain why the TPT700 mutants did not return to sensitivity once the TPT was withdrawn. Mutants of this type have the potential to become a major risk for the spread of drug resistance into an environment devoid of antileishmanial drugs. However, at this point, whether this phenomenon would be stable in vivo or in the vector remains to be evaluated.
This study represents the first whole-genome characterization of Leishmania parasites repeatedly exposed to a TOP1B inhibitor. Unlike current antileishmanial agents, TPT resistance did not have a major impact on leishmanial genomic organization or TOP1B expression levels, and did not lead to changes in gene dosage of known genes coding for efflux pumps; a phenomenon previously reported for CPT analogues in cancer cells. Of note, we found that these mutations could decrease the binding of TPT to the DNA-TOP1B binary complex, as well as lower the stabilization of the drug in the binding pocket of the leishmanial heterodimeric enzyme. Altogether, this work helps elucidate the previously unclear potential mechanisms of TPT resistance in Leishmania by mutations in the large subunit of TOP1B, and provides a valuable clue for the design of improved inhibitors to combat drug resistance. Due to the conserved nature of the mutated amino acids, this knowledge could also provide important means of overcoming resistance to TOP1B-specific drugs and developing diagnostic tools to detect TPT-resistant tumor cells. Finally, this study supports and expands the importance of including drug resistance assessments in drug discovery and drug repurposing cascades before proposing a molecule as a potential prototype for treatment of parasitic diseases.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its additional files. The sequence data for L. infantum TPT-resistant mutants is available at the NCBI BioProject (https://www.ncbi.nlm.nih.gov/bioproject/) under study accession PRJNA647847 and sample accessions SAMN15599759, SAMN15599760, and SAMN15599761, corresponding to TPT-resistant clones TPT700.1, TPT700.2 and TPT700.3, respectively.
- A600 :
Absorbance at 600 nm
ATP binding cassette
DNA copy number variations
- EC50 :
Half maximal effective concentration
US Food and Drug Administration
- hyg :
Hygromycin phosphotransferase gene
Homo sapiens DNA topoisomerase 1B
Leishmania donovani infantum DNA topoisomerase 1B
Multidrug resistance-associated protein
Neglected tropical disease
Open reading frame
Type 1B DNA topoisomerase
Fernández-Prada C, Douanne N, Minguez-Menendez A, Pena J, Tunes LG, Pires DEV et al. (2019) Repurposed Molecules A new hope in tackling neglected infectious diseases. In: Roy K, editor. In silico drug design. Academic Press; pp. 119–60.
Aronson N, Herwaldt BL, Libman M, Pearson R, Lopez-Velez R, Weina P, et al. Diagnosis and treatment of leishmaniasis: clinical practice guidelines by the Infectious Diseases Society of America (IDSA) and the American Society of Tropical Medicine and Hygiene (ASTMH). Am J Trop Med Hyg. 2017;96:24–45.
Bennis I, De Brouwere V, Belrhiti Z, Sahibi H, Boelaert M. Psychosocial burden of localised cutaneous leishmaniasis: a scoping review. BMC Public Health. 2018;18:358.
Heydarpour F, Sari AA, Mohebali M, Shirzadi M, Bokaie S. Incidence and disability-adjusted life years (DALYS) attributable to leishmaniasis in Iran, 2013. Ethiop J Health Sci. 2016;26:381–8.
Dantas-Torres F. Canine leishmaniosis in South America. Parasit Vectors. 2009;2(Suppl. 1):S1.
Dostalova A, Volf P. Leishmania development in sand flies: parasite–vector interactions overview. Parasit Vectors. 2012;5:276.
Souza Silva JA, Tunes LG, Coimbra RS, Ascher DB, Pires DEV, Monte-Neto RL. Unveiling six potent and highly selective antileishmanial agents via the open source compound collection ‘Pathogen Box’ against antimony-sensitive and -resistant Leishmania braziliensis. Biomed Pharmacother. 2021;133:111049.
Prada CF, Álvarez-Velilla R, Balaña-Fouce R, Prieto C, Calvo-Álvarez E, Escudero-Martínez JM, et al. Gimatecan and other camptothecin derivatives poison Leishmania DNA-topoisomerase IB leading to a strong leishmanicidal effect. Biochem Pharmacol. 2013;85:1433–40.
Wang JC. Interaction between DNA and an Escherichia coli protein ω. J Mol Biol. 1971;55:523-IN16.
Pommier Y, Sun Y, Huang SN, Nitiss JL. Roles of eukaryotic topoisomerases in transcription, replication and genomic stability. Nat Rev Mol Cell Biol. 2016;17:703–21.
Pommier Y. Topoisomerase I inhibitors: camptothecins and beyond. Nat Rev Cancer. 2006;6:789–802.
Ubhi T, Brown GW. Exploiting DNA replication stress for cancer treatment. Cancer Res. 2019;79:1730–9.
Balaña-Fouce R, Álvarez-Velilla R, Fernández-Prada C, García-Estrada C, Reguera RM. Trypanosomatids topoisomerase re-visited. New structural findings and role in drug discovery. Int J Parasitol Drugs Drug Resist. 2014;4:326–37.
Villa H, Otero Marcos AR, Reguera RM, Balana-Fouce R, Garcia-Estrada C, Perez-Pertejo Y, et al. A novel active DNA topoisomerase I in Leishmania donovani. J Biol Chem. 2003;278:3521–6.
Prada CF, Álvarez-Velilla R, Diaz-González R, Prieto C, Pérez-Pertejo Y, Balaña-Fouce R, et al. A pentapeptide signature motif plays a pivotal role in Leishmania DNA topoisomerase IB activity and camptothecin sensitivity. Biochim Biophys Acta. 2012;1820:2062–71.
Prada CF, Alvarez-Velilla R, Diaz-Gozalez R, Perez-Pertejo Y, Balana-Fouce R, Reguera RM. Identification and characterization of the regions involved in the nuclear translocation of the heterodimeric leishmanial DNA topoisomerase IB. PLoS ONE. 2013;8:e73565.
Balana-Fouce R, Prada CF, Requena JM, Cushman M, Pommier Y, Alvarez-Velilla R, et al. Indotecan (LMP400) and AM13-55: two novel indenoisoquinolines show potential for treating visceral leishmaniasis. Antimicrob Agents Chemother. 2012;56:5264–70.
D’Annessa I, Castelli S, Desideri A. Topoisomerase 1B as a target against leishmaniasis. Mini Rev Med Chem. 2015;15:203–10.
Kundig C, Leblanc E, Papadopoulou B, Ouellette M. Role of the locus and of the resistance gene on gene amplification frequency in methotrexate resistant Leishmania tarentolae. Nucleic Acids Res. 1999;27:3653–9.
Gazanion E, Fernandez-Prada C, Papadopoulou B, Leprohon P, Ouellette M. Cos-Seq for high-throughput identification of drug target and resistance mechanisms in the protozoan parasite Leishmania. Proc Natl Acad Sci USA. 2016;113(21):E3012–21.
El Fadili K, Messier N, Leprohon P, Roy G, Guimond C, Trudel N, et al. Role of the ABC transporter MRPA (PGPA) in antimony resistance in Leishmania infantum axenic and intracellular amastigotes. Antimicrob Agents Chemother. 2005;49:1988–93.
Perez-Victoria FJ, Gamarro F, Ouellette M, Castanys S. Functional cloning of the miltefosine transporter. A novel P-type phospholipid translocase from Leishmania involved in drug resistance. J Biol Chem. 2003;278:49965–71.
Douanne N, Wagner V, Roy G, Leprohon P, Ouellette M, Fernandez-Prada C. MRPA-independent mechanisms of antimony resistance in Leishmania infantum. Int J Parasitol Drugs Drug Resist. 2020;13:28–37.
Fernandez-Prada C, Vincent IM, Brotherton MC, Roberts M, Roy G, Rivas L, et al. Different mutations in a p-type ATPase Transporter in Leishmania parasites are associated with cross-resistance to two leading drugs by distinct mechanisms. PLoS Negl Trop Dis. 2016;10:e0005171.
Hefnawy A, Berg M, Dujardin JC, De Muylder G. Exploiting knowledge on Leishmania drug resistance to support the quest for new drugs. Trends Parasitol. 2017;33:162–74.
Van den Kerkhof M, Mabille D, Hendrickx S, Leprohon P, Mowbray CE, Braillard S, et al. Antileishmanial aminopyrazoles: studies into mechanisms and stability of experimental drug resistance. Antimicrob Agents Chemother. 2020;64:e00152-e220.
Bhattacharya A, Corbeil A, do Monte-Neto RL, Fernandez-Prada C. Of drugs and trypanosomatids: new tools and knowledge to reduce bottlenecks in drug discovery. Genes (Basel). 2020;11:722.
Aslett M, Aurrecoechea C, Berriman M, Brestelli J, Brunk BP, Carrington M, et al. TriTrypDB: a functional genomic resource for the Trypanosomatidae. Nucleic Acids Res. 2010;38(Database issue):D457–62.
Ritt JF, Raymond F, Leprohon P, Legare D, Corbeil J, Ouellette M. Gene amplification and point mutations in pyrimidine metabolic genes in 5-fluorouracil resistant Leishmania infantum. PLoS Negl Trop Dis. 2013;7:e2564.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
Chiang DY, Getz G, Jaffe DB, O’Kelly MJ, Zhao X, Carter SL, et al. High-resolution mapping of copy-number alterations with massively parallel sequencing. Nat Methods. 2009;6:99–103.
Papadopoulou B, Roy G, Ouellette M. A novel antifolate resistance gene on the amplified H circle of Leishmania. EMBO J. 1992;11:3601–8.
Fernandez-Prada C, Sharma M, Plourde M, Bresson E, Roy G, Leprohon P, et al. High-throughput Cos-Seq screen with intracellular Leishmania infantum for the discovery of novel drug-resistance mechanisms. Int J Parasitol Drugs Drug Resist Int J Parasitol Drugs Drug Resist. 2018;8:165–73.
Davies DR, Mushtaq A, Interthal H, Champoux JJ, Hol WG. The structure of the transition state of the heterodimeric topoisomerase I of Leishmania donovani as a vanadate complex with nicked DNA. J Mol Biol. 2006;357:1202–10.
Roy A, Chowdhury S, Sengupta S, Mandal M, Jaisankar P, D’Annessa I, et al. Development of derivatives of 3, 3′-diindolylmethane as potent Leishmania donovani bi-subunit topoisomerase IB poisons. PLoS ONE. 2011;6:e28493.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera—a visualization system for exploratory research and analysis. J Comput Chem. 2004;25:1605–12.
Staker BL, Hjerrild K, Feese MD, Behnke CA, Burgin AB Jr, Stewart L. The mechanism of topoisomerase I poisoning by a camptothecin analog. Proc Natl Acad Sci USA. 2002;99:15387–92.
Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983;79:926.
Maier JA, Martinez C, Kasavajhala K, Wickstrom L, Hauser KE, Simmerling C. ff14SB: improving the accuracy of protein side chain and backbone parameters from ff99SB. J Chem Theory Comput. 2015;11:3696–713.
Darden T, York D, Pedersen L. Particle mesh Ewald: an N⋅log(N) method for Ewald sums in large systems. J Chem Phys. 1993;98:10089.
Ceriotti M, Bussi G, Parrinello M. Langevin equation with colored noise for constant-temperature molecular dynamics simulations. Phys Rev Lett. 2009;102(2):020601.
Feller SE, Zhang Y, Pastor RW. Constant pressure molecular dynamics simulation: the Langevin piston method. J Chem Phys. 1995;103:4613.
Ryckaert J-P, Ciccotti G, Berendsen HJC. Numerical integration of the Cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys. 1977;23:327–41.
Miyamoto S, Kollman PA. Settle: an analytical version of the SHAKE and RATTLE algorithm for rigid water models. J Comput Chem. 1992;8:952–62.
Hess B, Kutzner C, van der Spoel D, Lindahl E. GROMACS 4: algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theory Comput. 2008;4:435–47.
Garcia-Estrada C, Prada CF, Fernandez-Rubio C, Rojo-Vazquez F, Balana-Fouce R. DNA topoisomerases in apicomplexan parasites: promising targets for drug discovery. Proc Biol Sci. 2010;277:1777–87.
Chong CR, Sullivan DJ Jr. New uses for old drugs. Nature. 2007;448:645–6.
Hendrickx S, Guerin PJ, Caljon G, Croft SL, Maes L. Evaluating drug resistance in visceral leishmaniasis: the challenges. Parasitology. 2018;145:453–63.
Tian Q, Zhang J, Chan SY, Tan TM, Duan W, Huang M, et al. Topotecan is a substrate for multidrug resistance associated protein 4. Curr Drug Metab. 2006;7:105–18.
Allen JD, Brinkhuis RF, Wijnholds J, Schinkel AH. The mouse Bcrp1/Mxr/Abcp gene: amplification and overexpression in cell lines selected for resistance to topotecan, mitoxantrone, or doxorubicin. Cancer Res. 1999;59:4237–41.
Legare D, Cayer S, Singh AK, Richard D, Papadopoulou B, Ouellette M. ABC proteins of Leishmania. J Bioenerg Biomembr. 2001;33:469–74.
Leprohon P, Legare D, Girard I, Papadopoulou B, Ouellette M. Modulation of Leishmania ABC protein gene expression through life stages and among drug-resistant parasites. Eukaryot Cell. 2006;5:1713–25.
BoseDasgupta S, Ganguly A, Roy A, Mukherjee T, Majumder HK. A novel ATP-binding cassette transporter, ABCG6 is involved in chemoresistance of Leishmania. Mol Biochem Parasitol. 2008;158(2):176–88.
Laffitte MN, Leprohon P, Papadopoulou B, Ouellette M. Plasticity of the Leishmania genome leading to gene copy number variations and drug resistance. F1000 Res. 2016;5:2350.
Vasudevan G, Ullman B, Landfear SM. Point mutations in a nucleoside transporter gene from Leishmania donovani confer drug resistance and alter substrate selectivity. Proc Natl Acad Sci USA. 2001;98:6092–7.
Monte-Neto R, Laffitte MC, Leprohon P, Reis P, Frezard F, Ouellette M. Intrachromosomal amplification, locus deletion and point mutation in the aquaglyceroporin AQP1 gene in antimony resistant Leishmania (Viannia) guyanensis. PLoS Negl Trop Dis. 2015;9:e0003476.
Gueiros-Filho FJ, Beverley SM. Selection against the dihydrofolate reductase-thymidylate synthase (DHFR-TS) locus as a probe of genetic alterations in Leishmania major. Mol Cell Biol. 1996;16(10):5655–63.
Bhattacharya A, Leprohon P, Bigot S, Padmanabhan PK, Mukherjee A, Roy G, et al. Coupling chemical mutagenesis to next generation sequencing for the identification of drug resistance mutations in Leishmania. Nat Commun. 2019;10:5627.
Patino LH, Imamura H, Cruz-Saavedra L, Pavia P, Muskus C, Mendez C, et al. Major changes in chromosomal somy, gene expression and gene dosage driven by Sb(III) in Leishmania braziliensis and Leishmania panamensis. Sci Rep. 2019;9:9485.
Marchand C, Antony S, Kohn KW, Cushman M, Ioanoviciu A, Staker BL, et al. A novel norindenoisoquinoline structure reveals a common interfacial inhibitor paradigm for ternary trapping of the topoisomerase I-DNA covalent complex. Mol Cancer Ther. 2006;5:287–95.
Diaz-Gonzalez R, Perez-Pertejo Y, Pommier Y, Balana-Fouce R, Reguera RM. Mutational study of the “catalytic tetrad” of DNA topoisomerase IB from the hemoflagellate Leishmania donovani: role of Asp-353 and Asn-221 in camptothecin resistance. Biochem Pharmacol. 2008;76:608–19.
Rubin E, Pantazis P, Bharti A, Toppmeyer D, Giovanella B, Kufe D. Identification of a mutant human topoisomerase I with intact catalytic activity and resistance to 9-nitro-camptothecin. J Biol Chem. 1994;269:2433–9.
Li XG, Haluska P Jr, Hsiang YH, Bharti AK, Kufe DW, Liu LF, et al. Involvement of amino acids 361 to 364 of human topoisomerase I in camptothecin resistance and enzyme catalysis. Biochem Pharmacol. 1997;53:1019–27.
Seol Y, Zhang H, Pommier Y, Neuman KC. A kinetic clutch governs religation by type IB topoisomerases and determines camptothecin sensitivity. Proc Natl Acad Sci USA. 2012;109:16125–30.
Orr HA. Fitness and its role in evolutionary genetics. Nat Rev Genet. 2009;10:531–9.
Jensen NF, Stenvang J, Beck MK, Hanakova B, Belling KC, Do KN, et al. Establishment and characterization of models of chemotherapy resistance in colorectal cancer: towards a predictive signature of chemoresistance. Mol Oncol. 2015;9:1169–85.
We would like to thank to Professor Marc Ouellette for his outstanding support with WGS experiments. We would like also to thank Dr. Aida Minguez-Menendez for her help with the conception and creation of the graphical abstract and figures.
This research was funded by the Fonds de Recherche du Québec-Nature et Technologies (FRQNT) through its New University Researchers Start-up Program, Grant Number 2019-NC-253631; the Canada foundation for Innovation (www.innovation.ca), Grant Number 37324; and by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant RGPIN-2017-04480 awarded to CFP. CRT, VW, and AC were supported by the FRQNT and the NSERC studentship programs. RMN is a CNPq (Brazilian National Council for Scientific and Technological Development) Research Fellow (#310640/2017-2).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no financial or non-financial competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Log2-transformed 16 × TPT-resistant mutants/WT reads ratios for non-overlapping 5 kb genomic windows for the 36 L. infantum chromosomes in the three TPT-resistant strains.
Homozygous and heterozygous SNPs in the different TPT-resistant mutants.
LinJ.34.3220 (top1B large subunit) and LinJ.04.0070 (top1B small subunit) RNA expression in Leishmania cells.
About this article
Cite this article
Rosa-Teijeiro, C., Wagner, V., Corbeil, A. et al. Three different mutations in the DNA topoisomerase 1B in Leishmania infantum contribute to resistance to antitumor drug topotecan. Parasites Vectors 14, 438 (2021). https://doi.org/10.1186/s13071-021-04947-4
- Protozoan parasites
- DNA topoisomerases
- Drug resistance