Methicillin-Resistant <i>Staphylococcus aureus</i>: Docking-Based Virtual Screening and Molecular Dynamics Simulations to Identify Potential Penicillin-Binding Protein 2a Inhibitors from Natural Flavonoids

Masumi, Motahareh; Noormohammadi, Fatemeh; Kianisaba, Fatemeh; Nouri, Fatemeh; Taheri, Mohammad; Taherkhani, Amir

doi:https://doi.org/10.1155/2022/9130700

International Journal of Microbiology

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Abbreviations Data Availability Ethical Approval Consent Conflicts of Interest Authors’ Contributions Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2022 | Article ID 9130700 | https://doi.org/10.1155/2022/9130700

Methicillin-Resistant Staphylococcus aureus: Docking-Based Virtual Screening and Molecular Dynamics Simulations to Identify Potential Penicillin-Binding Protein 2a Inhibitors from Natural Flavonoids

Motahareh Masumi,¹Fatemeh Noormohammadi,¹Fatemeh Kianisaba,¹Fatemeh Nouri,²Mohammad Taheri,³and Amir Taherkhani⁴

Academic Editor: Rossella Grande

Received02 Nov 2021

Revised04 Mar 2022

Accepted15 Apr 2022

Published04 May 2022

Abstract

Staphylococcus aureus (S. aureus) is responsible for several disorders including skin and soft tissue infections, bacteremia, pulmonary infections, septic arthritis, osteomyelitis, meningitis, gastroenteritis, toxic-shock syndrome, and urinary tract infections. Methicillin-resistant S. aureus (MRSA) contains penicillin-binding protein 2a (SauPBP2a) responsible for catalyzing the peptidoglycan production within the bacterial cell wall. The binding affinity of SauPBP2a to beta-lactam antibiotics is low, and thus, it is necessary to discover new effective SauPBP2a inhibitors to combat mortality and morbidity in patients affected by MRSA. The binding affinity of 46 natural flavonoids to the SauPBP2a active site was examined via molecular docking analysis. The stability of docked poses associated with the top-ranked flavonoids was tested by performing molecular dynamics (MD) in 10 nanoseconds (ns) computer simulations. Kaempferol 3-rutinoside-7-sophoroside and rutin demonstrated a considerable binding affinity to the SauPBP2a active site (Δ_binding < 11 kcal/mol). Their docked poses were found to be stable for 10 ns MD simulations. Kaempferol 3-rutinoside-7-sophoroside and rutin also exhibited salient binding affinity to the enzyme’s allosteric site. This study suggests that kaempferol 3-rutinoside-7-sophoroside and rutin may be considered as drug candidates for therapeutic aims in several human infections associated with MRSA. Nevertheless, in vitro and in vivo confirmations are warranted.

1. Introduction

Staphylococcus aureus (S. aureus), a Gram-positive bacteria, is the leading cause of severe invasive and toxin-mediated disorders including skin and soft tissue infections, bacteremia, pulmonary infections, septic arthritis, osteomyelitis, meningitis, gastroenteritis, toxic-shock syndrome, and urinary tract infections. Methicillin-resistant Staphylococcus aureus (MRSA) refers to any strains of S. aureus that carry the mecA gene encoding penicillin-binding protein 2a (PBP2a), responsible for catalyzing the peptidoglycan production within the bacterial cell wall. As PBP2a has a lower binding affinity to beta-lactam-containing antibiotics compared to other PBPs, MRSA continues to catalyze bacterial cell wall synthesis in the presence of penicillin-derived antibiotics such as methicillin, oxacillin, nafcillin, and cephalosporins [1–7]. Thus, PBP2a in S. aureus (SauPBP2a) has been considered a potential target for developing specific inhibitors to combat MRSA, which will be accelerated by analyzing the three-dimensional structure of the SauPBP2a active site [8]. Furthermore, the allosteric site of SauPBP2a has also been considered for negative regulation of the enzyme’s activity [9].

Organic compounds with antibacterial properties have sparked public interest in recent years for treating infections [10]. Flavonoids are one of the most important active compounds that naturally occur in photosynthesizing plants with several valuable properties in humans [11, 12]. Fruits and vegetables containing flavonoids have been widely used to treat human disorders [13]. Several flavonoids have been shown to possess antibacterial properties in previous research. In this regard, chrysin has been demonstrated to have bacteriostatic activity against Escherichia coli (E. coli) and Pseudomonas aeruginosa (P. aeruginosa).

Furthermore, baicalein has an inhibitory effect on the growth of S. aureus and Bacillus subtilis. Also, vitexin, saponarin, lucenin, apigenin, and luteolin have revealed bactericidal effects against E. coli, Proteus mirabilis, Enterobacter cloacae, Klebsiella pneumoniae, Proteus vulgaris, and P. aeruginosa [14]. In addition, kaempferol, myricetin, datiscetin, luteolin, and quercetin are effective against MRSA [15]. As a result, flavonoids have emerged as a primary research focus in the field of antibacterial medication discovery [13].

In the present study, we hypothesized that inhibition of PBP2a is a possible mechanism flavonoids use to combat MRSA growth. Thus, a docking-based virtual screening study was executed to examine the binding affinity of several natural flavonoids to the SauPBP2a active site. Molecular dynamics (MD) simulations were conducted to evaluate the stability of docked poses of top-ranked SauPBP2a active site inhibitors. The binding affinity of salient SauPBP2a active site inhibitors to the enzyme’s allosteric site was also tested.

