Next Article in Journal
Antimicrobial Evaluation of Two Polycyclic Polyprenylated Acylphloroglucinol Compounds: PPAP23 and PPAP53
Previous Article in Journal
Microbial and Metabolic Gut Profiling across Seven Malignancies Identifies Fecal Faecalibacillus intestinalis and Formic Acid as Commonly Altered in Cancer Patients
Previous Article in Special Issue
Integrative and Conjugative Elements and Prophage DNA as Carriers of Resistance Genes in Erysipelothrix rhusiopathiae Strains from Domestic Geese in Poland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

From Proteome to Potential Drugs: Integration of Subtractive Proteomics and Ensemble Docking for Drug Repurposing against Pseudomonas aeruginosa RND Superfamily Proteins

by
Gabriela Urra
1,
Elizabeth Valdés-Muñoz
2,
Reynier Suardiaz
3,
Erix W. Hernández-Rodríguez
1,4,
Jonathan M. Palma
5,
Sofía E. Ríos-Rozas
1,
Camila A. Flores-Morales
6,
Melissa Alegría-Arcos
7,
Osvaldo Yáñez
7,
Luis Morales-Quintana
8,
Vívian D’Afonseca
9,* and
Daniel Bustos
1,*
1
Laboratorio de Bioinformática y Química Computacional, Departamento de Medicina Traslacional, Facultad de Medicina, Universidad Católica del Maule, Talca 3480094, Chile
2
Doctorado en Biotecnología Traslacional, Facultad de Ciencias Agrarias y Forestales, Universidad Católica del Maule, Talca 3480094, Chile
3
Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
4
Unidad de Bioinformática Clínica, Centro Oncológico, Facultad de Medicina, Universidad Católica del Maule, Talca 3480094, Chile
5
Facultad de Ingeniería, Universidad de Talca, Curicó 3344158, Chile
6
Magíster en Ciencias de la Computación, Universidad Católica del Maule, Talca 3460000, Chile
7
Núcleo de Investigación en Data Science, Facultad de Ingeniería y Negocios, Universidad de las Américas, Santiago 7500000, Chile
8
Multidisciplinary Agroindustry Research Laboratory, Instituto de Ciencias Biomédicas, Facultad de Ciencias de la Salud, Universidad Autónoma de Chile, Cinco Pte. N° 1670, Talca 3467987, Chile
9
Departamento de Ciencias Preclínicas, Facultad de Medicina, Universidad Católica del Maule, Ave. San Miguel 3605, Talca 3466706, Chile
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2024, 25(15), 8027; https://doi.org/10.3390/ijms25158027
Submission received: 27 May 2024 / Revised: 8 July 2024 / Accepted: 17 July 2024 / Published: 23 July 2024
(This article belongs to the Special Issue Antibiotic Resistance: Appearance, Evolution, and Spread 2.0)

Abstract

:
Pseudomonas aeruginosa (P. aeruginosa) poses a significant threat as a nosocomial pathogen due to its robust resistance mechanisms and virulence factors. This study integrates subtractive proteomics and ensemble docking to identify and characterize essential proteins in P. aeruginosa, aiming to discover therapeutic targets and repurpose commercial existing drugs. Using subtractive proteomics, we refined the dataset to discard redundant proteins and minimize potential cross-interactions with human proteins and the microbiome proteins. We identified 12 key proteins, including a histidine kinase and members of the RND efflux pump family, known for their roles in antibiotic resistance, virulence, and antigenicity. Predictive modeling of the three-dimensional structures of these RND proteins and subsequent molecular ensemble-docking simulations led to the identification of MK-3207, R-428, and Suramin as promising inhibitor candidates. These compounds demonstrated high binding affinities and effective inhibition across multiple metrics. Further refinement using non-covalent interaction index methods provided deeper insights into the electronic effects in protein–ligand interactions, with Suramin exhibiting superior binding energies, suggesting its broad-spectrum inhibitory potential. Our findings confirm the critical role of RND efflux pumps in antibiotic resistance and suggest that MK-3207, R-428, and Suramin could be effectively repurposed to target these proteins. This approach highlights the potential of drug repurposing as a viable strategy to combat P. aeruginosa infections.

Graphical Abstract

1. Introduction

Each year, around 2 million individuals in the United States suffer infections caused by antibiotic-resistant bacteria, leading to approximately 23,000 deaths due to multi-drug resistant (MDR) infections. At the same time, Europe contends with a similar scenario of roughly 25,000 deaths annually [1]. The economic impact of these infections surpasses 2 billion dollars annually, affecting the healthcare system and citizen’s finances [2]. In 2024, the World Health Organization (WHO) designated seven bacteria as high-priority pathogens due to their capacity to induce deadly contagions. Among these are Salmonella Typhi, Shigella spp., non-typhoidal Salmonella, and Neisseria gonorrhoeae, all of them resistant to fluoroquinolone. Enterococcus faecium vancomycin-resistant, Staphylococcus aureus methicillin-resistant, and Pseudomonas aeruginosa carbapenem-resistant [3].
P. aeruginosa is an opportunistic gram-negative bacterium with MDR and is the cause of acute and chronic infections in immunocompromised individuals suffering from certain pathologies such as cancer, cystic fibrosis, trauma, chronic obstructive pulmonary disease, sepsis, burns, ventilator-associated pneumonia, and is also the main cause of nosocomial disease [4,5,6]. Concerning healthcare-associated P. aeruginosa infections, it is estimated that account for 16.2% of patient infections and contribute to 23% of all infections acquired in the Intensive Care Unit (ICU) [7].
This bacterium can adapt to diverse and adverse host environments through the secretion of virulence factors; moreover, these virulence factors are involved in infection, colonization, and antibiotic resistance [8]. Among the virulence factors, a finite number of polar pili expressed by P. aeruginosa are involved in host colonization. It also presents polar flagella which determines the motility of the bacterium [9,10], both flagella and pili are responsible for antibiotic resistance due to the formation of biofilms [11]. On the other hand, the resistance-nodulation-division (RND) efflux pumps play an important role in obtaining and expressing the phenotype that gives MDR to P. aeruginosa [12]. The function of these pumps is based on a tripartite complex formed by the RND protein, the membrane fusion protein (MFP) and the outer membrane protein (OMP), forming a channel through the cell envelope that transports toxic compounds out of the cell [13].
The early stage of drug discovery plays a pivotal role in identifying potential therapeutic targets and compounds, laying the groundwork for the development of effective treatments. Within the context of combating antibiotic resistance, this phase gains significance as it facilitates the identification of novel strategies to counteract the evolving mechanisms of resistance in bacteria such as P. aeruginosa. Two synergistic computational techniques, subtractive proteomics [14] and drug repurposing [15], emerge as promising approaches in this endeavor. Subtractive proteomics enables the systematic identification and prioritization of essential proteins within the microorganism proteome that are crucial for its survival and pathogenicity. By subtracting non-essential proteins, this technique streamlines the identification of potential drug targets, focusing efforts on those most likely to yield therapeutic benefits. Concurrently, drug repurposing employs existing pharmaceutical compounds for new therapeutic purposes. This strategy capitalizes on the knowledge of approved drugs, expediting the drug discovery process and potentially avoiding the lengthy and costly traditional drug development pipeline. By repurposing drugs to target specific proteins implicated in antibiotic resistance, this approach offers a pragmatic solution to address the urgent need for novel antimicrobial agents.

2. Results and Discussion

The study describes a pipeline that integrates various computational techniques to identify essential and druggable proteins within the P. aeruginosa proteome (PAO1—ATCC 15692). Specifically, these targets belonging to the RND efflux pump were three-dimensionally modeled and energy equilibrated through Molecular Dynamics (MD) simulations for repositioning 5709 approved drugs, using ensemble docking and promolecular densities to rank their protein-drug interactions. Every stage of the sequential protocol is depicted in Figure 1.

2.1. Subtractive Analysis from the P. aeruginosa Proteome

The whole proteome of P. aeruginosa was downloaded from the NCBI database, comprising 5564 protein sequences (Table 1). To streamline the dataset and ensure the inclusion of only one copy of highly similar proteins, 111 redundant sequences were eliminated. While the significance of small proteins (sequences with a length of up to 100 amino acids) remains a debated topic due to their potential roles in various cellular processes [16], we opted to remove them. The decision was made because the functionality of most small proteins remains largely unknown, mainly due to state of art of current experimental techniques and bioinformatics tools [17]. As a result, 333 small proteins with fewer than 100 amino acids were removed from the dataset.

2.2. Subtractive Analysis Removing Human and Human Microbiome Orthologs and Essentiality Analysis

In our pursuit of identifying pharmacologically relevant targets with high affinity for approved drugs, it is crucial to eliminate bacterial targets with significant homology to human proteins. This step is essential for mitigating the risk of unintended side effects resulting from the modulation of these human proteins by drugs intended for bacterial targets. By employing stringent criteria to identify and exclude homologous sequences between the proteomes of P. aeruginosa and humans, we ensure the specificity and safety of subsequent analyses. The remaining 4344 sequences were individually compared to the 79,684 proteins belonging to the human proteome by using BLAST+. There exist 776 host proteins nearly related to the pathogen, which were removed from the subsequent analysis (Table 1).
To discover novel protein inhibitors for P. aeruginosa, it is imperative to underscore the pivotal role of the human microbiome, as evidenced by recent scientific investigations [18,19,20,21]. These studies elucidate the intricate bidirectional interactions between the gut microbiota and non-antibiotic drugs, revealing that the microbiome not only influences the pharmacokinetics and pharmacodynamics of drugs but is also modulated by them. Furthermore, the microbiome’s capacity to enzymatically transform drug structures underscores its potential to influence drug bioavailability, bioactivity, and toxicity, thereby shaping individual responses to pharmacotherapy. In this sense, we compare the sequences of 299,611 proteins belonging to 73 bacteria from the gut microbiome to remove 984 of them with high homology to the P. aeruginosa proteome (Table 1). This ensures that any potential P. aeruginosa inhibitors identified in our pipeline do not induce inadequate pharmacokinetic conditions, such as poor absorption, distribution, metabolization, or excretion of the drugs, while also addressing pharmacodynamic concerns, including high toxicity in humans, to safeguard patient health.
The last step of this first filter involves comparing our dataset against the DEG database. DEG contains the essential genes for more than 66 bacteria, as well as other domains of life (Archaea and Eukaryotes). Essential genes represent fundamental components vital for the survival and sustenance of an organism, thus constituting a cornerstone of its existence. In bacteria, essential genes represent the minimal yet indispensable components of the genome, forming functional modules pivotal for key biological functions. Consequently, these genes and their protein products serve as prime candidates when seeking relevant pharmacological targets. Among the dataset comprising 3360 proteins obtained previously, 96.55% of these proteins have been categorized as non-essential, while only 116 proteins fulfill the critical subcellular requirements necessary for the survival of P. aeruginosa (Table 1).

2.3. Assessing the Attributes of Virulence, Antigenicity, Bacterial Resistance, and Druggability

In the P. aeruginosa pathogenicity context, a multitude of virulence factors (such as proteases, lipopolysaccharides, and outer membrane proteins among others) have been documented [22] to contribute to various pathogenic mechanisms underlying nosocomial infections in immunocompromised individuals. These factors intricately orchestrate diverse strategies by which the bacterium inflicts harm upon the host’s defenses. Another relevant characteristic of the pathogenicity of P. aeruginosa is that this bacterium exhibits a formidable array of intrinsic and acquired resistance mechanisms, contributing to its ability to withstand classical antibiotic treatments. This underscores the urgent imperative to explore novel strategies for combating these pathogenic infections [23]. Furthermore, the identification of antigenic-related proteins holds significant importance in the exploration of novel therapeutic avenues for treating infections caused by P. aeruginosa. In the past, several P. aeruginosa antigens have been a matter of vaccine and immunotherapy development [24].
The evaluation of the 116 previously identified proteins against the specified criteria revealed intriguing findings. Among them, 69 proteins demonstrated characteristics indicative of virulence factors, highlighting their potential role in the pathogenicity of P. aeruginosa. Additionally, 45 proteins exhibited features associated with conferring resistance, underscoring their contribution to the bacterium’s ability to withstand antibiotic treatments. Furthermore, 54 proteins displayed antigenic properties, suggesting their potential as targets for immune recognition and intervention. However, the intersection of these criteria revealed a select group of proteins that possess a particularly noteworthy profile. Remarkably, only 23 proteins met all three criteria simultaneously (Figure 2).
This subset of 23 proteins emerges as a compelling avenue for in-depth investigation and therapeutic targeting. Upon examination of their druggability within the DrugBank database, it was revealed that only 12 of these proteins are classified as pharmacological targets. Remarkably, each of these 12 proteins has either an FDA-approved drug or a candidate undergoing investigational phases in the drug development pipeline.

2.4. Characterization of Essential Proteins of P. aeruginosa

One of the proteins identified from the previous selection is a histidine kinase (HK; Table 2), which plays a critical role in coordinating resistance mechanisms and pathogenicity through its role in two-component systems (TCSs). These systems allow the bacterium to detect and respond to environmental cues essential for its cellular development, survival, and virulence [25]. HKs serve as integral components in sensing and transmitting signals associated with antibiotic exposure, enabling the bacterium to initiate downstream responses aimed at evading or neutralizing the antibiotic threat [26]. Specifically, this enzyme belongs to the transferase type and is situated in the cytoplasmic membrane (Figure S1A). It responds to stimuli originating from the periplasm by undergoing autophosphorylation of a histidine residue, which subsequently transmits a signal to an ATP-binding kinase domain.
Table 2 provides the functional characterization of these proteins using the KEGG database and the subcellular localization data from a consensus of two servers (PSORT and CELLO). Remarkably, excluding HK, the remaining 11 proteins belong to an efflux pump family comprising 3 different proteins: an outer porin-like protein called Outer Membrane Factor (OMF, in light green in Figure S1B) that spans the outer membrane in gram-negative bacteria, an inner membrane protein receiving the family name Resistance/Nodulation/cell Division (RND, in gray in Figure S1B), and a periplasmic protein named Membrane Fusion Protein (MFP, in dark green in Figure S1B) that connects the two previous proteins. These antiporter multidrug efflux pumps belong to the RND family [27,28]. The RND efflux pump systems are universally found across almost all biological kingdoms, utilizing common mechanisms that greatly enhance the diversity of resistance profiles [29].
The proposed functional mechanisms entail a coordinated rotation of the pump components, transitioning through three consecutive states: site access, substrate binding, and extrusion [29]. It has been reported that dynamic phenomena involved in the substrate binding and the subsequent extrusion differ from those observed in an inhibitor-bound state, where rotational movement is blocked, thus impeding the extrusion [30]. The main proteins forming the MFP group include MexA, MexX, MexC, and MexE; for the RND type, we found: MexB, MexY, MexD, and MexF; and OMF proteins comprise OprM, OprJ, and OprN proteins [31]. The MexAB-OprM, MexXY-OprM, MexCD-OprJ, and MexEF-OprN pumps, which belong to the RND superfamily, are linked to resistance against most antipseudomonal medications when overexpressed [32]. The MexCD-OprJ operon in P. aeruginosa is usually expressed at low levels, but it plays a crucial role in expelling various antibiotics and toxic substances, including quinolones, chloramphenicol, and tetracyclines [33]. The MexEF-OprN system typically remains inactive under standard conditions. However, its overexpression has been associated with resistance to chloramphenicol, quinolones, trimethoprim, and certain β-lactams [34]. Both the MexCD-OprJ and MexEF-OprN efflux systems play a role in quorum sensing (QS) by expelling signal molecules such as 4-hydroxy-2-heptylquinoline and kynurenine. This reduces the accumulation of QS signals within the cell, thereby influencing the regulation of QS-controlled virulence factors [35]. The MexXY efflux pump is one of the primary resistance mechanisms in P. aeruginosa but lacks a coding sequence for an OMF. However, it functions as a multidrug efflux pump in conjunction with OprM from the MexAB-OprM operon [36]. Like other RND pumps, MexXY is essential for resistance to aminoglycosides, fluoroquinolones, and certain β-lactams [37]. The MexAB-OprM operon, recognized as the first multidrug efflux pump in P. aeruginosa, is the major contributor to antibiotic resistance [38]. Increased expression of the MexAB-OprM efflux pump is linked to resistance against almost 50 antibiotics [29].
The identity percentages for OMF proteins (Figure S2A) and MFP proteins (Figure S2B) were determined for both canonical proteins and the 11 proteins studied. Among these, 7 proteins (Q9HXB9, Q9HU26, Q9I0V8, Q9HY88, Q9I006, Q9HWH3, and Q9I0Y7) were classified as OMF proteins, while the other 4 proteins (Q9I0V5, Q9I0Y9, Q9I3R2, and G3XD25) were categorized as the MFP type. In terms of OMF proteins, only two exhibit more than 35% identity with canonical proteins: Q9HWH3, with 37.10% identity to OprN, and Q9I0Y7, with 96.61% identity to OprN protein. Meanwhile, among the MFP proteins, Q9I0Y9 demonstrates 99.03% identity to MexE, and G3XD25 displays 71.00% identity to MexC. These identity percentages provide insights into the structural similarities between the studied proteins and their canonical counterparts, aiding in the elucidation of their functional roles and evolutionary relationships.

2.5. Predicting and Refining the Three-Dimensional Structure and Multimeric Assembly in RND Proteins

In the subsequent filtering step, the HK protein is discarded to narrow down the focus on conducting drug repositioning for the 11 proteins of the RND efflux pump. Upon conducting pairwise alignment of our sequences against the RCSB-Protein Data Bank (PDB) database [39], it was observed that only the Q9I0Y7 protein of the OMF type has been resolved through experimental methods. Specifically, it was resolved using the X-ray diffraction method, yielding a resolution of 1.69 Å. The corresponding PDB code for this resolved structure is 5ZAP [40] (Table S2). We used the AlphaFold server [41] to predict their 3D structures for the remaining proteins. AlphaFold is an AI system developed by DeepMind, capable of accurately predicting the 3D structures of proteins. It achieves this through the application of deep learning algorithms trained on extensive datasets comprising protein sequences and their corresponding 3D structures. However, its current version typically predicts monomeric forms of proteins and not the multimeric complexes that RND proteins form in their functional states. Table S2 in our study presents the sequences of the proteins under investigation, along with the corresponding PDB code and the percentage of identity of their closest homologs with resolved 3D structures and assembled multimeric units. This information served as input for constructing the multimers of each protein through structural alignments and replication of each subunit. Specifically, OMF type proteins are composed of three identical subunits, forming homotrimers. On the other hand, MFP proteins typically form homohexamers, as depicted in Figure 3. To mitigate steric impediments arising from multimerization, we conducted energy minimization and 50 ns NPT MD simulations for each protein. This process involved optimizing the atomic positions of the multimeric structures to relieve any potential clashes or overlaps and achieve more stable conformations.

2.6. Rescuing Relevant Protein Conformations and Repositioning Drugs through Ensemble-docking and Non-Covalent Interactions

