Next Article in Journal
MicroRNA as a Novel Biomarker in the Diagnosis of Head and Neck Cancer
Previous Article in Journal
Relationship between Serum Alkaline Phosphatase and Low Muscle Mass Index Among Korean Adults: A Nationwide Population-Based Study
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 11
Issue 6
10.3390/biom11060843
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nucleoside Analogues Are Potent Inducers of Pol V-mediated Mutagenesis

by
Balagra Kasim Sumabe
1,2,3,
Synnøve Brandt Ræder
1,
Lisa Marie Røst
4,
Animesh Sharma
5,
Eric S. Donkor
6,
Lydia Mosi
2,3,
Samuel Duodu
2,3,
Per Bruheim
4 and
Marit Otterlei
1,7,*
1
Department of Clinical and Molecular Medicine, Faculty of Medicine and Health Sciences, NTNU, Norwegian University of Science and Technology, NO-7489 Trondheim, Norway
2
West African Centre for Cell Biology of Infectious Pathogens (WACCBIP), University of Ghana, P.O. BOX LG 54 Accra, Ghana
3
Department of Biochemistry, Cell and Molecular Biology, University of Ghana, P.O. BOX LG 54 Accra, Ghana
4
Department of Biotechnology and Food Science, Faculty of Natural Sciences, NTNU Norwegian University of Science and Technology, NO-7481 Trondheim, Norway
5
Proteomics and Modomics Experimental Core Facility (PROMEC), NTNU Norwegian University of Science and Technology, NO-7481 Trondheim, Norway
6
Department of Medical Microbiology, University of Ghana Medical School, P.O. Box 4236 Accra, Ghana
7
Clinic of Laboratory medicine, St. Olav University Hospital, NO-7006 Trondheim, Norway
*
Author to whom correspondence should be addressed.
Biomolecules 2021, 11(6), 843; https://doi.org/10.3390/biom11060843
Submission received: 21 April 2021 / Revised: 1 June 2021 / Accepted: 2 June 2021 / Published: 5 June 2021
(This article belongs to the Section Molecular Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Drugs targeting DNA and RNA in mammalian cells or viruses can also affect bacteria present in the host and thereby induce the bacterial SOS system. This has the potential to increase mutagenesis and the development of antimicrobial resistance (AMR). Here, we have examined nucleoside analogues (NAs) commonly used in anti-viral and anti-cancer therapies for potential effects on mutagenesis in Escherichia coli, using the rifampicin mutagenicity assay. To further explore the mode of action of the NAs, we applied E. coli deletion mutants, a peptide inhibiting Pol V (APIM-peptide) and metabolome and proteome analyses. Five out of the thirteen NAs examined, including three nucleoside reverse transcriptase inhibitors (NRTIs) and two anti-cancer drugs, increased the mutation frequency in E. coli by more than 25-fold at doses that were within reported plasma concentration range (Pl.CR), but that did not affect bacterial growth. We show that the SOS response is induced and that the increase in mutation frequency is mediated by the TLS polymerase Pol V. Quantitative mass spectrometry-based metabolite profiling did not reveal large changes in nucleoside phosphate or other central carbon metabolite pools, which suggests that the SOS induction is an effect of increased replicative stress. Our results suggest that NAs/NRTIs can contribute to the development of AMR and that drugs inhibiting Pol V can reverse this mutagenesis.
Keywords:
NA; NRTIs; β-clamp; SOS; Pol V; TLS; AMR; MDR

1. Introduction

Antimicrobial resistance (AMR) is on the rise, posing serious global health threats, including high medical costs and mortality rates. In the United States, AMR is reported to account for nearly 3 million infections and 35 thousand deaths yearly [1]. In Europe, a higher prevalence of AMR infections is reported by countries in the southern and eastern parts of the continent, compared to countries in the north [2]. Exact data are not available for Africa because 40% of the countries do not have records on the incidence of AMR, but the fraction of infections with AMR is increasing [3]. Several Asian countries have a very high incidence of AMR [3] and, in the next 30 years, antibiotic resistance is estimated to cause an annual global mortality rate of 10 million deaths and a USD 100 trillion loss in the global economy, if not addressed. Most of these deaths and losses in the economy are predicted to be from Africa and Asia [4]. Following the outbreak of COVID-19 [5], the increased usage of antibiotics could further accelerate the development of antibiotic resistance [5,6,7].
Bacteria gain resistance to antibiotics via the development of endogenous mutations and/or transfer of resistance genes. Cellular stress, induced, for example, by antibiotics, activates SOS-dependent DNA translesion synthesis (TLS) and this is a major mechanism through which endogenous mutations lead to antibiotic resistance [8]. The ability of the TLS polymerases to introduce mutations is linked to their more spacious catalytic site, enabling a bypass of DNA lesions, and their lack of 3′-5’ exonuclease activity (proofreading activity). There are three TLS polymerases in Escherichia coli, Pol II, Pol IV and Pol V (gene product of polB, dinB and umuDC, respectively). In contrast to Pol II and Pol IV, which are also expressed at low levels in absence of SOS response, Pol V is strictly SOS-regulated and, therefore, vital for bacterial mutagenesis [9]. All three bacterial TLS polymerases have been implicated in the generation of mutations contributing to the development of resistance and virulence in bacteria. For example, all three polymerases contribute to the development of mutations conferring ciprofloxacin and rifampicin resistance [8] and both Pol IV and Pol V generate base substitution mutations in the ampD gene, conferring resistance to ampicillin [10]. Pol IV is also shown to be relevant for virulence in uropathogenic E. coli (UPEC) strains in infected mice [11]. All bacterial polymerases interact with a subunit of the DNA polymerase III (Pol III) holoenzyme, the β-clamp, a homodimer with structural similarities to the mammalian DNA sliding clamp, PCNA (proliferating cell nuclear antigen) [12].
Infections with resistant bacteria will certainly be more devastating in immunocompromised individuals, such as those having HIV/AIDS and cancer, compared to the normal population. Thus, with limited help from the immune system, efficient antibiotics become vital. It is therefore worrying that an increased number of multidrug-resistant (MDR) bacteria have been isolated from HIV-positive individuals compared to HIV-negative individuals in several studies [13,14,15,16]. Characterization of MDR bacteria isolated from HIV patients suggests that endogenous mutations contribute more to the development of resistance than uptake of plasmids carrying drug resistance [17].
Nucleoside analogues (NAs) are used not only in anti-cancer therapies, but also in anti-retroviral therapies as they effectively reduce viral load, and thereby improve the health of the patients and reduce the risk of viral transmission [18]. NAs mimic the natural building blocks of DNA, deoxyribonucleotide triphosphates (dNTPs), and are incorporated into DNA. However, elongations and/or base pairing are blocked by these drugs and this leads to termination of replication [19]. Some NAs, e.g., zidovudine, have a much higher affinity for viral reverse transcriptase than for human polymerases and these are therefore defined as nucleoside reverse transcriptase inhibitors (NRTI) [20].
Here, we have investigated the mutagenic activities of some NAs commonly used in anti-viral and anti-cancer therapies on E. coli MG1655. We show that NAs, particularly stavudine, didanosine, 5-fluoro-uracil (5-FU), emtricitabine and fludarabine, increase the mutation frequency by 50–100-fold in E. coli. Deletion mutant analysis showed that the increased mutagenesis was driven by the TLS polymerase Pol V, suggesting that the SOS response is induced in the bacteria upon NA treatment. Furthermore, the cell-penetrating APIM-peptide, previously shown to target the polymerase interaction site on the β-clamp and inhibit Pol V, significantly inhibited the NA-induced increase in mutation frequency. This supports a Pol V-dependent mutagenesis. Induction of SOS by stavudine was verified by proteomic analysis and metabolite profiling and showed that SOS was not induced by nucleotide starvation. Data presented here suggest that stavudine induced SOS by interfering with the bacterial replications machinery at doses far below the minimal inhibitory concentration (MIC).

2. Materials and Methods

2.1. Bacterial Species, Strains and Growth Measurements

E. coli K12 (MG1655) wild type (WT) and TLS polymerase deletion strains ΔumuDC (ΔPol V) from the Keio collection and ΔpolB/ΔdinB (ΔPol II/ΔPol IV) made in-house as previously described [21] were used. Briefly, polB was inactivated in E. coli MG1655 ΔPol IV using JW0059: polB. A phage stock (T4GT7) diluted in T4 buffer was added to the JW0059::polB culture, incubated and plated on LB agar plates. The lysate was harvested and added to 2 mL T4 buffer and 0.2 mL CHCl3 (T4GT7 lysate). The mixture was centrifuged and supernatant stored at 4 °C, over 20 μL CHCl3. An amount of 0.5 mL overnight culture of the recipient cells was re-suspended in 1 mL of T4 buffer and aliquoted. The T4GT7 lysate (10−1, 10−2, 10−3 and 10−4 dilutions) was added to the aliquots, incubated for 25 minutes at room temperature and plated on LB-agar plates with kanamycin. Colonies obtained were re-streaked twice to obtain a pure culture. Kanamycin resistance cassettes replacing the deleted polB genes were removed by an FLP-recombinase expressed by pCP20 plasmid.

2.2. Peptide

A cell-penetrating peptide targeting the β-clamp named APIM-peptide (Ac-MDRWLVK-GILQWRKI-RRRRRRRR-NH2) (Innovagen, Lund, Sweden) was used. This is the same peptide as RWLVK*, used in Nedal et al. [21], except it has three arginine residues less in its cell-penetrating tail.

2.3. Nucleoside Analogues (NAs)

NAs used in this study (listed in Table 1) were purchased from Glentham Life Sciences Ltd. (Wiltshire, UK).

2.4. Bacterial Cultivation

Luria Bertani (LB) broth Miller/Microbiological agar and Mueller Hinton broth/agar were used for bacterial cultivation in the rifampicin mutagenicity (RifR) and MIC assays, respectively. Bacterial cultures were incubated at 37 °C, with or without shaking at 250 rpm. Rifampicin (Tokyo chemicals industry Co Ltd., Tokyo, Japan) was used for the RifR assay.

2.5. MIC Assay

The MIC of the nucleoside analogues and APIM-peptide were determined following the protocol used for microtiter broth microdilution assays established by the clinical and laboratory standard institutes [22] with modifications as described [21]. MIC was determined after 24 h.

2.6. Rifampicin Mutagenicity Assay (RifR) Assay

This assay was conducted as described [21]. Briefly, pre-cultures of E. coli MG1655 (WT and mutant strains) were inoculated in LB broth (Miller) and grown at 37 °C, for 16 h. The pre-cultures were diluted in fresh LB broth (1:100) and 10 mL cultures were grown in 50 mL conical tubes for one hour at 37 °C (shaking incubator at 250 rpm), before they were treated with NAs and/or APIM-peptide. Untreated cultures were used as control. The cultures were further incubated for two hours at 37 °C, before the cultures were serially diluted in LB and plated on LB agar. Undiluted cultures (800 μL) were plated in LB soft agar containing rifampicin (100 µg/mL) on LB agar plates and incubated for 48 h, at 37 °C. The bacterial colony-forming units per millilitre (CFU/mL) were determined from plates with and without rifampicin and the frequency of rifampicin resistance per 108 CFU (RifR/108) was calculated.

2.7. Protein Extraction

When harvesting cells for the RifR assay, 2 × 2 mL of untreated cultures and culture treated with stavudine (Stav) (0.32 µg/mL), APIM-peptide (26 µg/mL), or the combination were pelleted for cell extracts (4000 rpm, 10 min) for each biological replica (n = 3). The pellets were resuspended in an OmniCleave™ storage buffer (Lucigen, Middleton, WI, USA) and incubated at room temperature in the presence of 200 U OmniCleave™ endonuclease (Lucigen), 20 µg/µL RNAse (Sigma-Aldrich, St. Louis, MO, USA), 10 U DNAse I (Sigma-Aldrich), 1× phosphatase inhibitor cocktails II and III (Sigma-Aldrich, St. Louis, MO, USA) and 1× cOmplete™ protease inhibitor cocktail (Roche Diagnostics GmbH, Basel, Switzerland). The cells were sonicated on ice at 1.5 output control and 15% duty cycle, for 30 s. The sonication was repeated until a clear phase was observed and cell debris was removed by centrifugation (1400 rpm, 15 min). The supernatants were snap-frozen and stored at −80 °C. The total protein concentration was measured using NanoDrop One C (Thermo Scientific, Waltham, MA, USA) and the Bio-Rad protein assay with a UV-1700 visible spectrophotometer (Pharmaspec, Redmond, WA, USA, 595 nm).

2.8. Multiplexed Inhibitor Beads (MIB) Assay and Mass Spectrometry (MS) Analysis of Proteome

Before MS analysis of the protein extract, kinases and other ATP/GTP binding proteins were enriched using the MIB assay as described [23]. A total of 100 µg of the protein extract was analysed per sample.

2.9. Metabolite Profiling

When harvesting cells for the RifR assay, 4 × 4 mL of untreated culture and culture treated with stavudine (Stav) (0.32 µg/mL), APIM-peptide (26 µg/mL), or the combination were also harvested for MS -based metabolite profiling of central carbon metabolite pools, including glycolytic, pentose phosphate pathway (PPP) and tricarboxylic acid cycle (TCA) intermediates, and the complete nucleoside phosphate pool. Complete sampling, processing and MS-based metabolite profiling protocols can be found in previous publications [24,25]. In short, cells were harvested by fast filtration using a vacuum filtering unit with controlled pressure 800 psi below the ambient pressure. After a quick water rinse, filters were transferred to tubes containing 50% ice-cold acetonitrile (VWR chemicals, city, France) for quenching of metabolism. The tubes were immediately snap-frozen and stored at −80 °C. Metabolite extraction was performed through three cycles of freeze-thaw, followed by centrifugation to remove cell debris before the freeze-drying of the cell-free metabolite containing supernatant. Dried supernatants were reconstituted in distilled water, filtered through 3 kDa spin filters (VWR) and analysed on a capillary ion chromatograph (capIC, Thermo Scientific) coupled to a TQ-XS triple quadrupole MS (Waters Corporation, Milford, PA, USA). Samples and analytical grade standards (Sigma-Aldrich, St. Louis, MO, USA) were spiked with yeast extract cultivated with U13C-glucose (Sigma-Aldrich, St. Louis, MO, USA) as the sole carbon source, as described in [25].

2.10. Data Analysis

Data from the MIB assay were analysed as described [23]. CapIC-MS/MS data were processed in the TargetLynx application manager of MassLynx 4.1 (Waters Corporation, Milford, PA, USA). Absolute quantification was performed by interpolation of calibration curves prepared from serial dilutions of corresponding analytical grade standards (Sigma-Aldrich), calculated by least squares regression with 1/x weighting. Response factors of the analytical standard and biological extracts were corrected by the corresponding response factor of the U13C-labeled isotopologue spiked into the samples, a strategy enabling the highest quantitative accuracy and precision [26,27]. Principal component analysis (PCA) was performed on auto-scaled metabolite concentrations normalized to the sum, applying MetaboAnalyst 4.0 [28]. Missing values were replaced by the mean concentration of each metabolite. Unpaired, two-tailed Student’s t-tests were used for comparing mutation frequencies; significance levels are described in figure legends.

2.11. Mass Spectrometry Data Analysis

Proteins were quantified by processing MS data using MaxQuant v.1.6.17.0 [29]. Open workflow [30] provided in FragPipe version 14 was used to inspect the raw files to determine optimal search criteria. Namely, the following search parameters were used: enzyme specified as trypsin with a maximum of two missed cleavages allowed; acetylation of protein N-terminal, oxidation of methionine, deamidation of asparagine/glutamine and phosphorylation of serine/threonine/tyrosine as dynamic post-translational modification. These were imported in MaxQuant, which uses m/z and retention time (RT) values to align each run against each other sample with 1 min window match-between-run function and 20 min overall sliding window, using a clustering-based technique. These were further queried against the E. coli (strain K12) proteome including isoforms downloaded from Uniprot (https://www.uniprot.org/proteomes/UP000000625; accessed on 24 August 2020) and MaxQuant’s internal contaminants database using Andromeda built into MaxQuant. Both protein and peptide identifications false discovery rate (FDR) was set to 1%; only unique peptides with high confidence were used for final protein group identification. Peak abundances were extracted by integrating the area under the peak curve. Each protein group abundance was normalized by the total abundance of all identified peptides for each run and protein by calculated median summing all unique and razor peptide-ion abundances for each protein, using label-free quantification (LFQ) algorithm [31] with minimum peptides ≥1. LFQ values for all samples were combined and log-transformed with base 2 and the transformed control values were subtracted. The resulting values reflecting the change relative to control for each condition were subjected to a two-sided non-parametric Wilcoxon signed rank test [32], as implemented in MATLAB R2020a (Math Works Inc., https://www.mathworks.com/, accessed on 18 March 2020), in order to check the consistency in the directionality of the change, namely, a negative sign reflecting decreased and a positive sign reflecting the increased expression of respective protein groups. The choice of this non-parametric test avoids the assumption of a certain type of null distribution as in the Student’s t-test by working over the Rank of the observation instead of the observation value itself. Further, it also makes it robust to outliers and extreme variations noticed in observed values. Differentially expressed (DE) protein groups were identified at p 0.25. The UniProt accession IDs of these DE were mapped to pathways (www.wikipathways.org/index.php/DownloadPathways version wikipathways-20201010-gmt-Homo_sapiens.gmt) using R (https://www.R-project.org/, accessed on 19 October 2020) libraries, org.Hs.eg.db and clusterProfiler (www.liebertpub.com/doi/10.1089/omi.2011.0118). Venn diagrams were built using the R package limma [33] and heatmap using pheatmap (https://cran.r-project.org/web/packages/pheatmap/index.html, accessed on 4 January 2020). The raw data with results have been deposited in the ProteomeXchange Consortium via the PRIDE [34] partner repository with the dataset identifier, project ID PXD025370 (Reviewer account Username: [email protected] Password: seTUpXp4).

3. Results and Discussion

3.1. NAs Induce Mutagenesis in Bacteria at Concentrations That Do Not Affect Bacterial Growth

Antibacterial activities of some anti-viral and anti-cancer NAs have been reported previously [19]. Since doses below the MIC of the different NAs must be used to test the mutagenic potential of these drugs, we first determined the MICs of the different NAs on E. coli MG1655. Among the NAs tested (Table 1), only stavudine, didanosine and 5-FU inhibited visible bacterial growth at concentrations below 256 µg/mL (maximal concentration examined). The lack of antibacterial activities of several of these NAs, e.g., cytarabine, vidarabine, fludarabine, abacavir, emtricitabine, lamivudine and tenofovir, are in agreement with published data [35,36].
To examine the mutagenic potential of the different NAs, E. coli MG1655 was treated with 0.5× MIC or 0.25× of the maximal concentration used in the MIC assay (256 µg/mL). The NRTIs stavudine, didanosine and emtricitabine induced a strong (45–125×) increase in mutation frequency (Figure 1A). Abacavir and lamivudine, two other NRTIs, also significantly increased the mutation frequency (~10–15×). 5-FU and fludarabine, two NAs commonly used in anticancer therapy, increased the mutation frequency by more than 50 times. The sub-MIC doses used in Figure 1A are in all cases except for 5-FU higher than plasma concentrations obtained in patients (see Table 1). We therefore examined doses that are in the human plasma concentrations range (Pl.CR) (Figure 1B) and found that in addition to 5-FU, stavudine, didanosine, emtricitabine and fludarabine were potent inducers of mutagenesis at clinically relevant doses. Thus, these drugs may increase bacterial mutagenesis and contribute to hard-to-treat bacterial infections when used clinically. It still remains elusive why these NAs are more mutagenic than the other NAs tested. Stavudine was the most potent inducer of mutagenesis. A more than 50× increase in mutation frequency was detected at 0.016 µg/mL, a dose 20-fold lower than the reported upper levels of Pl. CR (0.32 µg/mL). That NAs increase the mutation frequency in bacteria may, at least partly, explain why a high fraction of antibiotic-resistant bacteria is isolated from HIV-positive patients [13,14,15,16,35].

3.2. Pol V Is the TLS Polymerase Responsible for the Mutagenesis Induced by NAs

To determine which TLS polymerase is responsible for the induced mutagenesis, the mutation frequency in WT, a strain with deletion of Pol V (ΔPol V) and a Pol IV and Pol II double deletion strain (ΔPol II/ΔPol IV) of E. coli MG1655 were examined after treatment with stavudine. The results clearly indicate that the main TLS polymerase responsible for the increase in mutation frequency is Pol V and that the deletion of both the other two TLS polymerases in E. coli., Pol IV and Pol II, did not significantly reduce the mutation frequency (Figure 2A). Pol V is considered to be the main SOS inducible TLS polymerase in bacteria (reviewed in [51]).
We have previously shown that peptides containing the PCNA interacting motif APIM, bind to the β-clamp and block the function of Pol V [21]. Here we show that the APIM-peptide was able to reduce the mutation frequency induced by stavudine when these were combined (combination) (Figure 2B), further supporting a Pol V-dependent mutagenesis.

3.3. Stavudine Treatment Induces only Minor Changes in Central Carbon Metabolite Pools

Stavudine is a thymidine nucleoside analogue, which when incorporated into DNA by reverse transcriptase, terminates the replication of the virus. Pyrimidine synthesis is strictly regulated and some NAs, such as 5-FU, are also thymidylate synthase inhibitors. Thus, NAs may affect the nucleotide pools and this may both directly and indirectly stress the cells and increase the mutation frequency. To explore the effects of stavudine on the nucleotide pools, quantitative MS-based metabolite profiling of E. coli treated with stavudine and APIM-peptide alone and the combination of stavudine and APIM-peptide (combination) were performed. PCA did not reveal significant differences in the metabolite pool composition of the treated cells and the untreated control (95% confidence interval, Figure 3A). However, the combination treatment clustered closer to the APIM-peptide single treatment than the stavudine single treatment (Figure 3A), indicating a larger similarity between the two. This is also evident from inspecting fold changes on the level of individual metabolites; the central carbon metabolite pool in APIM-peptide and combination-treated E. coli were more similar and, more specifically, more down-regulated, compared to untreated control than stavudine single treatment (Figure 3B). Especially noteworthy is the slight increase in pyrimidine nucleoside pools in stavudine single treated cells, as opposed to down-regulated pools in both in the APIM-peptide- and the combination-treated group. Down-regulation of nucleoside pools after treatment with APIM-peptide has also been also detected in mammalian cancer cells [52]. However, the increased mutation frequency observed by stavudine as a single agent cannot be explained by cellular stress induced by the small increase in nucleotide pools detected.

3.4. Stavudine Treatment Activates the SOS System

When harvesting samples for mutation frequency and metabolite profiling, cells were also harvested for proteome analysis using the MIB assay [23]. This is an assay where activated signalling proteins, and proteins binding to these, are pulled down via their affinity for ATP/GTP using kinase inhibitors (also NAs) immobilized on sepharose beads. The proteins attached are trypsinized and identified by MS/MS. Approximately 1000 proteins were pulled down from each sample. During analysis, only proteins that were similarly changed relative to untreated control in all three biological replicas (significant change, Wilcoxon test [23]) were included in downstream analysis. More proteins were similarly changed in E. coli treated with stavudine and the combination, than in APIM-peptide-treated E. coli compared to either of the two groups (Figure 4A). This suggests that the APIM-peptide induces a different response than stavudine. When examining the change in proteins known to be related to the bacterial SOS response [53], more SOS proteins were found to change in stavudine and combination-treated cells, than in the APIM-peptide-treated cells (Figure 4B). Heatmaps of the log2 fold change in pull-down of proteins found to be significantly changed in at least one treatment group showed that most of the changed SOS regulated proteins were increased relative to control (Figure 4C). This was expected as the mutation frequencies increased. Interestingly, RecA was increased in both stavudine and the combination, but not in APIM-peptide-treated cells. The opposite was the case for LexA, which increased significantly only in APIM-peptide-treated cells. A trend towards a reduction in LexA was detected in stavudine and combination-treated bacteria, albeit not significant. These data support that stavudine is inducing an SOS response in bacteria at doses far below MIC. APIM-peptide is not inducing an SOS response, nor reducing the SOS response induced by stavudine when used in combination. This was expected as the APIM-peptide is believed to inhibit the interactions between the β-clamp and Pol V and thereby inhibit the increase in mutation frequency, but not to inhibit the SOS induction [21].
We next examined which pathways were changed upon treatment using information from stress responses previously described in ciprofloxacin-treated bacteria [54]. We found that stavudine treatment led to a reduced pull-down of several proteins involved in general metabolism, while APIM-peptide treatment led to an increased pull-down of more proteins involved in nucleotide metabolism (Figure 4D). The latter result is interesting as APIM-peptide also affected the nucleotide pools more than stavudine.
In summary, the most likely explanation for the strong increase in mutation frequency is that stavudine induces SOS at doses below MIC, likely via interference with normal bacterial replications leading to replication stress.

4. Conclusions

Our findings presented herein indicate that anti-viral and anti-cancer NAs are capable of inducing bacterial mutagenesis and contributing to high AMR in HIV and cancer patients at clinically relevant concentrations. Stavudine induces mutagenesis by activating the E. coli SOS system and, thereby, the TLS polymerase Pol V. Thus, the combination of NAs/NRTIs and antibiotic treatment in already vulnerable patients may drive the development of and selection for, antibiotic-resistant bacteria. APIM-peptide treatment significantly reduced the stavudine-induced mutagenesis suggesting that this peptide could be a good drug candidate in the fight against AMR.

Author Contributions

Conceptualization, B.K.S., E.S.D., L.M., S.D. and M.O.; methodology, L.M.R., P.B. and M.O.; validation, B.K.S., S.B.R., L.M.R., P.B. and M.O.; formal analysis, B.K.S., A.S., P.B., L.M.R. and M.O.; investigation, B.K.S., S.B.R. and L.M.R.; resources, M.O. and P.B.; writing—original draft preparation, B.K.S., S.B.R., P.B. and M.O.; writing—review and editing, B.K.S., L.M.R., A.S., P.B., M.O., L.M. and S.D.; visualization, B.K.S., S.B.R., A.S., P.B. and M.O.; supervision, M.O.; project administration, M.O.; funding acquisition, M.O. and P.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by WACCBIP-DELTAS Africa grant (DEL-15-007) and Welcome Trust 107755/Z/15/Z (awarded to Gordon Akanzuwine Awandare), the Joint Research Committee between St. Olavs and Faculty of Medicine and Health Science, NTNU, the program NTNU Health, at Norwegian University of Science and Technology (NTNU) and the Trond Mohn foundation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data is deposited in PRIDE, project ID PXD025370 (Reviewer account Username: [email protected] Password: seTUpXp4).

Acknowledgments

We would like to thank Siri Bachke for technical assistance and NorStore/Notur Projects NS/NN9036K for data storage/computation services.

Conflicts of Interest

No conflict of interest to declare.

References

  1. CDCP. Antibiotic Resistance Threats in the United States; US Department of Health and Human Service: Atlanta, GA, USA, 2019.
  2. European Centre for Disease Prevention and Control. Surveillance of Antimicrobial Resistance in Europe 2018; European Centre for Disease Prevention and Control: Stockholm, Sweden, 2019.
  3. Tadesse, B.T.; Ashley, E.A.; Ongarello, S.; Havumaki, J.; Wijegoonewardena, M.; González, I.J.; Dittrich, S. Antimicrobial resistance in Africa: A systematic review. BMC Infect. Dis. 2017, 17, 1–17. [Google Scholar] [CrossRef]
  4. O’Neill, J. Review on Antimicrobial Resistance. Antimicrobial Resistance: Tackling a Crisis for the Health and Wealth of Nations. 2014. Available online: https://amr-review.org (accessed on 4 January 2020).
  5. Getahun, H.; Smith, I.; Trivedi, K.; Paulin, S.; Balkhy, H.H. Tackling antimicrobial resistance in the COVID-19 pandemic. Bull. World Health Organ. 2020, 98, 442–442A. [Google Scholar] [CrossRef]
  6. World Health Organization. Clinical Management of COVID-19: Interim Guidance, 27 May 2020; World Health Organization: Geneva, Switzerland , 2020. [Google Scholar]
  7. Rawson, T.M.; Moore, L.S.; Zhu, N.; Ranganathan, N.; Skolimowska, K.; Gilchrist, M.; Satta, G.; Cooke, G.; Holmes, A. Bacterial and fungal co-infection in individuals with coronavirus: A rapid review to support COVID-19 antimicrobial prescribing. Clin. Infect. Dis. 2020, 71, 2459–2468. [Google Scholar]
  8. Cirz, R.T.; Chin, J.K.; Andes, D.R.; De Crécy-Lagard, V.; Craig, W.A.; Romesberg, F.E. Inhibition of Mutation and Combating the Evolution of Antibiotic Resistance. PLoS Biol. 2005, 3, e176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Merrikh, H.; Kohli, R.M. Targeting evolution to inhibit antibiotic resistance. FEBS J. 2020, 287, 4341–4353. [Google Scholar] [CrossRef] [PubMed]
  10. Petrosino, J.F.; Galhardo, R.S.; Morales, L.D.; Rosenberg, S.M. Stress-induced beta-lactam antibiotic resistance mutation and sequences of stationary-phase mutations in the Escherichia coli chromosome. J. Bacteriol. 2009, 191, 5881–5889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Gawel, D.; Seed, P.C. Urinary tract infection drives genome instability in uropathogenic Escherichia coli and necessitates translesion synthesis DNA polymerase IV for virulence. Virulence 2011, 2, 222–232. [Google Scholar] [CrossRef] [Green Version]
  12. Altieri, A.S.; Kelman, Z. DNA Sliding Clamps as Therapeutic Targets. Front. Mol. Biosci. 2018, 5, 87. [Google Scholar] [CrossRef]
  13. Campos, P.E.; Suarez, P.G.; Sánchez, J.; Zavala, D.; Arevalo, J.; Ticona, E.; Nolan, C.M.; Hooton, T.M.; Holmes, K.K. Multidrug-resistant Mycobacterium tuberculosis in HIV-Infected Persons, Peru. Emerg. Infect. Dis. 2003, 9, 1571–1578. [Google Scholar] [CrossRef] [Green Version]
  14. Marami, D.; Balakrishnan, S.; Seyoum, B. Prevalence, Antimicrobial Susceptibility Pattern of Bacterial Isolates, and Associated Factors of Urinary Tract Infections among HIV-Positive Patients at Hiwot Fana Specialized University Hospital, Eastern Ethiopia. Can. J. Infect. Dis. Med Microbiol. 2019, 2019, 1–8. [Google Scholar] [CrossRef] [PubMed]
  15. Romana, K.; Rieger, A.; Presterl, E.; Wrba, T.; Diab-Elschahawi, M. Epidemiology of Urinary Tract Infections in HIV Positive Patients at a Tertiary Care University Hospital in Central Europe (2011–2016). Infect. Disord. Drug Targets 2018, 18, 199–206. [Google Scholar]
  16. Ssengooba, W.; Lukoye, D.; Meehan, C.J.; Kateete, D.P.; Joloba, M.L.; De Jong, B.C.; Cobelens, F.G.; Van Leth, F. Tuberculosis resistance-conferring mutations with fitness cost among HIV-positive individuals in Uganda. Int. J. Tuberc. Lung Dis. 2017, 21, 531–536. [Google Scholar] [CrossRef]
  17. Adeyemi, F.M.; Ako-nai, K.A.; Adejuyigbe, E.; Ebhodaghe, B.I.; Osho, P.O.; Oyeniyi, T.T.; Kassim, O.O. Molecular characterization and antibiotic resistance profiles of bacterial isolates cultured from HIV seropositive patients. Arch. Clin. Microbiol. 2015, 6, e2. [Google Scholar]
  18. Jordheim, L.P.; Durantel, D.; Zoulim, F.; Dumontet, C. Advances in the development of nucleoside and nucleotide analogues for cancer and viral diseases. Nat. Rev. Drug Discov. 2013, 12, 447–464. [Google Scholar] [CrossRef]
  19. Thomson, J.M.; Lamont, I.L. Nucleoside Analogues as Antibacterial Agents. Front. Microbiol. 2019, 10, 952. [Google Scholar] [CrossRef] [Green Version]
  20. Furman, P.A.; Fyfe, J.A.; Clair, M.H.S.; Weinhold, K.; Rideout, J.L.; Freeman, G.A.; Lehrman, S.N.; Bolognesi, D.P.; Broder, S.; Mitsuya, H. Phosphorylation of 3′-azido-3′-deoxythymidine and selective interaction of the 5′-triphosphate with human immunodeficiency virus reverse transcriptase. Proc. Natl. Acad. Sci. USA 1986, 83, 8333–8337. [Google Scholar] [CrossRef] [Green Version]
  21. Nedal, A.; Ræder, S.B.; Dalhus, B.; Helgesen, E.; Forstrøm, R.J.; Lindland, K.; Sumabe, B.K.; Martinsen, J.H.; Kragelund, B.B.; Skarstad, K.; et al. Peptides containing the PCNA interacting motif APIM bind to the β-clamp and inhibit bacterial growth and mutagenesis. Nucleic Acids Res. 2020, 48, 5540–5554. [Google Scholar] [CrossRef]
  22. CLSI. Methods for Dilution Antimicrobial Susceptibility Tests for Bacteria That Grow Aerobically; Approved Standard—Ninth Edition; CLSI Document M07-A9; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2012. [Google Scholar]
  23. Petrovic, V.; Olaisen, C.; Sharma, A.; Nepal, A.; Bugge, S.; Sundby, E.; Hoff, B.H.; Slupphaug, G.; Otterlei, M. On-column trypsinization allows for re-use of matrix in modified multiplexed inhibitor beads assay. Anal. Biochem. 2017, 523, 10–16. [Google Scholar] [CrossRef] [PubMed]
  24. Røst, L.M.; Thorfinnsdottir, L.B.; Kumar, K.; Fuchino, K.; Langørgen, I.E.; Bartosova, Z.; Kristiansen, K.A.; Bruheim, P. Absolute Quantification of the Central Carbon Metabolome in Eight Commonly Applied Prokaryotic and Eukaryotic Model Systems. Metabolites 2020, 10, 74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Stafsnes, M.H.; Røst, L.M.; Bruheim, P. Improved phosphometabolome profiling applying isotope dilution strategy and capillary ion chromatography-tandem mass spectrometry. J. Chromatogr. B 2018, 1083, 278–283. [Google Scholar] [CrossRef]
  26. Røst, L.M.; Shafaei, A.; Fuchino, K.; Bruheim, P. Zwitterionic HILIC tandem mass spectrometry with isotope dilution for rapid, sensitive and robust quantification of pyridine nucleotides in biological extracts. J. Chromatogr. B 2020, 1144, 122078. [Google Scholar] [CrossRef]
  27. Lien, S.K.; Kvitvang, H.F.N.; Bruheim, P. Utilization of a deuterated derivatization agent to synthesize internal standards for gas chromatography–tandem mass spectrometry quantification of silylated metabolites. J. Chromatogr. A 2012, 1247, 118–124. [Google Scholar] [CrossRef]
  28. Chong, J.; Wishart, D.S.; Xia, J. Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr. Protoc. Bioinform. 2019, 68, e86. [Google Scholar] [CrossRef]
  29. Tyanova, S.; Temu, T.; Cox, J. The MaxQuant computational platform for mass spectrometry-based shotgun proteomics. Nat. Protoc. 2016, 11, 2301–2319. [Google Scholar] [CrossRef] [PubMed]
  30. Geiszler, D.J.; Kong, A.T.; Avtonomov, D.M.; Yu, F.; Leprevost, F.D.V.; Nesvizhskii, A.I. PTM-Shepherd: Analysis and Summarization of Post-Translational and Chemical Modifications From Open Search Results. Mol. Cell. Proteom. 2021, 20, 100018. [Google Scholar] [CrossRef]
  31. Cox, J.; Hein, M.; Luber, C.A.; Paron, I.; Nagaraj, N.; Mann, M. Accurate Proteome-wide Label-free Quantification by Delayed Normalization and Maximal Peptide Ratio Extraction, Termed MaxLFQ. Mol. Cell. Proteom. 2014, 13, 2513–2526. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Gibbons, J.D.; Chakraborti, S. Nonparametric Statistical Inference, 6th ed.; CRC Press: Boca Raton, FL, USA, 2021. [Google Scholar]
  33. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  34. Perez-Riverol, Y.; Csordas, A.; Bai, J.; Bernal-Llinares, M.; Hewapathirana, S.; Kundu, D.J.; Inuganti, A.; Griss, J.; Mayer, G.; Eisenacher, M.; et al. The PRIDE database and related tools and resources in 2019: Improving support for quantification data. Nucleic Acids Res. 2019, 47, D442–D450. [Google Scholar] [CrossRef] [PubMed]
  35. Sandrini, M.P.B.; Clausen, A.R.; On, S.L.W.; Aarestrup, F.M.; Munch-Petersen, B.; Piškur, J. Nucleoside analogues are activated by bacterial deoxyribonucleoside kinases in a species-specific manner. J. Antimicrob. Chemother. 2007, 60, 510–520. [Google Scholar] [CrossRef] [PubMed]
  36. Shilaih, M.; Angst, D.C.; Marzel, A.; Bonhoeffer, S.; Günthard, H.F.; Kouyos, R.D. Short communication Antibacterial effects of antiretrovirals, potential implications for microbiome studies in HIV. Antivir. Ther. 2018, 23, 91–94. [Google Scholar] [CrossRef]
  37. Hurst, M.; Noble, S. Stavudine. Drugs 1999, 58, 919–949. [Google Scholar] [CrossRef]
  38. Cruciani, M.; Liuzzi, G.; Chirianni, A.; Audagnotto, S.; Bonora, S.; Di Biagio, A.; Sinicco, A.; Bassetti, M.; Gatti, G. Penetration of didanosine in semen of HIV-1-infected men. J. Antimicrob. Chemother. 2006, 57, 1244–1247. [Google Scholar] [CrossRef] [Green Version]
  39. Bang, L.M.; Scott, L.J. Emtricitabine. Drugs 2003, 63, 2413–2424. [Google Scholar] [CrossRef]
  40. Calcagno, A.; De Requena, D.G.; Simiele, M.; D’Avolio, A.; Tettoni, M.C.; Salassa, B.; Orofino, G.; Bramato, C.; Libanore, V.; Motta, I.; et al. Tenofovir Plasma Concentrations According to Companion Drugs: A Cross-Sectional Study of HIV-Positive Patients with Normal Renal Function. Antimicrob. Agents Chemother. 2013, 57, 1840–1843. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Capparelli, E.V.; Letendre, S.L.; Ellis, R.J.; Patel, P.; Holland, D.; McCutchan, J.A. Population Pharmacokinetics of Abacavir in Plasma and Cerebrospinal Fluid. Antimicrob. Agents Chemother. 2005, 49, 2504–2506. [Google Scholar] [CrossRef] [Green Version]
  42. Yuen, G.J.; Lou, Y.; Bumgarner, N.F.; Bishop, J.P.; Smith, G.A.; Otto, V.R.; Hoelscher, D.D. Equivalent Steady-State Pharmacokinetics of Lamivudine in Plasma and Lamivudine Triphosphate within Cells following Administration of Lamivudine at 300 Milligrams Once Daily and 150 Milligrams Twice Daily. Antimicrob. Agents Chemother. 2004, 48, 176–182. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Garré, B.; Shebany, K.; Gryspeerdt, A.; Baert, K.; Van Der Meulen, K.; Nauwynck, H.; Deprez, P.; De Backer, P.; Croubels, S. Pharmacokinetics of Acyclovir after Intravenous Infusion of Acyclovir and after Oral Administration of Acyclovir and Its Prodrug Valacyclovir in Healthy Adult Horses. Antimicrob. Agents Chemother. 2007, 51, 4308–4314. [Google Scholar] [CrossRef] [Green Version]
  44. Lerner, A.M.; Bailey, E.J. Concentrations of Idoxuridine in Serum, Urine, and Cerebrospinal Fluid of Patients with Suspected Diagnoses of Herpesvirus hominis Encephalitis. J. Clin. Investig. 1972, 51, 45–49. [Google Scholar] [CrossRef] [Green Version]
  45. McEvoy, G.K. American Hospital Formulary Service-Drug Information 92; American Society of Hospital Pharmacists Inc.: Bethesda, MD, USA, 1993. [Google Scholar]
  46. National Center for Biotechnology Information. In PubChem Compound Summary for CID 21704, Vidarabine; National Center for Biotechnology Information: Bethesda, MD, USA, 2020.
  47. Robinson, D.M.; Scott, L.J.; Plosker, G.L. Entecavir. Drugs 2006, 66, 1605–1622. [Google Scholar] [CrossRef] [PubMed]
  48. Kosovec, J.E.; Egorin, M.J.; Gjurich, S.; Beumer, J.H. Quantitation of 5-fluorouracil (5-FU) in human plasma by liquid chromatography/electrospray ionization tandem mass spectrometry. Rapid Commun. Mass Spectrom. 2007, 22, 224–230. [Google Scholar] [CrossRef]
  49. Gandhi, V.; Plunkett, W. Cellular and Clinical Pharmacology of Fludarabine. Clin. Pharmacokinet. 2002, 41, 93–103. [Google Scholar] [CrossRef] [PubMed]
  50. Hande, K.R.; Stein, R.S.; McDonough, D.A.; Greco, F.A.; Wolff, S.N. Effects of high-dose cytarabine. Clin. Pharmacol. Ther. 1982, 31, 669–674. [Google Scholar] [CrossRef]
  51. Sutton, M.D.; Smith, B.T.; Godoy, V.G.; Walker, G.C. The SOS response: Recent insights into umuDC-dependent mutagenesis and DNA damage tolerance. Annu. Rev. Genet. 2000, 34, 479–497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Røst, L.M.; Olaisen, C.; Sharma, A.; Nedal, A.; Petrovic, V.; Kvitvang, H.F.N.; Bruheim, P.; Otterlei, M. PCNA has specific functions in regulation of metabolism in haematological cells. bioRxiv 2020. [Google Scholar] [CrossRef]
  53. Simmons, L.A.; Foti, J.J.; Cohen, S.E.; Walker, G.C. The SOS Regulatory Network. EcoSal Plus 2008, 3. [Google Scholar] [CrossRef] [Green Version]
  54. Cirz, R.T.; Jones, M.B.; Gingles, N.A.; Minogue, T.D.; Jarrahi, B.; Peterson, S.N.; Romesberg, F.E. Complete and SOS-Mediated Response of Staphylococcus aureus to the Antibiotic Ciprofloxacin. J. Bacteriol. 2007, 189, 531–539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Biomolecules 11 00843 g001 550
Figure 1. Nucleoside analogues (NAs) induce mutagenesis in Escherichia coli MG1655. Data represents the average increase in mutation frequency in E. coli MG1655 cells detected in the RifR assay after treatment, compared to untreated control. (A) Mutation frequencies determined using 0.5× MIC or 0.25× of maximal examined concentration (256 µg/mL) (Table 1) of the different NAs. The relative increase in mutation frequency is given from 5 independent biological replicas. Stavudine (Stav, 32 µg/mL), didanosine (Dida, 4 µg/mL), emitrcitabine (Emtri, 64 µg/mL), tenofovir (Teno, 64 µg/mL), abacavir (Abaca, 64 µg/mL), lamuvudine (Lamu, 64 µg/mL, acyclovir (Acyc, 64 µg/mL), idocuridine (Idox, 64 µg/mL), vidarabine (Vida, 64 µg/mL), entecavir (Ente, 64 µg/mL), 5-fluro uracil (5-FU, 2 µg/mL), fludarabine, (Fluda, 64 µg/mL) and cytarabine (Cyta, 64 µg/mL). (B) Mutation frequency at plasma concentrations of Stav (0.32 µg/mL and 0.016 µg/mL), Dida (0.5 µg/mL), Emtri (10 µg/mL) and Fluda (1.1 µg/mL). The relative increase in mutation frequency is given from a minimum of 3 independent biological replicas, except for 0.32 µg/mL Stav, which are from 14 independent biological replicas. Each biological replica consists of 2–3 technical replicates. Unpaired Student’s t-test, two-tailed, **** p < 0.0001, *** p < 0.001, ** p < 0.01 and * p < 0.05.
Figure 1. Nucleoside analogues (NAs) induce mutagenesis in Escherichia coli MG1655. Data represents the average increase in mutation frequency in E. coli MG1655 cells detected in the RifR assay after treatment, compared to untreated control. (A) Mutation frequencies determined using 0.5× MIC or 0.25× of maximal examined concentration (256 µg/mL) (Table 1) of the different NAs. The relative increase in mutation frequency is given from 5 independent biological replicas. Stavudine (Stav, 32 µg/mL), didanosine (Dida, 4 µg/mL), emitrcitabine (Emtri, 64 µg/mL), tenofovir (Teno, 64 µg/mL), abacavir (Abaca, 64 µg/mL), lamuvudine (Lamu, 64 µg/mL, acyclovir (Acyc, 64 µg/mL), idocuridine (Idox, 64 µg/mL), vidarabine (Vida, 64 µg/mL), entecavir (Ente, 64 µg/mL), 5-fluro uracil (5-FU, 2 µg/mL), fludarabine, (Fluda, 64 µg/mL) and cytarabine (Cyta, 64 µg/mL). (B) Mutation frequency at plasma concentrations of Stav (0.32 µg/mL and 0.016 µg/mL), Dida (0.5 µg/mL), Emtri (10 µg/mL) and Fluda (1.1 µg/mL). The relative increase in mutation frequency is given from a minimum of 3 independent biological replicas, except for 0.32 µg/mL Stav, which are from 14 independent biological replicas. Each biological replica consists of 2–3 technical replicates. Unpaired Student’s t-test, two-tailed, **** p < 0.0001, *** p < 0.001, ** p < 0.01 and * p < 0.05.
Biomolecules 11 00843 g001
Biomolecules 11 00843 g002 550
Figure 2. Stavudine-induced mutagenesis in Escherichia coli MG1655 is driven by the translesion synthesis (TLS) Pol V. (A) Data represent the average increase in mutation frequency using the RifR assay after treatment of stavudine (Stav, 0.32 µg/mL), compared to untreated control for E. coli MG1655 wild type (WT) and TLS polymerase deletion strains (ΔPol V and ΔPol II/Pol IV). (B) Data represent the average increase in mutation frequency using the RifR assay after treatment of stavudine (Stav, 0.32 µg/mL), APIM-peptide (26 µg/mL) and a combination (Stav + APIM-peptide), compared to untreated control for E. coli MG1655 wild type (WT). The mutation frequency from the stavudine-treated WT cells (WT and Stav) are from 14 independent biological replicas. The mutation frequencies from the remaining samples (ΔPol V and ΔPol II/Pol IV, APIM-peptide and combination) are from 5 independent biological replicas, where 3 technical replicates were plated from each biological replicate. Unpaired Student’s t-test, two-tailed, **** p < 0.0001 and * p < 0.05.
Figure 2. Stavudine-induced mutagenesis in Escherichia coli MG1655 is driven by the translesion synthesis (TLS) Pol V. (A) Data represent the average increase in mutation frequency using the RifR assay after treatment of stavudine (Stav, 0.32 µg/mL), compared to untreated control for E. coli MG1655 wild type (WT) and TLS polymerase deletion strains (ΔPol V and ΔPol II/Pol IV). (B) Data represent the average increase in mutation frequency using the RifR assay after treatment of stavudine (Stav, 0.32 µg/mL), APIM-peptide (26 µg/mL) and a combination (Stav + APIM-peptide), compared to untreated control for E. coli MG1655 wild type (WT). The mutation frequency from the stavudine-treated WT cells (WT and Stav) are from 14 independent biological replicas. The mutation frequencies from the remaining samples (ΔPol V and ΔPol II/Pol IV, APIM-peptide and combination) are from 5 independent biological replicas, where 3 technical replicates were plated from each biological replicate. Unpaired Student’s t-test, two-tailed, **** p < 0.0001 and * p < 0.05.
Biomolecules 11 00843 g002
Biomolecules 11 00843 g003 550
Figure 3. Stavudine treatment of Escherichia coli MG1655 leads to small changes in glycolytic, PPP and TCA intermediates and in nucleoside phosphate pools. (A) Scores plot (PC1 and PC2) from PCA of intracellular metabolite concentrations, as measured by capIC-MS/MS, for the four treatment groups: stavudine (Stav, 0.32 µg/mL), APIM-peptide (26 µg/mL) and combination (Stav + APIM-peptide) and untreated control harvested two hours after treatment, with 95% confidence regions indicated. (B) Log2 fold change of treated cultures vs. untreated control. Data are averages of four independent biological replicas normalized by sum.
Figure 3. Stavudine treatment of Escherichia coli MG1655 leads to small changes in glycolytic, PPP and TCA intermediates and in nucleoside phosphate pools. (A) Scores plot (PC1 and PC2) from PCA of intracellular metabolite concentrations, as measured by capIC-MS/MS, for the four treatment groups: stavudine (Stav, 0.32 µg/mL), APIM-peptide (26 µg/mL) and combination (Stav + APIM-peptide) and untreated control harvested two hours after treatment, with 95% confidence regions indicated. (B) Log2 fold change of treated cultures vs. untreated control. Data are averages of four independent biological replicas normalized by sum.
Biomolecules 11 00843 g003
Biomolecules 11 00843 g004 550
Figure 4. Stavudine induces SOS in Escherichia coli. Results of the MIB assay of extract from E. coli MG1655 treated with stavudine (Stav, 0.32 µg/mL), APIM-peptide (26 µg/mL) and the combination (Stav + APIM-peptide) for 2 h. Data from three independent biological replicas are shown. Protein concentrations are given as log2 fold change relative to untreated control. (A) Venn diagrams displaying the number of proteins significantly changed in pull-downs from treated cells relative to untreated control according to the Wilcoxon signed-rank test. (B) Venn diagram showing an overlap in the different treatment groups with proteins involved in SOS [53]. (C) Heatmaps displaying log2 fold change in protein concentration of SOS proteins relative to untreated control in the different treatment groups. The grey boxes represent no data. (D) Heatmaps displaying log2 fold change in proteins concentration of proteins involved in general metabolism, nucleotide metabolism, lipid biosynthesis and DNA metabolism relative to untreated control in the different treatment groups. * Proteins that are significantly changed according to the Wilcoxon signed-rank test (i.e., changed the same way in all three replicas).
Figure 4. Stavudine induces SOS in Escherichia coli. Results of the MIB assay of extract from E. coli MG1655 treated with stavudine (Stav, 0.32 µg/mL), APIM-peptide (26 µg/mL) and the combination (Stav + APIM-peptide) for 2 h. Data from three independent biological replicas are shown. Protein concentrations are given as log2 fold change relative to untreated control. (A) Venn diagrams displaying the number of proteins significantly changed in pull-downs from treated cells relative to untreated control according to the Wilcoxon signed-rank test. (B) Venn diagram showing an overlap in the different treatment groups with proteins involved in SOS [53]. (C) Heatmaps displaying log2 fold change in protein concentration of SOS proteins relative to untreated control in the different treatment groups. The grey boxes represent no data. (D) Heatmaps displaying log2 fold change in proteins concentration of proteins involved in general metabolism, nucleotide metabolism, lipid biosynthesis and DNA metabolism relative to untreated control in the different treatment groups. * Proteins that are significantly changed according to the Wilcoxon signed-rank test (i.e., changed the same way in all three replicas).
Biomolecules 11 00843 g004
Table
Table 1. Nucleoside analogues (NAs) and their minimal inhibitory concentrations (MICs) in Escherichia coli MG1655.
Table 1. Nucleoside analogues (NAs) and their minimal inhibitory concentrations (MICs) in Escherichia coli MG1655.
Nucleoside Analogue (NA)Reported Plasma Concentration (Pl.CR.) (µg/mL)MIC (µg/mL)StructureUse
Stavudine
(Stav)		≤0.9
[37]		64 Biomolecules 11 00843 i001HIV, NRTI. dTTP analogue
Didanosine
(Dida)		0.05–0.86
[38]		8 Biomolecules 11 00843 i002HIV NRTI. dATP analogue
Emtricitabine
(Emtri)		1.01–11.5
[39]		>256 Biomolecules 11 00843 i003HIV NRTI. dCTP analogue
Tenofovir
(Teno)		0.050–0.077
[40]		>256 Biomolecules 11 00843 i004HIV NRTI. dTTP analogue
Abacavir
(Abaca)		0.04–1.13
[41]		>256 Biomolecules 11 00843 i005HIV NRTI. dGTP analogue
Lamivudine
(Lamu)		0.02–0.11
[42]		>256 Biomolecules 11 00843 i006HIV-1 NRTI. dTTP analogue
Acyclovir
(Acyc)		10.0
[43]		>256 Biomolecules 11 00843 i007Herpesviruses. dGTP analogue
Idoxuridine
(Idox)		10.0–36.0
[44]		>256 Biomolecules 11 00843 i008Herpesviruses. dUTP analogue. Thymidylate phosphatase inhibitor.
Vidarabine
(Vida)		0.2–0.4
[45,46]		>256 Biomolecules 11 00843 i009Herpesviruses. dATP analogue
Entecavir
(Ente)		0.015–0.036
[47]		>256 Biomolecules 11 00843 i010Hepatitis virus. dGTP analogue
5-Fluoro Uracil
(5-FU)		0.025–8.0
[48]		4 Biomolecules 11 00843 i011Multiple cancers. Thymidylate synthase inhibitor
Fludarabine
(Fluda)		0.4–2.0
[49]		>256 Biomolecules 11 00843 i012Leukaemia and lymphoma. dATP analogue
Cytarabine
(Cyta)		~486.0
[50]		>256 Biomolecules 11 00843 i013Leukaemia. dCTP analogue
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sumabe, B.K.; Ræder, S.B.; Røst, L.M.; Sharma, A.; Donkor, E.S.; Mosi, L.; Duodu, S.; Bruheim, P.; Otterlei, M. Nucleoside Analogues Are Potent Inducers of Pol V-mediated Mutagenesis. Biomolecules 2021, 11, 843. https://doi.org/10.3390/biom11060843

AMA Style

Sumabe BK, Ræder SB, Røst LM, Sharma A, Donkor ES, Mosi L, Duodu S, Bruheim P, Otterlei M. Nucleoside Analogues Are Potent Inducers of Pol V-mediated Mutagenesis. Biomolecules. 2021; 11(6):843. https://doi.org/10.3390/biom11060843

Chicago/Turabian Style

Sumabe, Balagra Kasim, Synnøve Brandt Ræder, Lisa Marie Røst, Animesh Sharma, Eric S. Donkor, Lydia Mosi, Samuel Duodu, Per Bruheim, and Marit Otterlei. 2021. "Nucleoside Analogues Are Potent Inducers of Pol V-mediated Mutagenesis" Biomolecules 11, no. 6: 843. https://doi.org/10.3390/biom11060843

APA Style

Sumabe, B. K., Ræder, S. B., Røst, L. M., Sharma, A., Donkor, E. S., Mosi, L., Duodu, S., Bruheim, P., & Otterlei, M. (2021). Nucleoside Analogues Are Potent Inducers of Pol V-mediated Mutagenesis. Biomolecules, 11(6), 843. https://doi.org/10.3390/biom11060843

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop