Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Pharmacol., 30 January 2026

Sec. Pharmacology of Infectious Diseases

Volume 16 - 2025 | https://doi.org/10.3389/fphar.2025.1715533

This article is part of the Research TopicUnravelling the Secrets of Multidrug Transport ProteinsView all articles

Molecular features of prenylated (iso)flavonoids from Fabaceae in relation to their potential NorA inhibition in Staphylococcus aureus

  • 1Laboratory of Food Chemistry, Wageningen University and Research, Wageningen, Netherlands
  • 2Laboratory of Microbiology and Biotechnology, Faculty of Pharmacy, Universitas Indonesia, Depok, Indonesia
  • 3Biometris, Applied Statistics, Wageningen University and Research, Wageningen, Netherlands

Background: The overexpression of the NorA efflux pump is known to be an important factor in the antimicrobial resistance mechanism of Staphylococcus aureus. Therefore, NorA inhibition can help disarm this pathogen and tackle antimicrobial resistance. In this study, we aim to unravel the molecular properties of prenylated (iso)flavonoids from Fabaceae as potential NorA inhibitors and to propose leading compounds for future research.

Methods: A collection of 37 prenylated isoflavonoids and flavonoids (obtained by purification, chemically synthesized, or commercially purchased), collectively referred to as (iso)flavonoids, was evaluated for its activity against norA-overexpressing Staphylococcus aureus using the checkerboard assay with ciprofloxacin (NorA substrate) and erythromycin (primarily non-NorA substrate), in combination with ethidium accumulation assays. Moreover, a norA-knockout Staphylococcus aureus strain was used to corroborate the specificity of the observed effects. Subsequently, in silico binary QSAR and pharmacophore models were developed to elucidate the key molecular properties for potential NorA inhibition.

Results: Seven prenylated (iso)flavonoids, namely, 8-prenylnaringenin, 6-C,7-O-diprenylnaringenin, glabrene, neobavaisoflavone, wighteone, licoisoflavone A, and glycyrrhisoflavone, potentiated ciprofloxacin up to 8-fold in norA-overexpressing and up to 2-fold in norA-knockout strains at 10 μM, without any membrane permeabilization effects. Moreover, prenylated (iso)flavonoids potentiated erythromycin in norA-overexpressing Staphylococcus aureus only up to 2-fold. Binary QSAR models were generated using datasets from the checkerboard and ethidium accumulation assays with a total prediction accuracy of up to 90% for active and 88% for inactive compounds. Based on QSAR models, the polar surface area, the balance of hydrophobicity and hydrophilicity, and the overall hydrophobicity were correlated with antibiotic potentiation and efflux inhibition of prenylated (iso)flavonoids. Moreover, in our study, we revealed that fractional negative polar surface area and formal (negative) charges are key properties that differentiate prenylated (iso)flavonoids with antimicrobial activity from those that act as potential NorA inhibitors. A pharmacophore model provided the basis for further optimization of prenylated (iso)flavonoids, mainly neobavaisoflavone and wighteone, as potential NorA inhibitors.

Conclusion: In our study, we provide, for the first time, predictive QSAR models of prenylated (iso)flavonoids as potential NorA inhibitors and propose two potential leads based on this family of plant-derived compounds. Future research on the specificity and validation of prenylated (iso)flavonoids as NorA inhibitors is required.

1 Introduction

Staphylococcus aureus (S. aureus) is one of the most pathogenic Gram-positive bacteria responsible for infections in hospitals, community settings, and livestock farming (Silva et al., 2023). Methicillin-resistant S. aureus (MRSA) has become the most prevalent strain worldwide (Guo et al., 2020), and emerging resistance to fluoroquinolones in this strain is increasing (Alseqely et al., 2021). Unfortunately, the introduction of a new antibiotic to replace fluoroquinolones is often followed by the development of antimicrobial resistance (AMR) (Li et al., 2025). Therefore, discovering resistance-modifying agents (RMAs) to increase the potency of new and traditional antibiotics has become urgent (Kumar et al., 2023). AMR in S. aureus can be mediated by the presence of the NorA efflux pump from the major facilitator superfamily (MFS) (Yu et al., 2022). The norA efflux pump gene was found to be overexpressed in 43% of resistant S. aureus strains (Patel et al., 2010; Dashtbani-Roozbehani and Brown, 2021). Discovering efflux pump inhibitors (EPIs) against NorA can help tackle the resistance of S. aureus and restore the activity of existing antibiotics.

The Fabaceae family is listed as one of the largest plant families that contain more than 490 medicinal plant species (Abdelsalam et al., 2022). This plant family contains important secondary metabolites, such as phenolics, alkaloids, peptides, and terpenoids (Maroyi, 2023). Among the various secondary metabolites obtained from Fabaceae, prenylated (iso)flavonoids (i.e., prenylated isoflavonoids and flavonoids) are of particular interest. The addition of a prenyl group (a five-carbon isoprenoid substituent) to (iso)flavonoids enhances antimicrobial activity and the efflux pump-inhibitory activity against S. aureus (Kalli-Angel, 2021; Irianti et al., 2023). Prenyl substituents were found in some phytochemicals reported as NorA EPIs, such as nerol (Coelho et al., 2016), sophoraflavanone G (Sun et al., 2020), 5,7-dihydroxy-8-(2-methylbutanoyl)-6-[3,7-dimethylocta-2,6-dienyl]-4-phenyl-2H-chromen-2-one (Roy et al., 2013), 5,7-dihydroxy-6-(2-methyl-butanoyl)-8-(3-methylbut-2-enyl)-4-phenyl-2H-chromen-2-one (Roy et al., 2013), 5,7-dihydroxy-8-(2-methylbutanoyl)-6-(3-methylbut-2-enyl)-4-phenyl-2H-chromen-2-one (Roy et al., 2013), imperatorin (Joshi et al., 2014), and osthol (Joshi et al., 2014) (all the structures are shown in Supplementary Figure S1A). Among these reported phytochemicals with prenyl substituents, sophoraflavanone G is the only compound with a flavonoid backbone (Supplementary Figure S1A). Sophoraflavanone G showed a remarkable antibiotic potentiation effect, causing a 16-fold antibiotic MIC reduction at 2.4 µM (Sun et al., 2020). In addition to sophoraflavanone G, in our previous study revealed two other prenylated isoflavonoids, namely, neobavaisoflavone (isoflavone) and glabrene (isoflavene), as potential NorA EPIs in S. aureus, causing up to an 8-fold antibiotic MIC reduction in addition to their efflux inhibitory activity (Irianti et al., 2023). These findings further highlight the promising activity of prenylated (iso)flavonoids as potential NorA EPIs and their prospect to be developed further as potent NorA inhibitors.

Based on the characteristics of the top four EPI candidates (neobavaisoflavone, glabrene, glyceollin I, and glyceollin III) reported in our previous study, hydrophobicity and hydrophilic–lipophilic balance appear to play a crucial role in mediating the EPI activity (Irianti et al., 2023). Nevertheless, robust structure–activity relationships (SARs) of prenylated (iso)flavonoids as potential NorA EPIs remain to be established. To date, there have been limited quantitative structure–activity relationship (QSAR) studies of NorA EPIs involving phytochemical analogs. An example of a NorA EPI QSAR study using phytochemical analogs was performed by Nargotra et al., who used 25 analogs of piperine, which is a major constituent of Piper nigrum and Piper longum (Nargotra et al., 2009b). Piperine analogs potentiated ciprofloxacin by reducing the minimum inhibitory concentration (MIC) up to 8-fold in norA-overexpressing S. aureus, and the QSAR study highlighted the partial negative surface area as one of the key molecular features for NorA EPI activity (Nargotra et al., 2009b). Thus far, no QSAR model has been developed for prenylated (iso)flavonoids as potential NorA EPIs. Therefore, in this study, we aim to develop QSAR models for prenylated (iso)flavonoids as potential NorA inhibitors and propose new leading compounds for future research.