Ensemble-docking enhances the accuracy of virtual screening by considering multiple conformational states of the target protein, reflecting its dynamic nature and improving the identification of effective drug candidates [42]. In this study, we employed ensemble-docking combined with various docking energy-based metrics. Utilizing multiple metrics allows for a more robust selection process, as it mitigates the biases and limitations associated with relying on a single metric [43,44]. The three components of the OprM-MexAB pump, derived from the 6TA6 crystallographic structure, were employed as protein controls. Furthermore, a set of 8 known RND pump inhibitors served as ligand controls: P9D (also called ABI-PP or D13-9001) [45,46], cefiderocol (CEF) [29], doxorubicin (DM2) [46,47,48], minocycline (MIY) [45,46,48], oxacillin [29,46,49], MBX [45], rhodamine-6G [29,46,47], and erythromycin [29,45,46]. The components were simulated separately for 50 ns, and a root mean square deviation (RMSD, Figure S3) dendrogram was calculated for each, resulting in 5 clusters per component. Subsequently, a representative frame from each cluster was selected, and 5712 FDA-approved drugs were repositioned, along with the ligand controls. Given that most RND pump inhibitors have been reported to interact primarily with the internal component of the pump (MexB). We docked the ligand controls into the MexB binding site. Then, we computed the median binding energy for all conformations of MexB for each ligand control to identify the most affine inhibitor. This inhibitor was used as a benchmark in subsequent docking evaluations of the other proteins under study. Among the 8 inhibitors tested, P9D, CEF, and DM2 showed the strongest binding affinities to MexB, as depicted in Figure S4. Notably, P9D in its two tautomeric forms (P9D_1 and P9D_2) had the lowest median energy.
Remarkably, the binding energy trend for P9D, CEF, and DM2 persisted with the OprM and MexA controls (Figure S4), albeit with binding energies approximately 2 kcal/mol higher than those observed with MexB. Additionally, among the three components of the RND pump, the internal protein (MexB) consistently demonstrated superior binding energies with all the repositioned compounds (FDA dataset). Based on these findings, we selected P9D as the control ligand for the remaining evaluations, using its median binding energy in each protein as the cutoff criterion.
In our ensemble-docking protocol for drug repositioning with the 11 proteins under study, we observed consistent trends for both OMF and MFP type proteins. For OMF proteins, the repositioned compounds exhibited affinity energies comparable to those of the control protein (OprM), with a median energy around −6.0 kcal/mol (Figure 4A). Moreover, the data dispersion was consistent across all proteins. The median energy of the most effective known inhibitor, used as a cutoff, was approximately −7.0 kcal/mol for both OprM and the other seven OMF group proteins. Notably, for each protein, we identified over 100,000 repositioned drugs with binding affinity energies surpassing the median affinity energy of the control inhibitor (Figure 4B). Similarly, for MFP type proteins (Figure 4C), we observed a trend parallel to that of the control protein (MexA), with median energies around −6.0 kcal/mol and comparable data dispersions across all cases. Additionally, we identified over 100,000 repositioned drugs with binding affinity energies more favorable than the median energy of the control inhibitor (around −7.0 kcal/mol, Figure 4D).
The total energy population of our study is calculated by multiplying the number of proteins (p), frames (f), docking grids (g), docking poses (dp), and repositioned drugs (c). This calculation yields 5,997,600 docking values for the OMF proteins (7 × 5 × 3 × 10 × 5712) and 4,569,600 docking values for the MFP type proteins (4 × 5 × 4 × 10 × 5712), resulting in more than 10 million molecular docking evaluations in total. When selecting the best inhibitor candidates, we applied multiple criteria to the energetic populations of docked molecules [50]. For each repositioned drug, we considered all docking poses obtained from each protein conformation and each docking run (with varying grid positions). We calculated several metrics: the top-scoring energy (the most negative value), the median energy, and the arithmetic, geometric, and harmonic means. Each drug was then ranked based on these individual metrics. To determine the overall ranking, we calculated the mode of these rankings. For instance, if a compound consistently ranks first across several of the five metrics, its mode would be 1, indicating it as a top candidate. Additionally, we assessed the frequency of these rankings; if a compound ranks first in three out of the five metrics, its frequency is 3/5. The goal was to select candidates that ranked highest in at least 3 out of the 5 metrics and were present in at least 60% of the proteins evaluated (6 out of 11). This rigorous approach allowed us to robustly identify the most promising inhibitors.
We identified three compounds as the top candidates across the majority of the evaluated proteins:
  • MK-3207: This molecule is currently undergoing phase 2 clinical trials and is associated with migraine disorders [51]. It acts as an antagonist of the Calcitonin gene-related peptide type 1 receptor in humans. MK-3207 ranks first in the rankings for all MFP type proteins and in 6 out of 7 OMF proteins.
  • R-428 (Bemcentinib): This compound, in clinical phase 2, is being investigated for various pathologies, including myelodysplastic syndrome, melanoma, acute myeloid leukemia, and mesothelioma. It acts by inhibiting the Tyrosine-protein kinase receptor [52]. R-428 is the top-ranked compound across all metrics for 100% of the proteins evaluated.
  • Suramin: Currently advancing to phase 3 of clinical development, Suramin is associated with the treatment of non-small cell lung carcinoma, prostate adenocarcinoma, autism spectrum disorder, and acute kidney injury. This compound functions as an acidic fibroblast growth factor inhibitor [53,54]. Suramin appears as an inhibitor in all MFP proteins and 4 out of 7 OMF type proteins.
The final stage of our filtering process involved evaluating the three best inhibitor candidates for each protein studied using the non-covalent interaction index method (Figure 5 and Table S3). This method provides a more accurate quantification of all electronic effects involved in the protein–ligand interaction compared to traditional docking methods. Overall, we observed that the compound Suramin consistently exhibits better attractive energies across all evaluated proteins compared to MK-3207 and R-428.

3. Materials and Methods

3.1. Subtractive Proteomics