2. Materials and Methods

2.1. SauPBP2a and Preparation of Flavonoids’ Structures

The three-dimensional coordinate of SauPBP2a was obtained from the Protein Data Bank (PDB) database, which is available at https://www.rcsb.org. Two PDB codes including 1MWT (X-ray resolution = 2.45 Å [8, 16]) and 4CJN (X-ray resolution = 1.95 Å [17]) were selected for docking analyses with the active and allosteric sites of the enzyme, respectively. The 1MWT and 4CJN files included two polypeptide chains as A and B. Chain B was selected for further analysis, which contained 646 and 642 residues in the 1MWT and 4CJN files, respectively. The chain B also included penicillin G (PDB ID: PNM) and (E)-3-(2-(4-cyanostyryl)-4-oxoquinazolin-3(4H)-yl) benzoic acid (PDB ID: QNZ) in the 1MWT and 4CJN PDB files, respectively; these ligands were removed from the SauPBP2a structure before molecular docking and dynamics studies.

Energy minimizing (EM) of SauPBP2a with two different codes was carried out before docking operations using the Swiss-PDB Viewer version 4.1.0, which is available at https://spdbv.unil.ch/ [18]. A total of 46 natural flavonoids were docked with the active site of SauPBP2a to identify potential inhibitors of the enzyme’s catalytic domain. Next, molecular docking analysis was done using the known PBP inhibitors, including penicillin G and methicillin. Due to the low binding affinity of SauPBP2a to the beta-lactam-containing antibiotics, penicillin G and methicillin were considered control negative in this study. Furthermore, previous studies have reported that “oxadiazole” [19] and “ceftobiprole” [20] had exhibited high binding affinity to the SauPBP2a active site and hence were considered control positive in our work.

The two-dimensional (2D) coordinates of flavonoids, methicillin, and ceftobiprole were principally collected from the PubChem database (https://pubchem.ncbi.nlm.nih.gov) as a structure data file (SDF). Subsequently, the Cactus web server was used to achieve PDB formats from the SDF files, available at https://cactus.nci.nih.gov/translate/. The three-dimensional structure of penicillin G was collected from the 1MWT file, while the construction of oxadiazole was drawn by the ChemDraw software version 12.0.2.1076. The EM was administered for all small molecules using the HyperChem software version 8.0.10 [21].

2.2. Molecular Docking with the SauPBP2a Active Site

A windows-based program system was used with the following features: system type, 64 bit; processor, Intel Core i7; installed memory, 32 GB. Docking-based virtual screening was executed using AutoDock 4.0 (https://autodock.scripps.edu) [22]. The AutoDock software applies limited flexibility in the receptor and uses the Lamarckian genetic algorithm to predict the docked pose of the ligand [23]. The SauPBP2a catalytic site was considered the receptor for executing docking analyses with flavonoids, control negative, and control positive compounds. After analyzing the interactions between nitrocefin, penicillin G, methicillin, and the residues within the SauPBP2a catalytic domain from the study by Lim and Strynadka [8], a total of 14 main residues were identified inside the PBP2a active site including Ser337, Lys340, Ser403, Lys406, Tyr446, Ser462, Asn464, Thr500, Ser548, Gly549, Ser598, Gly599, Thr600, and Met681. Thus, these amino acids were considered as a receptor site for docking analysis to identify potential inhibitors of the enzyme’s active site. The grid box settings were as follows: X-dimension, 60; Y-dimension, 66; Z-dimension, 60; X-center, −34.264; Y-center, 44.65; Z-center, 66.778; spacing, 0.375 Å. The number of docked poses for each component was set at 50, where the docked model with the lowest ∆_binding involved in the largest cluster along with the largest number of models using a root mean square deviation (RMSD), tolerance of 2.0 Å, was considered as a binding affinity of the ligand to the receptor site.

2.3. Cross-Validation Study

Cross-validation of docking studies was executed for top-ranked SauPBP2a active site inhibitors based on the Δ_binding values achieved from the AutoDock tool (Δ_binding < 10 kcal/mol). This was performed using the Schrödinger Maestro docking software version 10.2 [24, 25]. In this regard, docking scores were calculated, and the prime MM-GBSA method was used to determine the relative binding energies of ligands. The Glide docking system was employed to evaluate the binding affinity of compounds to the SauPBP2a active site. The lowest Glide score (dock score) for each ligand was considered the best-docked model. More details regarding settings in the Schrödinger Maestro docking software have been reported in a previous study by Azadian et al. [26].

2.4. MD Simulations

To examine the stability of docked poses among top-ranked flavonoids and the SauPBP2a active sites, MD was simulated by applying the Discovery Studio Client software (version 16.1.0.15350). Only ligands that demonstrated considerable affinity to the enzyme’s active site using the AutoDock tool and Schrödinger docking software were considered for MD analysis. The advanced settings for computer simulations were as follows: MD time: 10 nanoseconds (10 ns), solvation model: explicit periodic boundary; cell shape: orthorhombic; minimum distance from the boundary: 10 Å; solvent: water; target temperature: 310 K; force field: CHARMm; charge distribution: point. In addition, the time evolution of root mean square deviation (RMSD) of backbone atoms and root mean square fluctuation (RMSF) of SauPBP2a complexed with the best inhibitors was also analyzed. The BIOVIA Discovery Studio Visualizer version 19.1.0.18287 was then utilized to demonstrate the interactions among top-ranked flavonoids and the residues inside the SauPBP2a catalytic site and to superimpose protein-ligand complexes before and after MD simulations.

2.5. Molecular Docking with the SauPBP2a Allosteric Site

This study considered the best SauPBP2a active site inhibitors to evaluate their binding affinity to the enzyme’s allosteric site. The grid box options at this stage followed the settings of the study by Ibrahim et al. [9]: X-dimension: 52; Y-dimension: 52; Z-dimension: 52; X-center: 9.658; Y-center: −1.662; Z-center: −70.269; spacing: 0.375 Å.

3. Results

3.1. Molecular Docking with the SauPBP2a Active Site

It was found that 9 of the flavonoids can potentially bind to the SauPBP2a catalytic site on the nanomolar scale (nM). A total of 37 compounds revealed an inhibitory effect on SauPBP2a active site at the micromolar (uM) scale. Four components, including kaempferol 3-rutinoside-7-sophoroside, rutin, amentoflavone, and quercetin, bounded to the SauPBP2a active site with the salient affinity of ∆_binding < 10 kcal/mol (Supplementary Table 1). Thus, these compounds found the top SauPBP2a active site inhibitors among the tested flavonoids based on the AutoDock tool and were considered for further analyses using the Schrödinger Maestro docking software. Ceftobiprole and oxadiazole exhibited a mild binding affinity to the enzyme’s active site with the ∆_binding values of 8.79 and 7.87 kcal/mol, respectively. Figure 1 shows the two-dimensional structures of kaempferol 3-rutinoside-7-sophoroside, rutin, amentoflavone, quercetin, and control negative and positive compounds in this study, which were achieved by the ChemDraw Ultra version 12.0.2.1076. Of note, the ∆_binding between penicillin G as well as methicillin and SauPBP2a active site was calculated to be 6.35 and 5.51 kcal/mol, respectively. Figure 2 shows the ∆_binding values between top-ranked flavonoids, control inhibitors, and the SauPBP2a active site. In comparison, the details of energies between top-ranked flavonoids, control inhibitors, and the enzyme’s catalytic domain are given in Table 1. All hydrogen, hydrophobic, electrostatic, and miscellaneous interactions between kaempferol 3-rutinoside-7-sophoroside, rutin, amentoflavone, quercetin, and residues inside the SauPBP2a catalytic site were analyzed before and after MD simulations (Table 2 and Figure 3). In this regard, kaempferol 3-rutinoside-7-sophoroside revealed the most H-bonds before and after MD analyses. The H-bonds with a distance of >5 Å were not considered significant interactions. Furthermore, Figure 4 shows the interaction modes between control compounds and the residues inside the enzyme’s active site. All possible interactions between kaempferol 3-rutinoside-7-sophoroside, rutin, amentoflavone, quercetin, and the residues within the SauPBP2a catalytic site are shown in Figure 5(a) as a single graph. This network was achieved using the Cytoscape software (3.8.0; https://www.cytoscape.org) [27]. Subsequently, the degree of each residue was calculated using the Cytoscape network analyzer tool, illustrating the number of interactions between each residue and the top-ranked flavonoids. Lys430 indicated the highest degree between the residues (Figure 5(b)).

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(a)

(b)

(c)

(d)

(e)

(f)

(g)

(h)

(a)

(b)

(c)

(d)

(a)

(b)

3.2. Cross-Validation and MM-GBSA Analyses

Schrödinger Maestro docking scores were calculated for kaempferol 3-rutinoside-7-sophoroside, rutin, amentoflavone, and quercetin. Kaempferol 3-rutinoside-7-sophoroside revealed the highest binding affinity to the active site of SauPBP2a, followed by rutin with the dock scores of 12.697 and 9.341 kcal/mol, respectively (Table 3). Thus, these two compounds were considered salient inhibitors of the SauPBP2a active site and were selected for further MD analysis. Meanwhile, MD simulations were performed based on the complexes achieved from the AutoDock software. Furthermore, the interaction modes between kaempferol 3-rutinoside-7-sophoroside, rutin, and residues within the sauPBP2a active site were also studied based on the complexes achieved from the Schrödinger Maestro docking tool (Table 2, Figures 3(g) and 3(h)). The prime MM-GBSA analysis showed each compound’s relative binding-free energy (Δ bind), with the results given in Table 4. The expanded formula is given as follows [28]: Δ (bind) = Δ (solv) + ΔE (MM) + Δ (SA).

3.3. MD Simulations

MD simulation was done to assess the stability of docked poses of kaempferol 3-rutinoside-7-sophoroside and rutin. Figure 6 shows the superimposed structures of SauPBP2a complexed with kaempferol 3-rutinoside-7-sophoroside and rutin before and after MD simulations. Furthermore, the RMSD of backbone atoms and RMSF for SauPBP2a complexed with kaempferol 3-rutinoside-7-sophoroside and rutin are shown in Figure 7. As shown in Figure 7, it seems that the helix structures are more stable than the other secondary structures in SauPBP2a.

(a)

(b)

(c)

(d)

(a)

(b)

3.4. Molecular Docking with the SauPBP2a Allosteric Site

The results showed that the potential inhibitors of SauPBP2a active site including kaempferol 3-rutinoside-7-sophoroside and rutin also had an excellent binding affinity to the allosteric site of the enzyme with the criteria of Δ_binding values as 14.35 and 9.68 kcal/mol, respectively (Table 1).

4. Discussion

Antibiotic resistance has emerged in response to the overuse of currently available antibiotics [10]. Thus, it is necessary to discover and design novel antibacterial compounds with new formulations to overcome serious infections [29]. MRSA infection is one of the most frequent causes of hospital-acquired conditions and is currently associated with poor prognosis and increased mortality/morbidity [30]. In the present study, molecular docking and dynamics simulations were performed to identify potential inhibitors of the SauPBP2a active site. However, targeting the PBP2a allosteric site is another alternative strategy to reduce SauPBP2a activity. As such, we also evaluated the binding affinity of our top-ranked PBP2a active site inhibitors including kaempferol 3-rutinoside-7-sophoroside and rutin to the allosteric site of PBP2a. We used the same protein Ibrahim et al. [9] used in their study (PDB ID: 4CJN; chain B). Interestingly, kaempferol 3-rutinoside-7-sophoroside and rutin exhibited high binding affinity to the PBP2a allosteric site with the Δ_binding values of 14.35 and 9.68 kcal/mol, respectively.

Kaempferol is a flavonoid compound with various health benefits such as anti-inflammatory, anticancer, liver-protective, antiobesity, antidiabetes, and heart-protective effects. It is found in fruits and vegetables, tea, beans, cabbage, broccoli, barriers, grapes, apples, and citrus fruits [31]. Previous research has shown that kaempferol and its many derivatives have antimicrobial properties. Qui et al. [32] demonstrated that two of the kaempferol derivatives named 3-O-[β-d-glucopyranosyl-(1 ⟶ 2)-α-l-rhamnopyranosyl-(1 ⟶ 6)]-β-d-glucopyranoside and kaempferol 3-O-[β-d-xylopyranosyl-(1 ⟶ 2)-α-l-rhamnopyranosyl-(1 ⟶ 6)]-β-d-glucopyranoside revealed considerable antibacterial actions against S. aureus, E. coli, Salmonella enteriditis, Bacillus thuringiensis, Aspergillus niger, and Rhizopus nigricans. These two compounds were extracted from Camellia oleifera using the continuous phase change extraction approach. Furthermore, Cruz et al. [10] demonstrated that the combination of kaempferol 7-O-β-D-(6″-O-cumaroyl)-glucopyranoside with the amikacin and gentamicin at a concentration of 128 μg/mL led to synergistic actions against S. aureus and E. coli and significantly diminished the minimum inhibitory concentration (MIC) from 16 to 4 and 88 μg/mL, respectively. One year later, Kannanoor et al. [33] reported that the silver nanoparticles (AgNPs) conjugated with kaempferol (K-AgNPs) demonstrated a considerable bactericidal activity against various bacterial strains including E. coli, Bacillus subtilis, and S. aureus. They used several experimental approaches, including live/dead bacterial ratio analysis, lactate dehydrogenase, and lipid peroxidation analyses. According to our results, it was estimated that kaempferol has a mild binding affinity to the catalytic site of SauPBP2a (∆_binding = 6.12 kcal/mol; Ki = 32.86 µM), while two of the kaempferol derivatives including kaempferol 3-rutinoside-7-sophoroside and kaempferol 3-rutinoside-4′-glucoside were found to bind to the SauPBP2a active site at the Ki values of 1.42 nM and 562.45 nM, respectively. This suggests that the sugar moiety has enhanced the binding affinity of kaempferol to SauPBP2a active site. Using the AutoDock and Schrödinger Maestro docking tools, kaempferol 3-rutinoside-7-sophoroside revealed a high binding affinity to the SauPBPa2 active site with the ∆_binding of 12.07 and 12.697 kcal/mol, respectively. Before MD, kaempferol 3-rutinoside-7-sophoroside indicated six hydrogen, one electrostatic, and one hydrophobic interaction with the Lys430, Thy444, Val448, Thr600, Glu602, and Ser643 inside the catalytic domain of SauPBPa2. It formed 12 hydrogen and 3 hydrophobic interactions with Ser403, Val443, Tyr444, Tyr446, Glu447, Val448, Asn464, Gln521, Gly640, and Ser643 within the SauPBP2a active site after 10 ns MD simulation.

Rutin is an organic flavonol with several pharmacological features such as antioxidant and antibacterial properties [34–37]. In a previous report, a combination of rutin and florfenicol exhibited enhanced antibacterial activity against multidrug-resistant (MDR) Aeromonas hydrophila, suggesting that the use of organic compounds in combination with traditional antibiotics may result in synergistic inhibition of the growth of MDR bacteria [29]. Recently, Rodríguez-Valdovinos and Salgado-Garciglia [38] designed a study to examine the antibacterial and antioxidant effects of Verbesina sphaerocephala (V. sphaerocephala) leaf and flower extracts, in which rutin is the main phenolic component. Rodríguez-Valdovinos and Salgado-Garciglia [38] reported that V. sphaerocephala extracts demonstrated significant antibacterial effects on S. aureus and E. coli. Our AutoDock results found that rutin can bind to the SauPBPa2 catalytic domain at the nanomolar concentration (Ki = 436.10 nM) with an estimated Gibbs free energy of 11.14 kcal/mol. Furthermore, the Schrödinger Maestro docking software estimated that the dock score between rutin and SauPBP2a active site was 9.341 kcal/mol. Before MD simulations, rutin formed 3 hydrogen and 1 hydrophobic interaction with Lys430, Tyr446, and Thr600 within the SauPBPa2 active site. This flavonoid showed 2 hydrogen and 3 hydrophobic interactions with Tyr444, Arg445, and Ser643 inside the SauPBP2a active site after 10 ns MD simulation.

Amentoflavone is a bioflavonoid metabolite with miscellaneous beneficial properties including anti-inflammatory, neuroprotective, antiviral, and anticancer effects [39–42]. It has been widely purified from Hypericum perforatum, Ginkgo biloba, and Selaginella tamariscina [43]. Hwang et al. [44] performed a study to evaluate the antibacterial effects of amentoflavone alone, isolated from Selaginella tamariscina, and in combination with other antibiotics including chloramphenicol, cefotaxime, and ampicillin. The authors reported synergistic results of combined amentoflavone with other traditional antibiotics on various Gram-positive and Gram-negative bacteria including S. aureus, E. coli, Enterococcus faecium, and Pseudomonas aeruginosa (P. aeruginosa). This was carried out by applying the Clinical and Laboratory Standards Institute guideline [45]. The authors suggested that the enhanced production of hydroxyl radicals by amentoflavone is potentially the main reason for the synergistic outcome. Our results found that amentoflavone can bind to the SauPBPa2 active site at the nanomolar concentration (Ki = 31.63 nM) and with an estimated free binding energy of 10.23 kcal/mol. Our docking results indicated that amentoflavone revealed 3 hydrogen, 4 hydrophobic, 2 electrostatic, and 2 miscellaneous interactions with Lys430, Thr444, Tyr446, Val448, Asn464, Thr600, Glu602, and Met641 within the SauPBPa2 active site.

Quercetin is predominantly found in vegetables and fruits such as apples, broccoli, berries, onions, kales, tee, and red grape. Quercetin has several valuable characteristics, including antioxidant, anti-inflammatory, anticarcinogenic, and antiviral properties. It also reduces the risk of infectious diseases and promotes mitochondrial biogenesis [46–51]. Wang et al. [52] performed in vivo and in vitro studies to examine the antibacterial effects of quercetin on cecal microbiota of Arbor Acre broiler chickens, as well as its mechanism of action. Wang et al. [52] demonstrated that quercetin significantly reduced the total amount of S. aureus (), P. aeruginosa (), Salmonella enterica serotype Typhimurium (), and E. coli () compared to the negative control. Furthermore, the authors provided the in vitro clue that quercetin damages the membrane and cell walls of S. aureus and E. coli, indicating that it could be used as an alternative antibacterial feed for animals. Based on our findings, quercetin revealed 1 hydrogen and 6 hydrophobic interactions with Lys430, Val443, Tyr446, Tyr519, and Met641 within the PBPa2 active site.

The stability of docked poses of kaempferol 3-rutinoside-7-sophoroside and rutin was studied by performing MD in 10 ns computer simulations. It was suggested that the docked poses of these 2 flavonoids were stable in 10 ns simulation. In addition, kaempferol 3-rutinoside-7-sophoroside and rutin also revealed considerable binding affinity to the allosteric site of the enzyme. Thus, these compounds may be considered effective inhibitors of SauPBP2a, though this needs in vitro and in vivo confirmations in future studies.

The present study had a principal limitation. MD is a computer simulation approach analyzing the physical movements of atoms and molecules in an environment similar to reality. Thus, a supercomputer with a powerful processor and random access memory (RAM) is an indispensable tool for executing these types of time-consuming procedures. Unfortunately, the supercomputer was unavailable for our research team to run MD simulations for more significant periods. Each MD simulation (10 ns) took about one month with the computer we had. We strongly believe that simulating MD for a greater period of times results in more reliable results.

5. Conclusion

46 flavonoids were docked to the SauPBP2a active site using the AutoDock tool to estimate the binding affinity of the selected herbal flavonoids to the target protein. A cross-validation study was performed for compounds with the criteria of Δ_binding < 10 kcal/mol. Next, MD simulation was executed for compounds that exhibited considerable binding affinity to the enzyme’s active site using AutoDock and Schrödinger Maestro docking tools. The results suggest that kaempferol 3-rutinoside-7-sophoroside and rutin have a high affinity of binding to the SauPBP2a active site. Their docked poses were stable in the 10 ns simulation. They also exhibited excellent binding affinity to the enzyme’s allosteric site. Thus, it may be speculated that kaempferol 3-rutinoside-7-sophoroside and rutin could be helpful in the therapeutic aims of MRSA infections. Nevertheless, further in vivo and in vitro inhibition experiments are warranted in the future.

Abbreviations

S. aureus:	Staphylococcus aureus
MRSA:	Methicillin-resistant Staphylococcus aureus
PBP2a:	Penicillin-binding protein 2a
SauPBP2a:	PBP2a in S. aureus
E. coli:	Escherichia coli
P. aeruginosa:	Pseudomonas aeruginosa
PDB:	Protein Data Bank
EM:	Energy minimizing
SDF:	Structure data file
RMSD:	Root mean square deviation
MD:	Molecular dynamics
nM:	Nanomolar
MIC:	Minimum inhibitory concentration
AgNPs:	Silver nanoparticles
MDR:	Multidrug-resistant
V. sphaerocephala:	Verbesina sphaerocephala
Ki:	Inhibition constant.

Data Availability

The datasets used and/or analyzed during the current study are available from the corresponding author upon request.

Ethical Approval

The current study was approved by the Ethics Committee of Hamadan University of Medical Sciences, Hamadan, Iran (ethics no. IR.UMSHA.REC.1400.345).

Not applicable.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

AT designed the study, performed MD simulations, conducted and interpreted network analysis, and wrote the manuscript. MM, FNoor, FK, and AT conducted docking operations. AT, FNouri, and MT analyzed and discussed the results. MT and FNouri edited the manuscript. All authors read and approved the final version of the manuscript.

Acknowledgments

The authors would like to thank the Deputy of Research and Technology and the Students Research Committee, Hamadan University of Medical Sciences, Hamadan, Iran, for their support. The authors would like to appreciate Dr. Zahra Azadian who provided guidance.

Supplementary Materials

Supplementary Table 1: estimated free binding energy and inhibition constant values between 46 flavonoids, control inhibitors, and the SauPBP2a active site (PBD ID: 1MWT; chain B) achieved from the AutoDock software. (Supplementary Materials)

References

J.-P. Rasigade and F. Vandenesch, “Staphylococcus aureus: a pathogen with still unresolved issues,” Infection, Genetics and Evolution, vol. 21, pp. 510–514, 2014.
View at: Publisher Site | Google Scholar
CFD Control and Prevention, “Outbreaks of community-associated methicillin-resistant Staphylococcus aureus skin infections–Los Angeles county, California, 2002-2003,” MMWR: Morbidity and Mortality Weekly Report, vol. 52, p. 88, 2003.
View at: Google Scholar
H. W. Boucher and G. R. Corey, “Epidemiology of methicillin‐resistant Staphylococcus aureus,” Clinical Infectious Diseases, vol. 46, no. S5, pp. S344–S349, 2008.
View at: Publisher Site | Google Scholar
T. A. Taylor and C. G. Unakal, Staphylococcus aureus, StatPearls Publishing, Treasure Island, FL, USA, 2021.
S. Y. C. Tong, J. S. Davis, E. Eichenberger, T. L. Holland, and V. G. Fowler Jr, “Staphylococcus aureus infections: epidemiology, pathophysiology, clinical manifestations, and management,” Clinical Microbiology Reviews, vol. 28, no. 3, pp. 603–661, 2015.
View at: Publisher Site | Google Scholar
F. R. DeLeo, B. A. Diep, and M. Otto, “Host defense and pathogenesis in Staphylococcus aureus infections,” Infectious Disease Clinics of North America, vol. 23, no. 1, pp. 17–34, 2009.
View at: Publisher Site | Google Scholar
V. Hemmadi and M. Biswas, “An overview of moonlighting proteins in Staphylococcus aureus infection,” Archives of Microbiology, vol. 203, no. 2, pp. 481–498, 2021.
View at: Publisher Site | Google Scholar
D. Lim and N. C. Strynadka, “Structural basis for the β lactam resistance of PBP2a from methicillin-resistant Staphylococcus aureus,” Nature Structural Biology, vol. 9, pp. 870–876, 2002.
View at: Google Scholar
M. A. A. Ibrahim, K. A. A. Abdeljawaad, A. H. M. Abdelrahman et al., “Non-β-lactam allosteric inhibitors target methicillin-resistant Staphylococcus aureus: an in silico drug discovery study,” Antibiotics, vol. 10, no. 8, p. 934, 2021.
View at: Publisher Site | Google Scholar
B. G. Cruz, H. S. Dos Santos, P. N. Bandeira et al., “Evaluation of antibacterial and enhancement of antibiotic action by the flavonoid kaempferol 7-O-β-D-(6″-O-cumaroyl)-glucopyranoside isolated from croton piauhiensis müll,” Microbial Pathogenesis, vol. 143, Article ID 104144, 2020.
View at: Publisher Site | Google Scholar
T. P. T. Cushnie and A. J. Lamb, “Detection of galangin-induced cytoplasmic membrane damage in Staphylococcus aureus by measuring potassium loss,” Journal of Ethnopharmacology, vol. 101, no. 1-3, pp. 243–248, 2005.
View at: Publisher Site | Google Scholar
K. Mori, H. Enokida, I. Kagara et al., “CpG hypermethylation of collagen type I α 2 contributes to proliferation and migration activity of human bladder cancer,” International Journal of Oncology, vol. 34, pp. 1593–1602, 2009.
View at: Google Scholar
M. U. Amin, M. Khurram, B. Khattak, and J. Khan, “Antibiotic additive and synergistic action of rutin, morin and quercetin against methicillin resistant Staphylococcus aureus,” BMC Complementary and Alternative Medicine, vol. 15, no. 1, p. 59, 2015.
View at: Publisher Site | Google Scholar
W. Bylka, I. Matlawska, and N. Pilewski, “Natural flavonoids as antimicrobial agents,” Jana, vol. 7, pp. 24–31, 2004.
View at: Google Scholar
H.-X. Xu and S. F. Lee, “Activity of plant flavonoids against antibiotic-resistant bacteria,” Phytotherapy Research, vol. 15, no. 1, pp. 39–43, 2001.
View at: Publisher Site | Google Scholar
H. M. Berman, J. Westbrook, Z. Feng et al., “The protein data bank,” Nucleic Acids Research, vol. 28, no. 1, pp. 235–242, 2000.
View at: Publisher Site | Google Scholar
R. Bouley, M. Kumarasiri, Z. Peng et al., “Discovery of antibiotic (E)-3-(3-carboxyphenyl)-2-(4-cyanostyryl)quinazolin-4(3H)-one,” Journal of the American Chemical Society, vol. 137, no. 5, pp. 1738–1741, 2015.
View at: Publisher Site | Google Scholar
N. Guex, M. C. Peitsch, and T. Schwede, “Automated comparative protein structure modeling with SWISS-MODEL and swiss-PdbViewer: a historical perspective,” Electrophoresis, vol. 30, no. S1, pp. S162–S173, 2009.
View at: Publisher Site | Google Scholar
P. I. O’Daniel, Z. Peng, H. Pi et al., “Discovery of a new class of non-β-lactam inhibitors of penicillin-binding proteins with gram-positive antibacterial activity,” Journal of the American Chemical Society, vol. 136, pp. 3664–3672, 2014.
View at: Google Scholar
A. L. Lovering, M. C. Gretes, S. S. Safadi et al., “Structural insights into the anti-methicillin-resistant Staphylococcus aureus (MRSA) activity of ceftobiprole,” Journal of Biological Chemistry, vol. 287, no. 38, pp. 32096–32102, 2012.
View at: Publisher Site | Google Scholar
D. Laxmi and S. Priyadarshy, “HyperChem 6.03,” Biotech Software and Internet Report, vol. 3, no. 1, pp. 5–9, 2002.
View at: Publisher Site | Google Scholar
G. M. Morris, R. Huey, and A. J. Olson, “Using autodock for ligand-receptor docking,” Current Protocols in Bioinformatics, vol. 24, no. 14, pp. 1–8, 2008.
View at: Google Scholar
G. M. Morris, R. Huey, W. Lindstrom et al., “AutoDock4 and AutoDockTools4: automated docking with selective receptor flexibility,” Journal of Computational Chemistry, vol. 30, no. 16, pp. 2785–2791, 2009.
View at: Publisher Site | Google Scholar
K. Zhu, K. W. Borrelli, J. R. Greenwood et al., “Docking covalent inhibitors: a parameter free approach to pose prediction and scoring,” Journal of Chemical Information and Modeling, vol. 54, no. 7, pp. 1932–1940, 2014.
View at: Publisher Site | Google Scholar
R. A. Friesner, R. B. Murphy, M. P. Repasky et al., “Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein−ligand complexes,” Journal of Medicinal Chemistry, vol. 49, no. 21, pp. 6177–6196, 2006.
View at: Publisher Site | Google Scholar
Z. Azadian, S. Hosseini, Z. P. Dizjikan et al., “Computational and in vitro validation of cardiogenic induction of quercetin on adipose-derived mesenchymal stromal cells through the inhibition of Wnt and non-Smad-dependent TGF-β pathways,” Journal of Cellular Biochemistry, vol. 123, 2021.
View at: Google Scholar
R. Saito, M. E. Smoot, K. Ono et al., “A travel guide to cytoscape plugins,” Nature Methods, vol. 9, no. 11, pp. 1069–1076, 2012.
View at: Publisher Site | Google Scholar
S. V. Pattar, S. A. Adhoni, C. M. Kamanavalli, and S. S. Kumbar, “In silico molecular docking studies and MM/GBSA analysis of coumarin-carbonodithioate hybrid derivatives divulge the anticancer potential against breast cancer,” Beni-Suef University Journal of Basic and Applied Sciences, vol. 9, pp. 1–10, 2020.
View at: Google Scholar
M. S. Deepika, R. Thangam, T. S. Vijayakumar et al., “Antibacterial synergy between rutin and florfenicol enhances therapeutic spectrum against drug resistant Aeromonas hydrophila,” Microbial Pathogenesis, vol. 135, Article ID 103612, 2019.
View at: Publisher Site | Google Scholar
A. H. Siddiqui and J. Koirala, Methicillin Resistant Staphylococcus aureus, StatPearls Publishing, Treasure Island, FL, USA, 2021.
J. Ren, Y. Lu, Y. Qian, B. Chen, T. Wu, and G. Ji, “Recent progress regarding kaempferol for the treatment of various diseases,” Experimental and Therapeutic Medicine, vol. 18, pp. 2759–2776, Oct 2019.
View at: Google Scholar
Y. Qiu, D. He, J. Yang, L. Ma, K. Zhu, and Y. Cao, “Kaempferol separated from Camellia oleifera meal by high-speed countercurrent chromatography for antibacterial application,” European Food Research and Technology, vol. 246, pp. 2383–2397, 2020.
View at: Google Scholar
M. Kannanoor, B. A. Lakshmi, and S. Kim, “Synthesis of silver nanoparticles conjugated with kaempferol and hydrocortisone and an evaluation of their antibacterial effects,” 3 Biotech, vol. 11, no. 7, p. 317, 2021.
View at: Publisher Site | Google Scholar
J. Yang, J. Guo, and J. Yuan, “In vitro antioxidant properties of rutin,” Lebensmittel-Wissenschaft und-Technologie-Food Science and Technology, vol. 41, no. 6, pp. 1060–1066, 2008.
View at: Publisher Site | Google Scholar
K. Gunathilake, K. Ranaweera, and H. Rupasinghe, “Analysis of rutin, β-carotene, and lutein content and evaluation of antioxidant activities of six edible leaves on free radicals and reactive oxygen species,” Journal of Food Biochemistry, vol. 42, Article ID e12579, 2018.
View at: Google Scholar
A. S. Silva and S. M. Nabavi, Nonvitamin and Nonmineral Nutritional Supplements, Academic Press, Cambridge, MA, USA, 2019.
B. Gullón, T. A. Lú-Chau, M. T. Moreira, J. M. Lema, and G. Eibes, “Rutin: a review on extraction, identification and purification methods, biological activities and approaches to enhance its bioavailability,” Trends in Food Science & Technology, vol. 67, pp. 220–235, 2017.
View at: Publisher Site | Google Scholar
K. Y. Rodríguez-Valdovinos and R. Salgado-Garciglia, “Quantitative analysis of rutin by HPTLC and in vitro antioxidant and antibacterial activities of phenolic-rich extracts from Verbesina sphaerocephala,” Plants, vol. 10, 2021.
View at: Publisher Site | Google Scholar
S. Rong, D. Wan, Y. Fan et al., “Amentoflavone affects epileptogenesis and exerts neuroprotective effects by inhibiting NLRP3 inflammasome,” Frontiers in Pharmacology, vol. 10, p. 856, 2019.
View at: Publisher Site | Google Scholar
J. Oh, H. S. Rho, Y. Yang et al., “Extracellular signal-regulated kinase is a direct target of the anti-inflammatory compound amentoflavone derived from Torreya nucifera,” Mediators of inflammation, vol. 2013, Article ID 761506, 2013.
View at: Publisher Site | Google Scholar
Q. Cao, L. Qin, F. Huang et al., “Amentoflavone protects dopaminergic neurons in MPTP-induced Parkinson’s disease model mice through PI3K/Akt and ERK signaling pathways,” Toxicology and Applied Pharmacology, vol. 319, pp. 80–90, 2017.
View at: Publisher Site | Google Scholar
C.-H. Chen, Y.-C. Huang, Y.-H. Lee et al., “Anticancer efficacy and mechanism of amentoflavone for sensitizing oral squamous cell carcinoma to cisplatin,” Anticancer Research, vol. 40, no. 12, pp. 6723–6732, Dec 2020.
View at: Publisher Site | Google Scholar
X. Shen, X. Niu, G. Li, X. Deng, and J. Wang, “Amentoflavone ameliorates Streptococcus suis-induced infection in vitro and in vivo,” Applied and Environmental Microbiology, vol. 84, 2018.
View at: Google Scholar
J. H. Hwang, H. Choi, E. R. Woo, and D. G. Lee, “Antibacterial effect of amentoflavone and its synergistic effect with antibiotics,” Journal of Microbiology and Biotechnology, vol. 23, no. 7, pp. 953–958, 2013.
View at: Publisher Site | Google Scholar
P. Wayne, Performance Standards for Antimicrobial Susceptibility Testing, Clinical and Laboratory Standards Institute, Wayne, PA, USA, 2011.
J. M. Davis, E. A. Murphy, and M. D. Carmichael, “Effects of the dietary flavonoid quercetin upon performance and health,” Current Sports Medicine Reports, vol. 8, no. 4, pp. 206–213, 2009.
View at: Publisher Site | Google Scholar
L. Aguirre, N. Arias, M. Teresa Macarulla, A. Gracia, and M. P. Portillo, “Beneficial effects of quercetin on obesity and diabetes,” The Open Nutraceuticals Journal, vol. 4, 2011.
View at: Google Scholar
Y. Li, J. Yao, C. Han et al., “Quercetin, inflammation and immunity,” Nutrients, vol. 8, no. 3, p. 167, 2016.
View at: Publisher Site | Google Scholar
D. Xu, M. J. Hu, Y. Q. Wang, and Y. L. Cui, “Antioxidant activities of quercetin and its complexes for medicinal application,” Molecules, vol. 24, 2019.
View at: Publisher Site | Google Scholar
A. F. Almeida, G. I. A. Borge, M. Piskula et al., “Bioavailability of quercetin in humans with a focus on interindividual variation,” Comprehensive Reviews in Food Science and Food Safety, vol. 17, no. 3, pp. 714–731, 2018.
View at: Publisher Site | Google Scholar
A. V. A. David, R. Arulmoli, and S. Parasuraman, “Overviews of biological importance of quercetin: a bioactive flavonoid,” Pharmacognosy Reviews, vol. 10, p. 84, 2016.
View at: Google Scholar
S. Wang, J. Yao, B. Zhou et al., “Bacteriostatic effect of quercetin as an antibiotic alternative in vivo and its antibacterial mechanism in vitro,” Journal of Food Protection, vol. 81, no. 1, pp. 68–78, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Motahareh Masumi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1428

Downloads

1038

Citations