In this research, we evaluated a total of 37 prenylated (iso)flavonoids, in contrast to our previous study, where we had a smaller dataset of only 11 compounds (Irianti et al., 2023). The activities of this extended collection of prenylated (iso)flavonoids from different subclasses (flavanones, isoflavans, isoflavene, isoflavones, pterocarpans, and pterocarpenes) with different numbers (mono- and di-prenylated) and different configuration of prenyl groups (ring, chain, and furan) were included (Figure 1). Then, the information obtained from the checkerboard and ethidium (Eth) accumulation assays was used to build binary classification QSAR models. A binary QSAR links the structures of the compounds through molecular descriptors with a “binary” expression of activity (i.e., 1 = active and 0 = inactive) (Gao et al., 1999). The developed binary QSAR models were used to systematically define and compare the important molecular features of prenylated (iso)flavonoids as potential NorA EPIs. Based on our previous work, we hypothesized that hydrophobicity, hydrophilic–lipophilic balance, and polar surface area are important molecular properties for prenylated (iso)flavonoids as potential NorA EPIs (Irianti et al., 2023). Moreover, because certain prenylated (iso)flavonoids exhibited both potential NorA EPI and antimicrobial activities (Irianti et al., 2023), we compared the key molecular properties associated with potential NorA EPI activity (identified in this study) to those linked to antimicrobial activity against MRSA, as determined by a previous QSAR model (Kalli et al., 2021). This comparison allowed us to identify which distinct properties defined prenylated (iso)flavonoids as potential NorA EPIs or as antimicrobials. Last, the 3D optimized structures of the candidate prenylated (iso)flavonoids were aligned with that of the known NorA EPI PQQ16P (Felicetti et al., 2018; Palazzotti et al., 2023) to highlight common pharmacophoric features and to propose possible structural improvements for prenylated (iso)flavonoids as potential NorA EPIs.

Figure 1
Chemical structures of flavonoids, isoflavonoids, and pterocarpans. Each section shows different compounds with structural formulas, emphasizing the diverse chemical modifications. Key examples include flavanones like bavachin, isoflavans like glabridin, and pterocarpans like glyceollidin. Each structure is annotated with names and numbered for reference.

Figure 1. Overview of prenylated (iso)flavonoids from different subclasses evaluated in this study. Prenyl groups are highlighted in red, and A, B, C, and D rings in (iso)flavonoids are labeled.

2 Materials and methods

2.1 Materials

Fluoroquinolone-resistant and norA-overexpressing strain of S. aureus (SA)-1199B and norA-knockout strain (SA-K1758) were obtained from BEI resources, as previously described (Irianti et al., 2023). Ethidium bromide (EtBr), propidium iodide (PI), phosphate buffer saline (PBS) (pH 7.4), reserpine, ciprofloxacin, and erythromycin were obtained from Sigma-Aldrich (United States). Dimethyl sulfoxide (DMSO) was purchased from Brunschwig Chemie B.V. (Netherlands), and ethyl acetate absolute was purchased from Biosolve B.V. (Netherlands). Peptone physiological salt solution (PPS) was obtained from Tritium Microbiologie (Netherlands). Tryptone soy agar (TSA) and tryptone soy broth (TSB) were purchased from Oxoid Ltd. (United Kingdom). Transparent 96-well plates (655161) and black 96-well plates (655900) were purchased from Greiner Bio-One transparent (Austria).

Prenylated isoflavonoids, namely, bavachin (1), isobavachin (2), 6-prenylnaringenin (3), isoxanthohumol (5), wighteone (19), lupiwighteone (20), isowighteone (21), licoisoflavone A (23), and glycyrrhisoflavone (25), were purchased from ChemFaces (China). 8-Prenylnaringenin (4), 7-O-prenylnaringenin (6), 6-C,7-O-diprenylnaringenin (7), 8-C,7-O-diprenylnaringenin (8), 3′-C,7-O-diprenylnaringenin (9), 7,4′-O-diprenylnaringenin (10), 7-O-prenylgenistein (26), and 8-C,7-O-diprenylgenistein (27) were previously synthesized and characterized (Ritsema et al., 2025). Glabridin was purchased from Wako (Japan). Neobavaisoflavone (18) was purchased from PhytoLab GmbH & Co. KG (Germany). Glabrene (17) was purchased from Arctom (United States). Licorisoflavan A (12), licoricidin (14), hispaglabridin A (15), glabrone (22), and licoisoflavone B (24) were previously purified from Glycyrrhiza spp. roots and characterized (van Dinteren et al., 2025). 4′-O-Methylglabridin (11) and hispaglabridin B (16) were previously purified from Glycyrrhiza spp. roots and characterized (Van de Schans et al., 2015). Pterocarpan and pterocarpene compounds, namely, glyceollidin II (28), glyceollin I (29), glyceollin II (30), glyceollin III (31), glyceollin IV (32), glyceollin VI (33), glyceofuran (34), dehydroglyceollin I (35), dehydroglyceollin II (36), and dehydroglyceollin III (37) were previously purified from soybeans (Glycine max (L.) Merrill) and characterized (Van De Schans et al., 2016). The identity and purity of all prenylated (iso)flavonoids were confirmed by reversed-phase ultrahigh-pressure liquid chromatography coupled to photo diode array and ion trap mass spectrometry detection (RP-UHPLC-PDA-IT-MSn), and proton nuclear magnetic resonance spectroscopy (1H NMR) was additionally performed for the synthesized or purified compounds (not the purchased ones). Prenylated (iso)flavonoids showed UV purities ≥80%, except licoricidin, glabrone, and glycyrrhisoflavone. The details of the compounds’ purities are provided in Supplementary Table S1.

2.2 Antibiotic potentiation (checkerboard) assay

The checkerboard assays were conducted with SA-1199B (with all 37 prenylated (iso)flavonoids) and SA-K1758 (only with the best candidate compounds, namely, glabrene, neobavaisoflavone, 8-prenylnaringenin, wighteone, licoisoflavone A, and glycyrrhisoflavone), as described previously with slight modifications (Irianti et al., 2025). One colony was grown in 3 mL TSB and further incubated overnight at 37 °C and 180 rpm. Stock solutions of phytochemicals (prepared in DMSO), along with ciprofloxacin or erythromycin (prepared in PBS and 0.1% v/v acetic acid), were diluted in TSB. Both phytochemicals and antibiotic (ciprofloxacin or erythromycin) TSB solutions were prepared at four times the final test concentrations. Phytochemicals were tested at a final concentration of 10, 20, and 40 µM. An overnight culture of bacteria was diluted in TSB to an OD600 of 0.012 (final OD600 of 0.006 in a 96-well plate). Next, 50 µL of each compound (phytochemical and antibiotic) was combined with 100 µL of the diluted culture in a 96-well plate, resulting in a final bacterial inoculum of 6.92 ± 0.29 log10 CFU/mL. Column 1 consisted of a serial dilution of antibiotic only, and row H consisted of a serial dilution of phytochemicals alone, which were subsequently used to determine the minimum inhibitory concentration (MIC) for each antibiotic and phytochemical. The 96-well plates (Greiner Bio-One transparent, Austria) were incubated at 200 rpm at 37 °C in a shaker incubator (Innova 42R, Eppendorf, Germany). Bacterial growth (OD600) was measured using a Tecan Infinite 200 Pro M nanoplate reader before and after 20 h of incubation. Experiments were conducted in three biological replicates. Final OD values were determined by subtracting the initial absorbance (ΔOD600). The concentration of compounds where samples showed ΔOD600 ≤ 0.18 was considered as the minimum inhibitory concentration (MIC). The fold-reduction (FR) was calculated by comparing the MIC of the antibiotic alone relative to the lowest MIC of the antibiotic in combination with phytochemicals (Equation 1).

Fold reduction FR=MICantibioticMICantibiotic+EPI.(1)

2.3 Eth accumulation assay

Eth accumulation assays were performed according to previous work (Irianti et al., 2023). Due to the limited availability of phytochemicals, Eth accumulation assay of all the prenylated (iso)flavonoids collection was performed only with SA-1199B. The OD600 of the overnight culture was adjusted to 0.6–0.8, with an average inoculum of 9.01 ± 0.13 log10 CFU/mL. Fold changes (FC) in increased Eth accumulation were calculated by comparing the final fluorescence unit (RFU) relative to the negative control at t = 60 min and presented as the means with a standard deviation (Irianti et al., 2023). All samples were measured in three biological replicates in technical triplicates.

2.4 Membrane permeability assay

Propidium iodide (PI) was used to evaluate the membrane permeabilization effect of all 37 prenylated (iso)flavonoids in SA-1199B, as previously reported with slight modifications (Irianti et al., 2023). After the washing step, the OD600 was adjusted to 0.6–0.8 with an average inoculum size of 9.18 ± 0.32 log10 CFU/mL. Experiments were performed in three biological replicates in technical triplicates. The permeabilization activity was represented as a percentage of propidium uptake and calculated with Equation 2.

Propidium uptake %=Exp.valuelowcontrolhighcontrollowcontrol×100.(2)

The experimental (Exp.) value represents the observed relative propidium uptake, which was calculated by taking the final fluorescence unit (RFU) at t = 60 min relative to the negative control. Low control denotes the relative propidium uptake in untreated cells (negative control), and high control denotes the relative propidium uptake in permeabilized cells, which were prepared by heating the cells for 10 min at 95 °C (positive control).

2.5 EPI definitions used in this study

Despite extensive research on EPI discovery, there is no consensus on the definition used to categorize EPI and non-EPI compounds. Therefore, in this study, we developed a classification method based on two assays to define and propose leading EPI candidates. EPI and non-EPI compounds were classified based on the maximum fold reduction in ciprofloxacin MIC (FRCIP) from three biological repetitions of the checkerboard assays and the fold change in increased ethidium accumulation (FC) obtained from Eth accumulation assays using SA-1199B. Due to the antimicrobial properties of most prenylated (iso)flavonoids (Kalli et al., 2021; Long et al., 2022), 10 µM was the chosen concentration to compare the EPI activity via checkerboard and Eth accumulation assays as this is a sub-inhibitory concentration for all the tested compounds. The threshold of FRCIP > 2 is considered a notable antibiotic potentiation effect, as a 2-fold reduction (FRCIP = 2) could still take place due to biological variation. Meanwhile, the threshold of FC > 1.2 was considered a notable Eth accumulation by comparing with the FC in Eth accumulation of the reported NorA EPI reserpine. The summary of EPI definitions used in this study is shown in Table 1. “EPIs” were defined as compounds that showed notable FRCIP (FRCIP > 2) and FC (FC > 1.2) at a concentration of 10 μM, which was equal to or lower than ¼ MIC (10 µM ≤ ¼ MIC). Antimicrobial EPIs were defined as compounds with notable FRCIP and FC at a concentration of 10 μM, which was equal to or higher than ½ MIC, and were labeled as “EPI (antimicrobial).” This distinction was made because the notable FRCIP and FC observed at or above ½ MIC might be affected by the antimicrobial properties of the compounds. Compounds that did not meet the FRCIP and/or FC thresholds at 10 µM were categorized as non-EPI.

Table 1
www.frontiersin.org

Table 1. Overview of the EPI and non-EPI categories used in this study.

2.6 Binary classification QSAR model

Binary classification models were developed according to previous work by Kalli et al. with some modifications (Kalli et al., 2022). Binary QSAR models for EPIs were built using FRCIP and the FC as response variables using the threshold values described in Section 2.5. In general, SMILES codes of the 37 prenylated (iso)flavonoids were obtained from PubChem or, when not available, the chemical structures were drawn manually using ChemDraw (PerkinElmer, version 18.1.2.18), and the SMILES codes were retrieved and imported to Molecular Operating Environment (MOE, v.2019.08, Chemical Computing Group). Partial charge correction followed by energy minimization was performed using MMFF94x as a force field and a root mean square (RMS) gradient of 0.01 kcal/mol/A2. Afterward, a conformational search was conducted using lowModeMD (RMS gradient 0.1 kcal/mol/A2, rejection limit = 100, conformation limit = 1, and iteration limit = 10,000 as in default setting). The best conformer was then used to calculate molecular descriptors (2D and i3D) in MOE. Subsequently, the molecular database was curated from descriptors that (i) had constant values or less than five different values (no variability), (ii) were highly intercorrelated (Rpearson > 0.99) (multicollinearity), and (iii) were complex in interpretation. As a result, 163 descriptors were present in the final dataset (the complete list of selected descriptors is provided in Supplementary Table S2).

Because of the relatively small dataset (37 molecules), binary QSAR models were constructed using the whole dataset, without splitting, using either FRCIP or FC as independent variables. The binary QSAR uses statistical probability estimation to predict the active and inactive compounds (Labute et al., 2002). The performance and validation of binary QSAR models was evaluated by assessing the total accuracy (A), accuracy on actives (A1), accuracy on inactives (A0), and the cross-validated accuracies [total (XA), active (XA1), and inactive (XA0)] (Kalli et al., 2022). The cross-validation of the binary QSAR model was obtained using a leave-one-out procedure (Gao et al., 1999). The binary QSAR models from 2–5 descriptors were built and evaluated for their accuracy. To check the multicollinearity of the selected models, variance inflation factors (VIFs) analysis was performed. The VIF value is calculated as 1/(1-r2), where r2 is the squared correlation coefficient between two molecular descriptors. A VIF<5 indicates the absence of collinearity (Nargotra et al., 2009a).

Within MOE, component limit and smoothing were optimized and adjusted to obtain the best binary QSAR model (Gao et al., 1999). The component limit specifies the number of principal components to be used in the model, and the smoothing parameter, also called observation error, is used to minimize the sensitivity of the derived model (Gao et al., 1999; Labute et al., 2002). To obtain the best model (highest accuracy), a component limit between 2 and 5 and a smooth factor between 0.08 and 0.25 were explored in generating the QSAR models, as previously described (Gao et al., 1999). The best model with the FRCIP dataset was obtained by using a component limit of three and a smooth factor of 0.15 or 0.16. Meanwhile, a component limit of four and a smooth factor of 0.13 generated the best QSAR model with the FC dataset. All the binary QSAR performances are shown in Supplementary Table S3.

2.7 Comparison of EPI with and without antimicrobial activity

Molecular descriptors were calculated as previously described in Section 2.6. Then, the key molecular descriptors defining good EPI activity based on our developed QSAR model were compared with the key molecular descriptors defining good antimicrobial activity against MRSA, obtained from a previous QSAR study for prenylated (iso)flavonoids (Kalli et al., 2021). For each molecular descriptor, the statistical difference between the mean values of the two groups (EPIs with antimicrobial activity and EPIs without antimicrobial activity) was calculated with GraphPad Prism 9 using a Student’s t-test with Welch’s correction (Ahad and Yahaya, 2014).

2.8 Structural alignment of prenylated (iso)flavonoids with PQQ16P

The structures of prenylated (iso)flavonoids were aligned with a 2-phenylquinolone derivative, namely, PQQ16P, a synthetic NorA EPI (Felicetti et al., 2018; Palazzotti et al., 2023), which has structural similarities with prenylated (iso)flavonoids. The alignment of prenylated (iso)flavonoids with PQQ16P was performed in MOE (v.2019.08, Chemical Computing Group). The compounds were energy-minimized as previously described in Section 2.6 and subsequently aligned with a flexible mode iteration limit of 100 and a failure limit of 10. The alignment was performed in three repetitions and evaluated based on the alignment score (S value), with lower values indicating better alignments.

2.9 Pharmacophore elucidation

A pharmacophore model was generated using the pharmacophore elucidation query module of MOE (active coverage 1.0, query spacing 0.6, feature limit 5, and query cluster 1.25). PQQ16P and the best-aligned prenylated (iso)flavonoids (S value < −80) were energy-minimized and used to extract common pharmacophoric features. The best pharmacophore model was selected based on the overlap score and the observed molecular features.

3 Results and discussion

3.1 Potential NorA EPI activity of prenylated (iso)flavonoids

The potential NorA EPI activities of 37 prenylated isoflavonoids were evaluated in norA-overexpressing S. aureus based on the antibiotic potentiation effect of the compounds and Eth accumulation assays (Irianti et al., 2023). In this study, the antibiotic (ciprofloxacin) potentiation effect (FRCIP) of prenylated (iso)flavonoids was determined as the ratio between the MIC of ciprofloxacin alone and its MIC in combination with prenylated (iso)flavonoids at 10 µM. Figure 2A shows the results of checkerboard assay of licoisoflavone A. The MIC of ciprofloxacin alone (48 μM, blue box) was divided by the MIC of ciprofloxacin in combination with licoisoflavone A (23) at 10 µM (6 μM, black box), which resulted in a fold-reduction of 8 (Figure 2A). In addition to that, efflux inhibition is also needed to select active EPIs, which is determined via FC increase in Eth accumulation after 1 h of measurement (Figure 2B). The FRCIP and FC values of all 37 prenylated (iso)flavonoids are described in Table 2.

Figure 2
Panel A shows a heat map of ΔOD₆₀₀ readings for varying concentrations of ciprofloxacin and licoisoflavone A, with a color gradient from purple to white. Highlighted boxes indicate specific values. Panel B displays a line graph of RFU over time, comparing a negative control and licoisoflavone A at 10 μM, with the latter showing higher fluorescence.

Figure 2. The representatives of (A) checkerboard results of licoisoflavone A in combination with ciprofloxacin, and (B) Eth accumulation result of licoisoflavone A. Both experiments were performed in a norA-overexpressing strain (SA-1199B).

Table 2
www.frontiersin.org

Table 2. Summary of the EPI activities of 37 prenylated isoflavonoids from different subclasses against norA-overexpressing S. aureus (SA-1199B)a.

The antimicrobial activities and membrane permeabilizing effects of prenylated (iso)flavonoids (Table 2; Supplementary Figure S2) were checked to avoid false-positive results, as some prenylated (iso)flavonoids are good antimicrobials and membrane permeabilizers (Araya-Cloutier et al., 2018; Kalli et al., 2021).

Based on the defined EPI criteria (as described in Section 2.5) (Table 1), in addition to our previously reported potential NorA EPIs glabrene (17) and neobavaisoflavone (18) (Irianti et al., 2023), five new prenylated (iso)flavonoid compounds, namely, 8-prenylnaringenin (4), 6-C,7-O-diperenylnaringenin (7), wighteone (19), licoisoflavone A (23), and glycyrrhisoflavone (25), exhibited efflux inhibition and potentiated ciprofloxacin at 10 µM without having antimicrobial activity. Compared to our previous study (Irianti et al., 2023), the FRCIP of glabrene in this study was two times higher (2-fold in the previous study vs. 4-fold in this study). This 2-fold difference in the FRCIP might be affected by biological variance or different experimental settings. Among these five new compounds, 8-prenylnaringenin (4), wighteone (19), and licoisoflavone A (23) exhibited better antibiotic potentiation activity (up to 8-fold) than glabrene (17) and neobavaisoflavone (18) at ¼ MIC. Moreover, these three compounds exhibited a better potentiation effect than reserpine (up to 4-fold) at the same concentration. Compared to sophoraflavanone G, a C10-prenylated flavonoid (Supplementary Figure S1A), the top three prenylated (iso)flavonoids demonstrated less antimicrobial activity (MIC of 40 µM). Although sophoraflavanone G potentiated fluoroquinolone antibiotic norfloxacin up to 16-fold at four times lower concentration than prenylated (iso)flavonoids, it showed very good antimicrobial activity (MIC of 9.4 µM) (Sun et al., 2020). The lower antimicrobial properties of our top three prenylated (iso)flavonoids make them interesting for further studies as potential NorA EPIs, as antimicrobial properties are less desired for EPI development (Sharma et al., 2019).

None of the compounds from the pterocarpan and pterocarpene subclasses exhibited promising efflux inhibition and antibiotic potentiation activities in norA-overexpressing S. aureus at 10 µM. In our previous work, we demonstrated that pterocarpans glyceollin I (29) and glyceollin III (31) showed notable activity as potential NorA EPIs, but only at ten times higher concentrations than the concentration used in the current study (Irianti et al., 2023). It should be taken into account that the classification of EPI and non-EPI (shown in Table 2) was based on the activity observed (FRCIP and FC) at the set threshold of 10 µM (as explained in Section 2.5; Table 1). Nonetheless, some “non-EPI” compounds could potentially still be active as potential EPIs at higher concentrations.

3.2 Specificity of potential NorA EPI candidates

In addition to remarkable ciprofloxacin potentiation and an increased level of Eth accumulation in the norA-overexpressing strain, a NorA EPI should neither show remarkable antibiotic potentiation with a primarily non-NorA substrate (such as erythromycin) nor in a norA-knockout strain (FRCIP ≤ 2), as shown in Figures 3A,B with the representative prenylated compound licoisoflavone A (23). The other top EPI candidates from the first screening were tested similarly, except compound 6-C,7-O-diprenylnaringenin (7), due to its limited availability.

Figure 3
Heatmaps showing the effect of Licoisoflavone A and antibiotics on ΔOD600. Panel A represents results with Erythromycin, showing varying shades of purple indicating different OD changes. Panel B represents results with Ciprofloxacin, using green shades. Concentrations of Licoisoflavone A and antibiotics are along the axes. Color bars indicate OD intensity from low to high.

Figure 3. Minimal antibiotic potentiation effect (represented by ΔOD600) shown in the checkerboard results of licoisoflavone A in combination with (A) erythromycin in a norA-overexpressing strain and (B) ciprofloxacin in a norA-knockout strain. The numbers represent the ΔOD600 before and after 20 h of growth. The color code indicates no growth (white) up to strong growth (dark color).

As shown in Table 3, all tested compounds showed minimal antibiotic potentiation effect (FRERY 1–2) in the presence of erythromycin, which is not the main substrate of NorA (Reynolds et al., 2003; Poole, 2007). Moreover, no remarkable potentiation with these prenylated isoflavonoids was seen in the norA-knockout strain in the presence of ciprofloxacin (FRCIP 1–2), which further supports their activity against NorA.

Table 3
www.frontiersin.org

Table 3. Checkerboard results of the most promising EPI candidates in the presence of erythromycin (ERY) and ciprofloxacin (CIP) in SA-1199B and SA-K1758, respectively.

The possibility that these prenylated (iso)flavonoids inhibit other pumps in S. aureus cannot be completely excluded, as minor Eth accumulation in the norA-knockout strain K-1758 was observed in our previous study (Irianti et al., 2023). Therefore, to validate the specificity of prenylated (iso)flavonoids for NorA, future research should involve the use of a “clean” complemented NorA-upregulated mutant, as previously suggested (Laws et al., 2021).

3.3 Effect of structural changes in prenylated (iso)flavonoids on their potential NorA EPI activity

To analyze the effect of structural changes on the potential NorA EPI activity of prenylated (iso)flavonoids, SARs were explored using some prenylated (iso)flavonoids in this study. The SARs were established by analyzing structural differences between the top three compounds, which showed FRCIP up to 8-fold, in comparison with their inactive EPI analogs (at 10 µM), which are further referred as non-EPI. The three pairs of EPI and non-EPI analogs demonstrated that the number of -OH groups and prenyl locations influenced the EPI activity of prenylated flavanones and isoflavones (Figure 4; Supplementary Figure S3). Based on SARs (Figure 4), some main conclusions can be derived:

1. The presence of -OH groups influenced the potential NorA EPI activity of prenylated flavanones and isoflavones. In some cases, hydroxylation increased the NorA EPI activity: from isoxanthohumol (5, non-EPI) to 8-prenylnaringenin (4, EPI) and from isowighteone (21, non-EPI) to licoisoflavone A (23, EPI) and to glycyrrhisoflavone (25, EPI) (Figure 4; Supplementary Figure S3). In other cases, the removal of an -OH group resulted in an increased EPI activity, as observed from isowighteone (21, non-EPI) to neobavaisoflavone (18, EPI) (Supplementary Figure S3).

2. The position of prenylation was also associated with the potential NorA EPI activities in prenylated flavanones and isoflavones. As shown in Figure 4, C6 prenylation improved the NorA EPI activities of isoflavones, as observed in isowighteone (21, non-EPI) vs. wighteone (19, EPI). The same trend was also observed for lupiwighteone (20, non-EPI) vs. wighteone (19, EPI) (Supplementary Figure S3). In contrast, C8 prenylation was favorable for EPI activity in flavanones, as observed in 8-prenylnaringenin (4, EPI) vs. 6-prenylnaringenin (3, non-EPI) (Supplementary Figure S3). The effect of prenyl location on EPI activity appears to be comparable to the effect of prenyl location observed in antibacterial activity against S. aureus (Sato et al., 2006; Zhou and Wan, 2015).

Figure 4
Graph comparing paired compounds on two y-axes: left for FR Cip MIC and right for FC Eth accumulation. Symbols represent compounds: black circle, white circle, black square, striped circle, and white square. Horizontal lines indicate reference levels. Chemical structures are displayed at the bottom associated with compound numbers.

Figure 4. Structure–activity relationship (SAR) analysis of the top-three potential EPIs in comparison with non-EPI analogs. The structural differences between EPI and non-EPI analogs are shown in gray. The dashed lines indicate the minimum of antibiotic potentiation effect (FRCIP) of 2-fold (black) and the minimum Eth accumulation (FC) of 1-fold (red).

The role of -OH groups in the efflux inhibition of (unprenylated) flavonoids in SA-1199B has been previously reported, wherein the additional -OH group in epigallocatechin gallate resulted in decreased efflux inhibition relative to epicatechin gallate (Gibbons et al., 2004; Waditzer and Bucar, 2021). Moreover, the role of prenyl configuration was reported as one of the key molecular properties involved in the antimicrobial activity of prenylated (iso)flavonoids (Araya-Cloutier et al., 2018). Within our dataset of flavanone and isoflavone compounds, changes in prenyl configuration appeared to affect the potential NorA EPI activity, but no clear SAR could be proposed (Supplementary Figure S3). Overall, our findings highlighted that primarily, the number of -OH groups and the prenyl position are linked to the potential EPI activity of prenylated (iso)flavonoids from the flavanone and isoflavone subclasses.

3.4 (Quantitative) structure–activity relationship and molecular properties of prenylated (iso)flavonoids as potential NorA EPIs

Using the EPI activities (checkerboard and Eth accumulation assays) and molecular structures of the collection of prenylated (iso)flavonoids tested (Table 2), robust binary classification QSAR models were built to obtain overall insight into the molecular features underlying potential NorA EPI activity by using the maximum FR and FC as dependent variables. The best model was chosen based on all accuracies (A, A1, A0, XA, XA1, and XA0) that were higher than 75%. Two descriptor models were chosen due to their balanced statistical performance, with total accuracy and cross-validated accuracy ≥75% (Supplementary Figures S4A,B), as well as to avoid overfitting (Kalli et al., 2022).

3.4.1 QSAR binary model from the FR dataset

The binary classification QSAR model from the FRCIP dataset was generated to know the molecular properties of prenylated (iso)flavonoids in potentiating ciprofloxacin (Supplementary Table S3A). Table 4 summarizes the best binary classification model of two molecular descriptors, PEOE_VSA_FPPOS and vsurf_HL1, that explained the ciprofloxacin potentiation activity of prenylated isoflavonoids. This best model obtained with the FR dataset had a total accuracy of 89% and a cross-validated accuracy of 78%. The first descriptor PEOE_VSA_FPPOS refers to the fractional positive polar van der Waals surface area. It is calculated as the sum of the surface area (vi) for which the partial charge (qi) is greater than 0.2, divided by the total surface area (Molecular Operating Environment (MOE), 2024). The second descriptor is vsurf_HL1, which represents the hydrophilic–lipophilic balance (Molecular Operating Environment (MOE), 2024). The vsurf_HL1 can be linked to the interaction of molecules with the hydrophilic and hydrophobic part of the biological target through their surface properties, such as hydrophobicity, hydrogen bond, shape, and electrostatic properties (Moorthy et al., 2014). The binary FRCIP QSAR model with PEOE_VSA_FPPOS and vsurf_HL1 descriptors indicated that both the fractional positive polar van der Waals surface area and balanced hydrophilicity and hydrophobicity influence their antibiotic’s potentiating activities. Within 37 prenylated (iso)flavonoids, it was observed that potential NorA EPIs (without antimicrobial activity) showed PEOE_VSA_FPPOS in a range of 0.05–0.16, and their vsurf_HL1 values were in a range of 0.03–0.14 (Supplementary Figure S5A; Supplementary Table S5).

Table 4
www.frontiersin.org

Table 4. Accuracy and cross-validated accuracy of the best binary QSAR model developed with the FRCIP dataset in this study.

3.4.2 QSAR binary model from the FC dataset

To gain insight into the molecular properties associated with Eth accumulation, the binary QSAR model was built from the FC dataset. Similar to the QSAR binary from the FRCIP dataset, the best binary classification model from the FC dataset was selected based on its overall performance accuracy (Supplementary Table S3B). Table 5 describes the best binary classification model obtained with two molecular descriptors, logP(o/w) and PEOE_VSA_FPNEG, that explained the FC in Eth accumulation. The best model obtained from the FC dataset had total and cross-validated accuracies of 84% and 78%, respectively. Descriptor logP(o/w) refers to the log of the octanol/water partition coefficient (Molecular Operating Environment (MOE), 2024), which is related to hydrophobicity. Descriptor PEOE_VSA_FPNEG represents the fractional negative polar van der Waals surface area, which is the sum of the vi such that the qi is less than −0.2, divided by the total surface area (Molecular Operating Environment (MOE), 2024). The active EPIs (without antimicrobial activity) showed logP(o/w) between 2 and 5 and PEOE_VSA_FPNEG in a range of 0.07–0.14 (Supplementary Figure S5B; Supplementary Table S5).

Table 5
www.frontiersin.org

Table 5. Accuracy and cross-validated accuracy of the best binary QSAR model developed with the FC dataset in this study.

It is interesting to highlight that descriptors PEOE_VSA_FPPOS (Table 4) and PEOE_VSA_FPNEG (Table 5) are positively correlated (Rpearson = 0.87), which indicates that these descriptors explain the same variation in the data. Overall, the obtained binary classification models can explain the FRCIP and FC variation in the datasets with high levels of accuracy. Moreover, the selected descriptors highlight that the fractional polar van der Waals surface area, hydrophilicity–hydrophobicity balance, and overall hydrophobicity of prenylated (iso)flavonoids play a role in efflux inhibition and antibiotic potentiation activity. The importance of hydrophobicity and hydrophilic–lipophilic balance is consistent with our previous study (Irianti et al., 2023) with a smaller subset of prenylated (iso)flavonoids. Moreover, the role of the fractional polar van der Waals surface area observed in this study was in line with the result of the QSAR study with piperine analogs (Nargotra et al., 2009b), where the partial charge surface area is involved in NorA EPI activity. Yet, external validation of the developed QSAR models with a different set of prenylated (iso)flavonoid compounds should be performed to completely validate the proposed binary models.

3.4.3 Distinct molecular properties of prenylated (iso)flavonoids as potential NorA EPIs and as antimicrobials against S. aureus

Previous studies have demonstrated antimicrobial properties of prenylated (iso)flavonoids against S. aureus and MRSA (Kalli et al., 2021; Long et al., 2022). Furthermore, the interplay between antimicrobial and potential NorA EPI activity of some prenylated (iso)flavonoids was already observed in our previous study (Irianti et al., 2023) and in this current study. It was found that some prenylated (iso)flavonoids were active as EPIs and had antimicrobial activity, such as 7-O-prenylnaringenin (6), 4′-O-methylglabridin (11), glabridin (13), licoricidin (14), hispaglabridin A (15), licoisoflavone B (24), and dehydroglyceollin I (35) (Table 2). To explore if there are distinct molecular properties that determine the activity of prenylated (iso)flavonoids as antimicrobials and as potential EPIs, we compared the values of molecular descriptors of active EPIs with and without antimicrobial activity. For this comparison, we employed the molecular descriptors obtained from this study (prenylated (iso)flavonoids as potential NorA EPIs) and descriptors from the previous QSAR study for prenylated (iso)flavonoids as antimicrobials against MRSA, consisting of vsurf_IW7 (hydrophilic integy moment), h_pavgQ (average total charge), vsurf_CW3 (ratio between the hydrophilic surface and the total molecular surface at −1.0 kcal/mol), PEOE_VSA_PPOS (total positive polar van der Waals surface area), E_vdw (van der Waals surface energy), vsurf_DD12 (contact distance of lowest and the second lowest hydrophobic energy points of a molecule), PEOE_VSA+2 (sum of van der Waals surface area for each atom of which the partial charge is between +0.100e and +0.149e), and vsurf_D4 (hydrophobic volume) (Kalli et al., 2021) (Supplementary Figure S6; Supplementary Table S6).

Interestingly, there was a significant difference between potential EPIs with and without antimicrobial activity in terms of their formal charge at pH 7 (h_pavgQ) (Nguyen et al., 2021) and fractional negative polar van der Waals surface area (PEOE_VSA_FPNEG) (p < 0.05) (Figure 5). No significant differences were found with other studied descriptors (Supplementary Figure S6). Based on this comparison, potential EPIs with antimicrobial properties showed h_pavgQ values close to 0, except for compounds 7-O-prenylnaringenin (6) and licoisoflavone B (24), as shown in Figure 5A. Potential EPI compounds without antimicrobial properties demonstrated h_pavgQ values between −0.1 and −0.64 (Figure 5; Supplementary Table S6), which indicated that they have more negative formal charges than those with antimicrobial properties (h_pavgQ ∼ 0). In addition, potential EPIs with antimicrobial properties demonstrated a lower negative polar surface area fraction than those without antimicrobial properties (Figure 5B). Based on this study, we hypothesize the following: (i) formal charges and fractional negative polar surface area are key molecular properties that differentiate potential EPI and antimicrobial prenylated (iso)flavonoids; and (ii) a higher fractional negative polar surface area and more negative formal charges help to reduce the antimicrobial action of potential EPIs for this family of compounds. Overall, this comparison provides insight into the key molecular properties that can be useful for selecting prenylated (iso) flavonoids as potential NorA EPIs without having antimicrobial properties as an EPI is preferably a different chemical entity from an antimicrobial to prevent the risk of resistance (Sharma et al., 2019).

Figure 5
Box plots labeled A and B compare two groups: blue for EPI with antimicrobial activity and red for EPI without. Plot A shows a significant difference in h_pavgQ values between groups, while Plot B indicates a significant difference in PEOE_VSA_FPNEG values. Asterisks denote statistical significance.

Figure 5. Two molecular descriptors, (A) h_pavgQ and (B) PEOE_VSA_FPNEG, that showed significant difference (p < 0.05) values between potential EPIs with and without antimicrobial activity. The descriptors’ values were obtained from the calculation of molecular descriptors (Supplementary Table S6). Compounds 6 and 24 that showed formal charges >0 are indicated in the figure.

3.5 Prospects to improve the activity of prenylated (iso)flavonoids as NorA EPIs

To gain insights into possible structural modifications to improve the potency of prenylated (iso)flavonoids as potential NorA EPI, the seven best prenylated (iso)flavonoids were compared with a known NorA EPI that showed better activity than prenylated (iso)flavonoids. The synthetic NorA EPI PQQ16P was previously reported, and it exhibited 8-fold ciprofloxacin MIC reduction at 1.85 µM, i.e., five times lower concentration than the tested prenylated (iso)flavonoids and reserpine (Felicetti et al., 2018; Palazzotti et al., 2023). In addition, PQQ16P had low antimicrobial activity (MIC >200 µM) and low cytotoxicity (Palazzotti et al., 2023). Interestingly, PQQ16P has some similarities with prenylated (iso)flavonoids regarding the structural features, such as (i) the presence of aromatic rings, (ii) the presence of oxygen and nitrogen (as hydrogen bond acceptors) in the aromatic and heterocyclic ring, and (iii) the propoxy (O-alkyl) group in PQQ16P resembles the O-prenyl group present in some prenylated (iso)flavonoids used in this study (PQQ16P structure is shown in Supplementary Figure S1).

To visualize and evaluate these structural similarities, the seven best prenylated (iso)flavonoids, namely, 8-prenylnaringenin (4), 6-C,7-O-diperenylnaringenin (7), glabrene (17), neobavaisoflavone (18), wighteone (19), licoisoflavone A (23), and glycyrrhisoflavone (25), were aligned or superposed with PQQ16P (Figure 6A). Based on the alignment scores (the lower the S value, the better the alignment), neobavaisoflavone (18) and wighteone (23) showed the best alignment (with S values of −81.50 and −80.06, respectively). The rest of the prenylated compounds demonstrated S scores between −60 and −80, except for glabrene (17), which showed S values between −40 and −60 (Figure 6A; Supplementary Table S7).

Figure 6
Panel A shows a graph with S values plotted against numeric labels, divided into 4-fold and 8-fold sections. Panels B and C display molecular structures with 3D mesh spheres in different colors, connected by dashed lines indicating measurements.

Figure 6. (A) Alignment scores (S value, y-axis) of prenylated (iso)flavonoids (indicated in numbers, x-axis) with PQQ16P from three repetitions. The 4-fold and 8-fold antibiotic potentiation (FRCIP) are indicated in the figure. The S values (calculated in MOE) lower than −80 and higher than −80 are represented by gray and black dots, respectively. (B) Pharmacophore elucidation of neobavaisoflavone (dark gray) and wighteone (pink) with the known EPI PQQ16P (yellow). The color of the spheres represents the following features: green spheres represent the hydrophobic features, and orange and blue spheres represent aromatic rings and hydrogen-bond acceptors, respectively. (C) Pharmacophore model showing the distances (in numbers) between the molecular features in Ångström.

To further highlight the similarities in the structural features, pharmacophoric features were extracted from the best-aligned prenylated (iso)flavonoids (neobavaisoflavone and wighteone) and PQQ16P using pharmacophore query in MOE. Based on the pharmacophore analysis, neobavaisoflavone and wighteone shared similarities with PQQ16P, consisting of two hydrophobic features (green), one aromatic ring feature (orange), and two hydrogen bond acceptor features (blue) (Figures 6B,C), which were previously reported as essential chemical features for NorA recognition and binding (Palazzotti et al., 2023).

As observed in Figure 6B, the -OH groups in the A-ring of neobavaisoflavone and the B-ring of wighteone overlay with the propoxy (O-alkyl) group of PQQ16P, whereas the -OH groups in the B-ring of neobavaisoflavone and the A-ring of wighteone align with the O-methyl (-OCH3) group of PQQ16P. Moreover, the prenyl group in both prenylated isoflavonoids and the piperidine group of PQQ16P were pointed in a similar direction. It is important to consider that the positive ionizable piperidine group in PQQ16P is mainly responsible for electrostatic interaction with negatively charged residues in the NorA-binding pocket (e.g., Glu222 and Asp307) (Palazzotti et al., 2023). The lack of a positive ionizable group in prenylated (iso)flavonoids makes this type of interaction impossible.

Reflecting on these findings, prenylated (iso)flavonoids, mainly neobavaisoflavone and wighteone, are a good starting point for further structural optimization. Based on the structural alignment, possible modifications might be considered, such as O-prenylation, O-methylation, or the addition of a positive ionizable group, as observed in PQQ16P. Overall, our structural alignment and pharmacophore analysis with PQQ16P can further expand the possibility of improving the activity of prenylated (iso)flavonoids as potential NorA as EPIs.

4 Conclusion

In this study, the evaluation and QSAR analysis of prenylated (iso)flavonoids were conducted to understand their molecular features as potential NorA EPIs. Binary QSAR models were developed with a total prediction accuracy of up to 90% for active and 88% for inactive compounds. These models indicate that the fractional polar surface area, the balance of hydrophobicity–hydrophilicity, and the overall hydrophobicity were the main properties related to the potential NorA EPI activity. Furthermore, differences in negative formal charge and fractional negative polar surface area appeared to be important in distinguishing between the antimicrobial and potential NorA EPI activity of prenylated (iso)flavonoids. Compounds 8-prenylnaringenin (4), 6-C,7-O-diprenylnaringenin (7), wighteone (19), licoisoflavone A (23), and glycyrrhisoflavones (25) showed promising activity at 10 μM, in addition to the previously reported candidates neobavaisoflavone and glabrene. Among these compounds, 8-prenylnaringenin, wighteone, and licoisoflavone A at a concentration of 10 µM (≤¼ MIC) showed the best potentiation effect, causing up to 8-fold ciprofloxacin MIC reduction.

