Next Article in Journal
Antimicrobial Activity of Eucalyptus globulus, Azadirachta indica, Glycyrrhiza glabra, Rheum palmatum Extracts and Rhein against Porphyromonas gingivalis
Previous Article in Journal
Antibiotic Therapy in the Treatment of COVID-19 Pneumonia: Who and When?
Previous Article in Special Issue
Green and Chemical Silver Nanoparticles and Pomegranate Formulations to Heal Infected Wounds in Diabetic Rats
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Antibiotics
Volume 11
Issue 2
10.3390/antibiotics11020185
antibiotics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Drug Repurposing Targeting Pseudomonas aeruginosa MvfR Using Docking, Virtual Screening, Molecular Dynamics, and Free-Energy Calculations

by
Tatiana F. Vieira
1,2,
Rita P. Magalhães
1,2,
Manuel Simões
3 and
Sérgio F. Sousa
1,2,*
1
UCIBIO/REQUIMTE, BioSIM, Departamento de Medicina, Faculdade de Medicina da Universidade do Porto, Alameda Prof. Hernâni Monteiro, 4200-319 Porto, Portugal
2
Associate Laboratory i4HB, Institute for Health and Bioeconomy, Faculdade de Medicina, Universidade do Porto, 4200-319 Porto, Portugal
3
LEPABE Laboratory for Process Engineering, Environment, Biotechnology and Energy, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, s/n, 4200-465 Porto, Portugal
*
Author to whom correspondence should be addressed.
Antibiotics 2022, 11(2), 185; https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics11020185
Submission received: 14 December 2021 / Revised: 26 January 2022 / Accepted: 27 January 2022 / Published: 31 January 2022
(This article belongs to the Special Issue 10th Anniversary of Antibiotics - Recent Advances in Antibiotic Biosynthesis)
Browse Figures
Versions Notes

Abstract

:
Pseudomonas aeruginosa is an opportunistic Gram-negative bacterium responsible for acute and chronic infections in planktonic state or in biofilms. The sessile structures are known to confer physical stability, increase virulence, and work as a protective armor against antimicrobial compounds. P. aeruginosa can control the expression of genes, population density, and biofilm formation through a process called quorum sensing (QS), a rather complex and hierarchical system of communication. A recent strategy to try and overcome bacterial resistance is to target QS proteins. In this study, a combined multi-level computational approach was applied to find possible inhibitors against P. aeruginosa QS regulator protein MvfR, also known as PqsR, using a database of approved FDA drugs, as a repurposing strategy. Fifteen compounds were identified as highly promising putative MvfR inhibitors. On those 15 MvfR ligand complexes, molecular dynamic simulations and MM/GBSA free-energy calculations were performed to confirm the docking predictions and elucidate on the mode of interaction. Ultimately, the five compounds that presented better binding free energies of association than the reference molecules (a known antagonist, M64 and a natural inducer, 2-nonyl-4-hydroxyquinoline) were highlighted as very promising MvfR inhibitors.

1. Introduction

Pseudomonas aeruginosa is a highly adaptable Gram-negative bacterium, responsible for acute and chronic infections that rapidly evolve to multi-drug resistance to antibiotics, and is becoming extremely difficult to eradicate, leading to high morbidity and mortality rates [1,2]. It can be found in planktonic state or as aggregates, called biofilms. This association of microorganisms is the ultimate way of protection to adverse conditions [3] and involves a self-produced matrix of extracellular polymeric substances that confers stability and protection from external stress conditions [4,5,6,7].
P. aeruginosa has a mechanism of intracellular communication that controls the population density, expression of genes, and biofilm formation called quorum sensing (QS). It involves the production, detection, and response to extracellular signaling molecules, called autoinducers (AIs) [2,8,9]. QS in P. aeruginosa is rather complex and hierarchical. Until now, for P. aeruginosa, there are four types of signaling systems described. Two are based on acyl homoserine lactones (LasR and RhlR): one that uses quinolone as signaling molecule (PQS) and one whose mechanism and targets are still unknown (IQS) [10,11]. All these systems communicate with each other to control the direct and indirect expression of several virulence genes [12]. Details on the intricate web of interactions between these systems have already been defined extensively by several other groups [1,13,14,15,16].
A recent strategy to try to overcome bacterial resistance is to target QS proteins, inhibiting their activity. In this way, the bacterial virulence is eliminated, but the bacteria do not die, leading to lower acquired and adaptive resistance [1,17,18]. Targeting the PQS system may be advantageous because it is specific to P. aeruginosa. By not interfering with other QS systems, such as las and rhl, there is minimization of potential side effects, as bacteria that are essential to human health are not affected [19].
Pseudomonas aeruginosa QS regulator (PqsR), also called a multiple virulence factor regulator (MvfR), is a LysR-like transcriptional regulator protein (LTTRs) [20] common in many species of bacteria. It has a highly conserved DNA binding domain and a poorly conserved ligand binding domain, also called a co-inducer or a ligand-binding domain (CBD). The natural AIs are: 2-Heptyl-3-hydroxy-4(1H)-quinolone (also called the Pseudomonas Quinolone Signal, PQS) and, its precursor, 2-heptyl-4-hydroxyquinoline (HHQ). The synthesis and regulation of these AIs are controlled by the pqsABCDE operon which, in turn, is regulated by MvfR, creating an auto-regulatory loop. The expression of MvfR, on the other hand, is controlled by LasR [1,14,20,21,22].
This protein is not only an important target due to the direct and indirect activation of a variety of virulence genes, but it has been proven that it is also responsible for the formation and development of biofilms as well as antibiotic-tolerant persister (AT/P) cells [23]. Thus, this is an excellent target to control and treat P. aeruginosa infections.
The binding pocket of PqsR/MvfR has two subdomains, defined as Pocked A and Pocket B. The bicyclic ring of the natural autoinducer fits perfectly into pocket B and the aliphatic chains interact with residues on pocket A, as represented in Figure 1. It is a mainly hydrophobic pocket, which is one of the reasons why it is quite a challenging and high-interest target [24,25].
A review paper [26] was published describing the recent advances in the development of MvfR antagonists, leading the way to the rational design of new and specific drugs to treat P. aeruginosa infections. Analyzing all the compounds described, it is clear that there are common features and it is even possible to design a pharmacophore model [26]. To sum up, so far, the classes of molecules that have been tested and present inhibitory activity against MvfR are quinazolinones, quinolinones, benzamide–benzimidazoles, and hydroxybenzamides [20,23,24,27,28,29,30]. However, it is still unknown if these compounds will work on their own or if it is still necessary to use antibiotics, as there are still no clinical trials being made on this subject.
In this work, docking and virtual screening experiments were applied to discover new inhibitors for PqsR/MvfR using the ZINC FDA-approved database as a starting point. De novo synthesis is a very time-consuming and expensive process, and the search for new antibiotics has led to disappointing results with the main restricting factor being the lack of translation of in vitro activity in the live bacterial cells [31]. Drug repurposing for antimicrobial purposes aims to generate new clinical uses for already-established drugs with possible antibiotic effects or potentiate it when used in combination with an antibiotic [31,32].

2. Results and Discussion

