- Department of Biological Sciences, Illinois Institute of Technology, Chicago, IL, United States
Pseudomonas aeruginosa is a common cause of urinary tract infections by strains that are often multidrug resistant, representing a major challenge to the world’s health care system. This microorganism has a highly adaptable metabolism that allows it to colonize many environments, including the urinary tract. In this work, we have characterized the metabolic strategies used by stationary phase P. aeruginosa cells cultivated in urine-like media to understand the adaptations used by this microorganism to survive and produce disease. Our proteomics results show that cells rely on the Entner-Duodoroff pathway, pentose phosphate pathway, the Krebs cycle/ glyoxylate shunt and the aerobic oxidative phosphorylation to survive in urine-like media and other conditions. A deep characterization of the oxidative phosphorylation showed that the respiratory rate of stationary phase cells is increased 3–4 times compared to cells in the logarithmic phase of growth, indicating that the aerobic metabolism plays critical roles in the stationary phase of cells grown in urine like media. Moreover, the data show that respiratory complex III, succinate dehydrogenase and the NADH dehydrogenase NQR have important functions and could be used as targets to develop new antibiotics against this bacterium.
1 Introduction
Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen, ubiquitously distributed in natural and clinical environments (Bodey et al., 1983; Wiehlmann et al., 2007; Selezska et al., 2012). This bacterium relies on its highly flexible metabolism to survive under diverse and often harsh environmental conditions, which include antibiotic exposure (Livermore, 2002; Laborda et al., 2021). Pseudomonas aeruginosa has intrinsic antibiotic resistance mechanisms and determinants for disinfectant resistance (Murray et al., 2015). Due to its natural high adaptability, P. aeruginosa is one of the most common pathogens in hospitals (de Abreu et al., 2014; Reynolds and Kollef, 2021). The WHO and CDC list multidrug-resistant P. aeruginosa as one of the greatest threats to the global health care system and have called for the urgent need to develop new antibiotics against this pathogen (WHO, 2017; CDC, 2019). In 2017, multidrug-resistant P. aeruginosa caused more than 30,000 infections in the US with a mortality rate of 10%, producing estimated economic losses of $800 million (CDC, 2019). Moreover, P. aeruginosa is particularly relevant as the second most common cause of catheter-associated urinary tract infections (CAUTIs) (Sievert et al., 2013; Weiner et al., 2016; Weiner-Lastinger et al., 2020), causing more than 14,000 cases, 14% of all CAUTIs between 2015 to 2017 (CDC, 2019).
Pseudomonas aeruginosa has an unparalleled ability to adapt to different environmental conditions, due to its large genome, coding for a variety of virulence factors and for a greatly adaptable metabolism, and regulation capacity (Stover et al., 2000; Mikkelsen et al., 2011). Pseudomonas aeruginosa contains one of the most complex and highly branched respiratory chains (Figure 1A) that includes 17 predicted dehydrogenases, five aerobic and two anaerobic terminal oxidases, providing the flexibility for ATP generation that allows its survival under a variety of environments (Williams et al., 2006). Pseudomonas aeruginosa contains three seemingly redundant NADH dehydrogenases that catalyze electron transfer form NADH to ubiquinone: Complex I, NDH-2 and NQR (Torres et al., 2019; Liang et al., 2020). Complex I couples NADH oxidation to proton pumping and contains 14–15 subunits (nuoA-N) (Yagi, 1991; Baradaran et al., 2013). The non-proton pumping NDH-2 is coded by a single gene and contains FAD as its sole cofactor (Yagi, 1991; Kerscher et al., 2008). NQR is encoded by nqrA-F genes and contains six confirmed cofactors (Juárez and Barquera, 2012; Reyes-Prieto et al., 2014; Raba et al., 2018), and in this microorganism it functions as a proton pump (Raba et al., 2018), rather than a sodium pump as in other cases (Reyes-Prieto et al., 2014; Dibrov et al., 2017b). Succinate dehydrogenase (SDH), lactate dehydrogenase, malate dehydrogenase, glucose dehydrogenase and other dehydrogenases provide electrons to the respiratory chain to adapt to different environmental conditions (Kretzschmar et al., 2002; Frimmersdorf et al., 2010; Wang et al., 2018; McGill et al., 2021). In addition to the diversity of dehydrogenases, P. aeruginosa contains five aerobic terminal oxidases, including three COX-type oxidases: one aa3 -type oxidase (homologous to mitochondrial complex IV) and two cbb3 -type oxidases (Jo et al., 2014), which are connected to the ubiquinol pool through complex III. Pseudomonas aeruginosa also contains two ubiquinol oxidases: a bo3-type oxidase (CYO), which contains 5 subunits (Arai, 2011), and a cyanide-insensitive bd oxidase (CIO), containing two subunits (Arai et al., 2014; Jo et al., 2014).
Figure 1. Pseudomonas aeruginosa aerobic respiratory chain. Respiratory enzymes annotated in P. aeruginosa genome (A). Main components of the respiratory chain in LB (B). mAUM (C). NQR, ion-pumping NADH: ubiquinone oxidoreductase; CI, complex I; NDH-2, NADH: ubiquinone oxidoreductase II; SDH, succinate dehydrogenase; LLD, L-lactate dehydrogenase; MQO, malate: quinone oxidoreductase; GCD, glucose dehydrogenase; CIII, complex III (the bc1 complex); CIO, cyanide-insensitive bd oxidase; CYO, bo3-type oxidase; COX, three cytochrome c oxidases; LOG, logarithmic phase; STA, stationary phase.
The main aim of the study is to understand the metabolic adaptations of P. aeruginosa to the urine. In a previous report, we characterized the function of P. aeruginosa respiratory enzymes in the logarithmic phase of growth of cells cultured in modified artificial urine media (mAUM), which was recently developed by our group and offers a reliable and stable medium to study the bacterial adaptations to urine (Liang et al., 2020). In this report, we have expanded the previous studies to the stationary phase of growth, as this condition might represent a relevant stage in CAUTIs and UTIs in general, since the majority of biofilm colonizers are static organisms (Paula et al., 2020). Our results offer a clear picture of the role and relevance of the metabolic pathways and respiratory enzymes employed by this microorganism to survive in urine-like media (Figure 1C). Moreover, the data reveal that NQR, SDH and complex III play critical roles in cell physiology and represent potential antibiotic targets against CAUTIs produced by P. aeruginosa (Figures 1B,C).
2 Materials and methods
2.1 Bacterial growth and membrane purification
Pseudomonas aeruginosa PAO1 was grown in LB media and in freshly-prepared mAUM (Table 1; Supplementary Figure S1; Liang et al., 2020) in 2 L flasks. Cultures were initiated by diluting overnight cultures (1:1000) in 1 L of media. Cells were grown at 37°C in baffled flasks and agitated at 250 rpm. Cultures were collected at the logarithmic (7 h) and stationary phases of growth (24 h in LB broth and 12 h in mAUM). Cells were harvested and washed twice in KHE buffer (150 mM KCl, 20 mM HEPES, 1 mM EDTA, pH 7.5). Bacteria cells were disrupted at 16,000 psi using the Avestin Emulsiflex-C5. The cell suspension was centrifuged at 10,000 × g to eliminate cell debris. The membrane fraction was obtained by ultracentrifugation at 100,000 × g. The membrane pellet was washed with KHE buffer and stored at −80°C.
Wild type PAO1 was obtained from Dr. John C. Alverdy, University of Chicago. Growth of wild-type P. aeruginosa was compared to that of the Δbc1, ΔcioB, ΔccoN-1 and ΔccoN-2 mutants. The mutants were obtained from the P. aeruginosa PA14 Transposon Insertion Mutants Library (Liberati et al., 2006), kindly provided by Dr. Martin Schuster and Dr. Claudia Hässe from Oregon State University. Cultures of PAO1 wild type strains and the mutants of PA14 were grown in LB and mAUM at an inoculation ratio of 1:40 in 4 mL and incubated at 37°C under vigorous agitation, 250 rpm. Growth curves were analyzed using a logistic function to calculate the yield of the growth, the growth rate and the time of lag phase (Kahm et al., 2010).
2.2 Respiratory activity
Respiratory activities were measured in P. aeruginosa membranes (0.2 mg/mL) using a Clark-type oxygen electrode (YSI 5300) in a customized oximetric chamber (Liang et al., 2018, 2020). Experiments were carried out in KHE buffer at 37°C with different substrates, including 200 μM NADH, 20 mM succinate, 10 mM L-malate, 10 mM lactate (racemic mixture) and 10 mM D-glucose. Ubiquinol oxidase was measured in the presence of 50 μM ubiquinone-1 and 500 μM DTT (1,4-Dithiothreitol), which rapidly reduces ubiquinone (Chu et al., 2006; Furt et al., 2010; Gong et al., 2018). Cytochrome c oxidase activity was measured using 1 mM ascorbate and 100 μM TMPD (N,N,N′,N′-tetramethyl-para-phenylene-diamine) (Kimelberg and Nicholls, 1969). Cyanide titrations (0.1 μM to 10 mM) of the terminal oxidases were measured using NADH as a substrate with freshly prepared KCN.
2.3 NADH dehydrogenase activity
NADH: ubiquinone oxidoreductase activity (Qred) was determined in membranes in the presence of NADH or deamino-NADH with or without 1 μM rotenone to inhibit complex I, in KHE buffer containing 0.05% DDM (β-D-dodecyl maltoside), 50 mM ubiquinone-1 and 10 mM KCN to inhibit terminal oxidases (Matsushita et al., 1987; Managò et al., 2015; Liang et al., 2020). The activity was followed spectrophotometrically at 282 nm as described previously (Juárez et al., 2011), using an Agilent Cary 8,454 UV–vis spectrophotometer. The activities reported here represent initial rates, collected within two minutes after the reaction was initiated. The Qred % was normalized comparing the individual dehydrogenase activities to the total activity by all dehydrogenases.
2.4 Blue native PAGE
Blue native polyacrylamide gel electrophoresis (BN-PAGE) was used to separate native and functional P. aeruginosa respiratory enzymes (complex I, NDH-2 and NQR) based on their molecular weights (Wittig et al., 2006). Pseudomonas aeruginosa membrane samples were washed and resuspended in buffer containing 50 mM Bis-Tris, 500 mM aminocaproic acid, pH 7.0, and were solubilized at 4°C for 40 min, using b-D-dodecylmaltoside at a detergent: protein ratio of 2, as described previously (Liang et al., 2020). Solubilized membranes (50 μg of protein) were loaded into a 4 to 16% gradient acrylamide-bisacrylamide gel. NADH dehydrogenase activity was assayed in-gel with MTT as electron acceptor using NADH as the substrate as described previously (Liang et al., 2020). The reaction buffer contained 100 mM Tris, 140 μM NADH, 50 μM MTT (3-(4,5-Dimethylthiazol-2-yl)-2,5-Diphenyltetrazolium Bromide), pH 7.4 (Liang et al., 2020). The Coomassie-blue was then washed off the gel by 0.1% SDS after the MTT staining. The reduction of MTT by the NADH dehydrogenase produces a purple color on the location of the enzyme bands.
2.5 RNA isolation and cDNA synthesis
Bacteria were collected at the stationary phase in LB and mAUM, and were suspended in Purelink RNA Mini Kit lysis buffer, containing 10 mg/mL lysozyme, 10% SDS and 1% 2-mercaptoethanol. The lysate was homogenized by passing it through a 21-gauge needle 5 times. The supernatant was collected after centrifugation at 12,000 × g for 2 min at room temperature. Total RNA was extracted using the Purelink RNA Mini Kit, according to the manufacturer’s instructions (ThermoFisher Scientific, Waltham, MA). RNA concentration was measured by AccuBlue Broad Range RNA Quantitation Kit (Biotium, Fremon, CA) on a Qubit 3.0 fluorometer. cDNA was synthesized with The ProtoScript II First Strand cDNA Synthesis Kit (New England Biolabs Ipswich, MA), using the random primer mix.
2.6 Real time quantitative PCR
RT-qPCR was performed following the guidelines for minimum information of quantitative real-time PCR for publication (Bustin et al., 2009). The rpoS gene was selected as reference, which shows transcriptional stability under carbon starvation (Alqarni et al., 2016). The forward and reverse primers (Table 2) were designed using the Primer3 software (Kõressaar et al., 2018). 100 ng cDNA were used as template in each amplification. Amplification and quantification were done by the fluorescent dye-based double-strand DNA detective PCR technology, using the iTaq Universal SYBR Green Supermix kit (Bio-Rad, Hercules, CA). The real-time PCR reaction was run in a Bio-Rad CFX Connect Real-Time PCR Detection System following the manufacturer’s instructions (Bio-Rad). The relative quantification of PCR expression was analyzed by the 2-ΔΔCT method that compared the differences between the cycle threshold of target gene and reference gene (rpoS) under the control condition (LB) and experimental condition (mAUM) (Livak and Schmittgen, 2001).
2.7 Proteomics
Bacterial membrane samples (20 μg of protein) were prepared using the Filter Aided Sample Preparation protocol, employing 30 K NMWL centrifugal filter units. Samples were reduced with 20 mM dithiothreitol, followed by 20 min alkylation with 50 mM iodoacetamide. The alkylated proteins underwent three washes with 200 μL of 8 M Urea in 0.1 M Ammonium Bicarbonate (ABC) and were subsequently equilibrated three times with 0.1 M ABC prior to trypsin digestion. Proteins were digested using Trypsin at an enzyme: protein ratio of 1:50 in 50 μL of 0.1 M ABC buffer at 37°C overnight. The resulting digested peptides were collected by centrifugation at 14,000 g for 20 min, followed by two elutions with 50 μL of 50 mM ABC and one elution with 50 μL of 0.5 M NaCl. The eluted peptides were desalted using C18 columns and were dried. The digested proteins were then reconstituted in a 30 μL solution of 5% acetonitrile and 0.1% formic acid buffer for LC–MS analysis.
Five microliters of digested peptides were analyzed using a Q Exactive HF mass spectrometer coupled to an UltiMate 3,000 RSLC nanosystem, featuring a Nanospray Flex Ion Source (Thermo Fisher Scientific). The digested peptides were initially loaded into a Waters nanoEase M/Z C18 (100 Å, 5 μm, 180 μm × 20 m) trap column and subsequently separated using a 75 μm × 150 mm Waters BEH C18 (130 Å, 1.7 μm, 75 μm × 15 cm) column, at a flow rate of 300 nL/min. Solvent A consisted of 0.1% formic acid (FA) in water, while solvent B comprised 0.1% FA and 80% acetonitrile in water. The liquid chromatography gradient profile started with 5% of solvent B from 0 to 3 min, increased to 8% B at 3.2 min, ramped from 8 to 35% B within 110 min, elevated from 35 to 95% B in 115 min, underwent a 95% wash until 119.8 min, and finally equilibrated back to 5% B until 130 min.
Full MS scans were acquired using the Q Exactive mass spectrometer over the range of 350 to 1,400 m/z, operating at a resolution of 60,000 (at 200 m/z) from 10 min to 130 min. The AGC (automatic gain control) target value was set at 3×106 for the full scan. The top 15 most intense peaks with charge states 2, 3, 4, and 5 were subjected to fragmentation in the HCD (higher-energy collisional dissociation) collision cell with a normalized collision energy of 28%. These peaks were subsequently excluded for the next 30 s within a mass window of 1.4 m/z. Tandem mass spectra were collected in the mass analyzer at a resolution of 30,000, with an AGC target value set at 1×105. The ion selection threshold was established at 1×104 counts, while the maximum allowed ion injection time was 50 msec for both full scans and fragment ion scans. Mass spectrometry data were deposited in the MassIVE database under the accession number PXD047660.
Spectra data were searched against the Uniprot P. aeruginosa database using Mascot daemon (2.6.0, updated on 08/11/20). Results were entered into Scaffold Q + S software (v5.1.2, Proteome Software, Portland, OR) for compilation, normalization, and comparison of spectral counts, etc. The filtering criteria of protein identification were a 1% false discovery rate (FDR, statistical measure of the certainty of protein identification) of protein and peptide with 1 minimum peptide count. Peptide sequences and their assigned UniProt accessions were parsed from the peptide report with parse_scaffold_peptide_summary.py v0.1.0. Protein sequences and enzyme classification (EC) numbers for each UniProt accession were downloaded utilizing the Consortium’s API as implemented in grab_UniProt_data.py v0.3.0. The KEGG (Kanehisa et al., 2023) organism code for P. aeruginosa PAO1 (pae) was obtained from the KEGG organism catalog.1 KEGG pathways were reconstructed from mass spec data using the MassSpec2KEGG pipeline2 as described below. UniProt accessions were converted to KEGG orthologs (KOs) leveraging the KEGG API as implemented in UniProt_to_KEGG.py v0.5.0 as follows. UniProt accessions were mapped to NCBI accessions by performing DIAMOND (Buchfink et al., 2021) BLASTP homology searches against the P. aeruginosa PAO1 reference proteome (NCBI accession number GCF_000006765.1). NCBI accessions were converted to P. aeruginosa PAO1 specific KEGG IDs with the P. aeruginosa PAO1 specific KEGG ID/NCBI Protein ID conversion file from the KEGG database.
3 Results
3.1 Proteomic analysis
To characterize the metabolic strategies used by P. aeruginosa when grown in LB or mAUM, we carried out an analysis of the membrane and membrane-associated proteins found in the logarithmic and the stationary phases in these media. Of the 5,570 possible ORFs, a total of 2,184 proteins were identified in the four conditions, employing a stringent filtering threshold of 1% FDR for both peptide and protein recognition. In mAUM, 1,201 proteins were identified in mid-log growth phase and 1,226 proteins were identified in the stationary phase. In LB medium, 1,403 proteins were detected in the mid-log phase and 1,512 proteins in the stationary phase (Supplementary Data 1). Protein expression patterns are shown in a Venn diagram (Figure 2A), which illustrates the proteins identified across these conditions. Our data reveals 708 membrane or membrane-associated proteins that are recognized as core proteins, common among all growth conditions. Figure 2B illustrates the role of the core proteins, evaluated from the analysis of KEGG pathways, GO terms, and the Pseudomonas Community Annotation project (Winsor et al., 2016). As can be observed in Table 3, aerobic oxidative phosphorylation, the Entner-Duodoroff pathway, the pentose phosphate pathway, and the Krebs cycle/ glyoxylate shunt are core pathways in P. aeruginosa.
Figure 2. Proteomic data from P. aeruginosa PAO1 membrane isolated samples grown in LB and mAUM media. (A) Venn diagram of the four conditions depicting the numbers and overlaps of differentially identified proteins. (B) Bar graph showing the distribution of the protein function among the 708 core proteins, including those of the central carbon metabolism, indicated in orange. Hypothetical proteins are not included.
Table 3. Metabolic proteins related to Entner-Duodoroff pathway, gluconeogenesis, fermentation, pentose phosphate pathway, Krebs cycle, and glyoxylate cycle.
The respiratory proteins identified across all growth conditions are detailed in Table 4, comprising a total of 12 enzymes. Our data shows that P. aeruginosa expresses the three NADH dehydrogenases: NQR, complex I, and NDH-2, in all conditions. SDH, complex III, CYO oxidase, along with the two cbb3 oxidases are also core proteins. Additionally, lactate dehydrogenase (lld, PA4771) and malate dehydrogenase (mqo) are consistently found across all growth conditions, as well as F1-F0 ATP synthase. Interestingly, the A-type cytochrome c oxidase (coxB, PA0105) was identified during mAUM growth in the stationary phase (Supplementary Data 2). This enzyme is considered more energy efficient compared to other oxidases (Osamura et al., 2017). Its participation in these conditions could be beneficial considering the limited nutrients in the media. The cyanide-insensitive CIO oxidase peptides (cioA, PA3930) were identified in LB stationary phase as well as in AUM growth (Supplementary Data 2), aligning with increased cyanide tolerance in P. aeruginosa cultivated in AUM. Correspondingly, the early stages of LB growth provide enough nutrients for P. aeruginosa which will not trigger the production of cyanide and the expression of CIO. An intriguing discovery is the presence of glucose dehydrogenase (gcd, PA2290, Supplementary Data 2) across all four growth conditions, despite the absence of glucose in mAUM media. This suggests an activated gluconeogenesis pathway that potentially generates glucose as an energy source. Overall, the proteomics data show that P. aeruginosa carries a very active aerobic respiratory metabolism in urine-like media, which was further characterized in this report.
Table 4. Respiratory proteins identified among all growth conditions tested for Pseudomonas aeruginosa PAO1.
3.2 Dehydrogenase and oxidase activities
In a previous report, we characterized the respiratory chain of P. aeruginosa in the logarithmic phase of growth in conditions resembling urine, in mAUM media (Liang et al., 2020). However, in CAUTIs the growth rate of bacteria is likely low (Chambers et al., 1999), and the physiology of the cells might be similar to the stationary phase of growth in this and other bacteria (Trautner and Darouiche, 2004; Jacobsen et al., 2008). Thus, we aimed to characterize the P. aeruginosa respiratory chain in the stationary phase. Pseudomonas aeruginosa PAO1 plasma membranes of cells grown in mAUM were obtained in the stationary phase, where most of the cells are no longer dividing or dividing at a reduced rate and before the death phase has started. The respiratory activity was measured as the oxygen consumption rate in the presence of different substrates such as NADH, succinate, lactate, malate, glucose, and ethanol, to evaluate different dehydrogenases. In LB media, the activity of NADH dehydrogenases decreased by half in the stationary stage compared to the rapid growing stage (Figure 3A). On the other hand, in mAUM the activity of the NADH dehydrogenases increased more than threefold in the stationary stage compared to the rapid growth phase (Figure 3B). Succinate oxidation also increased in LB and mAUM in the stationary phase, with the highest increase in mAUM (>five-fold). Lactate and malate oxidation were not significantly modified in LB but increased four times in mAUM in the stationary phase. As reported before, glucose and ethanol oxidation did not contribute substantially to the respiratory activity (Liang et al., 2020).
Figure 3. Pseudomonas aeruginosa respiratory activity in LB and mAUM membranes. Oxygen consumption rate (OCR) was measured in bacterial membranes grown in LB (A) and mAUM (B) collected in logarithmic (black bars) and stationary (white bars) phases. The data from the logarithmic phase was obtained from a previous publication (Liang et al., 2020). Respiratory activity was measured in KHE buffer with the following substrates: 200 μM NADH, 10 mM succinate, 10 mM lactate, 10 mM glucose, 0.3% ethanol, 50 μM ubiquinol or 100 μM TMPD+ plus 5 mM ascorbic acid. n > 3 independent experiments for each group. Data represented as mean ± SD. *p < 0.05 and **p < 0.01 determined by t-test of the logarithmic versus stationary phase.
In addition to the dehydrogenases, the activity of the terminal aerobic oxidases was investigated using ubiquinol and TMPD-ascorbate as substrates. Ubiquinol oxidation reports the activity of CIO and CYO oxidases, as well as the activity of the pathway composed by complex III and COX oxidases, such as cytochrome aa3 and cbb3 oxidases (Liang et al., 2020), and TMPD-ascorbate oxidation reports the activity of cytochrome C oxidases directly (Liang et al., 2020). In cells in the stationary phase grown in LB media, the activity of aerobic ubiquinol oxidases doubled, while the activity of cytochrome C oxidases decreased by >3-fold (Figure 3A). In the stationary phase, the cells grown in mAUM show a 3-fold activation of the oxidases (Figure 3B). Overall, in urine-like media, P. aeruginosa undergoes deep metabolic changes in the stationary phase, which include a 3–4 -fold activation of the respiratory activity.
3.3 Participation of NADH dehydrogenases in respiratory activity
The participation of NADH dehydrogenases in the respiratory activity of stationary phase P. aeruginosa cells was also characterized. Although the three NADH dehydrogenases coded in the genome catalyze the same chemical reaction: electron transfer from NADH to ubiquinone, they can be differentiated by their substrate selectivity and inhibitor sensitivity. Complex I is able to use as substrates NADH and deamino-NADH equally well (Matsushita et al., 1987), and it is specifically inhibited by rotenone (Degli Esposti, 1998; Zickermann et al., 2000). NQR also uses NADH and deamino-NADH, but it is insensitive to rotenone (Zhou et al., 1999). On the other hand, NDH-2 exclusively uses NADH as substrate and is insensitive to rotenone (Matsushita et al., 1987; Zambrano and Kolter, 1993). Our results show that in all conditions tested, the NADH and deamino-NADH oxidation were very similar and were largely insensitive to rotenone, indicating that NQR is the main functional NADH dehydrogenase, sustaining 75–95% of the electron transfer, with a modest participation of complex I and a negligible function of NDH-2 (Figure 4A).
Figure 4. Participation of NADH dehydrogenases in P. aeruginosa metabolism in LB and mAUM media. (A) Normalized activity of the three NADH dehydrogenases in logarithmic (LOG) and stationary (STA) phases of cells grown in LB and mAUM media. The data from the logarithmic phase was obtained from a previous publication (Liang et al., 2020). (B) In-gel activity staining of NADH dehydrogenases in BN-PAGE. (C) Relative expression of the three NADH dehydrogenases, as determined by quantitative PCR in stationary phase of cells grown in LB (black bars) and mAUM (white bars).
3.4 Blue native gel electrophoresis
In addition to the substrate specificity, the three NADH dehydrogenases can be differentiated by their sizes. Complex I, NDH-2 and NQR have an estimated molecular weight of 530, 46, and 210 kDa, respectively (Liang et al., 2020), and can be separated electrophoretically using Blue Native PAGE, as described previously (Liang et al., 2020). The complexes can be stained in-gel following the precipitation of the electron acceptor MTT using NADH as substrate (Berridge et al., 2005). As shown in Figure 4B, the gel contains a prominent band of 210 kDa, which corresponds to NQR, as confirmed previously by our group using mass spectroscopy (Liang et al., 2020). This band increases 3.3-fold in the mAUM stationary phase membranes, consistent with the increase in respiratory activity discussed above. Another band of 110 kDa is also visible but as described before (data not shown), it does not contain any NADH dehydrogenase peptides (Liang et al., 2020), and does not show an increased intensity in the log state cells. Bands of 370 and 300 kDa were also evident in the gel, which may correspond to complex I subcomplexes. However, they have low intensity and required much longer to be stained, consistent with the low activity or complex I in these samples.
3.5 NADH dehydrogenase mRNA expression
To understand the genetic regulation of the respiratory chain of P. aeruginosa grown in LB and mAUM, we carried out a quantification of mRNAs of the three NADH dehydrogenases through qPCR, using the rpoS gene as a reference gene, which presents stability under carbon-limiting conditions (Alqarni et al., 2016) to normalize the results. Figure 4A shows that although there are changes in the activities of the three dehydrogenases depending on the media and the growth stage, there are no significant differences in the expression of these genes at the mRNA level (Figure 4C), indicating that other levels of genetic regulation may play major roles in the control of the respiratory activity in vivo.
3.6 Contribution of terminal oxidases to respiratory activity
As mentioned above, P. aeruginosa contains five different aerobic oxidases: CIO oxidase, CYO oxidase, and the pathway composed of complex III and cytochrome oxidases aa3, cbb3-1 and cbb3-2. Although these enzymes catalyze the same reaction: electron transfer from ubiquinol to oxygen, they can be differentiated by their sensitivity to cyanide. In a previous study we have shown that by conducting a KCN titration, the contribution of each oxidase can be unambiguously determined, since the IC50 for cytochrome oxidases (<1 μM), CYO oxidase (10–30 μM), and CIO oxidase (>5 mM) for cyanide are vastly different (Liang et al., 2020). KCN titrations were carried out in membranes obtained in the stationary phase of growth and the data was fitted to a three component Michaelis–Menten inhibition equation (Figures 5A,B) to obtain the contributions of each oxidase, as reported previously (Liang et al., 2020). In LB media we have reported that the main oxidase in the logarithmic phase is CYO oxidase. However, when the cells enter the stationary phase, the respiratory chain changes dramatically, with close to 80% of the activity carried by the cyanide-insensitive CIO oxidase (Figure 5C). Even though a three-fold increase in the respiratory activity was induced in the stationary phase of growth in cells grown in mAUM, the contributions of the different oxidases did not change significantly, with complex III and the COX oxidases contributing to 50% of the respiratory activity, compared to close to 60% in the log phase (Figure 5C).
Figure 5. Participation of terminal oxidases in P. aeruginosa metabolism in LB and mAUM media. KCN titration of the NADH-dependent respiratory activity of LB (A) and mAUM (B) membranes collected from logarithmic (black) and stationary (gray) phase cells. The data from the logarithmic phase was obtained from a previous publication (Liang et al., 2020). (C) Relative contribution of the oxidases to the respiratory activity of cells grown in LB (left) and mAUM (right) under logarithmic and stationary phases.
3.7 Studies on mutants
To corroborate the results of the functional analysis of the respiratory chain, we carried out an examination of the growth of mutants lacking genes of the respiratory complexes, using the PA14 library. The genes encoding the enzymes for the bc1 complex, cioB, ccoN-1, and ccoN-2 are identical between the PA14 and PAO1 strains, according to the PA14 Transposon Insertion Mutants Library,3 which provides the comparison data for each gene between the mutant of PA14 strain and its corresponding homologous gene in PAO1. A protein BLAST also reveals 100% identity for these enzymes in PAO1 and PA14 strains. As previously reported by our group (Liang et al., 2020) and others (Torres et al., 2019), the mutants of single NADH dehydrogenases do not produce noticeable effects on the growth in LB or mAUM, indicating that in the absence of any of them, the other two NADH dehydrogenases can compensate and rescue the growth. In contrast, the mutants of complex III (Δbc1) and CIO oxidase (ΔcioB) show significant defects in the growth. The growth of the Δbc1 mutant in LB and mAUM shows a 2-fold decrease in the growth rate and a modest decrease in the yield in LB media (Figures 6A,B), which indicates that complex III plays significant roles in both media. The growth of the ΔcioB mutant in LB showed a reduction in the yield (Figure 6C), which is consistent with the large role that this enzyme plays in the stationary phase in this medium, but noticeable effects were not found in mAUM (Figure 6D). The cbb3-1 (ΔccoN) and cbb3-2 (ΔccoN-2) mutants do not show significant changes in the growth profile (Figures 6E–H).
Figure 6. Growth of wild-type and mutant P. aeruginosa cells in LB and mAUM. Pseudomonas aeruginosa cells were grown in LB (A,C,E,G) and mAUM (B,D,F,H). The growth of wild-type cells (black line) is compared to the mutants (gray line) Δbc1 (A and B), ΔcioB (C,D), ΔccoN-1 (E,F) and ΔccoN-2 (G,H). Data represented as mean and SD.
4 Discussion
Pseudomonas aeruginosa is one of the main causal agents of multidrug-resistant infections and CAUTIs in the world, representing a significant health problem (Mittal et al., 2009; Weiner-Lastinger et al., 2020). Pseudomonas aeruginosa has a highly adaptive metabolism that allows it to survive in different environments, which likely contributes to the evolution of antibiotic resistance. In order to obtain new therapeutic treatments against P. aeruginosa UTIs and CAUTIs, it is necessary to understand its metabolic adaptations, in particular the roles played by the aerobic respiratory chain, of this notorious multi-drug resistant and ubiquitous pathogen in urine-like media. Previous studies demonstrate that the urine environment is aerobic, with oxygen’s concentration close to 100 μM, even in patients with urinary tract infections (Giannakopoulos et al., 1997), indicating that the aerobic respiratory chain has an important role in the physiology of this bacterium. Our previous characterization of the respiratory chain showed a very clear picture of the respiratory enzymes used by P. aeruginosa in the actively multiplying state. Moreover, our results clearly showed that transcriptomic and mutant studies must be supplemented with functional studies, as the results of gene expression profiles and mutant growth do not necessarily correlate with the functional architecture of the pathway (Liang et al., 2020). In this study, we have characterized the respiratory chain of P. aeruginosa in the stationary phase of growth in urine-like media, to understand the metabolism of this pathogen when CAUTIs or UTIs have been established. The stationary phase is particularly relevant in the clinic, as non-replicating cells might resemble the physiological state of cells in the urinary tract, as bacterial growth is prevented by the chemical composition of urine itself (Chambers et al., 1999; Bono et al., 2022) and by antibacterial coatings on the catheter (Andersen and Flores-Mireles, 2020). Moreover, long-term infection by this pathogen often involves plateauing in bacterial growth and biofilm formation (Trautner and Darouiche, 2004).
4.1 Metabolic activation in AUM
Urine is a complex fluid which could composed of more than 3,000 substance and varies among individuals and makes the study of an in vitro experimental condition difficult (Sarigul et al., 2019), especially because urine-like media is prone to salt precipitation (Mittal et al., 2009) and to drastic changes in pH, due to the chemical decomposition of urea (Levinsky and Berliner, 1959). AUM is made with the major components of human urine to generate a practical and reproducible formulation. In our previous study, we showed that mAUM, which was developed in our group, offers higher stability and reproducibility compared to other urine-like published media (Liang et al., 2020), and thus was selected to study and understand the physiology of bacteria in the urine.
Urine contains low nutrients, which limits energy production, explaining the reduced growth compared to the nutrient rich condition (Liang et al., 2020). Proteomics data show that P. aeruginosa has a very active aerobic metabolism and uses essentially the same pathways in LB and in mAUM media: the Entner-Duodoroff pathway, the pentose phosphate pathway, the Krebs cycle/ glyoxylate shunt and oxidative phosphorylation. Even though all respiratory enzymes were found as core proteins, present in all conditions evaluated, the regulation of the activity of some of these enzymes is tightly controlled by the cell and several of them are activated in urine-like media. Our data show that the respiratory activity of stationary phase cells grown in mAUM exhibits a significant increase of approximately 3 times compared to the activity in the logarithmic phase. In contrast with the expectation that the cells might become quiescent and perhaps decrease their metabolic rate due to the lack of nutrients, which is observed at a low extent in LB media (Figure 3), we noticed that in urine-like media the activity of all respiratory enzymes is higher. This metabolic activation could be related to the need to maximize the energy production from the limited resources available, in particular from lactate which is one of the main components of urine (Scheijen et al., 2012). Although major metabolic changes were found, our qPCR results on NQR did not point to changes in the transcriptional activity, indicating that other levels of protein or enzymatic regulations are involved in the respiratory activation of the cells.
4.2 Role of NADH dehydrogenases
In this study, we are showing that the NADH-dependent oxygen consumption rate increases significantly in the stationary phase in cells grown in mAUM. Through a functional characterization of the respiratory activity, we were able to determine that NQR is the main NADH dehydrogenase in all conditions tested, with a minor role by complex I and a negligible participation of NDH-2. These results are similar to the behavior that we have reported previously in the logarithmic phase of growth (Liang et al., 2020), which also show that NQR is the main NADH dehydrogenase. The reasons behind this metabolic strategy have been unclear until now, and we can start to elucidate the advantages of this pathway architecture. It is likely that NDH-2 plays a minor role in the respiratory activity since is not electrogenic and it is not coupled to the synthesis of ATP (Kerscher et al., 2008; Heikal et al., 2014). On the other hand, complex I and NQR pump protons with different stoichiometries, of 2 (Wikström and Hummer, 2012) and 1 H+ (Bogachev et al., 1997) per electron delivered to ubiquinone, respectively. While complex I appears to be more energy efficient, there are several advantages of NQR that could explain why P. aeruginosa uses it as the main NADH dehydrogenase. For instance, complex I contains two 2Fe-2S and six 4Fe-2S clusters (Yagi, 1991; Yagi et al., 1998; Kerscher et al., 2008), while NQR contains only two 2Fe-2S centers (Kishikawa et al., 2022). The requirement of minimum iron use of NQR supports the major role that it plays in pathogens (Reyes-Prieto et al., 2014), as iron is normally a limiting factor in pathogen’s growth (Cassat and Skaar, 2013). Moreover, the molecular weight of NQR is only half compared to complex I, which would require less energy and faster assembly. In addition, complex I is one of the main sources of reactive oxygen species (McLennan and Degli Esposti, 2000). While the rates of NADH oxidation and ubiquinone reduction are similar in NQR (Juárez et al., 2009), indicating that most electrons are delivered to the physiological substrate, in complex I the oxidation of NADH is 20–30 times higher compared to ubiquinone reduction (Fedor et al., 2017), which would produce significantly more oxygen radicals, and also would decrease the effective energy yield of the complex. Another important advantage of NQR is related to the resistance against self-poisoning. HQNO is a molecule produced by P. aeruginosa for quorum sensing and biofilm formation, which eliminates competing bacteria by inhibiting their respiratory activity (Déziel et al., 2004; Raba et al., 2018). Our previous study shows that P. aeruginosa NQR is only partially sensitive to HQNO inhibition, which allows the bacteria to survive autopoisoning and outcompete other bacteria (Raba et al., 2019). Even though NQR is the preferred NADH dehydrogenase by the bacteria, our previous studies (Liang et al., 2020) and other reports (Torres et al., 2019) indicate that its genetic inactivation does not produce significant changes in the growth parameters, which is puzzling considering the major role that it plays. However, by analyzing the contribution of all respiratory chain enzymes we can start to understand this and other phenotypes (see below).
4.3 Role of terminal oxidases
The role of aerobic terminal oxidases was also studied, as they should play a major part in the aerobic urinary medium. In LB medium, there is a major change in the metabolic strategy used between the two phases of growth, from the use of CYO as the main oxidase in the logarithmic phase, to the use of CIO in the stationary phase. In rich media, P. aeruginosa as well as many other bacteria such as P. putida, Escherichia coli, and Rhizobium etli (Iuchi et al., 1990; Ugidos et al., 2008; Lunak and Noel, 2015) use CYO oxidases, which could be a response to high oxygen concentration in the initial part of growth, as described before (Liang et al., 2020). As the cell growth rate decreases, the concentration of oxygen might decrease, and the cells start producing cyanide (Cody et al., 2009), which can reach up to 500 μM (Lenney and Gilchrist, 2011), triggering the switch to the cyanide-resistant CIO oxidase that allows the continuous operation of the electron transfer chain, supporting ATP synthesis (Kawakami et al., 2010; Yan et al., 2019). On the other hand, in mAUM the cells maintain their metabolic profile, with complex III and COX oxidases participating with 50–60% of the electron flow. However, an increase in CIO activity in the stationary phase (10% approximately) was noticeable, which could be related to an increase in cyanide production. Previous studies indicate that P. aeruginosa adapts its metabolic machinery to optimize the energy yield and survive environmental changes, especially in nutrient barren conditions (Frimmersdorf et al., 2010), such as the urine. Because of the constraints of the mAUM components, ATP yield is likely to be a major factor influencing the metabolic strategy utilized. The use of the pathway composed of complex III and COXs appears as optimal, as these enzymes pump 6 (4 and 2, respectively) protons per oxygen atom reduced (Brand, 1994), while both CYO (Graf et al., 2019) and CIO (D’mello et al., 1996) oxidases only pump four and two protons per oxygen reduced, respectively.
In contrast with the behavior of the NADH dehydrogenase mutants. The growth is significantly affected in mutants of the terminal oxidases. For instance, our data show that the CIO mutant is not affected in mAUM, consistent with its minor role in this condition, and shows a decrease in the yield in LB media, reflecting its importance in the stationary phase of growth. On the other hand, the growth rate of the of the Δbc1 mutant was significantly reduced in mAUM, consistent with the main role that we are describing above. Interestingly, this mutant also shows significant changes in the growth rate in LB media, even though it seems to have a small role compared to other oxidases in this condition. Our results also show that the mutants of cbb3 oxidases do not produce a significant decrease in the growth, indicating that the two other COX oxidases can compensate the lack of any of the other two, but also highlighting the role of complex III in energy generation. These studies were performed with the more virulent strain PA14, carrying two pathogenicity islands that are absent in PAO1, due to accessibility of mutants. However, the genomes of PAO1 and PA14 strains are very similar (Mikkelsen et al., 2011) and the genes and protein products for all the respiratory complexes mutants that we tested, are identical between the two strains. Therefore, we expect a similar growth behavior for both strains. Moreover, genome-wide genetic studies performed for P. aeruginosa PA14 and several clinical isolates grown in different conditions (including urine) reveal that the essential core genes across the isolates and growth conditions include those that encode for the respiratory enzymes described in this study. Thus, different metabolic strategies are not expected among strains (Poulsen et al., 2019).
4.4 Contribution of the respiratory complexes to Pseudomonas aeruginosa energy metabolism
In order to understand the role that each of the respiratory complexes play in the energy metabolism of P. aeruginosa, we carried out an analysis of their relative contribution to the respiratory activity and to ATP synthesis (Figure 7). To estimate the contribution to the respiratory activity we assumed that the dehydrogenases activities are rate limiting, since the terminal oxidases appear to be much more active in comparison (Figure 3). We also assumed that the physiologic concentrations of their substrates are close to the concentrations used in our experiments, and thus they all can participate with rates close to those observed individually. Finally, we assumed that substrate channeling between individual complexes, as well as microdomains, which might divert preferentially electron flow to some complexes, is negligible. With this information and considering the relative participation of NADH dehydrogenases and terminal oxidases to the respiratory rate, we were able to estimate the relative contribution of all complexes to the respiratory metabolism (Supplementary Tables S1, S2). The data obtained are presented in Figure 7A, showing that the main dehydrogenases in all conditions are NQR and SDH, carrying out about 40% of all the respiratory activity. Depending on the media, the main oxidases are CYO, CIO, or complex III plus COX oxidases. To fully understand the role of these enzymes in energy metabolism, we also estimated their roles in the generation of ATP, using the known ion pumping (H+/e-) stoichiometry for NQR [1, (Bogachev et al., 1997)], Complex I [2, (Wikström and Hummer, 2012)], Complex III and COX [4 and 2, (Brand, 1994)], CYO [4, (Graf et al., 2019)] and CIO [3, (D’mello et al., 1996)]. We would like to mention that other stoichiometries have been reported for P. aeruginosa respiratory complexes (Arai et al., 2014), but the ratios reported do not consider the contribution of NADH dehydrogenases and thus are likely overestimated. Figure 7B, shows the relative contributions of all measured respiratory enzymes to the proton pumping activity. In all conditions tested, the vast majority (>80%) of the ATP synthesis is carried out by the terminal oxidases, with NQR being the most important dehydrogenase, but playing a comparatively smaller role. These results allow us to understand the lack of effects of NQR elimination on the growth of the bacteria, as its contribution can be compensated by other enzymes, not only by the NADH dehydrogenases. The data also show that complex III and COX oxidases play significant roles in all conditions (>20%), but especially in mAUM, with a 50% participation to the overall ATP production. This effect is very well reflected in the phenotype of the Δbc1 mutant, showing important decreases in the growth rate and yield, which can be compensated, but only partially by other oxidases. This analysis finally provides a framework to understand the metabolism used by this important pathogen, which has remained unknown despite several decades of work by different groups.
Figure 7. Relative contribution of respiratory enzymes to the respiratory rate (A) and to the ion pumping rate of P. aeruginosa cells grown in different conditions (B). (A) Normalized respiratory rate of major respiratory enzymes and (B) normalized ion pumping rate of respiratory enzymes in logarithmic (LOG) and stationary (STA) phases of cells grown in LB and mAUM.
4.5 Metabolism guided antibiotic design
The results of this study provide a way to understand the role of respiratory enzymes in the metabolism of P. aeruginosa. Understanding the significance of the aerobic respiration is crucial as they play a pivotal role in maintaining bacterial growth, particularly in the context of urinary tract infections. In particular, we are describing that Complex III, NQR and SDH play important roles in cell physiology and may be targeted with new antibiotics. Complex III appears as a suitable target, as previous reports show that antiparasitic drugs targeting the microbial complex III can be obtained (Stickles et al., 2015). NQR also plays an important role, as the main NADH dehydrogenase, which could be targeted by inhibitors, as recently described for Chlamydia trachomatis (Dibrov et al., 2017a; Liang et al., 2018). However, our data shows that the elimination of the NQR gene does not produce significant effects on the growth of the bacteria. Nonetheless, it is likely that its disruption can produce a dysregulation of the metabolism and affect pathogenicity. Indeed, it has been shown that NQR deletion in V. cholerae decreases the expression of the cholera toxin (Minato et al., 2014). Finally, SDH is another potential target that could be used to design new drugs. Indeed, a recent report (Mogi et al., 2009) has shown that the inhibitor siccanin is able to specifically inhibit P. aeruginosa SDH, with no significant toxicity to human cells. Further characterizations are required to identify and validate these and other pharmacologic targets, which are urgently needed to treat CAUTIs and multidrug resistant infections by this microorganism.
Data availability statement
The Mass spectrometry data presented in the study are deposited in the MassIVE repository, accession number PXD047660. http://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=7184a7a931fe4a08b7f4bdb03c4cbc84.
Author contributions
YH: Formal analysis, Investigation, Validation, Visualization, Writing – original draft, Writing – review & editing. MY: Investigation, Validation, Visualization, Writing – original draft. AJ: Formal analysis, Investigation, Software, Writing – original draft. KT: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing. OJ: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by NIH grants NIGMS – R15GM131292 (OJ) and NIAID – R01AI151152 (OJ and KT).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2024.1347466/full#supplementary-material
Footnotes
1. ^https://www.kegg.jp/kegg/catalog/org_list.html
References
Alqarni, B., Colley, B., Klebensberger, J., McDougald, D., and Rice, S. A. (2016). Expression stability of 13 housekeeping genes during carbon starvation of Pseudomonas aeruginosa. J. Microbiol. Methods 127, 182–187. doi: 10.1016/j.mimet.2016.06.008
Andersen, M. J., and Flores-Mireles, A. L. (2020). Urinary catheter coating modifications: the race against catheter-associated infections. Coatings 10:23. doi: 10.3390/coatings10010023
Arai, H. (2011). Regulation and function of versatile aerobic and anaerobic respiratory metabolism in Pseudomonas aeruginosa. Front. Microbiol. 2:103. doi: 10.3389/fmicb.2011.00103
Arai, H., Kawakami, T., Osamura, T., Hirai, T., Sakai, Y., and Ishii, M. (2014). Enzymatic characterization and in vivo function of five terminal oxidases in Pseudomonas aeruginosa. J. Bacteriol. 196, 4206–4215. doi: 10.1128/JB.02176-14
Baradaran, R., Berrisford, J. M., Minhas, G. S., and Sazanov, L. A. (2013). Crystal structure of the entire respiratory complex I. Nature 494, 443–448. doi: 10.1038/nature11871
Berridge, M. V., Herst, P. M., and Tan, A. S. (2005). Tetrazolium dyes as tools in cell biology: new insights into their cellular reduction. Biotechnol. Annu. Rev. 11, 127–152. doi: 10.1016/S1387-2656(05)11004-7
Bodey, G. P., Bolivar, R., Fainstein, V., and Jadeja, L. (1983). Infections caused by Pseudomonas aeruginosa. Rev. Infect. Dis. 5, 279–313. doi: 10.1093/clinids/5.2.279
Bogachev, A. V., Murtazina, R. A., and Skulachev, V. P. (1997). The Na+/e- stoichiometry of the Na+−motive nadh:quinone oxidoreductase in Vibrio alginolyticus. FEBS Lett. 409, 475–477. doi: 10.1016/s0014-5793(97)00536-x
Bono, M. J., Leslie, S. W., and Reygaert, W. C. (2022). Urinary tract infection. StatPearls. StatPearls Publishing: Treasure Island, FL.
Brand, D. (1994). The stoichiometry of proton pumping and ATP synthesis in mitochondria. Biochemist 16, 20–24.
Buchfink, B., Reuter, K., and Drost, H.-G. (2021). Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 18, 366–368. doi: 10.1038/s41592-021-01101-x
Bustin, S. A., Benes, V., Garson, J. A., Hellemans, J., Huggett, J., Kubista, M., et al. (2009). The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 55, 611–622. doi: 10.1373/clinchem.2008.112797
Cassat, J. E., and Skaar, E. P. (2013). Iron in infection and immunity. Cell Host Microbe 13, 509–519. doi: 10.1016/j.chom.2013.04.010
CDC. (2019). Antibiotic resistance threats in the United States, 2019. Atlanta, GA: Centers for Disease Control and Prevention.
Chambers, S. T., Peddie, B. A., Randall, K., and Lever, M. (1999). Inhibitors of bacterial growth in urine: what is the role of betaines? Int. J. Antimicrob. Agents 11, 293–296. doi: 10.1016/S0924-8579(99)00033-3
Chu, P.-H., Huang, T.-Y., Williams, J., and Stafford, D. W. (2006). Purified vitamin K epoxide reductase alone is sufficient for conversion of vitamin K epoxide to vitamin K and vitamin K to vitamin KH2. Proc. Natl. Acad. Sci. USA 103, 19308–19313. doi: 10.1073/pnas.0609401103
Cody, W. L., Pritchett, C. L., Jones, A. K., Carterson, A. J., Jackson, D., Frisk, A., et al. (2009). Pseudomonas aeruginosa AlgR controls cyanide production in an AlgZ-dependent manner. J. Bacteriol. 191, 2993–3002. doi: 10.1128/JB.01156-08
de Abreu, P. M., Farias, P. G., Paiva, G. S., Almeida, A. M., and Morais, P. V. (2014). Persistence of microbial communities including Pseudomonas aeruginosa in a hospital environment: a potential health hazard. BMC Microbiol. 14:118. doi: 10.1186/1471-2180-14-118
Degli Esposti, M. (1998). Inhibitors of NADH–ubiquinone reductase: an overview. Biochimica et Biophysica Acta 1364, 222–235. doi: 10.1016/S0005-2728(98)00029-2
Déziel, E., Lépine, F., Milot, S., He, J., Mindrinos, M. N., Tompkins, R. G., et al. (2004). Analysis of Pseudomonas aeruginosa 4-hydroxy-2-alkylquinolines (HAQs) reveals a role for 4-hydroxy-2-heptylquinoline in cell-to-cell communication. Proc. Natl. Acad. Sci. 101, 1339–1344. doi: 10.1073/pnas.0307694100
Dibrov, P., Dibrov, E., Maddaford, T. G., Kenneth, M., Nelson, J., Resch, C., et al. (2017a). Development of a novel rationally designed antibiotic to inhibit a nontraditional bacterial target. Can. J. Physiol. Pharmacol. 95, 595–603. doi: 10.1139/cjpp-2016-0505
Dibrov, P., Dibrov, E., and Pierce, G. N. (2017b). Na+-NQR (Na+−translocating NADH:ubiquinone oxidoreductase) as a novel target for antibiotics. FEMS Microbiol. Rev. 41, 653–671. doi: 10.1093/femsre/fux032
D’mello, R., Hill, S., and Poole, R. K. (1996). The cytochrome bd quinol oxidase in Escherichia coli has an extremely high oxygen affinity and two oxygen-binding haems: implications for regulation of activity in vivo by oxygen inhibition. Microbiology 142, 755–763. doi: 10.1099/00221287-142-4-755
Fedor, J. G., Jones, A. J. Y., Di Luca, A., Kaila, V. R. I., and Hirst, J. (2017). Correlating kinetic and structural data on ubiquinone binding and reduction by respiratory complex I. Proc. Natl. Acad. Sci. USA 114, 12737–12742. doi: 10.1073/pnas.1714074114
Frimmersdorf, E., Horatzek, S., Pelnikevich, A., Wiehlmann, L., and Schomburg, D. (2010). How Pseudomonas aeruginosa adapts to various environments: a metabolomic approach. Environ. Microbiol. 12, 1734–1747. doi: 10.1111/j.1462-2920.2010.02253.x
Furt, F., Van Oostende, C., Widhalm, J. R., Dale, M. A., Wertz, J., and Basset, G. J. C. (2010). A bimodular oxidoreductase mediates the specific reduction of phylloquinone (vitamin K1) in chloroplasts. Plant J. 64, 38–46. doi: 10.1111/j.1365-313X.2010.04305.x
Giannakopoulos, X., Evangelou, A., Kalfakakou, V., Grammeniatis, E., Papandropoulos, I., and Charalambopoulos, K. (1997). Human bladder urine oxygen content: implications for urinary tract diseases. Int. Urol. Nephrol. 29, 393–401. doi: 10.1007/BF02551103
Gong, H., Li, J., Xu, A., Tang, Y., Ji, W., Gao, R., et al. (2018). An electron transfer path connects subunits of a mycobacterial respiratory supercomplex. Science 362:eaat8923. doi: 10.1126/science.aat8923
Graf, S., Brzezinski, P., and von Ballmoos, C. (2019). The proton pumping Bo oxidase from Vitreoscilla. Sci. Rep. 9:4766. doi: 10.1038/s41598-019-40723-2
Heikal, A., Nakatani, Y., Dunn, E., Weimar, M. R., Day, C. L., Baker, E. N., et al. (2014). Structure of the bacterial type II NADH dehydrogenase: a monotopic membrane protein with an essential role in energy generation. Mol. Microbiol. 91, 950–964. doi: 10.1111/mmi.12507
Iuchi, S., Chepuri, V., Fu, H. A., Gennis, R. B., and Lin, E. C. (1990). Requirement for terminal cytochromes in generation of the aerobic signal for the arc regulatory system in Escherichia coli: study utilizing deletions and lac fusions of cyo and cyd. J. Bacteriol. 172, 6020–6025. doi: 10.1128/jb.172.10.6020-6025.1990
Jacobsen, S., Stickler, D. J., Mobley, H., and Shirtliff, M. (2008). Complicated catheter-associated urinary tract infections due to Escherichia coli and Proteus mirabilis. Clin. Microbiol. Rev. 21, 26–59. doi: 10.1128/CMR.00019-07
Jo, J., Price-Whelan, A., and Dietrich, L. E. P. (2014). An aerobic exercise: defining the roles of Pseudomonas aeruginosa terminal oxidases. J. Bacteriol. 196, 4206–4215. doi: 10.1128/JB.02336-14
Juárez, O., Athearn, K., Gillespie, P., and Barquera, B. (2009). Acid residues in the transmembrane helices of the Na+−pumping NADH:quinone oxidoreductase (Na+-NQR) from Vibrio cholerae involved in sodium translocation. Biochemistry 48, 9516–9524. doi: 10.1021/bi900845y
Juárez, O., and Barquera, B. (2012). Insights into the mechanism of electron transfer and sodium translocation of the Na+−pumping NADH:quinone oxidoreductase. Biochimica et Biophysica Acta 1817, 1823–1832. doi: 10.1016/j.bbabio.2012.03.017
Juárez, O., Shea, M. E., Makhatadze, G. I., and Barquera, B. (2011). The role and specificity of the catalytic and regulatory cation-binding sites of the Na+−pumping NADH:quinone oxidoreductase from Vibrio cholerae. J. Biol. Chem. 286, 26383–26390. doi: 10.1074/jbc.M111.257873
Kahm, M., Hasenbrink, G., Lichtenberg-Fraté, H., Ludwig, J., and Kschischo, M. (2010). Grofit: fitting biological growth curves with R. J. Stat. Softw. 33, 1–21. doi: 10.18637/jss.v033.i07
Kanehisa, M., Furumichi, M., Sato, Y., Kawashima, M., and Ishiguro-Watanabe, M. (2023). KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 51, D587–D592. doi: 10.1093/nar/gkac963
Kawakami, T., Kuroki, M., Ishii, M., Igarashi, Y., and Arai, H. (2010). Differential expression of multiple terminal oxidases for aerobic respiration in Pseudomonas aeruginosa. Environ. Microbiol. 12, 1399–1412. doi: 10.1111/j.1462-2920.2009.02109.x
Kerscher, S., Dröse, S., Zickermann, V., and Brandt, U. (2008). “The three families of respiratory NADH dehydrogenases” in Bioenergetics: Energy conservation and conversion, results and problems in cell differentiation. eds. G. Schäfer and H. S. Penefsky (Berlin, Heidelberg: Springer), 185–222.
Kimelberg, H. K., and Nicholls, P. (1969). Kinetic studies on the interaction of TMPD with cytochrome c and cytochrome c oxidase. Arch. Biochem. Biophys. 133, 327–335. doi: 10.1016/0003-9861(69)90461-5
Kishikawa, J., Ishikawa, M., Masuya, T., Murai, M., Kitazumi, Y., Butler, N. L., et al. (2022). Cryo-EM structures of Na+−pumping NADH-ubiquinone oxidoreductase from Vibrio cholerae. Nat. Commun. 13:4082. doi: 10.1038/s41467-022-31718-1
Kõressaar, T., Lepamets, M., Kaplinski, L., Raime, K., Andreson, R., and Remm, M. (2018). Primer3_masker: integrating masking of template sequence with primer design software. Bioinformatics 34, 1937–1938. doi: 10.1093/bioinformatics/bty036
Kretzschmar, U., Rückert, A., Jeoung, J.-H., and Görisch, H. (2002). Malate:quinone oxidoreductase is essential for growth on ethanol or acetate in Pseudomonas aeruginosa. Microbiology (Reading) 148, 3839–3847. doi: 10.1099/00221287-148-12-3839
Laborda, P., Sanz-García, F., Hernando-Amado, S., and Martínez, J. L. (2021). Pseudomonas aeruginosa: an antibiotic resilient pathogen with environmental origin. Curr. Opin. Microbiol. 64, 125–132. doi: 10.1016/J.MIB.2021.09.010
Lenney, W., and Gilchrist, F. J. (2011). Pseudomonas aeruginosa and cyanide production. Eur. Respir. J. 37, 482–483. doi: 10.1183/09031936.00122810
Levinsky, N. G., and Berliner, R. W. (1959). Changes in composition of the urine in ureter and bladder at low urine flow. Am. J. Physiol. Legacy Content 196, 549–553. doi: 10.1152/ajplegacy.1959.196.3.549
Liang, P., Fang, X., Hu, Y., Yuan, M., Raba, D. A., Ding, J., et al. (2020). The aerobic respiratory chain of Pseudomonas aeruginosa cultured in artificial urine media: role of NQR and terminal oxidases. PLoS One 15:e0231965. doi: 10.1371/journal.pone.0231965
Liang, P., Rosas-Lemus, M., Patel, D., Fang, X., Tuz, K., and Juárez, O. (2018). Dynamic energy dependency of Chlamydia trachomatis on host cell metabolism during intracellular growth: role of sodium-based energetics in chlamydial ATP generation. J. Biol. Chem. 293, 510–522. doi: 10.1074/jbc.M117.797209
Liberati, N. T., Urbach, J. M., Miyata, S., Lee, D. G., Drenkard, E., Wu, G., et al. (2006). An ordered, nonredundant library of Pseudomonas aeruginosa strain PA14 transposon insertion mutants. Proc. Natl. Acad. Sci. USA 103, 2833–2838. doi: 10.1073/pnas.0511100103
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Livermore, D. M. (2002). Multiple mechanisms of antimicrobial resistance in Pseudomonas aeruginosa: our worst nightmare? Clin. Infect. Dis. 34, 634–640. doi: 10.1086/338782
Lunak, Z. R., and Noel, K. D. (2015). A quinol oxidase, encoded by cyoABCD, is utilized to adapt to lower O2 concentrations in Rhizobium etli CFN42. Microbiology (Reading) 161, 203–212. doi: 10.1099/mic.0.083386-0
Managò, A., Becker, K. A., Carpinteiro, A., Wilker, B., Soddemann, M., Seitz, A. P., et al. (2015). Pseudomonas aeruginosa pyocyanin induces neutrophil death via mitochondrial reactive oxygen species and mitochondrial acid sphingomyelinase. Antioxid. Redox Signal. 22, 1097–1110. doi: 10.1089/ars.2014.5979
Matsushita, K., Ohnishi, T., and Kaback, H. R. (1987). NADH-ubiquinone oxidoreductases of the Escherichia coli aerobic respiratory chain. Biochemistry 26, 7732–7737. doi: 10.1021/bi00398a029
McGill, S. L., Yung, Y., Hunt, K. A., Henson, M. A., Hanley, L., and Carlson, R. P. (2021). Pseudomonas aeruginosa reverse diauxie is a multidimensional, optimized, resource utilization strategy. Sci. Rep. 11:1457. doi: 10.1038/s41598-020-80522-8
McLennan, H. R., and Degli Esposti, M. (2000). The contribution of mitochondrial respiratory complexes to the production of reactive oxygen species. J. Bioenerg. Biomembr. 32, 153–162. doi: 10.1023/a:1005507913372
Mikkelsen, H., McMullan, R., and Filloux, A. (2011). The Pseudomonas aeruginosa reference strain PA14 displays increased virulence due to a mutation in ladS. PLoS One 6:e29113. doi: 10.1371/journal.pone.0029113
Minato, Y., Fassio, S. R., Reddekopp, R. L., and Häse, C. C. (2014). Inhibition of the sodium-translocating NADH-ubiquinone oxidoreductase [Na+-NQR] decreases cholera toxin production in Vibrio cholerae O1 at the late exponential growth phase. Microb. Pathog. 66, 36–39. doi: 10.1016/j.micpath.2013.12.002
Mittal, R., Aggarwal, S., Sharma, S., Chhibber, S., and Harjai, K. (2009). Urinary tract infections caused by Pseudomonas aeruginosa: a minireview. J. Infect. Public Health 2, 101–111. doi: 10.1016/j.jiph.2009.08.003
Mogi, T., Kawakami, T., Arai, H., Igarashi, Y., Matsushita, K., Mori, M., et al. (2009). Siccanin rediscovered as a species-selective succinate dehydrogenase inhibitor. J. Biochem. 146, 383–387. doi: 10.1093/jb/mvp085
Murray, J. L., Kwon, T., Marcotte, E. M., and Whiteley, M. (2015). Intrinsic antimicrobial resistance determinants in the superbug Pseudomonas aeruginosa. MBio 6, e01603–e01615. doi: 10.1128/mbio.01603-15
Osamura, T., Kawakami, T., Kido, R., Ishii, M., and Arai, H. (2017). Specific expression and function of the A-type cytochrome c oxidase under starvation conditions in Pseudomonas aeruginosa. PLoS One 12:e0177957. doi: 10.1371/journal.pone.0177957
Paula, A. J., Hwang, G., and Koo, H. (2020). Dynamics of bacterial population growth in biofilms resemble spatial and structural aspects of urbanization. Nat. Commun. 11:1354. doi: 10.1038/s41467-020-15165-4
Poulsen, B. E., Yang, R., Clatworthy, A. E., White, T., Osmulski, S. J., Li, L., et al. (2019). Defining the core essential genome of Pseudomonas aeruginosa. Proc. Natl. Acad. Sci. 116, 10072–10080. doi: 10.1073/pnas.1900570116
Raba, D. A., Rosas-Lemus, M., Menzer, W. M., Li, C., Fang, X., Liang, P., et al. (2018). Characterization of the Pseudomonas aeruginosa NQR complex, a bacterial proton pump with roles in autopoisoning resistance. J. Biol. Chem. 293, 15664–15677. doi: 10.1074/jbc.RA118.003194
Raba, D. A., Yuan, M., Fang, X., Menzer, W. M., Xie, B., Liang, P., et al. (2019). Role of subunit D in ubiquinone-binding site of Vibrio cholerae NQR: pocket flexibility and inhibitor resistance. ACS Omega 4, 19324–19331. doi: 10.1021/acsomega.9b02707
Reyes-Prieto, A., Barquera, B., and Juárez, O. (2014). Origin and evolution of the sodium -pumping NADH: ubiquinone oxidoreductase. PLoS One 9:e96696. doi: 10.1371/journal.pone.0096696
Reynolds, D., and Kollef, M. (2021). The epidemiology and pathogenesis and treatment of Pseudomonas aeruginosa infections: an update. Drugs 81, 2117–2131. doi: 10.1007/s40265-021-01635-6
Sarigul, N., Korkmaz, F., and Kurultak, İ. (2019). A new artificial urine protocol to better imitate human urine. Sci. Rep. 9:20159. doi: 10.1038/s41598-019-56693-4
Scheijen, J. L. J. M., Hanssen, N. M. J., van de Waarenburg, M. P. H., Jonkers, D. M. A. E., Stehouwer, C. D. A., and Schalkwijk, C. G. (2012). L(+) and D(−) lactate are increased in plasma and urine samples of type 2 diabetes as measured by a simultaneous quantification of L(+) and D(−) lactate by reversed-phase liquid chromatography tandem mass spectrometry. J. Diabetes Res. 2012:e234812, 1–10. doi: 10.1155/2012/234812
Selezska, K., Kazmierczak, M., Müsken, M., Garbe, J., Schobert, M., Häussler, S., et al. (2012). Pseudomonas aeruginosa population structure revisited under environmental focus: impact of water quality and phage pressure. Environ. Microbiol. 14, 1952–1967. doi: 10.1111/j.1462-2920.2012.02719.x
Sievert, D. M., Ricks, P., Edwards, J. R., Schneider, A., Patel, J., Srinivasan, A., et al. (2013). Antimicrobial-resistant pathogens associated with healthcare-associated infections: summary of data reported to the National Healthcare Safety Network at the Centers for Disease Control and Prevention, 2009–2010. Infect. Control Hosp. Epidemiol. 34, 1–14. doi: 10.1086/668770
Stickles, A. M., Ting, L.-M., Morrisey, J. M., Li, Y., Mather, M. W., Meermeier, E., et al. (2015). Inhibition of cytochrome bc1 as a strategy for single-dose, multi-stage antimalarial therapy. Am. J. Trop. Med. Hyg. 92, 1195–1201. doi: 10.4269/ajtmh.14-0553
Stover, C. K., Pham, X. Q., Erwin, A. L., Mizoguchi, S. D., Warrener, P., Hickey, M. J., et al. (2000). Complete genome sequence of Pseudomonas aeruginosa PAO1, an opportunistic pathogen. Nature 406, 959–964. doi: 10.1038/35023079
Torres, A., Kasturiarachi, N., DuPont, M., Cooper, V. S., Bomberger, J., and Zemke, A. (2019). NADH dehydrogenases in Pseudomonas aeruginosa growth and virulence. Front. Microbiol. 10:75. doi: 10.3389/fmicb.2019.00075
Trautner, B. W., and Darouiche, R. O. (2004). Role of biofilm in catheter-associated urinary tract infection. Am. J. Infect. Control 32, 177–183. doi: 10.1016/j.ajic.2003.08.005
Ugidos, A., Morales, G., Rial, E., Williams, H. D., and Rojo, F. (2008). The coordinate regulation of multiple terminal oxidases by the Pseudomonas putida ANR global regulator. Environ. Microbiol. 10, 1690–1702. doi: 10.1111/j.1462-2920.2008.01586.x
Wang, Y., Xiao, D., Liu, Q., Zhang, Y., Hu, C., Sun, J., et al. (2018). Two NAD-independent l-lactate dehydrogenases drive l-lactate utilization in Pseudomonas aeruginosa PAO1. Environ. Microbiol. Rep. 10, 569–575. doi: 10.1111/1758-2229.12666
Weiner, L. M., Webb, A. K., Limbago, B., Dudeck, M. A., Patel, J., Kallen, A. J., et al. (2016). Antimicrobial-resistant pathogens associated with healthcare-associated infections: summary of data reported to the National Healthcare Safety Network at the Centers for Disease Control and Prevention, 2011–2014. Infect. Control Hosp. Epidemiol. 37, 1288–1301. doi: 10.1017/ice.2016.174
Weiner-Lastinger, L. M., Abner, S., Edwards, J. R., Kallen, A. J., Karlsson, M., Magill, S. S., et al. (2020). Antimicrobial-resistant pathogens associated with adult healthcare-associated infections: summary of data reported to the National Healthcare Safety Network, 2015–2017. Infect. Control Hosp. Epidemiol. 41, 1–18. doi: 10.1017/ice.2019.296
WHO (2017). WHO publishes list of bacteria for which new antibiotics are urgently needed [WWW Document]. Available at: https://www.who.int/news/item/27-02-2017-who-publishes-list-of-bacteria-for-which-new-antibiotics-are-urgently-needed (Accessed January 21, 2022).
Wiehlmann, L., Wagner, G., Cramer, N., Siebert, B., Gudowius, P., Morales, G., et al. (2007). Population structure of Pseudomonas aeruginosa. Proc. Natl. Acad. Sci. USA 104, 8101–8106. doi: 10.1073/pnas.0609213104
Wikström, M., and Hummer, G. (2012). Stoichiometry of proton translocation by respiratory complex I and its mechanistic implications. Proc. Natl. Acad. Sci. 109, 4431–4436. doi: 10.1073/pnas.1120949109
Williams, H. D., Zlosnik, J. E. A., and Ryall, B. (2006). Oxygen, cyanide and energy generation in the cystic fibrosis pathogen Pseudomonas aeruginosa. Adv. Microb. Physiol. 52, 1–71. doi: 10.1016/S0065-2911(06)52001-6
Winsor, G. L., Griffiths, E. J., Lo, R., Dhillon, B. K., Shay, J. A., and Brinkman, F. S. L. (2016). Enhanced annotations and features for comparing thousands of Pseudomonas genomes in the Pseudomonas genome database. Nucleic Acids Res. 44, D646–D653. doi: 10.1093/nar/gkv1227
Wittig, I., Braun, H.-P., and Schägger, H. (2006). Blue native PAGE. Nat. Protoc. 1, 418–428. doi: 10.1038/nprot.2006.62
Yagi, T. (1991). Bacterial NADH-quinone oxidoreductases. J. Bioenerg. Biomembr. 23, 211–225. doi: 10.1007/BF00762218
Yagi, T., Yano, T., Di Bernardo, S., and Matsuno-Yagi, A. (1998). Procaryotic complex I (NDH-1), an overview. Biochimica et Biophysica Acta Bioenergetics 1364, 125–133. doi: 10.1016/S0005-2728(98)00023-1
Yan, H., Asfahl, K. L., Li, N., Sun, F., Xiao, J., Shen, D., et al. (2019). Conditional quorum-sensing induction of a cyanide-insensitive terminal oxidase stabilizes cooperating populations of Pseudomonas aeruginosa. Nat. Commun. 10:4999. doi: 10.1038/s41467-019-13013-8
Zambrano, M. M., and Kolter, R. (1993). Escherichia coli mutants lacking NADH dehydrogenase I have a competitive disadvantage in stationary phase. J. Bacteriol. 175, 5642–5647. doi: 10.1128/jb.175.17.5642-5647.1993
Zhou, W., Bertsova, Y. V., Feng, B., Tsatsos, P., Verkhovskaya, M. L., Gennis, R. B., et al. (1999). Sequencing and preliminary characterization of the Na+−translocating NADH:ubiquinone oxidoreductase from Vibrio harveyi. Biochemistry 38, 16246–16252. doi: 10.1021/bi991664s
Keywords: Pseudomonas aeruginosa, aerobic metabolism, respiratory chain, NQR, bc1 complex, SDH, urinary tract infection, metabolic activation
Citation: Hu Y, Yuan M, Julian A, Tuz K and Juárez O (2024) Identification of complex III, NQR, and SDH as primary bioenergetic enzymes during the stationary phase of Pseudomonas aeruginosa cultured in urine-like conditions. Front. Microbiol. 15:1347466. doi: 10.3389/fmicb.2024.1347466
Edited by:
Davide Zannoni, University of Bologna, ItalyReviewed by:
Martine Pestel-Caron, Université de Rouen, FranceRaymond J. Turner, University of Calgary, Canada
Viviana Teresa Orlandi, University of Insubria, Italy
Copyright © 2024 Hu, Yuan, Julian, Tuz and Juárez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Oscar Juárez, ojuarez@iit.edu