The proteomic datasets for the P. aeruginosa PAO1 reference strain (ATCC: 15692; UniProt ID: UP000002438) with 5564 predicted proteins and Homo sapiens (UniProt ID: UP000005640) with 79,684 proteins were retrieved from the UniProt repository [55] (www.uniprot.org). To enhance the quality of the datasets, we implemented an additional step in our pipeline aimed at eliminating redundant information, specifically duplicated proteins. To achieve this, we employed CD-Hit software v.4.8.1 [56] (https://sites.google.com/view/cd-hit, accessed on 30 December 2023) and OrthoMCL software v.2.0 [57] (https://orthomcl.org/orthomcl/app, accessed on 30 December 2023). These tools were applied to the P. aeruginosa proteome, ensuring the removal of paralogous sequences by retaining only one copy of each protein. The criteria for eliminating homologous proteins included an amino acid identity greater than 60% and an e-value exceeding 10−5. Additionally, protein sequences shorter than 100 amino acids were excluded from the final dataset. This refined approach enhanced the accuracy and reliability of the proteomic dataset generated in this work, providing a concise protein subset for subsequent analyses.

3.2. Search for Essential Protein Sequences of P. aeruginosa with Low Similarity to the Human Proteome and Human Bacterial Microbiome

Since the objective of our study is to repurpose drugs for proteins identified through the subtractive proteomics process, to ensure the safety and efficacy of these repurposed drugs, it is crucial to avoid potential side effects that may arise from the drugs binding to homologous human proteins. Therefore, we excluded proteins with high homology to the human proteome and the proteomes in the human microbiota. The screening process used BLAST+ protein software v.2.15.1 [58] (https://blast.ncbi.nlm.nih.gov/, accessed on 10 January 2024) installed locally. This software was employed to identify homologous sequences between the proteomes from P. aeruginosa and humans, subsequently excluding them from the dataset. To ensure the specificity of the analysis, a cutoff criterion was set, requiring a bit score value greater than 100 and an e-value exceeding 10−5.
The human gut, known for harboring a diverse community of approximately 1000 microorganisms, plays a crucial role in maintaining health. In consideration of the microbial landscape, a strategic filtering process was implemented to prevent potential cross-interactions between molecular targets of P. aeruginosa and the vast array of microorganisms constituting the human bacterial microbiome. To execute this, we utilized the BLAST+ protein software. The objective was to filter out homologous sequences from the P. aeruginosa proteome by comparing them against the extensive repertoire of the human bacterial microbiome. Our reference dataset consisted of 299,611 proteins curated from 73 bacteria commonly found in the gut (Table S1), as reported by Shanmugham and Pan [59]. It eliminated homologous sequences between the P. aeruginosa proteome and the human microbiome available, employing stringent cutoff criteria—a bit score value exceeding 100 and an e-value greater than 10−5. In the final step of the first filter of our protocol, we searched for essential proteins in P. aeruginosa through analysis using the DEG Database [60] (Database of Essential Genes, http://origin.tubic.org/deg/public/index.php/search/bacteria, accessed on 10 January 2024). It used P. aeruginosa PAO1 as the reference genome (RefSeq: NC_002516), which encompasses 117 essential genes; the analysis involved a comparison with our protein dataset obtained from the preceding steps. A cutoff was set at a bit score value greater than 100 and an e-value exceeding 10−10.

3.3. Characterizing Pharmacological Targets in P. aeruginosa: Integrating Virulence, Resistance, and Antigenicity Analyses

To identify suitable pharmacological targets among the 116 essential proteins of P. aeruginosa, we conducted a comprehensive analysis focusing on three critical characteristics: virulence, resistance, and antigenicity. This approach ensures that the selected targets are not only essential for bacterial survival but also relevant for therapeutic intervention.
(a)
Virulence Factors: We used the Virulence Factor Database (VFDB) [61] (http://www.mgc.ac.cn/VFs/, accessed on 10 January 2024) to identify proteins involved in virulence. Virulence factors are crucial for the pathogenicity of P. aeruginosa, making them ideal targets for drug development. We applied stringent criteria with a cutoff set at a bit score value greater than 100 and an e-value exceeding 10−5;
(b)
Resistance Factors: Resistance factors were identified using the BacMet2 database, including both experimental and predicted data [62] (http://bacmet.biomedicine.gu.se/, accessed on 10 January 2024). Resistance mechanisms enable P. aeruginosa to withstand antibiotic treatments, so targeting these factors can enhance the efficacy of therapeutic strategies. We employed the same stringent cutoff criteria (bit score > 100, e-value < 10−5) to accurately detect resistance-related proteins;
(c)
Antigenicity Features: The antigenicity of essential proteins was predicted using VaxiJen [63] (https://www.ddg-pharmfac.net/vaxijen/VaxiJen/VaxiJen.html, accessed on 10 January 2024). This tool assesses the potential of proteins to elicit an immune response, which is vital for the development of vaccines and immunotherapies. Identifying antigenic proteins ensures that our targets can be used to stimulate protective immunity against P. aeruginosa infections.

3.4. Identifying Pharmacological Candidates for P. aeruginosa Essential Proteins

The 23 sequences that meet the criteria of virulence, bacterial resistance, and antigenicity were further assessed for their druggability using the DrugBank Database [64] (https://go.drugbank.com/, accessed on 30 January 2024). This database served as a comprehensive resource for assessing whether each essential protein identified in our pipeline is recognized as a viable druggable target in P. aeruginosa or other microorganisms by the FDA.
The primary focus of our search was on drugs capable of inhibiting the activity of specific essential proteins in P. aeruginosa. To validate the uniqueness of each target, the protein sequences were submitted to the DrugBank web server for comparison with existing druggable targets. The filters applied included considering only approved targets and specifying protein types as targets, enzymes, carriers, or transporters.

3.5. Characterizing the Biological Classification of P. aeruginosa Essential Proteins

The remaining 11 sequences capable of meeting all the previously evaluated criteria were analyzed for biological characterization:
(a)
Subcellular localization: The PSORT v.3.0 [65] (https://www.psort.org/psortb/, accessed on 2 March 2024) and CELLO2GO [66] (http://cello.life.nctu.edu.tw/cello2go/, https://www.psort.org/psortb/, accessed on 2 March 2024) databases were employed to determine the subcellular localization of the proteins.
(b)
Functional classification: The KEGG [67] (Kyoto Encyclopedia of Genes and Genomes, https://www.genome.jp/kegg/pathway.html, https://www.psort.org/psortb/, accessed on 2 March 2024) database was used to classify the proteins into functional biological categories.

3.6. Predicting 3D Structures and Supramolecular Assembly of Essential Proteins in P. aeruginosa

The filtration process identified 11 essential proteins necessitating the prediction of their three-dimensional (3D) structures to facilitate inhibitor identification. Utilizing their sequences, a BLAST search against the PDB database was conducted. For proteins lacking experimental structural data (NMR, X-ray, or Cryo-EM) in the PDB, the AlphaFold server was employed to predict their structures based on amino acid sequences. Once the structures were determined, an extensive review of the literature was conducted to elucidate their assembly patterns, whether as monomers, dimers, trimers, or higher-order supramolecular complexes. Additionally, BLAST searches against the PDB database were performed to identify structurally resolved proteins closely related to our proteins of interest. These homologous structures served as templates for three-dimensional alignment, aiding in the construction of supramolecular assemblies by replicating and aligning each subunit based on the templates to form functional multimeric complexes.

3.7. Refining 3D Structures through Molecular Dynamics Simulations

To enhance the accuracy and stability of the predicted 3D structures, MD simulations were conducted. This computational approach enables the exploration of protein dynamics over time, facilitating the attainment of energetically favorable conformations. In our simulations, proteins were immersed in a periodic boundary box of pre-equilibrated Single Point Charge (SPC) water molecules and neutralized with sodium chloride, ensuring an appropriate physiological environment. For membrane proteins, they were embedded within a patch of pre-equilibrated 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) membrane at 300 K before the addition of water molecules. Each system underwent thermodynamic equilibration following a protocol similar to that described by Bustos et al. [68], implemented using the Maestro/Schrodinger suite [69]. Subsequently, a 50 ns production simulation was conducted under an NPT ensemble (T = 300 K and P = 1 atm), using the OPLS3e force field [70]. Notably, no energetic constraints were applied to the protein during the simulation, allowing it to energetically relax solely based on environmental constraints.

3.8. Identification of Relevant Protein Conformations and Drug Repurposing through Ensemble-docking and Non-Covalent Interaction Index

To explore diverse protein conformations and facilitate drug repurposing, each MD trajectory was analyzed using the Bio3D library [71] in R to generate dendrograms based on RMSD. The cutree() function was employed to partition the trajectories into 5 clusters. From each cluster, 5 representative frames were randomly selected and converted into the PDBQT format with the Open Babel tool [72].
Subsequently, a collection of over 7000 FDA-approved, clinical trials, and pre-clinical compounds from various pathologies was obtained from the CLUE.io database [73]. These compounds underwent preprocessing using the Ligprep tool in the Maestro/Schrodinger suite to ensure proper ionization, stereochemistry, and 3D structure optimization. The resulting set comprised 5712 compounds after excluding cyclic peptides and accounting for different tautomeric and protonation states at pH 7.0. These compounds were further converted to PDBQT format using the Open Babel tool [72]. For each MFP protein, four different grids (32 Å × 32 Å × 32 Å) were generated in each frame, positioned at distinct internal regions. Similarly, for OMF proteins, three grids per frame were generated. These grids served as docking sites for the 5712 compounds, with 10 poses collected per compound, grid, and protein. This exhaustive docking procedure generated more than 2.5 million docking outputs for subsequent analysis.
To serve as a control, we utilized the tridimensional structure of MexAB-OprM from P. aeruginosa, obtained via electron microscopy at a resolution of 3.20 Å (PDB id: 6TA6 [74]). Additionally, we selected eight known RND pump inhibitors: P9D, CEF, DM2, MIY, oxacillin, MBX, rhodamine-6G, and erythromycin as ligand controls. The tripartite pump was separated into its three components (MexA, MexB, and OprM), and we conducted the same MD simulations, frame selections, and ensemble-docking protocol previously described for the studied proteins.
We computed the ranking of the eight ligand controls in MexB based on their median binding affinities. Once the top three ligand controls were identified, these were used as cutoffs in the docking evaluations for the 11 proteins under study, as well as for the MexA and OprM components. The top three median docking energies of these selected inhibitors served as benchmarks to evaluate the binding affinities of the 5712 drugs in the proteins studied.
The dockings were run with the VINA GPU implementation [75]. AutoDock Vina employs an empirical scoring function to estimate the binding affinity between ligands and proteins. This scoring function is inspired by the X-Score function [76] and includes terms for steric interactions, hydrophobic interactions, and hydrogen bonding [77]. Specifically, it combines attractive Gaussian and repulsive parabolic functions for steric interactions, linear functions for hydrophobic interactions, and piecewise linear functions for hydrogen bonding. Each term is weighted and summed to predict the overall binding affinity. The search algorithm used by Vina is an Iterated Local Search global optimization algorithm, which involves generating random initial conformations, performing mutations followed by local optimization.
Selection of the most promising inhibitor candidates was conducted using custom R scripts, which integrated various metrics such as minima (top scoring), arithmetic, geometric, and harmonic means, and medians. This approach, inspired by the methodology described by Yáñez et al. [50], aimed to identify energetically favorable interactions between compounds and proteins across the large population of data, ensuring robust inhibitor selection tailored to each protein under investigation.
To complement the extensive exploration carried out on molecular docking, approximately 4000 integration calculations were performed on the NCI regions for protein–ligand interactions [78]. These calculations provided detailed insight into the non-covalent forces involved in protein–ligand binding, including hydrogen bonds, hydrophobic interactions, and other types of intermolecular interactions. To perform the NCI calculation, the promolecular approximation was employed, which provides a description of the approximate electron density constructed by linearly combining the densities of the individual atoms/fragments that make up the molecular system, without considering their mutual polarization effects. The use of promolecular density allows NCI analysis to be applied to very large systems, such as proteins or nucleic acids, where a full quantum calculation is too computationally expensive. The promolecular approach provides a reasonable approximation of the true electron density at a much lower computational cost. This approach allowed an efficient and reliable evaluation of the characteristics of the non-covalent interactions in the protein–ligand complex, thus facilitating the interpretation of the molecular docking results and providing valuable information for the rational design of new bioactive compounds.
This analysis consists of deriving analytical expressions for the NCI equations and integrating them in the NCI regions to obtain closed-form solutions for the binding energies. Details of the approach used for this work are given below:
The NCI equations relate properties of the electron density (ρ) and its derivatives like the reduced density gradient (s) and the density Hessian eigenvalues (λi) to identify and characterize non-covalent interaction regions. Analytical expressions can be derived for quantities like the volume (V) and integrated density (q) of the attractive (λ2 < 0) and repulsive (λ2 > 0) NCI regions in terms of ρ, s, λi, etc. [79]
q a t t n = Ω ( N C I ) ρ n r   d r          λ 2 r < 0
q r e p n = Ω ( N C I ) ρ n r   d r          λ 2 r > 0
q b i n d n = q a t t n q r e p n
  • These analytical expressions for q can then be integrated over the NCI regions to obtain closed-form solutions relating them to the intermolecular separation;
  • The binding energy can be approximated as a functional of q a t t n q r e p n , and fit to reproduce known binding curves, e.g., E b i n d = f ( q a t t n q r e p n ) where n is an optimized exponent;
  • The key advantage is obtaining analytical, transferable expressions for E b i n d in terms of NCI region properties that are valid over the entire potential energy surface.
This fully analytical approach aims to directly relate binding energies to real-space properties of the electron density within the NCI regions through derived closed-form expressions. All these calculations have been carried out with the NCIplot V4.2 software (https://github.com/juliacontrerasgarcia/NCIPLOT-4.2, accessed on 10 January 2024) [80,81].
The mechanisms of action of the control inhibitors used in our study are well-documented. For instance, MBX2319 targets the AcrB efflux pump in E. coli by binding to the “hydrophobic trap” within the deep substrate binding pocket, inhibiting the pump’s function by preventing necessary conformational changes and hindering substrate binding [45]. D13-9001 (P9D) specifically inhibits the MexAB-OprM efflux pump in P. aeruginosa by binding to a unique site near the deep substrate binding pocket, preventing conformational changes required for pump function [46]. Both AcrB and MexB use a peristaltic mechanism powered by the proton motive force to transport substrates, and inhibitors like MBX2319 and P9D disrupt their function by binding to specific sites, enhancing antibiotic activity [48]. Our study, however, focuses on the upper two-thirds of the RND efflux pumps. Given this focus, it is challenging to directly identify the mechanisms of action, as the well-characterized mechanisms of the control inhibitors occur predominantly in the lower third of the pumps. Despite this, the relevance of our study lies in its emphasis on the less studied upper regions of the efflux pumps, opening new research avenues. We speculate that our identified molecules may obstruct the conduction pore or inhibit the assembly of the pump. Although this research is in its early stages, our identified inhibitor candidates exhibit high affinity across all proteins studied, presenting two significant advantages: (1) These compounds demonstrate high multi-target or broad-spectrum activity within the RND family, despite the studied proteins sharing a low percentage of identity. (2) They bind with high affinity not only to one component of the pump but to all three components, suggesting multiple potential modulatory sites within the channel.

4. Conclusions

In this study, we employed a comprehensive approach to identify and characterize essential proteins in P. aeruginosa (PAO1 strain), elucidating their potential as pharmacological targets.
Our analysis revealed 12 essential proteins, including a histidine kinase and 11 members of the RND efflux pump family (OMF and MFP proteins), crucial for coordinating resistance mechanisms and pathogenicity. By predicting their three-dimensional structures and multimeric assemblies, we facilitated the identification of potential inhibitor candidates through molecular docking simulations. Notably, we identified MK-3207, R-428, and Suramin as promising inhibitor candidates across multiple metrics, demonstrating their potential efficacy against a range of pathologies.
Furthermore, our study leveraged non-covalent interaction index methods to refine the selection of inhibitor candidates, providing insights into the precise electronic effects involved in protein–ligand interactions. Suramin consistently exhibited superior attractive energies across all evaluated proteins compared to MK-3207 and R-428, suggesting its potential as a broad-spectrum inhibitor.
The newly identified drug candidates, selected through a repurposing strategy, are already FDA-approved for other indications. This provides an advantage as these drugs have established pharmacokinetic and pharmacodynamic profiles, ensuring their efficacy in human patients. Candidates such as MK-3207, R-428 (Bemcentinib), and Suramin have shown promising binding affinities against the RND efflux pumps of P. aeruginosa, suggesting their potential effectiveness in combating multidrug-resistant strains by targeting and modulating efflux pump functions. Additionally, the safety profiles of these drugs are well-documented, reducing the risk of unforeseen adverse effects. Our study further ensures specificity by excluding targets with high homology to human proteins, minimizing potential off-target interactions and cross-reactivity. Overall, our findings contribute to the understanding of essential proteins in P. aeruginosa and provide a framework for the rational design of novel therapeutic interventions targeting bacterial virulence and resistance mechanisms.
The accuracy of the subtractive proteomics approach heavily depends on the comprehensiveness and accuracy of the databases used (e.g., VFDB, BacMet2, VaxiJen, and DEG). Any gaps or inaccuracies in these databases could affect the identification of virulence, resistance, and antigenic factors. Ensemble docking relies on a limited number of conformational states (five in our study). This may not fully capture the dynamic nature of the proteins, potentially overlooking some critical conformations that could affect binding affinity predictions. Given these limitations, experimental validation of our findings would be highly beneficial. Conducting in vitro and in vivo assays would confirm the computational predictions and provide more reliable insights into the efficacy of the identified drug candidates. Such validation steps are crucial to bridge the gap between computational predictions and practical therapeutic applications, ensuring that the proposed drug candidates are both effective and safe for use.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms25158027/s1.

Author Contributions

Conceptualization, D.B. and V.D.; methodology, G.U., S.E.R.-R. and O.Y.; software, E.W.H.-R.; validation, E.V.-M., J.M.P. and M.A.-A.; formal analysis, G.U. and C.A.F.-M.; investigation, R.S., L.M.-Q., V.D. and D.B.; resources, E.W.H.-R., R.S. and D.B.; data curation, E.V.-M. and O.Y.; writing—original draft preparation, V.D. and G.U.; writing—review and editing, R.S. and D.B.; visualization, S.E.R.-R. and G.U.; supervision, V.D. and D.B.; project administration, D.B.; funding acquisition, D.B. All authors have read and agreed to the published version of the manuscript.

Funding

D.B. would like to offer thanks to ANID FONDECYT de Iniciación en Investigación #11220444, and FOVI230136. The APC was partially funded by FOVI230136. R.S. would like to thank the Royal Society of Chemistry for research funding (Grants R19-3409, R20-6912 and R21-6448709305) and MCIN/AEI/10.13039/501100011033 (Grants PID2020-113147GA-I00 and PID2021-122839NB-I00). L.M.-Q. thanks ANILLO #ATE220014. E.W.H-R. thanks ANID FONDECYT de Iniciación en Investigación #11230033 and ANID FONDECYT de Postdoctorado #3170107. J.M.P. thanks ANID FONDECYT Iniciación en Investigación finance code #11201049 and ANID FONDECYT Regular #1241305.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Provenzani, A.; Hospodar, A.R.; Meyer, A.L.; Vinci, D.L.; Hwang, E.Y.; Butrus, C.M.; Polidori, P. Multidrug-resistant gram-negative organisms: A review of recently approved antibiotics and novel pipeline agents. Int. J. Clin. Pharm. 2020, 42, 1016–1025. [Google Scholar] [CrossRef] [PubMed]
  2. Thorpe, K.E.; Joski, P.; Johnston, K.J. Antibiotic-resistant infection treatment costs have doubled since 2002, now exceeding $2 billion annually. Health Aff. 2018, 37, 662–669. [Google Scholar] [CrossRef] [PubMed]
  3. World Health Organization. WHO Bacterial Priority Pathogens List, 2024; World Health Organization: Geneva, Switzerland, 2024. [Google Scholar]
  4. Cendra, M.d.M.; Torrents, E. Pseudomonas aeruginosa biofilms and their partners in crime. Biotechnol. Adv. 2021, 49, 107734. [Google Scholar] [CrossRef] [PubMed]
  5. Rossi, E.; La Rosa, R.; Bartell, J.A.; Marvig, R.L.; Haagensen, J.A.J.; Sommer, L.M.; Molin, S.; Johansen, H.K. Pseudomonas aeruginosa adaptation and evolution in patients with cystic fibrosis. Nat. Rev. Microbiol. 2021, 19, 331–342. [Google Scholar] [CrossRef] [PubMed]
  6. Bustos, D.; Hernández-Rodríguez, E.W.; Poblete, H.; Alzate-Morales, J.; Challier, C.; Boetsch, C.; Vergara-Jaque, A.; Beassoni, P. Structural Insights into the Inhibition Site in the Phosphorylcholine Phosphatase Enzyme of Pseudomonas aeruginosa. J. Chem. Inf. Model. 2022, 62, 3067–3078. [Google Scholar] [CrossRef] [PubMed]
  7. Vincent, J.-L.; Sakr, Y.; Singer, M.; Martin-Loeches, I.; Machado, F.R.; Marshall, J.C.; Finfer, S.; Pelosi, P.; Brazzi, L.; Aditianingsih, D.; et al. Prevalence and Outcomes of Infection among Patients in Intensive Care Units in 2017. JAMA J. Am. Med. Assoc. 2020, 323, 1478–1487. [Google Scholar] [CrossRef] [PubMed]
  8. Vidaillac, C.; Chotirmall, S.H. Pseudomonas aeruginosa in bronchiectasis: Infection, inflammation, and therapies. Expert Rev. Respir. Med. 2021, 15, 649–662. [Google Scholar] [CrossRef] [PubMed]
  9. Adamo, R.; Sokol, S.; Soong, G.; Gomez, M.I.; Prince, A. Pseudomonas aeruginosa flagella activate airway epithelial cells through asialoGM1 and toll-like receptor 2 as well as toll-like receptor 5. Am. J. Respir. Cell Mol. Biol. 2004, 30, 627–634. [Google Scholar] [CrossRef] [PubMed]
  10. Soong, G.; Reddy, B.; Sokol, S.; Adamo, R.; Prince, A. TLR2 is mobilized into an apical lipid raft receptor complex to signal infection in airway epithelial cells. J. Clin. Investig. 2004, 113, 1482–1489. [Google Scholar] [CrossRef]
  11. Ozer, E.; Yaniv, K.; Chetrit, E.; Boyarski, A.; Meijler, M.M.; Berkovich, R.; Kushmaro, A.; Alfonta, L. An inside look at a biofilm: Pseudomonas aeruginosa flagella biotracking. Sci. Adv. 2021, 7, eabg8581. [Google Scholar] [CrossRef]
  12. Colclough, A.L.; Alav, I.; Whittle, E.E.; Pugh, H.L.; Darby, E.M.; Legood, S.W.; McNeil, H.E.; Blair, J.M. RND efflux pumps in Gram-negative bacteria; regulation, structure and role in antibiotic resistance. Future Microbiol. 2020, 15, 143–157. [Google Scholar] [CrossRef] [PubMed]
  13. Bialvaei, A.Z.; Rahbar, M.; Hamidi-Farahani, R.; Asgari, A.; Esmailkhani, A.; Dashti, Y.M.; Soleiman-Meigooni, S. Expression of RND efflux pumps mediated antibiotic resistance in Pseudomonas aeruginosa clinical strains. Microb. Pathog. 2021, 153, 104789. [Google Scholar] [CrossRef] [PubMed]
  14. Maurya, S.; Akhtar, S.; Siddiqui, M.H.; Khan, M.K.A. Subtractive Proteomics for Identification of Drug Targets in Bacterial Pathogens: A Review. Int. J. Eng. Res. 2020, V9, 262–273. [Google Scholar] [CrossRef]
  15. Farha, M.A.; Brown, E.D. Drug repurposing for antimicrobial discovery. Nat. Microbiol. 2019, 4, 565–577. [Google Scholar] [CrossRef] [PubMed]
  16. Fuchs, S.; Engelmann, S. Small proteins in bacteria—Big challenges in prediction and identification. Proteomics 2023, 23, e2200421. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, F.; Xiao, J.; Pan, L.; Yang, M.; Zhang, G.; Jin, S.; Yu, J. A Systematic Survey of Mini-Proteins in Bacteria and Archaea. PLoS ONE 2008, 3, e4027. [Google Scholar] [CrossRef] [PubMed]
  18. Steiner, H.E.; Patterson, H.K.; Giles, J.B.; Karnes, J.H. Bringing pharmacomicrobiomics to the clinic through well-designed studies. Clin. Transl. Sci. 2022, 15, 2303–2315. [Google Scholar] [CrossRef]
  19. Tarasiuk, A.; Fichna, J. Gut microbiota: What is its place in pharmacology? Expert Rev. Clin. Pharmacol. 2019, 12, 921–930. [Google Scholar] [CrossRef]
  20. Aziz, R.K.; Hegazy, S.M.; Yasser, R.; Rizkallah, M.R.; ElRakaiby, M.T. Drug pharmacomicrobiomics and toxicomicrobiomics: From scattered reports to systematic studies of drug–microbiome interactions. Expert Opin. Drug Metab. Toxicol. 2018, 14, 1043–1055. [Google Scholar] [CrossRef]
  21. Weersma, R.K.; Zhernakova, A.; Fu, J. Interaction between drugs and the gut microbiome. Gut 2020, 69, 1510–1519. [Google Scholar] [CrossRef]
  22. Liao, C.; Huang, X.; Wang, Q.; Yao, D.; Lu, W. Virulence Factors of Pseudomonas aeruginosa and Antivirulence Strategies to Combat Its Drug Resistance. Front. Cell. Infect. Microbiol. 2022, 12, 926758. [Google Scholar] [CrossRef] [PubMed]
  23. Qin, S.; Xiao, W.; Zhou, C.; Pu, Q.; Deng, X.; Lan, L.; Liang, H.; Song, X.; Wu, M. Pseudomonas aeruginosa: Pathogenesis, virulence factors, antibiotic resistance, interaction with host, technology advances and emerging therapeutics. Signal Transduct. Target. Ther. 2022, 7, 199. [Google Scholar] [CrossRef] [PubMed]
  24. Stanislavsky, E.S.; Lam, J.S. Pseudomonas aeruginosa antigens as potential vaccines. FEMS Microbiol. Rev. 1997, 21, 243–277. [Google Scholar] [CrossRef] [PubMed]
  25. Zschiedrich, C.P.; Keidel, V.; Szurmant, H. Molecular Mechanisms of Two-Component Signal Transduction. J. Mol. Biol. 2016, 428, 3752–3775. [Google Scholar] [CrossRef] [PubMed]
  26. Fadel, F.; Bassim, V.; Francis, V.I.; Porter, S.L.; Botzanowski, T.; Legrand, P.; Perez, M.M.; Bourne, Y.; Cianférani, S.; Vincent, F. Insights into the atypical autokinase activity of the Pseudomonas aeruginosa GacS histidine kinase and its interaction with RetS. Structure 2022, 30, 1285–1297.e5. [Google Scholar] [CrossRef] [PubMed]
  27. Johnson, J.M.; Church, G.M. Alignment and Structure Prediction of Divergent Protein Families: Periplasmic and Outer Membrane Proteins of Bacterial Efflux Pumps. J. Mol. Biol. 1999, 287, 695–715. [Google Scholar] [CrossRef] [PubMed]
  28. Abadi, M.S.S.; Gholipour, A.; Hadi, N. The highly conserved domain of RND multidrug efflux pumps in pathogenic Gram-negative bacteria. Cell. Mol. Biol. 2018, 64, 79–83. [Google Scholar] [CrossRef]
  29. Lorusso, A.B.; Carrara, J.A.; Barroso, C.D.N.; Tuon, F.F.; Faoro, H. Role of Efflux Pumps on Antimicrobial Resistance in Pseudomonas aeruginosa. Int. J. Mol. Sci. 2022, 23, 15779. [Google Scholar] [CrossRef]
  30. Jamshidi, S.; Sutton, J.M.; Rahman, K.M. Mapping the dynamic functions and structural features of AcrB Efflux pump transporter using accelerated molecular dynamics simulations. Sci. Rep. 2018, 8, 10470. [Google Scholar] [CrossRef]
  31. López-Causapé, C.; Sommer, L.M.; Cabot, G.; Rubio, R.; Ocampo-Sosa, A.A.; Johansen, H.K.; Figuerola, J.; Cantón, R.; Kidd, T.J.; Molin, S.; et al. Evolution of the Pseudomonas aeruginosa mutational resistome in an international Cystic Fibrosis clone. Sci. Rep. 2017, 7, 5555. [Google Scholar] [CrossRef]
  32. Oliveira, W.K.; Ferrarini, M.; Morello, L.G.; Faoro, H. Resistome analysis of bloodstream infection bacterial genomes reveals a specific set of proteins involved in antibiotic resistance and drug efflux. NAR Genom. Bioinform. 2020, 2, lqaa055. [Google Scholar] [CrossRef] [PubMed]
  33. Alcalde-Rico, M.; Olivares-Pacheco, J.; Alvarez-Ortega, C.; Cámara, M.; Martínez, J.L. Role of the multidrug resistance efflux pump MexCD-OprJ in the Pseudomonas aeruginosa quorum sensing response. Front. Microbiol. 2018, 9, 2752. [Google Scholar] [CrossRef] [PubMed]
  34. Linares, J.F.; López, J.A.; Camafeita, E.; Albar, J.P.; Rojo, F.; Martínez, J.L. Overexpression of the multidrug efflux pumps MexCD-OprJ and MexEF-OprN is associated with a reduction of type III secretion in Pseudomonas aeruginosa. J. Bacteriol. 2005, 187, 1384–1391. [Google Scholar] [CrossRef] [PubMed]
  35. Kristensen, R.; Andersen, J.B.; Rybtke, M.; Jansen, C.U.; Fritz, B.G.; Kiilerich, R.O.; Uhd, J.; Bjarnsholt, T.; Qvortrup, K.; Tolker-Nielsen, T.; et al. Inhibition of Pseudomonas aeruginosa quorum sensing by chemical induction of the MexEF-oprN efflux pump. Antimicrob. Agents Chemother. 2024, 68, e0138723. [Google Scholar] [CrossRef] [PubMed]
  36. Mine, T.; Morita, Y.; Kataoka, A.; Mizushima, T.; Tsuchiya, T. Expression in Escherichia coli of a New Multidrug Efflux Pump, MexXY, from Pseudomonas aeruginosa. Antimicrob. Agents Chemother. 1999, 43, 415–417. [Google Scholar] [CrossRef] [PubMed]
  37. Seupt, A.; Schniederjans, M.; Tomasch, J.; Häussler, S. Expression of the MexXY Aminoglycoside Efflux Pump and Presence of an Aminoglycoside-Modifying Enzyme in Clinical Pseudomonas aeruginosa Isolates Are Highly Correlated. Antimicrob. Agents Chemother. 2020, 65, e01166-20. [Google Scholar] [CrossRef] [PubMed]
  38. Poole, K.; Krebes, K.; McNally, C.; Neshat, S. Multiple Antibiotic Resistance in Pseudomonas aeruginosa: Evidence for Involvement of an Efflux Operon. J. Bacteriol. 1993, 175, 7363–7372. [Google Scholar] [CrossRef] [PubMed]
  39. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef]
  40. Yonehara, R.; Yamashita, E.; Nakagawa, A. Crystal structures of OprN and OprJ, outer membrane factors of multidrug tripartite efflux pumps of Pseudomonas aeruginosa. Proteins Struct. Funct. Bioinform. 2016, 84, 759–769. [Google Scholar] [CrossRef]
  41. Jumper, J.; Evans, R.; Pritzel, A.; Green, T.; Figurnov, M.; Ronneberger, O.; Tunyasuvunakool, K.; Bates, R.; Žídek, A.; Potapenko, A.; et al. Highly accurate protein structure prediction with AlphaFold. Nature 2021, 596, 583–589. [Google Scholar] [CrossRef]
  42. Greenidge, P.A.; Kramer, C.; Mozziconacci, J.-C.; Sherman, W. Improving docking results via reranking of ensembles of ligand poses in multiple X-ray protein conformations with MM-GBSA. J. Chem. Inf. Model. 2014, 54, 2697–2717. [Google Scholar] [CrossRef] [PubMed]
  43. Palacio-Rodríguez, K.; Lans, I.; Cavasotto, C.N.; Cossio, P. Exponential consensus ranking improves the outcome in docking and receptor ensemble docking. Sci. Rep. 2019, 9, 5142. [Google Scholar] [CrossRef] [PubMed]
  44. Bajusz, D.; Rácz, A.; Héberger, K. Comparison of data fusion methods as consensus scores for ensemble docking. Molecules 2019, 24, 2690. [Google Scholar] [CrossRef] [PubMed]
  45. Aron, Z.; Opperman, T.J. Optimization of a novel series of pyranopyridine RND efflux pump inhibitors. Curr. Opin. Microbiol. 2016, 33, 1–6. [Google Scholar] [CrossRef] [PubMed]
  46. Nakashima, R.; Sakurai, K.; Yamasaki, S.; Hayashi, K.; Nagata, C.; Hoshino, K.; Onodera, Y.; Nishino, K.; Yamaguchi, A. Structural basis for the inhibition of bacterial multidrug exporters. Nature 2013, 500, 102–106. [Google Scholar] [CrossRef] [PubMed]
  47. Nishino, K.; Yamasaki, S.; Nakashima, R.; Zwama, M.; Hayashi-Nishino, M. Function and Inhibitory Mechanisms of Multidrug Efflux Pumps. Front. Microbiol. 2021, 12, 737288. [Google Scholar] [CrossRef] [PubMed]
  48. Murakami, S.; Nakashima, R.; Yamashita, E.; Matsumoto, T.; Yamaguchi, A. Crystal structures of a multidrug transporter reveal a functionally rotating mechanism. Nature 2006, 443, 173–179. [Google Scholar] [CrossRef] [PubMed]
  49. Zwama, M.; Nishino, K. Ever-adapting rnd efflux pumps in gram-negative multidrug-resistant pathogens: A race against time. Antibiotics 2021, 10, 774. [Google Scholar] [CrossRef] [PubMed]
  50. Yáñez, O.; Alegría-Arcos, M.; Suardiaz, R.; Morales-Quintana, L.; Castro, R.I.; Palma-Olate, J.; Galarza, C.; Catagua-González, Á.; Rojas-Pérez, V.; Urra, G.; et al. Calcium-Alginate-Chitosan Nanoparticle as a Potential Solution for Pesticide Removal, a Computational Approach. Polymers 2023, 15, 3020. [Google Scholar] [CrossRef]
  51. Salvatore, C.A.; Moore, E.L.; Calamari, A.; Cook, J.J.; Michener, M.S.; O’Malley, S.; Miller, P.J.; Sur, C.; Williams, D.L.; Zeng, Z.; et al. Pharmacological properties of MK-3207, a potent and orally active calcitonin gene-related peptide receptor antagonist. J. Pharmacol. Exp. Ther. 2010, 333, 152–160. [Google Scholar] [CrossRef]
  52. Holland, S.J.; Pan, A.; Franci, C.; Hu, Y.; Chang, B.; Li, W.; Duan, M.; Torneros, A.; Yu, J.; Heckrodt, T.J.; et al. R428, a selective small molecule inhibitor of Axl kinase, blocks tumor spread and prolongs survival in models of metastatic breast cancer. Cancer Res. 2010, 70, 1544–1554. [Google Scholar] [CrossRef]
  53. Wu, Z.-S.; Liu, C.F.; Fu, B.; Chou, R.-H.; Yu, C. Suramin blocks interaction between human FGF1 and FGFR2 D2 domain and reduces downstream signaling activity. Biochem. Biophys. Res. Commun. 2016, 477, 861–867. [Google Scholar] [CrossRef] [PubMed]
  54. De Clercq, E. Suramin in the treatment of AIDS: Mechanism of action. Antivir. Res. 1987, 7, 1–10. [Google Scholar] [CrossRef] [PubMed]
  55. Bateman, A.; Martin, M.-J.; Orchard, S.; Magrane, M.; Ahmad, S.; Alpi, E.; Bowler-Barnett, E.H.; Britto, R.; Bye-A-Jee, H.; Cukura, A.; et al. UniProt: The Universal Protein Knowledgebase in 2023. Nucleic Acids Res. 2023, 51, D523–D531. [Google Scholar] [CrossRef]
  56. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [PubMed]
  57. Chen, F. OrthoMCL-DB: Querying a comprehensive multi-species collection of ortholog groups. Nucleic Acids Res. 2006, 34, D363–D368. [Google Scholar] [CrossRef] [PubMed]
  58. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef]
  59. Shanmugham, B.; Pan, A. Identification and Characterization of Potential Therapeutic Candidates in Emerging Human Pathogen Mycobacterium abscessus: A Novel Hierarchical In Silico Approach. PLoS ONE 2013, 8, e59126. [Google Scholar] [CrossRef]
  60. Zhang, R.; Ou, H.-Y.; Zhang, C.-T. DEG: A database of essential genes. Nucleic Acids Res. 2004, 32, D271–D272. [Google Scholar] [CrossRef]
  61. Chen, L.; Yang, J.; Yu, J.; Yao, Z.; Sun, L.; Shen, Y.; Jin, Q. VFDB: A reference database for bacterial virulence factors. Nucleic Acids Res. 2005, 33, D325–D328. [Google Scholar] [CrossRef]
  62. Pal, C.; Bengtsson-Palme, J.; Rensing, C.; Kristiansson, E.; Larsson, D.G.J. BacMet: Antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 2014, 42, D737–D743. [Google Scholar] [CrossRef] [PubMed]
  63. Doytchinova, I.A.; Flower, D.R. VaxiJen: A server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinform. 2007, 8, 4. [Google Scholar] [CrossRef]
  64. Wishart, D.S.; Feunang, Y.D.; Guo, A.C.; Lo, E.J.; Marcu, A.; Grant, J.R.; Sajed, T.; Johnson, D.; Li, C.; Sayeeda, Z.; et al. DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res. 2018, 46, D1074–D1082. [Google Scholar] [CrossRef] [PubMed]
  65. Yu, N.Y.; Wagner, J.R.; Laird, M.R.; Melli, G.; Rey, S.; Lo, R.; Dao, P.; Sahinalp, S.C.; Ester, M.; Foster, L.J.; et al. PSORTb 3.0: Improved protein subcellular localization prediction with refined localization subcategories and predictive capabilities for all prokaryotes. Bioinformatics 2010, 26, 1608–1615. [Google Scholar] [CrossRef] [PubMed]
  66. Yu, C.-S.; Cheng, C.-W.; Su, W.-C.; Chang, K.-C.; Huang, S.-W.; Hwang, J.-K.; Lu, C.-H. CELLO2GO: A Web Server for Protein subCELlular LOcalization Prediction with Functional Gene Ontology Annotation. PLoS ONE 2014, 9, e99368. [Google Scholar] [CrossRef] [PubMed]
  67. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  68. Bustos, D.; Hernández-Rodríguez, E.W.; Castro, R.I.; Morales-Quintana, L. Structural Effects of pH Variation and Calcium Amount on the Microencapsulation of Glutathione in Alginate Polymers. Biomed Res. Int. 2022, 2022, 5576090. [Google Scholar] [CrossRef] [PubMed]
  69. Maestro, S. Schrödinger Release 2021-1; Schrödinger LLC: New York, NY, USA, 2021. [Google Scholar]
  70. Roos, K.; Wu, C.; Damm, W.; Reboul, M.; Stevenson, J.M.; Lu, C.; Dahlgren, M.K.; Mondal, S.; Chen, W.; Wang, L.; et al. OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules. J. Chem. Theory Comput. 2019, 15, 1863–1874. [Google Scholar] [CrossRef]
  71. Grant, B.J.; Skjærven, L.; Yao, X.Q. The Bio3D packages for structural bioinformatics. Protein Sci. 2021, 30, 20–30. [Google Scholar] [CrossRef] [PubMed]
  72. O’Boyle, N.M.; Banck, M.; James, C.A.; Morley, C.; Vandermeersch, T.; Hutchison, G.R. Open Babel: An Open chemical toolbox. J. Cheminform. 2011, 3, 33. [Google Scholar] [CrossRef]
  73. Corsello, S.; Bittker, J.A.; Liu, Z.; Gould, J.; McCarren, P.; Hirschman, J.E.; Johnston, S.E.; Vrcic, A.; Wong, B.; Khan, M.; et al. The Drug Repurposing Hub: A next-generation drug library and information resource. Nat. Med. 2017, 23, 405–408. [Google Scholar] [CrossRef]
  74. Glavier, M.; Puvanendran, D.; Salvador, D.; Decossas, M.; Phan, G.; Garnier, C.; Frezza, E.; Cece, Q.; Schoehn, G.; Picard, M.; et al. Antibiotic export by MexB multidrug efflux transporter is allosterically controlled by a MexA-OprM chaperone-like complex. Nat. Commun. 2020, 11, 4948. [Google Scholar] [CrossRef] [PubMed]
  75. Ding, J.; Tang, S.; Mei, Z.; Wang, L.; Huang, Q.; Hu, H.; Ling, M.; Wu, J. Vina-GPU 2.0: Further Accelerating AutoDock Vina and Its Derivatives with Graphics Processing Units. J. Chem. Inf. Model. 2023, 63, 1982–1998. [Google Scholar] [CrossRef]
  76. Wang, R.; Lai, L.; Wang, S. Further development and validation of empirical scoring functions for structure-based binding affinity prediction. J. Comput. Aided Mol. Des. 2002, 16, 11–26. [Google Scholar] [CrossRef] [PubMed]
  77. Trott, O.; Olson, A. Autodock vina: Improving the speed and accuracy of docking. J. Comput. Chem. 2019, 31, 455–461. [Google Scholar] [CrossRef]
  78. Johnson, E.R.; Keinan, S.; Mori-Sánchez, P.; Contreras-García, J.; Cohen, A.J.; Yang, W. Revealing noncovalent interactions. J. Am. Chem. Soc. 2010, 132, 6498–6506. [Google Scholar] [CrossRef] [PubMed]
  79. Contreras-García, J.; Yang, W.; Johnson, E.R. Analysis of hydrogen-bond interaction potentials from the electron density: Integration of noncovalent interaction regions. J. Phys. Chem. A 2011, 115, 12983–12990. [Google Scholar] [CrossRef] [PubMed]
  80. Contreras-García, J.; Johnson, E.R.; Keinan, S.; Chaudret, R.; Piquemal, J.-P.; Beratan, D.N.; Yang, W. NCIPLOT: A program for plotting noncovalent interaction regions. J. Chem. Theory Comput. 2011, 7, 625–632. [Google Scholar] [CrossRef]
  81. Boto, R.; Peccati, F.; Laplaza, R.; Quan, C.; Carbone, A.; Piquemal, J.-P.; Maday, Y.; Contreras-García, J. NCIPLOT4: Fast, Robust, and Quantitative Analysis of Noncovalent Interactions. J. Chem. Theory Comput. 2020, 16, 4150–4158. [Google Scholar] [CrossRef]
Figure 1. Schematic workflow: outlines each step involved in our computational methodology to first identify pharmacologically relevant proteins in P. aeruginosa, and subsequently predict potential drugs capable of modulating proteins in the RND family efflux pumps. The workflow includes stages of subtractive proteomics, identification of essential proteins, 3D modeling, and drug repurposing predictions through ensemble-docking and non-covalent interactions index.
Figure 1. Schematic workflow: outlines each step involved in our computational methodology to first identify pharmacologically relevant proteins in P. aeruginosa, and subsequently predict potential drugs capable of modulating proteins in the RND family efflux pumps. The workflow includes stages of subtractive proteomics, identification of essential proteins, 3D modeling, and drug repurposing predictions through ensemble-docking and non-covalent interactions index.
Ijms 25 08027 g001
Figure 2. Venn diagrams depicting characteristics of P. aeruginosa proteins. (A) Evaluation of 116 most relevant proteins based on their virulence (V), antigenicity (A), and bacterial resistance (B) properties. (B) Druggability assessment of the 23 proteins that satisfy the combined criteria of virulence, antigenicity, and bacterial resistance (VAB).
Figure 2. Venn diagrams depicting characteristics of P. aeruginosa proteins. (A) Evaluation of 116 most relevant proteins based on their virulence (V), antigenicity (A), and bacterial resistance (B) properties. (B) Druggability assessment of the 23 proteins that satisfy the combined criteria of virulence, antigenicity, and bacterial resistance (VAB).
Ijms 25 08027 g002
Figure 3. Three-dimensional structures of the 11 proteins were modeled in the study. (A) Outer Membrane Factors (OMF)-type proteins. (B) Membrane Fusion Protein (MFP)-type proteins. The membrane model represents the outer membrane of Gram-negative bacteria. Each color represents a distinct monomer of the protein, with three monomers in OMF and four in MFP proteins.
Figure 3. Three-dimensional structures of the 11 proteins were modeled in the study. (A) Outer Membrane Factors (OMF)-type proteins. (B) Membrane Fusion Protein (MFP)-type proteins. The membrane model represents the outer membrane of Gram-negative bacteria. Each color represents a distinct monomer of the protein, with three monomers in OMF and four in MFP proteins.
Ijms 25 08027 g003
Figure 4. Population of docking energies. (A) The boxplots depict the distribution of docking poses for the 5712 drugs repositioned against each of the OMF-type proteins. Within each boxplot, the red dot represents the median docking energy for the control inhibitor (P9D), serving as a reference. Additionally, the green boxplot represents the control protein (OprM) used in the study for comparative analysis. (B) This section presents the descriptive statistical analysis of the docking energy values obtained for each protein. It includes the first quartile (Q1), median, third quartile (Q3), average, and standard deviation (S.D.) of the docking energies. Furthermore, it highlights the median docking energy value for the control inhibitor (P9D), akin to the red dot in section (A). Additionally, it provides the number of compounds exhibiting docking energy values better than the control used as cutoff. (C,D) These sections mirror sections A and B but focus on MFP-type proteins, respectively, using the MexA protein as control (the green boxplot). Notate that OprM and MexA data are the same showed in Figure S4.
Figure 4. Population of docking energies. (A) The boxplots depict the distribution of docking poses for the 5712 drugs repositioned against each of the OMF-type proteins. Within each boxplot, the red dot represents the median docking energy for the control inhibitor (P9D), serving as a reference. Additionally, the green boxplot represents the control protein (OprM) used in the study for comparative analysis. (B) This section presents the descriptive statistical analysis of the docking energy values obtained for each protein. It includes the first quartile (Q1), median, third quartile (Q3), average, and standard deviation (S.D.) of the docking energies. Furthermore, it highlights the median docking energy value for the control inhibitor (P9D), akin to the red dot in section (A). Additionally, it provides the number of compounds exhibiting docking energy values better than the control used as cutoff. (C,D) These sections mirror sections A and B but focus on MFP-type proteins, respectively, using the MexA protein as control (the green boxplot). Notate that OprM and MexA data are the same showed in Figure S4.
Ijms 25 08027 g004
Figure 5. Promolecular energy distribution calculations for the top 3 drug candidates. (A) Energy distribution was calculated from the integral of the Non-Covalent Interactions (NCI) indices for all docking poses obtained for the top 3 candidates: MK-3207 (pink), R-428 (yellow), and Suramin (blue) against all OMF-type proteins. The NCI integral provides an energy value, with more negative values indicating stronger protein–ligand interactions. (B) Descriptive statistics for each OMF-type protein–ligand pair, including the first quartile (Q1), median, third quartile (Q3), average, and standard deviation (S.D.) of the affinity energy. (C) Energy distribution similar to (A), for the top 3 candidates against all MFP-type proteins. (D) Descriptive statistics for each MFP-type protein–ligand pair, including Q1, median, Q3, average, and S.D. of the affinity energy.
Figure 5. Promolecular energy distribution calculations for the top 3 drug candidates. (A) Energy distribution was calculated from the integral of the Non-Covalent Interactions (NCI) indices for all docking poses obtained for the top 3 candidates: MK-3207 (pink), R-428 (yellow), and Suramin (blue) against all OMF-type proteins. The NCI integral provides an energy value, with more negative values indicating stronger protein–ligand interactions. (B) Descriptive statistics for each OMF-type protein–ligand pair, including the first quartile (Q1), median, third quartile (Q3), average, and standard deviation (S.D.) of the affinity energy. (C) Energy distribution similar to (A), for the top 3 candidates against all MFP-type proteins. (D) Descriptive statistics for each MFP-type protein–ligand pair, including Q1, median, Q3, average, and S.D. of the affinity energy.
Ijms 25 08027 g005
Table 1. Subtractive proteomics analysis using P. aeruginosa proteins for finding relevant pharmacological targets. This table outlines the sequential filtering steps applied to the P. aeruginosa proteome to isolate essential proteins suitable for drug targeting. Each step details the number of proteins selected and discarded, refining the dataset to identify 116 essential proteins as potential pharmacological targets.
Table 1. Subtractive proteomics analysis using P. aeruginosa proteins for finding relevant pharmacological targets. This table outlines the sequential filtering steps applied to the P. aeruginosa proteome to isolate essential proteins suitable for drug targeting. Each step details the number of proteins selected and discarded, refining the dataset to identify 116 essential proteins as potential pharmacological targets.
ProtocolNumber of Proteins SelectedNumber of Proteins Discarded
Total of proteins
(P. aeruginsa proteome)
55640
Remove paralogs sequences5453111
Remove sequences lower than 100 amino acids5120333
Remove sequences homologous to human proteom4344776
Remove sequences orthologous to microbiota in humans3360984
Identification of proteins essential for bacterial survival1163244
Table 2. Functional characterization of P. aeruginosa essential proteins identified in our pipeline, categorized into Histidine Kinase, MFP-type proteins, and OMF-type proteins from the RND family.
Table 2. Functional characterization of P. aeruginosa essential proteins identified in our pipeline, categorized into Histidine Kinase, MFP-type proteins, and OMF-type proteins from the RND family.
Uniprot IDKEGG OrthologyKEGG Additional InformationSubcellular Localization Consensus
Q9HTZ0Two-component system
(Histidine Kinase)
Enzyme
(Transferases)
Signal transductionCytoplasmatic Membrane
Q9I0Y7Outer Membrane Factor
(Multidrug Efflux System)
TransportersAntimicrobial resistance genesOuter Membrane
Q9HXB9 *
Q9I0V8 *
Q9I006 *
Q9HU26 *
Q9HWH3
Q9HY88
Q9I0V5Membrane Fussion Protein
(Multidrug Efflux System)
TransportersAntimicrobial resistance genesCytoplasmatic Membrane
Q9I0Y9
Q9I3R2
G3XD25
* Hypothetical protein into PAO1 genome. The subcellular localization is a consensus of PSORT and CELLO servers.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Urra, G.; Valdés-Muñoz, E.; Suardiaz, R.; Hernández-Rodríguez, E.W.; Palma, J.M.; Ríos-Rozas, S.E.; Flores-Morales, C.A.; Alegría-Arcos, M.; Yáñez, O.; Morales-Quintana, L.; et al. From Proteome to Potential Drugs: Integration of Subtractive Proteomics and Ensemble Docking for Drug Repurposing against Pseudomonas aeruginosa RND Superfamily Proteins. Int. J. Mol. Sci. 2024, 25, 8027. https://doi.org/10.3390/ijms25158027

AMA Style

Urra G, Valdés-Muñoz E, Suardiaz R, Hernández-Rodríguez EW, Palma JM, Ríos-Rozas SE, Flores-Morales CA, Alegría-Arcos M, Yáñez O, Morales-Quintana L, et al. From Proteome to Potential Drugs: Integration of Subtractive Proteomics and Ensemble Docking for Drug Repurposing against Pseudomonas aeruginosa RND Superfamily Proteins. International Journal of Molecular Sciences. 2024; 25(15):8027. https://doi.org/10.3390/ijms25158027

Chicago/Turabian Style

Urra, Gabriela, Elizabeth Valdés-Muñoz, Reynier Suardiaz, Erix W. Hernández-Rodríguez, Jonathan M. Palma, Sofía E. Ríos-Rozas, Camila A. Flores-Morales, Melissa Alegría-Arcos, Osvaldo Yáñez, Luis Morales-Quintana, and et al. 2024. "From Proteome to Potential Drugs: Integration of Subtractive Proteomics and Ensemble Docking for Drug Repurposing against Pseudomonas aeruginosa RND Superfamily Proteins" International Journal of Molecular Sciences 25, no. 15: 8027. https://doi.org/10.3390/ijms25158027

APA Style

Urra, G., Valdés-Muñoz, E., Suardiaz, R., Hernández-Rodríguez, E. W., Palma, J. M., Ríos-Rozas, S. E., Flores-Morales, C. A., Alegría-Arcos, M., Yáñez, O., Morales-Quintana, L., D’Afonseca, V., & Bustos, D. (2024). From Proteome to Potential Drugs: Integration of Subtractive Proteomics and Ensemble Docking for Drug Repurposing against Pseudomonas aeruginosa RND Superfamily Proteins. International Journal of Molecular Sciences, 25(15), 8027. https://doi.org/10.3390/ijms25158027

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