The results presented in Table 1 show that the scoring functions (SFs) that can more effectively reproduce the binding poses of the crystallographic ligands in MvfR are GOLD’s CHEMPLP and ASP, with an average RMSD of 1.78 and 1.81 Å, respectively. ChemScore and VINA also gave reasonable results with average RMSD values of 3.00 Å and 3.57 Å, respectively. The worst results were obtained with LeDock and GoldScore with average RMSD values of 4.51 Å and 3.27 Å. Within individual scoring values, there were also significant variations among the different PDB–ligand combinations, arising from the diversity of molecules considered. QZN, NNQ, HLH, QAE, and Q25 have long aliphatic tails that are very flexible and hydrophobic, and many of the scoring functions could not reproduce the exact crystallographic pose and, hence, the high RMSD values. M64 is less flexible and bulkier, and most scoring functions were able to reproduce the original pose. Only ChemScore failed this test with a RMSD of 3.48 Å.
The same tendency is seen when we analyze the remaining data. It is more challenging for all the scoring functions to reproduce the crystallographic pose in ligands with long, aliphatic chains.
Upon further analysis, it was clear that the different SF showed a clear tendency in RMSD values and there were two targets that stood out: 4JVI and 6B8A. In these cases, the cross-docking scores (Supplementary Table S1) are consistently higher than any other target. 4JVI seems to accommodate ligands with aliphatic chains better, while 6B8A showed good scores for different types of ligand structures. The main difference between the two structures lies in the position of residue Ile186, as seen on Figure 2—a feature that can be important to accommodate a variety of ligands in a virtual screening (VS) run. At the cross-docking stage with structure 4JVI, the SFs were not able to dock M64 in the same pose as the crystallographic due to the position of Ile186, evidenced by the RMSD values.
Based on the results presented, two main structures were chosen (PDB codes: 4JVI and 6B8A) to continue the experiments, for the validation of the virtual screening (VS) protocol. Furthermore, at this stage, LeDock and GoldScore SFs were eliminated due to poor RMSD results.
Docking of the actives vs decoys dataset (2290 compounds) was performed against only the two PDB structures selected previously (4JVI and 6B8A). At the end of the calculations, each compound was ranked according to their binding score and the ROC curves were plotted (Figure 3). The curves show similar results for both targets, with GoldScore and LeDock being the worst at discriminating between binders and non-binders early on.
Table 2 summarizes the results obtained for these two targets for all scoring functions tested. Vina, CHEMPLP, ChemScore, and ASP provided good discriminatory ability between binders and non-binders, with ChemScore giving the lowest enrichment factor (EF) 1% for both targets. For 4JVI, Vina, CHEMPLP, and ASP provided the best EF 1% but then presented TG below 0.25. For 6B8A, Vina and ASP consistently provided the best results across most of the metrics studied. ASP exhibited an EF of 10%, an AUC above 0.5, and a TG of 0.25, indicating good performance and reproducibility of the protocol. Ultimately, ASP was the SF chosen to proceed to the VS of the FDA-approved database.
After performing the virtual screening protocol for the ZINC FDA-approved database, the molecules present in the top 1% were analyzed in detail, corresponding to 32 compounds. The top 15 molecules identified in the VS stage are listed in Table 3, along with their respective structure and docking score. A brief description of the pharmaceutical use of each compound is also provided.
Mellini and co-workers performed a similar VS study using a library of 1467 FDA-approved drugs against MvfR [25]. The five top hits were selected (Conivaptan, Ergotamine, Eltrombopag, Pimozide, and Dutasteride) and tested experimentally in in vitro assays. Only pimozide presented inhibitory activity toward MvfR. Eltrombopag was also one of the compounds present in the top 1% of our VS results (rank position nº 19) but, since this, a different protocol and a different SF were used; thus, it is only logical that different hits were obtained. Among the collection of compounds selected, the protocol was able to place two antibiotics in the top 15 results which reinforces the virtual screening methodology even further.
The subsequent stage of the study involved performing the 100 ns of MD simulations for each complex with 6B8A, followed by MM/GBSA calculations. As reference, a MD simulation of NNQ in complex with 6B8A was also performed. NNQ is one of the natural agonists of MvfR. The first antagonist with activity against MvfR and M64 was also used as reference.
The stability of the complexes was accessed by RMSD calculations for the Cα of the protein and ligands as well (supporting information for more detail). As presented in Table 4, all the MvfR–ligand complexes displayed low protein and ligand RMSD values throughout the simulations meaning that they are indeed stable.
Smaller average SASA values and a higher percentage of average buried ligand area indicate that the ligand is less exposed to the solvent, more shielded by the pocket, also suggesting greater stability. Of all the compounds simulated, isavuconazonium and methotrexate are these more exposed to the solvent (522.5 ± 47.4 Å2; 42.7 ± 0.1% and 312.6 ± 69.2 Å2; 53.8 ± 0.1%, respectively). The reference compound M64 exhibits the lowest SASA and higher percentage of potential ligand SASA buried (90.4 ± 19.3 Å2; 86.3 ± 0.03%), mainly because of the π–π interaction with Tyr258 and the ring portion of M64. The position of this residue and the position of Ile186 seem to be important for the “closing” or “opening” of pocket A, as seen on Figure 4.
Of the 15 FDA compounds studied, eight presented better binding free energies than the reference, M64 (−39.0 ± 0.1 kcal/mol). These were venetoclax (−70.1 ± 0.3 kcal/mol), indocyanine green (−58.6 ± 0.3 kcal/mol), nilotinib (−48.1 ± 0.2 kcal/mol), cabozantinib (−44.6 ± 0.2 kcal/mol), montelukast (−43.2 ± 0.2 kcal/mol), cefoperazone (−42.4 ± 0.4 kcal/mol), valrubicin (−41.1 ± 0.2 kcal/mol), and lomitapide (−40.6 ± 0.3 kcal/mol). The compounds highlighted as of potential therapeutic interest to target QS consider the SASA, % of buried ligand, and binding free-energy values, narrowing the selection to the top five (venetoclax, indocyanine green, nilotinib, cabozantinib, and montelukast).
Figure 5 compares the results from the MM/GBSA calculations and the GOLD ASP scores. The graph indicates that all the compounds selected present a higher affinity toward MvfR than the reference ligand, with venetoclax (−70.1 ± 5.7 kcal/mol for MM/GBSA), showing the highest affinity of all the molecules in the set. Isavuconazonium (−25.0 ± 0.1 kcal/mol), methotrexate (−22.8 ± 0.3 kcal/mol), and cafsulodin (−27.4 ± 0.3 kcal/mol) were the compounds with the lowest binding affinity, as previously mentioned. There is a clear tendency in terms of LogP values, with the most promising drugs showing higher LogP values. This can be an important aspect to consider because these compounds need to cross the bacterial cell wall to reach their specific target. The molecular weight of all the molecules selected is above 500 g/mol.
When the free-energy values obtained using MM/GBSA were decomposed, the contribution of each residue was analyzed. There are three main residues that contribute to an increased affinity in all MvfR–ligand complexes: Ile186, Ile236, and Ile263. The contributions of Val170, Leu207, Leu208, and Tyr258 does not exist for all the complexes, but their presence benefits the interaction for nilotinib, valrubicin, venetoclax, isavuconazonium, montelukast, and cefoperazone (Figure 6). Lapatinib also interacts strongly with Val211 and cefoperazone with Arg209 (data not shown in Figure 6).
When comparing these residues with the main residues interacting with NNQ (blue line), there was an increase in the number of interacting residues (in the NNQ case the interacting residues are only Leu208 and Ile236). This means that the association between these 15 selected FDA molecules is stronger than the interaction with its autoinducer. This might be explained by the fact that NNQs has a big aliphatic chain and the interactions in the binding pocket are mainly due to hydrophobic effects.
In Figure 7 and Figure 8, the interaction maps of the five selected FDA-approved compounds are depicted. Venetoclax is stabilized through hydrophobic interactions with residues Leu207, Tyr258, and Ile236. It can form hydrogen bond interactions with His184. It is a big molecule and occupies the total area of pocket B and pocket A, and also induces a shift in the position of Tyr258, opening pocket A even more.
Indocyanine green is mainly stabilized through hydrophobic interactions with residues Ile186, Leu207, and Ile236. It can form π–π interactions with Tyr258 and the oxygen atoms of the sulfate group form a salt bridge with Arg209. Indocyanine green causes a shift in the position of Tyr258 and Ile186, inducing the opening of pocket A.
Nilotinib is stabilized by hydrophobic interactions with residues Val170, Ile186, Leu208, Ile236, and Tyr258. Tyr258 also forms a hydrogen bond with nilotinib. It also induces a shift in the position of Tyr258 and Ile186, causing an opening of the binding pocket.
Cabozantinib is stabilized by hydrophobic interactions with residues Val170, Ile208, Ile236, and Tyr258. It induces a big shift in the position of Tyr258, forming a π–π interaction and closing pocket A at its extremity.
Montelukast is stabilized through hydrophobic interactions with residues Ile186, Leu189, Leu207, Ile236, and Tyr258. This is the only compound that induces a significative shift in the position of Ile186, closing pocket A, but leaving one of the aliphatic arms of montelukast outside.

3. Materials and Methods

3.1. Structure Identification and Analysis