Based on our findings, there appears to be a window of opportunity for prenylated (iso)flavonoids as potential NorA EPIs without having antimicrobial activity. The structural alignment and the similar pharmacophoric features of prenylated (iso)flavonoids and the known NorA EPI PQQ16P further indicated that prenylated (iso)flavonoids, mainly neobavaisoflavone and wighteone, might be amenable for further enhancing their activity as potential NorA EPIs. In conclusion, this research provides new insight into NorA EPI discovery by proposing new potential EPIs and unraveling the key molecular features of prenylated (iso)flavonoids as potential NorA inhibitors.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

MI: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Writing – original draft, Writing – review and editing. JR: Data curation, Investigation, Writing – review and editing. JH: Writing – review and editing. J-PV: Supervision, Writing – review and editing. CA-C: Conceptualization, Funding acquisition, Resources, Supervision, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by NWO (Nederlandse Organisatie voor Wetenschappelijk Onderzoek, Netherlands) through the grant number VI.Veni.192.095.

Acknowledgements

The authors thank Khoa Nguyen for his support during data processing.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2025.1715533/full#supplementary-material

References

Abdelsalam, N. R., Hasan, M. E., Javed, T., Rabie, S. M., El-Wakeel, H.E.-D. M., Zaitoun, A. F., et al. (2022). Endorsement and phylogenetic analysis of some fabaceae plants based on DNA barcoding. Mol. Biol. Rep. 49, 5645–5657. doi:10.1007/s11033-022-07574-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahad, N. A., and Yahaya, S. S. S. (2014). Sensitivity analysis of Welch’s t-test. AIP Conference proceedings 1605, 888–893. doi:10.1063/1.4887707

CrossRef Full Text | Google Scholar

Alseqely, M., Newton-Foot, M., Khalil, A., El-Nakeeb, M., Whitelaw, A., and Abouelfetouh, A. (2021). Association between fluoroquinolone resistance and MRSA genotype in Alexandria, Egypt. Sci. Rep. 11, 4253. doi:10.1038/s41598-021-83578-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Araya-Cloutier, C., Vincken, J. P., Van De Schans, M. G. M., Hageman, J., Schaftenaar, G., Den Bestens, H. M. W., et al. (2018). QSAR-based molecular signatures of prenylated (iso)flavonoids underlying antimicrobial potency against and membrane-disruption in gram positive and gram negative bacteria. Sci. Rep. 8, 1–14. doi:10.1038/s41598-018-27545-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Coelho, M. L., Ferreira, J. H., De Siqueira Junior, J. P., Kaatz, G. W., Barreto, H. M., and De Carvalho Melo Cavalcante, A. A. (2016). Inhibition of the NorA multi-drug transporter by oxygenated monoterpenes. Microb. Pathog. 99, 173–177. doi:10.1016/j.micpath.2016.08.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Dashtbani-Roozbehani, A., and Brown, M. H. (2021). Efflux pump mediated antimicrobial resistance by Staphylococci in health-related environments: challenges and the quest for inhibition. Antibiotics (Basel) 10, 10. doi:10.3390/antibiotics10121502

PubMed Abstract | CrossRef Full Text | Google Scholar

Felicetti, T., Cannalire, R., Pietrella, D., Latacz, G., Lubelska, A., Manfroni, G., et al. (2018). 2-Phenylquinoline S. aureus NorA efflux pump inhibitors: evaluation of the importance of methoxy group introduction. J. Med. Chem. 61, 7827–7848. doi:10.1021/acs.jcim.3c00516

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, H., Williams, C., Labute, P., and Bajorath, J. (1999). Binary quantitative structure-activity relationship (QSAR) analysis of estrogen receptor ligands. J. Chem. Inf. Comput. Sci. 39, 164–168. doi:10.1021/ci980140g

PubMed Abstract | CrossRef Full Text | Google Scholar

Gibbons, S., Moser, E., and Kaatz, G. W. (2004). Catechin gallates inhibit multidrug resistance (MDR) in Staphylococcus aureus. Planta Med. 70, 1240–1242. doi:10.1055/s-2004-835860

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Y., Song, G., Sun, M., Wang, J., and Wang, Y. (2020). Prevalence and therapies of antibiotic-resistance in Staphylococcus aureus. Front. Cell Infect. Microbiol. 10, 107. doi:10.3389/fcimb.2020.00107

PubMed Abstract | CrossRef Full Text | Google Scholar

Irianti, M. I., Vincken, J.-P., Van Dinteren, S., Ter Beest, E., Pos, K. M., and Araya-Cloutier, C. (2023). Prenylated isoflavonoids from fabaceae against the NorA efflux pump in Staphylococcus aureus. Sci. Rep. 13, 22548. doi:10.1038/s41598-023-48992-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Irianti, M. I., Malloci, G., Ruggerone, P., Lodinsky, E. V., Vincken, J.-P., Pos, K. M., et al. (2025). Indole phytochemical camalexin as a promising scaffold for AcrB efflux pump inhibitors against Escherichia coli. Biomed. Pharmacother. 182, 117779. doi:10.1016/j.biopha.2024.117779

PubMed Abstract | CrossRef Full Text | Google Scholar

Joshi, P., Singh, S., Wani, A., Sharma, S., Jain, S. K., Singh, B., et al. (2014). Osthol and curcumin as inhibitors of human Pgp and multidrug efflux pumps of Staphylococcus aureus: reversing the resistance against frontline antibacterial drugs. MedChemComm 5, 1540–1547. doi:10.1039/C4MD00196F

CrossRef Full Text | Google Scholar

Kalli, S., Araya-Cloutier, C., Hageman, J., and Vincken, J.-P. (2021). Insights into the molecular properties underlying antibacterial activity of prenylated (iso)flavonoids against MRSA. Sci. Rep. 11, 14180. doi:10.1038/s41598-021-92964-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalli, S., Araya-Cloutier, C., Chapman, J., Sanders, J.-W., and Vincken, J.-P. (2022). Prenylated (iso)flavonoids as antifungal agents against the food spoiler Zygosaccharomyces parabailii. Food Control 132, 108434. doi:10.1016/j.foodcont.2021.108434

CrossRef Full Text | Google Scholar

Kalli-Angel, S. (2021). Prenylated (Iso)flavonoids as antimicrobial agents: production activity and mode of action. PhD Thesis. Wageningen, Netherlands: Wageningen University and Research.

CrossRef Full Text | Google Scholar

Kumar, V., Yasmeen, N., Pandey, A., Ahmad Chaudhary, A., Alawam, A. S., Ahmad Rudayni, H., et al. (2023). Antibiotic adjuvants: synergistic tool to combat multi-drug resistant pathogens. Front. Cell Infect. Microbiol. 13, 1293633. doi:10.3389/fcimb.2023.1293633

PubMed Abstract | CrossRef Full Text | Google Scholar

Labute, P., Nilar, S., and Williams, C. (2002). A probabilistic approach to high throughput drug discovery. Comb. Chem. High Throughput Screen. 5, 135–145. doi:10.2174/1386207024607329

PubMed Abstract | CrossRef Full Text | Google Scholar

Laws, M., Hind, C., Rahman, K. M., Sutton, J. M., and Wand, M. E. (2021). Whole genome sequencing of Staphylococcus aureus SA-1199B reveals previously unreported mutations. Int. J. Antimicrob. Agents 57, 106225. doi:10.1016/j.ijantimicag.2020.106225

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Cheng, F., Wei, X., Bai, Y., Wang, Q., Li, B., et al. (2025). Methicillin-resistant Staphylococcus aureus (MRSA): resistance, prevalence, and coping strategies. Antibiotics (Basel) 14, 14. doi:10.3390/antibiotics14080771

PubMed Abstract | CrossRef Full Text | Google Scholar

Long, G.-Q., Li, X.-M., Wang, D.-D., Yang, Y.-C., Zheng, Y.-Y., Liu, X.-L., et al. (2022). Prenylated flavonoids and derivatives isolated from the root barks of Sophora flavescens as potent antibacterial agents against Staphylococcus aureus. Ind. Crops Prod. 189, 115834. doi:10.1016/j.indcrop.2022.115834