The Protein Databank [45] and the Biofilms Structural Database [46] were explored to find molecular structures of P. aeruginosa PqsR/MvfR. A total of 12 X-ray structures were identified. Details are summarized in Table 5.
Structures 4JVC, 4JVD, and 4JVI (2013) [20] were the first to present the MvfR ligand binding domain in the apo form and in complex with its native agonist 2-nonyl-4-hydroxyquinoline or NNQ (4JVC and 4JVD, respectively). The authors have shown that the structure of this binding pocket is unusually large, in which a native AQ agonist is stabilized entirely by hydrophobic interactions [20]. They also presented a structure with an analogous to NNQ in the binding pocket (4JVI), i.e., the compound 3-amino-7-chloro-2-nonylquinazolin-4-one (QZN).
The PDB structure 6B8A (2018) presents the X-ray structure of the MvfR ligand-binding domain (LBD) in complex with one potent benzamide–benzimidazole (BB) inhibitor, i.e., 2-[(5-nitro-1H-benzimidazol-2-yl)sulfanyl]-N-(4-phenoxyphenyl)acetamide (or M64) [23]. The authors were able to demonstrate that this is a competitive inhibitor that binds in the same hydrophobic MvfR pocket as the natural inducers. The subtle conformational difference in the first three (4JVC, 4JVD, 4JVI) and fourth (6B8A) structures is mainly due to rearrangements in residues 181 to 191, in the presence of the natural inducers versus the inhibitor M64. This inhibitor, combined with sub-therapeutic doses of ciprofloxacin, was the first to show in vivo activity when tested in the treatment of a mouse lung infection model [28].
6Q7U, 6Q7V, and 6Q7W are structures of the LBD of MvfR in complex with 2-heptyl-4-quinolone or HHQ (6Q7U) and 2-aminopyridine derivatives, N-(4-(4-fluorophenyl)methyl)-6-(trifluoromethyl)pyridine-2,4-diamine or compound 11 (6Q7V), and N-(4-[3-(4-fluorophenyl)propyl]-6-(trifluoromethyl)pyridine-2,4-diamine or compound 20 (6Q7W). These derivatives were described by the Hartmann group as novel lead-like structure for the design of MvfR antagonists [22].
6TPR presented the MvfR LBD structure bound to another possible antagonist found through the optimization of a bacterial cell-based reporter assay hit. Soukarieh and co-workers were able to confirm that 2-[(5-methyl-[1,2,4]triazino[5,6-b]indol-3-yl)sulfanyl]-N-(4-pyridin-2-yloxyphenyl)ethanamide (or NV5) did reduce the biosynthesis of the alkylquinolone AIs even though it was unable to potentiate the effect of ciprofloxacin when combined [47].
6Z07, 6Z17, 6Z5K, and 6YZ3 are structures of the LBD of MvfR bound to four different thiazole-containing quinazolinones capable of inhibiting it. Grossman and his team were able to understand that the four ligands occupied the binding pocket in a similar manner to the AIs and that increasing the length of the aliphatic chain improved potency. For this reason, ligands 6-chloro-3((2-pentylthiazol-4-yl)methyl)quinazolin-4(3H)-one (or QAE) bound to 6Z5K and 6-chloro-3((2-hexylthiazol-4-yl)methyl)quinazolin-4(3H)-one (or Q25) bound to 6YZ3 were the ones that attenuated the production of the MvfR-regulated virulence factor pyocyanin while also showing low cytotoxicity in the in vitro assays [48].
A schematic representation of the workflow used is presented in Figure 9.
Initially, all the PDB structures of the proteins were analyzed, aligned, and treated using Pymol 2.3.0. Crystallographic waters were removed, and ligands were extracted and saved in separate files. All the structures were aligned so that the docking coordinates and conditions were the same for each structure in the re-docking and cross-docking experiments. After that, Gasteiger charges and polar hydrogens were assigned to the proteins using AutoDockTools [49].
The OpenBabel software [50] was used to prepare the ligands for docking. The chemical structures of the crystallographic ligands used in the validation stage are listed in Figure 10.

3.2. Protein–Ligand Docking Protocol Validation

Three docking software alternatives and six scoring functions (SFs) have been used for this study: Autodock Vina [51], LeDock, and Gold [52] (CHEMPLP, GoldScore, ChemScore, and ASP). The goal in testing all these different SFs is to find out which one works best for this specific target. The characteristics of the binding pocket and ligands in study have a big impact on the docking results [53,54].
To ensure reproducibility, the docking conditions were kept consistent for every software and every target. The parameters adjusted were the size of the docking box, the binding site definition, the number of runs, and search efficiency.
To validate the protocol, re-docking was performed by removing the crystallographic ligands and docking them again. This allows the user to evaluate the ability of the docking software to reproduce the geometry and orientation of the crystallographic pose. This was evaluated through the calculation of the root mean square deviation (RMSD) between the heavy atoms of the crystallographic and docked poses. Cross docking was also performed as a measure of validation. The goal here is to evaluate the robustness of different target structures and how they can accommodate different ligands. For cross docking, all the X-ray ligand structures isolated from the twelve PDB structures of this target were “docked” into the different X-ray structures. The RMSDs in both cases were calculated using DockRMSD [55].
The scoring functions used employ different metrics and scales. A good result is the one that presents a good score (depending on the SF used as different SFs use different metrics, it can be a high positive score—GOLD CHEMPLP, ASP, ChemScore, and GoldScore, or a more negative score—Vina, LeDock) and a RMSD below 2 Å.

3.3. Virtual Screening Protocol Optimization

The virtual screening protocol must be validated with benchmark datasets to ensure that it provides reliable results in the subsequent VS stage. This active versus decoy protocol is essential to validate the virtual screening conditions and evaluate the ability of each scoring function to discriminate between real binders and non-binders. A perfect scoring function would be capable of finding all the true binders early on and rank them higher than the decoys; however, that is not always the case due to the simplifications of the scoring functions. Virtual screening is meant to be fast as we want to screen large databases of compounds and that comes with a cost in accuracy [56]. So, special care was dedicated to improving the protocol’s ability in maximizing the number of known binders among the top solutions. For that, a PqsR/MvfR-specific virtual screening training library was prepared to evaluate and optimize the ability of the protocol in discriminating between known binders and non-binders.
After an initial query in the ChEMBL [57] and BindingDB [58] databases, 29 ligands were found to have reported experimental activity against MvfR. These 29 MVfR antagonists were first described by the Hartmann group in 2012 and 2013 using a ligand-based drug design approach [29] and fragment identification approach using surface plasmon resonance screening, respectively [24]. After going through the literature, an additional list of 11 manually curated active ligands was created [30,59,60,61]. Only the ligands that presented activity toward PqsR/MvfR were selected (i.e., compounds with IC50 values), raising the total of active molecules in the test set to 40.
Using the DUD-E [62] decoy generator, based on the ligands previously mentioned, a set of 50 decoys for each ligand was created. Decoys are molecules that resemble the ligands in their physical properties but are chemically and topologically different so that they are likely non-binders. This approach is useful to validate our scoring functions and build a benchmarking dataset. The total number of decoys generated was 2250, as 5 of the active molecules are present in two different protonated forms. As previously mentioned, both active molecules and decoys were prepared for docking using OpenBabel [50] and converted into pdbqt or mol2, depending on the specific docking software considered. A database of 2295 compounds was finally created and made ready for docking.
The six scoring functions were evaluated, and the metrics were calculated using a web-based application, Screening Explorer [63], as well as Excel. The metrics used for the evaluation of the VS results were the enrichment factor at 1%, receiver-operating characteristic (ROC) curves and the respective area under the curve (AUC), Boltzmann-enhanced discrimination of ROC (BEDROC), the robust initial enhancement (RIE), and total gain (TG). TG quantifies the discrimination of actives over decoys attributable to score variations. TG values over 0.25, combined with an AUC over 0.5, indicate a good performance and reproducibility from the VS protocol [63].

3.4. Virtual Screening for Drug Repurposing

After careful validation of the docking and virtual screening protocols, virtual screening (VS) experiments with the PqsR/MvfR target 6B8A were conducted using the ZINC FDA database of compounds to evaluate a possible drug repurposing strategy. For this stage, only the best scoring function and the X-ray structure that yield better discrimination of actives and decoys in the validation stages were selected.
The FDA-approved drugs library used was a subset of the ZINC [64] library, a free database of commercially available compounds for virtual screening. ZINC contains over 230 million purchasable compounds. At the time of the VS experiments, the FDA-approved drugs dataset had 3207 compounds that were all docked against the target.
The top 15 compounds identified in the VS were selected to be further studied by molecular dynamics simulations and free-energy calculations.

3.5. Molecular Dynamics Simulations

Molecular dynamic simulations are useful to validate the docking and VS predictions and provide insight into the protein–ligand interactions. It allows the study of the physical evolution of the system through time and is a valuable tool in the interpolation between theory and experiments [65,66,67]. The fifteen MvfR–ligand complexes were treated with the Leap module of Amber18 [68]. The MvfR protein was treated with the ff14SB force field [69]. All the FDA compounds were parameterized with Gaussian16 [70] using ANTECHAMBER, with RESP HF/6-31G(d) charges calculated with Gaussian16 and the general AMBER force field (GAFF) [71]. As control, two of the crystallographic ligands were used, i.e., NNQ as the natural agonist and M64 as an antagonist.
Sodium counter-ions were added to neutralize the overall charge in the system. The protein–ligand complexes were embedded into a box of TIP3P water molecules, whose edges were placed at least 12 Å away from each complex atom. Periodic boundaries were applied, and the long-range electrostatic interactions were calculated using the particle mesh Ewald summation method. The cut-off value for the short-range electrostatic and Lennard–Jones interactions was set at 10.0 Å. The hydrogen bonds were constrained using the SHAKE algorithm. A time step of 2 fs was used.
Four minimization steps were performed to remove clashed and applied in the following order: one water molecule (2500 steps); two hydrogens atoms (2500 steps); three side chains of all the amino acid residues (2500 steps); and four full systems (10,000 steps). After minimization, a molecular dynamics equilibration procedure was applied and divided into two stages: in the first stage (50 ps), the systems were gradually heated to 298 K using a Langevin thermostat at constant volume (NVT ensemble); in the second stage (50 ps), the density of the systems was further equilibrated at 298 K. Finally, the production phase was run for a total of 100 ns in an NPT ensemble at a pressure of 1 bar and a temperature of 298 K. This overall procedure has been previously used with success in the treatment of several biomolecular systems [72,73,74,75].
The molecular dynamic (MD) trajectories were analyzed using VMD [76] and the cpptraj tool [77].

3.6. Free-Energy Calculations

To estimate the binding free energies of the ligands toward MvfR, the molecular mechanics—generalized born surface area (MM-GBSA) method [78] was applied using the MM/PBSA.py script available in AMBER [79]. The calculations considered the last 90 ns of the MD simulation of every complex, using an interval of 10 frames, representing a total of 1800 frames per complex. For the MM/GBSA calculations, a salt concentration of 0.100 mol dm−3 was used.
The free-energy decomposition option was used to obtain information about the local interactions of the complex. Using per-residue decomposition, the contribution of each residue to the total free energy was estimated. This approach has been used with success in the study of several other systems, including quorum sensing inhibitors [72,73,80,81,82].
For the MvfR–FDA complexes with the top five binding free energies, two additional MD simulation replicas were performed per complex. MM-GBSA binding free energies were calculated, together with the additional MD simulations properties (average RMSd, SASA, number of hydrogen bonds), confirming the overall tendencies obtained (Suppelmentary Table S2).

4. Conclusions

A docking protocol was optimized for P. aeruginosa MfvR using the crystallographic ligands as validation tools in the reproducibility of the pose generated by the docking software. The VS protocol was adjusted based on known MfvR active ligands to obtain the best discriminatory ability between real binders and non-binders, and it was applied to a database of 3207 FDA-approved compounds.
The results obtained using the optimized VS protocol were further analyzed by MD simulations followed by free-energy calculations using the MM/GBSA method. This resulted in confirming five FDA-approved compounds with a high probability of exhibiting activity as P. aeruginosa QS inhibitors: venetoclax, indocyanine green, nilotinib, cabozantinib, and montelukast. Throughout the MD simulations, it was possible to observe shifts in the position of two residues that might be key in the activation mechanism of MvfR, i.e., Ile186 and Tyr258. π–π interactions with Tyr258 seem to have a crucial role in protein ligand affinity. In addition to the identification of this new potential application for drug repurposing of these molecules and demonstrating their mode of interaction, the protocol here described and validated can now be applied to large libraries of drug-like compounds to highlight new promising candidates for P. aeruginosa QS inhibitors. This study also highlighted the key role in molecular recognition played by residues Ile186 and Tyr258, a feature that can be taken into consideration in future rational drug design and optimization efforts targeting this important QS protein to develop new anti-biofilm molecules.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/antibiotics11020185/s1, Table S1: Cross-docking scores obtained with the six scoring functions tested. Table S2. Summary of the MD simulations results on the replicas performed for the top 5 MvfR–FDA complexes.

Author Contributions

Conceptualization, S.F.S.; methodology, T.F.V. and R.P.M.; software, T.F.V. and R.P.M.; validation, T.F.V.; formal analysis, T.F.V.; investigation, T.F.V.; resources, T.F.V.; data curation, T.F.V.; writing—original draft preparation, T.F.V.; writing—review and editing, S.F.S. and M.S.; visualization, T.F.V.; supervision, S.F.S. and M.S.; project administration, S.F.S.; funding acquisition, S.F.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by national funds from Fundação para a Ciência e a Tecnologia (grant numbers: SFRH/BD/137844/2018, 2020.09087.BD, UIDP/04378/2020, UIDB/04378/2020, UIDB/00511/2020, PTDC/BII-BTI/30219/201-POCI-01-0145-FEDER-030219 and 2020.01423. CEECIND/CP1596/CT0003). Some of the calculations were produced with the support of the INCD funded by the FCT and FEDER under project 01/SAICT/2016 number 022153, and projects CPCA/A00/7140/2020 and CPCA/A00/7145/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The database of compounds used in this project is available at https://zinc.docking.org/substances/home/ (accessed on 1 January 2021).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Allegretta, G.; Maurer, C.K.; Eberhard, J.; Maura, D.; Hartmann, R.W.; Rahme, L.; Empting, M. In-depth Profiling of MvfR-Regulated Small Molecules in Pseudomonas aeruginosa after Quorum Sensing Inhibitor Treatment. Front. Microbiol. 2017, 24, 924. [Google Scholar] [CrossRef] [PubMed]
  2. Williams, P.; Cámara, M. Quorum sensing and environmental adaptation in Pseudomonas aeruginosa: A tale of regulatory networks and multifunctional signal molecules. Curr. Opin. Microbiol. 2009, 12, 182–191. [Google Scholar] [CrossRef]
  3. Kamaruzzaman, N.F.; Tan, L.P.; Mat Yazid, K.A.; Saeed, S.I.; Hamdan, R.H.; Choong, S.S.; Wong, W.K.; Chivu, A.; Gibson, A.J. Targeting the Bacterial Protective Armour; Challenges and Novel Strategies in the Treatment of Microbial Biofilm. Materials 2018, 11, 1705. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jamal, M.; Tasneem, U.; Hussain, T.; Andleeb, S. Bacterial Biofilm: Its Composition, Formation and Role in Human Infections. J. Microbiol. Biotechnol. 2015, 4, 1–14. [Google Scholar]
  5. Dufour, D.; Leung, V.; Lévesque, C.M. Bacterial biofilm: Structure, function, and antimicrobial resistance. Endod. Top. 2010, 22, 2–16. [Google Scholar] [CrossRef]
  6. Jamal, M.; Ahmad, W.; Andleeb, S.; Jalil, F.; Imran, M.; Nawaz, M.A.; Hussain, T.; Ali, M.; Rafiq, M.; Kamil, M.A. Bacterial biofilm and associated infections. J. Chin. Med. Assoc. 2018, 81, 7–11. [Google Scholar] [CrossRef] [PubMed]
  7. Satpathy, S.; Sen, S.K.; Pattanaik, S.; Raut, S. Review on bacterial biofilm: An universal cause of contamination. Biocatal. Agric. Biotechnol. 2016, 7, 56–66. [Google Scholar] [CrossRef]
  8. Moradali, M.F.; Ghods, S.; Rehm, B.H.A. Pseudomonas aeruginosa Lifestyle: A Paradigm for Adaptation, Survival, and Persistence. Front. Cell. Infect. Microbiol. 2017, 7, 1–29. [Google Scholar] [CrossRef] [Green Version]
  9. Ghosh, R.; Das, A.; Mallik, S. Inhibition of Quorum Sensing in Pseudomonas aeruginosa: A Review. Indian J. Pharm. Sci. 2019, 81, 797–806. [Google Scholar] [CrossRef]
  10. Li, S.; Chen, S.; Fan, J.; Cao, Z.; Ouyang, W.; Tong, N.; Hu, X.; Hu, J.; Li, P.; Feng, Z.; et al. Anti-biofilm effect of novel thiazole acid analogs against Pseudomonas aeruginosa through IQS pathways. Eur. J. Med. Chem. 2018, 145, 64–73. [Google Scholar] [CrossRef]
  11. Abelyan, N.; Grabski, H.; Tiratsuyan, S. In silico Screening of Flavones and its Derivatives as Potential Inhibitors of Quorum-Sensing Regulator LasR of Pseudomonas aeruginosa. Mol. Biol. 2020, 54, 153–163. [Google Scholar] [CrossRef]
  12. Lazdunski, A.M.; Ventre, I.; Sturgis, J.N. Regulatory circuits and communication in gram-negative bacteria. Nat. Rev. Microbiol. 2004, 2, 581–592. [Google Scholar] [CrossRef] [PubMed]
  13. Maura, D.; Hazan, R.; Kitao, T.; Ballok, A.E.; Rahme, L.G. Evidence for Direct Control of Virulence and Defense Gene Circuits by the Pseudomonas aeruginosa Quorum Sensing Regulator, MvfR. Sci. Rep. 2016, 6, 34083. [Google Scholar] [CrossRef]
  14. Schütz, C.; Empting, M. Targeting the Pseudomonas quinolone signal quorum sensing system for the discovery of novel anti-infective pathoblockers. Beilstein J. Org. Chem. 2018, 14, 2627–2645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Venturi, V. Regulation of quorum sensing in Pseudomonas. FEMS Microbiol. Rev. 2006, 30, 274–291. [Google Scholar] [CrossRef] [Green Version]
  16. Miller, M.B.; Bassler, B.L. Quorum Sensing in Bacteria. Annu. Rev. Microbiol. 2001, 55, 165–199. [Google Scholar] [CrossRef] [Green Version]
  17. DIckey, S.W.; Cheung, G.Y.C.; Otto, M. Different drugs for bad bugs: Antivirulence strategies in the age of antibiotic resistance. Nat. Rev. Drug Discov. 2017, 16, 457–471. [Google Scholar] [CrossRef]
  18. Mühlen, S.; Dersch, P. Anti-virulence strategies to target bacterial infections. In Current Topics in Microbiology and Immunology; Springer: Cham, Switzerland, 2015; Volume 398, pp. 147–183. [Google Scholar]
  19. Thomann, A.; Brengel, C.; Börger, C.; Kail, D.; Steinbach, A.; Empting, M.; Hartmann, R.W. Structure-Activity Relationships of 2-Sufonylpyrimidines as Quorum-Sensing Inhibitors to Tackle Biofilm Formation and eDNA Release of Pseudomonas aeruginosa. ChemMedChem 2016, 11, 2522–2533. [Google Scholar] [CrossRef]
  20. Ilangovan, A.; Fletcher, M.; Rampioni, G.; Pustelny, C.; Rumbaugh, K.; Heeb, S.; Cámara, M.; Truman, A.; Chhabra, S.R.; Emsley, J.; et al. Structural Basis for Native Agonist and Synthetic Inhibitor Recognition by the Pseudomonas aeruginosa Quorum Sensing Regulator PqsR (MvfR). PLoS Pathog. 2013, 9, e1003508. [Google Scholar] [CrossRef] [Green Version]
  21. Lin, J.; Cheng, J.; Wang, Y.; Shen, X. The Pseudomonas Quinolone Signal (PQS): Not Just for Quorum Sensing Anymore. Front. Cell. Infect. Microbiol. 2018, 8, 230. [Google Scholar] [CrossRef]
  22. Zender, M.; Witzgall, F.; Kiefer, A.; Kirsch, B.; Maurer, C.K.; Kany, A.M.; Xu, N.; Schmelz, S.; Börger, C.; Blankenfeldt, W.; et al. Flexible Fragment Growing Boosts Potency of Quorum-Sensing Inhibitors against Pseudomonas aeruginosa Virulence. ChemMedChem 2019, 14, 188–194. [Google Scholar] [CrossRef] [Green Version]
  23. Kitao, T.; Lepine, F.; Babloudi, S.; Walte, F.; Steinbacher, S.; Maskos, K.; Blaesse, M.; Negri, M.; Pucci, M.; Zahler, B.; et al. Molecular Insights into Function and Competitive Inhibition of Pseudomonas aeruginosa Multiple Virulence Factor Regulator. mBio 2018, 9, 611–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Zender, M.; Klein, T.; Henn, C.; Kirsch, B.; Maurer, C.K.; Kail, D.; Ritter, C.; Dolezal, O.; Steinbach, A.; Hartmann, R.W. Discovery and Biophysical Characterization of 2-Amino-oxadiazoles as Novel Antagonists of PqsR, an Important Regulator of Pseudomonas aeruginosa Virulence. J. Med. Chem. 2013, 56, 6761–6774. [Google Scholar] [CrossRef] [PubMed]
  25. Mellini, M.; Di Muzio, E.; D’Angelo, F.; Baldelli, V.; Ferrillo, S.; Visca, P.; Leoni, L.; Polticelli, F.; Rampioni, G. In silico Selection and Experimental Validation of FDA-Approved Drugs as Anti-quorum Sensing Agents. Front. Microbiol. 2019, 10, 2355. [Google Scholar] [CrossRef] [PubMed]
  26. Soheili, V.; Tajani, A.S.; Ghodsi, R.; Bazzaz, B.S.F. Anti-PqsR compounds as next-generation antibacterial agents against Pseudomonas aeruginosa: A review. Eur. J. Med. Chem. 2019, 172, 26–35. [Google Scholar] [CrossRef] [PubMed]
  27. Klein, T.; Henn, C.; De Jong, J.C.; Zimmer, C.; Kirsch, B.; Maurer, C.K.; Pistorius, D.; Müller, R.; Steinbach, A.; Hartmann, R.W. Identification of small-molecule antagonists of the Pseudomonas aeruginosa transcriptional regulator PqsR: Biophysically guided hit discovery and optimization. ACS Chem. Biol. 2012, 7, 1496–1501. [Google Scholar] [CrossRef]
  28. Starkey, M.; Lepine, F.; Maura, D.; Bandyopadhaya, A.; Lesic, B.; He, J.; Kitao, T.; Righi, V.; Milot, S.; Tzika, A.; et al. Identification of Anti-virulence Compounds That Disrupt Quorum-Sensing Regulated Acute and Persistent Pathogenicity. PLoS Pathog. 2014, 10, e1004321. [Google Scholar] [CrossRef]
  29. Lu, C.; Kirsch, B.; Zimmer, C.; de Jong, J.C.; Henn, C.; Maurer, C.K.; Müsken, M.; Häussler, S.; Steinbach, A.; Hartmann, R.W. Discovery of antagonists of PqsR, a key player in 2-alkyl-4-quinolone-dependent quorum sensing in Pseudomonas aeruginosa. Chem. Biol. 2012, 19, 381–390. [Google Scholar] [CrossRef] [Green Version]
  30. Soukarieh, F.; Vico Oton, E.; Dubern, J.F.; Gomes, J.; Halliday, N.; de Pilar Crespo, M.; Ramírez-Prada, J.; Insuasty, B.; Abonia, R.; Quiroga, J.; et al. In Silico and in Vitro-Guided Identification of Inhibitors of Alkylquinolone-Dependent Quorum Sensing in Pseudomonas aeruginosa. Molecules 2018, 23, 257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Boyd, N.K.; Teng, C.; Frei, C.R. Brief Overview of Approaches and Challenges in New Antibiotic Development: A Focus on Drug Repurposing. Front. Cell. Infect. Microbiol. 2021, 11, 442. [Google Scholar] [CrossRef]
  32. Farha, M.A.; Brown, E.D. Drug repurposing for antimicrobial discovery. Nat. Microbiol. 2019, 4, 565–577. [Google Scholar] [CrossRef] [PubMed]
  33. Choulis, N.H. Miscellaneous drugs, materials, medical devices, and techniques. In Side Effects of Drugs Annual; Elsevier: Amsterdam, The Netherlands, 2009; pp. 757–769. [Google Scholar]
  34. Alonso, R.; Cuevas, A.; Mata, P. Lomitapide: A review of its clinical use, efficacy, and tolerability. Core Evid. 2019, 14, 19–30. [Google Scholar] [CrossRef] [Green Version]
  35. Opdam, F.L.; Guchelaar, H.; Beijnen, J.H.; Schellens, J.H.M. Lapatinib for Advanced or Metastatic Breast Cancer. Oncologist 2012, 17, 536–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Lee, A.T.J.; Jones, R.L.; Huang, P.H. Pazopanib in advanced soft tissue sarcomas. Signal Transduct. Target Ther. 2019, 4, 16. [Google Scholar] [CrossRef] [Green Version]
  37. Longo, N.; Diaz, G.A.; Lichter-Konecki, U.; Schulze, A.; Inbar-Feigenberg, M.; Conway, R.L.; Bannick, A.A.; McCandless, S.E.; Zori, R.; Hainline, B.; et al. Glycerol phenylbutyrate efficacy and safety from an open label study in pediatric patients under 2 months of age with urea cycle disorders. Mol. Genet. Metab. 2021, 132, 19–26. [Google Scholar] [CrossRef] [PubMed]
  38. Abdelaziz, A.; Vaishampayan, U. Cabozantinib for the treatment of kidney cancer. Expert Rev. Anticancer Ther. 2017, 17, 577–584. [Google Scholar] [CrossRef] [PubMed]
  39. Li, Q.; Cheng, L.; Shen, K.; Jin, H.; Li, H.; Cheng, Y.; Ma, X. Efficacy and Safety of Bcl-2 Inhibitor Venetoclax in Hematological Malignancy: A Systematic Review and Meta-Analysis of Clinical Trials. Front. Pharmacol. 2019, 10, 697. [Google Scholar] [CrossRef]
  40. McManus, D.S. Antifungal drugs. In Side Effects of Drugs Annual; Elsevier: Amsterdam, The Netherlands, 2016; pp. 243–253. [Google Scholar]
  41. Bedoui, Y.; Guillot, X.; Sélambarom, J.; Guiraud, P.; Giry, C.; Jaffar-Bandjee, M.C.; Ralandison, S.; Gasque, P. Methotrexate an Old Drug with New Tricks. Int. J. Mol. Sci. 2019, 20, 5023. [Google Scholar] [CrossRef] [Green Version]
  42. Smith, B.R. Cefsulodin and ceftazidime, two antipseudomonal cephalosporins. Clin. Pharm. 1984, 3, 373–385. [Google Scholar] [CrossRef]
  43. Paggiaro, P.; Bacci, E. Montelukast in asthma: A review of its efficacy and place in therapy. Ther. Adv. Chronic Dis. 2011, 2, 47–58. [Google Scholar] [CrossRef]
  44. Castle, S.S. Cefoperazone. In XPharm: The Comprehensive Pharmacology Reference; Elsevier: Amsterdam, The Netherlands, 2007; pp. 1–5. [Google Scholar]
  45. Berman, H.M.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] [PubMed] [Green Version]
  46. Magalhães, R.P.; Vieira, T.F.; Fernandes, H.S.; Melo, A.; Simões, M.; Sousa, S.F. The Biofilms Structural Database. Trends Biotechnol. 2020, 38, 937–940. [Google Scholar] [CrossRef]
  47. Soukarieh, F.; Liu, R.; Romero, M.; Roberston, S.N.; Richardson, W.; Lucanto, S.; Oton, E.V.; Qudus, N.R.; Mashabi, A.; Grossman, S.; et al. Hit Identification of New Potent PqsR Antagonists as Inhibitors of Quorum Sensing in Planktonic and Biofilm Grown Pseudomonas aeruginosa. Front. Chem. 2020, 8, 204. [Google Scholar] [CrossRef] [PubMed]
  48. Grossman, S.; Soukarieh, F.; Richardson, W.; Liu, R.; Mashabi, A.; Emsley, J.; Williams, P.; Cámara, M.; Stocks, M.J. Novel quinazolinone inhibitors of the Pseudomonas aeruginosa quorum sensing transcriptional regulator PqsR. Eur. J. Med. Chem. 2020, 208, 112778. [Google Scholar] [CrossRef] [PubMed]
  49. Morris, G.M.; Huey, R.; Lindstrom, W.; Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 2009, 30, 2785–2791. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. 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, 757–769. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2009, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
  52. Jones, G.; Willett, P.; Glen, R.C.; Leach, A.R.; Taylor, R. Development and validation of a genetic algorithm for flexible docking. J. Mol. Biol. 1997, 267, 727–748. [Google Scholar] [CrossRef] [Green Version]
  53. Vieira, T.F.; Sousa, S.F. Comparing AutoDock and Vina in Ligand/Decoy Discrimination for Virtual Screening. Appl. Sci. 2019, 9, 4538. [Google Scholar] [CrossRef] [Green Version]
  54. Vieira, T.F.; Magalhães, R.P.; Sousa, S.F. Tailoring Specialized Scoring Functions for More Efficient Virtual Screening. Front. Drug Chem. Clin. Res. 2019, 2, 1–4. [Google Scholar]
  55. Bell, E.W.; Zhang, Y. DockRMSD: An open-source tool for atom mapping and RMSD calculation of symmetric molecules through graph isomorphism. J. Cheminform. 2019, 11, 40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Cerqueira, N.M.F.S.A.; Gesto, D.; Oliveira, E.F.; Santos-Martins, D.; Brás, N.F.; Sousa, S.F.; Fernandes, P.A.; Ramos, M.J. Receptor-based virtual screening protocol for drug discovery. Arch. Biochem. Biophys. 2015, 582, 56–67. [Google Scholar] [CrossRef] [PubMed]
  57. Gaulton, A.; Hersey, A.; Nowotka, M.; Bento, A.P.; Chambers, J.; Mendez, D.; Mutowo, P.; Atkinson, F.; Bellis, L.J.; Cibrián-Uhalte, E.; et al. The ChEMBL database in 2017. Nucleic Acids Res. 2017, 45, D945–D954. [Google Scholar] [CrossRef]
  58. Gilson, M.K.; Liu, T.; Baitaluk, M.; Nicola, G.; Hwang, L.; Chong, J. BindingDB in 2015: A public database for medicinal chemistry, computational chemistry and systems pharmacology. Nucleic Acids Res. 2016, 44, D1045–D1053. [Google Scholar] [CrossRef] [PubMed]
  59. Lu, C.; Kirsch, B.; Maurer, C.K.; de Jong, J.C.; Braunshausen, A.; Steinbach, A.; Hartmann, R.W. Optimization of anti-virulence PqsR antagonists regarding aqueous solubility and biological properties resulting in new insights in structure–activity relationships. Eur. J. Med. Chem. 2014, 79, 173–183. [Google Scholar] [CrossRef]
  60. Hossain, M.A.; Sattenapally, N.; Parikh, H.I.; Li, W.; Rumbaugh, K.P.; German, N.A. Design, synthesis, and evaluation of compounds capable of reducing Pseudomonas aeruginosa virulence. Eur. J. Med. Chem. 2020, 185, 111800. [Google Scholar] [CrossRef]
  61. Lesic, B.; Lépine, F.; Déziel, E.; Zhang, J.; Zhang, Q.; Padfield, K.; Castonguay, M.H.; Milot, S.; Stachel, S.; Tzika, A.A.; et al. Inhibitors of Pathogen Intercellular Signals as Selective Anti-Infective Compounds. PLoS Pathog. 2007, 3, e126. [Google Scholar] [CrossRef]
  62. Mysinger, M.M.; Carchia, M.; Irwin, J.J.; Shoichet, B.K. Directory of Useful Decoys, Enhanced (DUD-E): Better Ligands and Decoys for Better Benchmarking. J. Med. Chem. 2012, 55, 6582–6594. [Google Scholar] [CrossRef]
  63. Empereur-Mot, C.; Zagury, J.F.; Montes, M. Screening Explorer–An Interactive Tool for the Analysis of Screening Results. J. Chem. Inf. Model. 2016, 56, 2281–2286. [Google Scholar] [CrossRef]
  64. Sterling, T.; Irwin, J.J. ZINC 15—Ligand Discovery for Everyone. J. Chem. Inf. Model. 2015, 55, 2324–2337. [Google Scholar] [CrossRef]
  65. Aminpour, M.; Montemagno, C.; Tuszynski, J.A. An Overview of Molecular Modeling for Drug Discovery with Specific Illustrative Examples of Applications. Molecules 2019, 24, 1693. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Durrant, J.D.; McCammon, J.A. Molecular dynamics simulations and drug discovery. BMC Biol. 2011, 9, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Liu, X.; Shi, D.; Zhou, S.; Liu, H.; Liu, H.; Yao, X. Molecular dynamics simulations and novel drug discovery. Expert Opin. Drug Discov. 2018, 13, 23–37. [Google Scholar]
  68. Case, D.A.; Cheatham, T.E.; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.; Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comput. Chem. 2005, 26, 1668–1688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [Green Version]
  70. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, Revision A.02; Published Online 2016; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  71. Wang, J.; Wolf, R.M.; Caldwell, J.W.; Kollman, P.A.; Case, D.A. Development and Testing of a General Amber Force Field. J. Comp. Chem. 2004, 25, 1157–1174. [Google Scholar] [CrossRef]
  72. Martins, F.G.; Melo, A.; Sousa, S.F. Identification of New Potential Inhibitors of Quorum Sensing through a Specialized Multi-Level Computational Approach. Molecules 2021, 26, 2600. [Google Scholar] [CrossRef]
  73. Vieira, T.F.; Martins, F.G.; Moreira, J.P.; Barbosa, T.; Sousa, S.F. In Silico Identification of Possible Inhibitors for Protein Kinase B (PknB) of Mycobacterium tuberculosis. Molecules 2021, 26, 6162. [Google Scholar] [CrossRef]
  74. Pereira, R.B.; Pinto, N.F.S.; Fernandes, M.J.G.; Vieira, T.F.; Rodrigues, A.R.O.; Pereira, D.M.; Sousa, S.F.; Castanheira, E.M.S.; Fortes, A.G.; Gonçalves, M.S.T. Amino Alcohols from Eugenol as Potential Semisynthetic Insecticides: Chemical, Biological, and Computational Insights. Molecules 2021, 26, 6616. [Google Scholar] [CrossRef]
  75. Natal, C.M.; Fernandes, M.J.G.; Pinto, N.F.S.; Pereira, R.B.; Vieira, T.F.; Rodrigues, A.R.O.; Pereira, D.M.; Sousa, S.F.; Fortes, A.G.; Castanheira, E.M.S.; et al. New carvacrol and thymol derivatives as potential insecticides: Synthesis, biological activity, computational studies and nanoencapsulation. RSC Adv. 2021, 11, 34024–34035. [Google Scholar] [CrossRef]
  76. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  77. Roe, D.R.; Cheatham, T.E. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef] [PubMed]
  78. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef] [PubMed]
  79. Miller, B.R.; McGee, T.D.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [PubMed]
  80. Sousa, S.F.; Coimbra, J.T.S.; Paramos, D.; Pinto, R.; Guimarães, R.S.; Teixeira, V.; Fernandes, P.A.; Ramos, M.J. Molecular dynamics analysis of a series of 22 potential farnesyltransferase substrates containing a CaaX-motif. J. Mol. Model. 2013, 19, 673–688. [Google Scholar] [CrossRef] [PubMed]
  81. Sousa, S.F.; Fernandes, P.A.; Ramos, M.J. Molecular dynamics simulations on the critical states of the farnesyltransferase enzyme. Bioorg. Med. Chem. 2009, 17, 3369–3378. [Google Scholar] [CrossRef] [PubMed]
  82. Martins, S.A.; Perez, M.A.S.; Moreira, I.S.; Sousa, S.F.; Ramos, M.J.; Fernandes, P.A. Computational Alanine Scanning Mutagenesis: MM-PBSA vs TI. J. Chem. Theory Comput. 2013, 9, 1311–1319. [Google Scholar] [CrossRef] [PubMed]
Antibiotics 11 00185 g001 550
Figure 1. MvfR pocket A and pocket B depiction. Protein bound to its natural inducer NNQ. (PDB:4JVI).
Figure 1. MvfR pocket A and pocket B depiction. Protein bound to its natural inducer NNQ. (PDB:4JVI).
Antibiotics 11 00185 g001
Antibiotics 11 00185 g002 550
Figure 2. Cartoon representation of the active site residues for 4JVI (cyan), 6B8A (yellow), NNQ ligand (white) and M64 (yellow). The main difference in these three structures is the position of residue Ile186 that in 6B8A structure tends to close the binding site, accommodating bulkier ligands.
Figure 2. Cartoon representation of the active site residues for 4JVI (cyan), 6B8A (yellow), NNQ ligand (white) and M64 (yellow). The main difference in these three structures is the position of residue Ile186 that in 6B8A structure tends to close the binding site, accommodating bulkier ligands.
Antibiotics 11 00185 g002
Antibiotics 11 00185 g003 550
Figure 3. ROC curves for 4JVI and 6B8A.
Figure 3. ROC curves for 4JVI and 6B8A.
Antibiotics 11 00185 g003
Antibiotics 11 00185 g004 550
Figure 4. 6B8A–M64 interaction map.
Figure 4. 6B8A–M64 interaction map.
Antibiotics 11 00185 g004
Antibiotics 11 00185 g005 550
Figure 5. Graph comparing the free-energy calculation results using MM/GBSA with the GOLD ASP scores. The reference ligands are NNQ, one of MvfR auto-inducer molecules, and M64, an antagonist. Blue-colored ligands indicate a higher LogP than those colored in pink. The size of the marker is a representation of the molecular weight.
Figure 5. Graph comparing the free-energy calculation results using MM/GBSA with the GOLD ASP scores. The reference ligands are NNQ, one of MvfR auto-inducer molecules, and M64, an antagonist. Blue-colored ligands indicate a higher LogP than those colored in pink. The size of the marker is a representation of the molecular weight.
Antibiotics 11 00185 g005
Antibiotics 11 00185 g006 550
Figure 6. Per-residue decomposition of the free-energy calculations using MM/GBSA for the fifteenth selected FDA compounds.
Figure 6. Per-residue decomposition of the free-energy calculations using MM/GBSA for the fifteenth selected FDA compounds.
Antibiotics 11 00185 g006
Antibiotics 11 00185 g007 550
Figure 7. Interaction maps for venetoclax, indocyanine green, and nilotinib.
Figure 7. Interaction maps for venetoclax, indocyanine green, and nilotinib.
Antibiotics 11 00185 g007
Antibiotics 11 00185 g008 550
Figure 8. Interaction maps of cabozantinib and montelukast.
Figure 8. Interaction maps of cabozantinib and montelukast.
Antibiotics 11 00185 g008
Antibiotics 11 00185 g009 550
Figure 9. Workflow of the current study.
Figure 9. Workflow of the current study.
Antibiotics 11 00185 g009
Antibiotics 11 00185 g010 550
Figure 10. Chemical structures of the six crystallographic ligands.
Figure 10. Chemical structures of the six crystallographic ligands.
Antibiotics 11 00185 g010
Table
Table 1. RMSD of the re-docking for the other targets and their respective crystallographic ligand.
Table 1. RMSD of the re-docking for the other targets and their respective crystallographic ligand.
Redocking RMSD (Å)
PDB CodeLigandVinaLeDockCHEMPLPGoldScoreChemScoreASPAverage per Target
4JVDNNQ6.673.512.163.201.182.683.23
4JVIQZN1.593.071.332.933.181.852.33
6B8AM640.340.580.460.633.481.921.24
6Q7UHLH7.265.803.713.142.642.114.11
6Q7VC115.776.163.495.733.761.774.98
6Q7WC203.544.501.995.153.171.713.34
6TPRNV59.155.231.541.164.591.353.84
6Z07Q4E1.221.410.920.871.311.321.18
6Z17Q4W1.187.251.594.082.321.963.06
6Z5KQAE1.434.461.178.133.951.603.46
6YZ3Q251.117.671.250.993.461.562.67
Average by SF 3.574.511.783.273.001.81
Table
Table 2. Available X-ray structures of MvfR on the PDB and BSD.
Table 2. Available X-ray structures of MvfR on the PDB and BSD.
4JVI6B8A
EF 1%AUC%TGRIEBEDROCEF 1%AUCTGRIEBEDROC
Vina10.230.680.212.340.1410.350.630.192.130.13
CHEMPLP10.400.550.082.220.135.200.540.072.070.12
ChemScore5.200.490.0051.680.102.600.520.021.650.10
ASP10.400.660.212.390.1410.390.660.252.330.14
Table
Table 3. Top 15 hits of the FDA-approved drugs database.
Table 3. Top 15 hits of the FDA-approved drugs database.
Drug Name & CodeDescriptionStructureASP Score
NilotinibBcr-Abl tyrosine kinase inhibitor (TKI) used in the treatment of chronic myelogenous leukemia (CML) Antibiotics 11 00185 i00154.96
Indocyanine GreenDye used in medical diagnosis. It has been used to measure cardiac output, liver function, and in ophthalmic angiography [33]. Antibiotics 11 00185 i00250.55
LomitapideUsed to treat patients with Homozygous familial hypercholesterolaemia (HoFH). It is an inhibitor of MTP, an enzyme responsible for the synthesis of low-density lipoproteins in the liver [34]. Antibiotics 11 00185 i00350.01
ValrubicinChemotherapy drug used to treat carcinoma in situ bladder tumors. Antibiotics 11 00185 i00449.89
LapatinibInhibitor of tyrosine kinase domains of epidermal growth factor receptor and human epidermal growth factor receptor (HER)-2. Used to treat metastatic HER-2 + breast cancer [35]. Antibiotics 11 00185 i00549.84
PazopanibMultitarget tyrosine kinase inhibitor approved for the treatment of multiple histological subtypes of soft tissue sarcoma (STS) [36]. Antibiotics 11 00185 i00649.69
RavictiUsed for the treatment of patients with urea cycle disorders (UCDs) [37]. Antibiotics 11 00185 i00748.49
CabozantinibTyrosine kinase inhibitor that targets pathways that have been linked to tumor growth. Used for the treatment of metastatic renal cell carcinoma [38]. Antibiotics 11 00185 i00848.32
VenetoclaxInhibitor of B-cell leukemia/lymphoma-2 protein [39]. Antibiotics 11 00185 i00948.05
IsavuconazoniumProdrug used as antifungal for the treatment of invasive aspergillosis and invasive mucormycosis [40]. Antibiotics 11 00185 i01048.05
EmendNK1 antagonist to prevent chemotherapy-induced nausea and vomiting. Antibiotics 11 00185 i01146.77
Amethopterin or MethotrexateAnalog and antagonist of folic acid, is commonly used in the treatment of a wide range of malignant and non-malignant diseases [41]. Antibiotics 11 00185 i01246.10
CefsulodinBroad-spectrum beta-lactamase stable cephalosporin with excellent activity against gram-negative bacilli, including P. aeruginosa [42]. Antibiotics 11 00185 i01346.01
MontelukastLeukotriene receptor antagonist (LTRAs) used for asthma treatment [43]. Antibiotics 11 00185 i01445.27
CefoperazoneIs a parenteral, third-generation cephalosporin that can be given intravenously or intramuscularly [44]. Antibiotics 11 00185 i01545.05
Table
Table 4. Average protein RMSD values (Å), average ligand RMSD (Å), average MvfR–ligand complex SASA (Å2), percentage of SASA for the buried ligand, and average number of ligand hydrogen bonds, obtained from the MD simulations of MvfR–FDA complexes. ΔG binding energy determined using MM/GBSA and per-residue decomposition, calculated for the last 90 ns of the simulation. Compounds are ordered by increasing ΔGbind (kcal/mol).
Table 4. Average protein RMSD values (Å), average ligand RMSD (Å), average MvfR–ligand complex SASA (Å2), percentage of SASA for the buried ligand, and average number of ligand hydrogen bonds, obtained from the MD simulations of MvfR–FDA complexes. ΔG binding energy determined using MM/GBSA and per-residue decomposition, calculated for the last 90 ns of the simulation. Compounds are ordered by increasing ΔGbind (kcal/mol).
IDAverage Protein RMSd (Å)Average Ligand RMSDSASA (Å2)Percentage of Potential Ligand SASA Buried (%)Average Number H BondsΔGbind (kcal/mol)Main Contributors (kcal/mol)
Venetoclax2.3 ± 0.22.7 ± 0.4344.8 ± 27.070.4 ± 0.020.5 ± 0.6−70.1 ± 0.3TYR258 (−4.9), ILE186 (−4.3),
ILE236 (−2.9)
Indocyanine Green2.2 ± 0.22.3 ± 0.2298.9 ± 28.372.3 ± 0.031.1 ± 0.8−58.6 ± 0.3ILE186 (−4.2), ILE236 (−3.1),
TYR258 (−3.1)
Nilotinib2.0 ± 0.21.7 ± 0.4105.2 ± 23.286.8 ± 0.030.1 ± 0.3−48.1 ± 0.2ILE186 (−2.4), LEU208 (−2.1), ILE236 (−3.4)
Cabozantinib2.3 ± 0.31.6 ± 0.4135.3 ± 40.582.3 ± 0.10.03 ± 0.2−44.6 ± 0.2ILE236 (−3.2), ILE186 (−2.4), LEU208 (−2.3)
Montelukast2.2 ± 0.22.2 ± 0.4228.8 ± 35.673.8 ± 0.040.04 ± 0.2−43.2 ± 0.2ILE236 (−2.8), LEU207 (−2.4), ILE263 (−2.0)
Cefoperazone2.6 ± 0,32.9 ± 1.2289.5 ± 41.466.7 ± 0.041.0 ± 0.9−42.4 ± 0.4ARG209 (−3.2), ILE236 (−2.5), LEU208 (−2.3)
Valrubicin2.4 ± 0.43.3 ± 0.3276.7 ± 31.470.6 ± 0.030.04 ± 0.2−41.1 ± 0.2LEU207 (−4.2), ILE236 (−2.9), TYR258 (−3.9)
Lomitapide2.3 ± 0.33.8 ± 0.4323.4 ± 81.665.6 ± 0.10.1 ± 0.3−40.6 ± 0.3ILE186 (−2.3), LEU208 (−1.9), ILE236 (−2.5)
M64 (antagonist)2.2 ± 0.21.2 ± 0.290.4 ± 19.386.3 ± 0.030.1 ± 0.2−39.0 ± 0.1IlE236 (−3.2), ILE186 (−1.7), TYR258 (−1.7)
Emend2.3 ± 0.41.5 ± 0.2130.2 ± 25.880.7 ± 0.040.4 ± 0.5−38.9 ± 0.2ILE236 (−3.5), LEU207 (−109), TYR258 (−1.7)
Pazopanib2.3 ± 0.21.9 ± 0.5203.6 ± 40.569.4 ± 0.10.7 ± 0.8−37.7 ± 0.3LEU208 (−3.7), ILE236 (−3.1), SER196 (−2.0)
lapatinib2.5 ± 0.42.8 ± 0.9249.0 ± 56.769.6 ± 0.10.5 ± 0.7−35.4 ± 0.3LEU207 (−2.1), ILE236 (−2.9), ILE263 (−1.8)
Ravicti2.1 ± 0.23.3 ± 0.7263.8 ± 50.069.3 ± 0.10.03 ± 0.2−34.5 ± 0.3LEU207 (−2.6), ILE236 (−2.6), ILE263 (−1.7)
Cefsulodin2.3 ± 0.21.6 ± 0.4235.4 ± 39.966.5 ± 0.10.9 ± 0.8−27.4 ± 0.3LEU207 (−3.5), ILE236 (−2.2), LEU208 (−2.1)
NNQ (natural inducer)2.3 ± 0.31.6 ± 0.4154.8 ± 53.771.6 ± 0.10.1 ± 0.3−26.1 ± 0.3LEU208 (−1.7), ILE236 (−1.5)
Isavuconazonium2.4 ± 0.22.7 ± 0.5522.5 ± 47.442.7 ± 0.10.04 ± 0.2−25.0 ± 0.1TYR258 (−3.4), ILE186 (−2.6), LEU189 (−1.1)
Methotrexate2.8 ± 0.42.4 ± 0.4312.6 ± 69.253.8 ± 0.10.7 ± 0.8−22.8 ± 0.3ILE186 (−4.1), TYR258 (−3.4), ARG209 (−2.7)
Table
Table 5. Available X-ray structures of MvfR on the PDB and BSD.
Table 5. Available X-ray structures of MvfR on the PDB and BSD.
PDB CodeProteinResolution (Å)LigandStrainReferences
4JVCLigand-Binding Domain2.50 UCBPP-PA14[20]
4JVDLigand-Binding Domain2.95NNQ
4JVILigand-Binding Domain2.90QZN
6B8ALigand-Binding Domain2.65M64PAO1[23]
6Q7ULigand-Binding Domain3.15HLHPAO1[22]
6Q7VLigand-Binding Domain2.56HLK
6Q7WLigand-Binding Domain2.82HLQ
6TPRLigand-Binding Domain3.20NV5UCBPP-PA14[47]
6Z07Ligand-Binding Domain2.95Q4EUCBPP-PA14[48]
6Z17Ligand-Binding Domain3.15Q4W
6Z5KLigand-Binding Domain3.20QAE
6YZ3Ligand-Binding Domain3.00Q25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vieira, T.F.; Magalhães, R.P.; Simões, M.; Sousa, S.F. Drug Repurposing Targeting Pseudomonas aeruginosa MvfR Using Docking, Virtual Screening, Molecular Dynamics, and Free-Energy Calculations. Antibiotics 2022, 11, 185. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics11020185

AMA Style

Vieira TF, Magalhães RP, Simões M, Sousa SF. Drug Repurposing Targeting Pseudomonas aeruginosa MvfR Using Docking, Virtual Screening, Molecular Dynamics, and Free-Energy Calculations. Antibiotics. 2022; 11(2):185. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics11020185

Chicago/Turabian Style

Vieira, Tatiana F., Rita P. Magalhães, Manuel Simões, and Sérgio F. Sousa. 2022. "Drug Repurposing Targeting Pseudomonas aeruginosa MvfR Using Docking, Virtual Screening, Molecular Dynamics, and Free-Energy Calculations" Antibiotics 11, no. 2: 185. https://0-doi-org.brum.beds.ac.uk/10.3390/antibiotics11020185

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