CrossRef Full Text | Google Scholar

Maroyi, A. (2023). Medicinal uses of the fabaceae family in Zimbabwe: a review. Plants (Basel) 12, 12. doi:10.3390/plants12061255

PubMed Abstract | CrossRef Full Text | Google Scholar

Molecular Operating Environment (MOE) (2024). 2024.06 Chemical Computing Group ULC. Montreal, Quebec.

Google Scholar

Moorthy, N. H. N., Ramos, M. J., and Fernandes, P. A. (2014). Predictive QSAR models development and validation for human ether-a-go-go related gene (hERG) blockers using newer tools. J. Enzyme Inhibition Med. Chem. 29, 317–324. doi:10.3109/14756366.2013.779264

PubMed Abstract | CrossRef Full Text | Google Scholar

Nargotra, A., Koul, S., Sharma, S., Khan, I. A., Kumar, A., Thota, N., et al. (2009a). Quantitative structure-activity relationship (QSAR) of aryl alkenyl amides/imines for bacterial efflux pump inhibitors. Eur. J. Med. Chem. 44, 229–238. doi:10.1016/j.ejmech.2008.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Nargotra, A., Sharma, S., Koul, J. L., Sangwan, P. L., Khan, I. A., Kumar, A., et al. (2009b). Quantitative structure activity relationship (QSAR) of piperine analogs for bacterial NorA efflux pump inhibitors. Eur. J. Med. Chem. 44, 4128–4135. doi:10.1016/j.ejmech.2009.05.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, P. T. V., Van Dat, T., Mizukami, S., Nguyen, D. L. H., Mosaddeque, F., Kim, S. N., et al. (2021). 2D-quantitative structure–activity relationships model using PLS method for anti-malarial activities of anti-haemozoin compounds. Malar. J. 20, 1–15. doi:10.1186/s12936-021-03775-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Palazzotti, D., Felicetti, T., Sabatini, S., Moro, S., Barreca, M. L., Sturlese, M., et al. (2023). Fighting antimicrobial resistance: insights on how the Staphylococcus aureus NorA efflux pump recognizes 2-phenylquinoline inhibitors by supervised molecular dynamics (SuMD) and molecular docking simulations. J. Chem. Inf. Model. 63, 4875–4887.

PubMed Abstract | CrossRef Full Text | Google Scholar

Patel, D., Kosmidis, C., Seo, S. M., and Kaatz, G. W. (2010). Ethidium bromide MIC screening for enhanced efflux pump gene expression or efflux activity in Staphylococcus aureus. Antimicrob. Agents Chemother. 54, 5070–5073. doi:10.1128/aac.01058-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Poole, K. (2007). Efflux pumps as antimicrobial resistance mechanisms. Ann. Med. 39, 162–176. doi:10.1080/07853890701195262

PubMed Abstract | CrossRef Full Text | Google Scholar

Reynolds, E., Ross, J. I., and Cove, J. H. (2003). Msr(A) and related macrolide/streptogramin resistance determinants: incomplete transporters? Int. J. Antimicrob. Agents 22, 228–236. doi:10.1016/s0924-8579(03)00218-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ritsema, J. H., Kramer, N. I., De Bruijn, W. J., Van Dinteren, S., Franssen, M. C., Vincken, J.-P., et al. (2025). Effective membrane permeabilization of methicillin-resistant Staphylococcus aureus by prenylated phenolics. J. Nat. Prod. 88, 2065–2075. doi:10.1021/acs.jnatprod.5c00540

PubMed Abstract | CrossRef Full Text | Google Scholar

Roy, S. K., Kumari, N., Pahwa, S., Agrahari, U. C., Bhutani, K. K., Jachak, S. M., et al. (2013). NorA efflux pump inhibitory activity of coumarins from Mesua ferrea. Fitoterapia 90, 140–150. doi:10.1016/j.fitote.2013.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, M., Tanaka, H., Tani, N., Nagayama, M., and Yamaguchi, R. (2006). Different antibacterial actions of isoflavones isolated from Erythrina poeppigiana against methicillin-resistant Staphylococcus aureus. Lett. Appl. Microbiol. 43, 243–248. doi:10.1111/j.1472-765x.2006.01963.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, A., Gupta, V. K., and Pathania, R. (2019). Efflux pump inhibitors for bacterial pathogens: from bench to bedside. Indian J. Med. Res. 149, 129–145. doi:10.4103/ijmr.IJMR_2079_17

PubMed Abstract | CrossRef Full Text | Google Scholar

Silva, V., Araújo, S., Monteiro, A., Eira, J., Pereira, J. E., Maltez, L., et al. (2023). Staphylococcus aureus and MRSA in livestock: antimicrobial resistance and genetic lineages. Microorganisms 11, 124. doi:10.3390/microorganisms11010124

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Z. L., Sun, S. C., He, J. M., Lan, J. E., Gibbons, S., and Mu, Q. (2020). Synergism of sophoraflavanone G with norfloxacin against effluxing antibiotic-resistant Staphylococcus aureus. Int. J. Antimicrob. Agents 56, 106098. doi:10.1016/j.ijantimicag.2020.106098

PubMed Abstract | CrossRef Full Text | Google Scholar

Van De Schans, M. G., Ritschel, T., Bovee, T. F., Sanders, M. G., De Waard, P., Gruppen, H., et al. (2015). Involvement of a hydrophobic pocket and helix 11 in determining the modes of action of prenylated flavonoids and isoflavonoids in the human estrogen receptor. Chembiochem 16, 2668–2677. doi:10.1002/cbic.201500343

PubMed Abstract | CrossRef Full Text | Google Scholar

Van De Schans, M. G., Vincken, J.-P., De Waard, P., Hamers, A. R., Bovee, T. F., and Gruppen, H. (2016). Glyceollins and dehydroglyceollins isolated from soybean act as SERMs and ER subtype-selective phytoestrogens. J. Steroid Biochem. Mol. Biol. 156, 53–63. doi:10.1016/j.jsbmb.2015.11.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Dinteren, S., Ritsema, J. H., Sanders, M. G., Meijerink, J., Vincken, J. P., and Araya-Cloutier, C. (2025). Unraveling the molecular drivers of antibacterial prenylated (iso)flavonoids and chalcones against Streptococcus mutans. Sci. Rep. 15, 14776. doi:10.1038/s41598-025-98782-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Waditzer, M., and Bucar, F. (2021). Flavonoids as inhibitors of bacterial efflux pumps. Molecules 26, 6904. doi:10.3390/molecules26226904

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, X. H., Hao, Z. H., Liu, P. L., Liu, M. M., Zhao, L. L., and Zhao, X. (2022). Increased expression of efflux pump norA drives the rapid evolutionary trajectory from tolerance to resistance against ciprofloxacin in Staphylococcus aureus. Antimicrob. Agents Chemother. 66, e0059422. doi:10.1128/aac.00594-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, B., and Wan, C.-X. (2015). Phenolic constituents from the aerial parts of Glycyrrhiza inflata and their antibacterial activities. J. Asian Nat. Prod. Res. 17, 256–261. doi:10.1080/10286020.2014.966095

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: efflux pump inhibitors, Fabaceae, fluoroquinolone, NorA, prenylated (iso)flavonoids, QSAR, Staphylococcus aureus

Citation: Irianti MI, Ritsema JH, Hageman J, Vincken J-P and Araya-Cloutier C (2026) Molecular features of prenylated (iso)flavonoids from Fabaceae in relation to their potential NorA inhibition in Staphylococcus aureus. Front. Pharmacol. 16:1715533. doi: 10.3389/fphar.2025.1715533

Received: 29 September 2025; Accepted: 15 December 2025;
Published: 30 January 2026.

Edited by:

Rajesh Vishwanathan, Indian Institute of Science Education and Research Tirupati, India

Reviewed by:

Sarit Prabha, Maulana Azad National Institute of Technology, India
Abd. Wahid Rizaldi Akili, Padjadjaran University, Indonesia

Copyright © 2026 Irianti, Ritsema, Hageman, Vincken and Araya-Cloutier. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Carla Araya-Cloutier, Y2FybGEuYXJheWFjbG91dGllckB3dXIubmw=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.