1. Introduction
Antimicrobial resistance (AMR) has been seriously threatening public health and global economic development, and the COVID-19 pandemic has worsened this issue [
1]. Therefore, there is an urgent need to discover new antimicrobial agents, especially those against Gram-negative bacteria [
2,
3,
4,
5]. Due to the unique structure of the cell envelope, the development of antibacterial agents against Gram-negative bacteria has always been a challenge [
3,
4]. Fortunately, via the structural modification of tetracyclines, quinolones, aminoglycosides, etc., several antibiotics have been approved for marketing, such as plazomicin and delafloxacin [
4,
5,
6]. Simultaneously, some natural antibacterial products have been also undergoing clinical research [
5,
6]. Eventually, the resistant mutants will emerge after they spend some time circulating in practice, especially in hospitals [
7]. Meanwhile, some plant secondary metabolites have been gradually gaining attention since they not only have antimicrobial activities but also show a good safety profile [
8,
9]. Among these, plant flavonoids are increasingly gaining widespread recognition [
10,
11,
12,
13,
14,
15].
Flavonoids are an important class of secondary metabolites widely distributed in various plants, and many of them have different degrees of inhibitory activity against numerous pathogenic bacteria. Above all, some of them can enhance the antimicrobial activities of antimicrobial agents and/or even reverse AMR [
16,
17]. Various antibacterial mechanisms, involving the inhibition of DNA, proteins, and cell envelope biosynthesis, and the damage to the cell membrane, have been reported for plant flavonoids [
16,
17]. Simultaneously, many structure–activity relationships have been summarized, while some of them contradict one another [
17,
18,
19]. Recently, the antimicrobial quantitative relationship of plant flavonoids against Gram-positive bacteria has been established that [
20], through the statistical analyses and validation of large sample sets. Based on this relationship, the minimum inhibitory concentration (MIC) of plant flavonoids against Gram-positive bacteria can be predicted from their lipophilic physicochemical parameter LogP (Log
10 of partition coefficient) or LogD
7.40 (Log
10 of distribution coefficient at pH 7.40), which is calculated from their chemical structures using the software ACD/Labs. Moreover, it has been indicated that the cell membrane is the major action site of plant flavonoids against Gram-positive bacteria, which includes damage to phospholipid bilayers, respiratory chain inhibition, etc. [
20,
21]. Also, experiments interfering with menaquinone-4 or menaquinones extracted from
Staphylococcus aureus further suggest that the quinone pool is a key target [
22].
Considering the difficulties and biological obstacles in discovering natural products with strong antibacterial activity toward Gram-negative bacteria, it was spontaneously wondered from above whether the antimicrobial activities of plant flavonoids against Gram-negative bacteria are also related to their lipophilicities, and whether the antimicrobial quantitative relationship could also be established. In addition, do plant flavonoids have similar antibacterial effects and mechanisms against both Gram-negative and Gram-positive bacteria? The aim to explore these above will enable the prediction of Gram-negative bacteria MIC values and speed up the discovery of novel antimicrobial compounds. Therefore, the correlation between the lipophilic physicochemical parameter (LogP or LogD7.40) and the MIC, the antibacterial quantitative structure-activity relationship, and the possible mechanisms of plant flavonoids against Gram-negative bacteria were further explored by following a procedure similar to those used to explore Gram-positive species.
3. Discussion
To predict the antimicrobial activities of plant flavonoids against Gram-negative bacteria and to discover compounds with effective inhibitory activity against Gram-negative bacteria, correlation and regression analyses were performed individually using 37 tested and 46 reported data pairs comprising the LogP (or LogD7.40) values and MICs of plant flavonoids against Gram-positive bacteria. The results indicate that, similar to those against Gram-positive bacteria, the antibacterial activities of plant flavonoids against Gram-negative species are related to their lipophilicities. However, the antimicrobial quantitative relationship could not be established for plant flavonoids against Gram-negative bacteria, and so their antimicrobial activities against Gram-negative bacteria could not be predicted from their LogP values. Moreover, different from mainly acting on the cell membrane of Gram-positive bacteria, it was further proposed that plant flavonoids have multiple mechanisms against Gram-negative bacteria, while the cell membrane is also an important action site among them.
Lipophilicity, generally expressed as LogP, is a very important descriptor indicating membrane permeation. LogD, which refers to a pH-dependent mixture of all electrical species at any given pH value, is considered a better descriptor reflecting the actual partitioning and lipophilicity. LogD
7.40 is the LogD at the pH value of 7.40 which is approximately equal to that in human blood or the media of the MIC test. Both lipophilic parameters are important for drug effects. In our previous study [
21], various physicochemical parameters of plant flavonoids were investigated to analyze the correlation between them and the MIC of plant flavonoids against Gram-positive bacteria, and the results indicate that there is a good correlation only between the parameter LogP or LogD
7.40 and the MIC of plant flavonoids against Gram-positive species. Based on this, both physicochemical parameters were selected for the correlation and regression analyses and to conveniently compare the results from plant flavonoids against Gram-negative and Gram-positive bacteria.
As the structural skeleton of plant flavonoids is similar to that of quinolone antimicrobial agents, 14 flavonoids had been designed and synthesized in a previous study [
39]. An evaluation of the DNA gyrase inhibitory and antibacterial activities showed that the antibacterial activities of these flavonoids did not parallel the potency at the enzyme level, and the inhibition of DNA gyrase was only one possible reason for this [
39]. Therefore, the study concluded that some mechanisms other than DNA gyrase inhibition may play a role in the antibacterial activity of flavonoids. As inferred in our previous study [
20,
21], the antibacterial activities of plant flavonoids against Gram-positive bacteria mainly depend on their lipophilicities and not on the inhibition of DNA gyrase as the cell membrane is a major site where plant flavonoids act on these bacteria. Inspired by the multiple action mechanisms of plant flavonoids against Gram-negative bacteria, the LogP values of those designed flavonoids were calculated to explore the correlation between their enzyme inhibitory activities and their LogP values, which ranged from 2.07 to 3.76. Surprisingly, the results (
Tables S2 and S3) demonstrate a significant correlation between the enzyme inhibitory activities and LogP values of the plant flavonoids, and the character of the obtained regression curve (
Figure S3) was very similar to that of the curve (
Figure 2a or
Figure 2c) with LogP values ranging from 2.07 to 3.76. Furthermore, the antibacterial activities of these flavonoids against Gram-negative bacteria, especially
E. coliss, completely paralleled the potency at the enzyme level. Although no specific MIC values were provided for those flavonoids against some pathogenic strains of Gram-negative bacteria, their relative antibacterial activities could be roughly obtained from the corresponding MIC values in combination with polymyxin B nonapeptide (PMBN) [
35]. As a result, it was further suggested that DNA gyrase is another important target of plant flavonoids against Gram-negative bacteria, which is also supported by other studies [
40,
41,
42].
According to this, the antibacterial activities of plant flavonoids against Gram-positive bacteria should also present a concave feature when the LogP values are less than about 3.0 (
Figure 3b). Unfortunately, the LogP value of 3.0 is close to the boundary one, which could possibly cause the curve character to be inconsistent with actual one when the LogP values are less than about 3.0 (
Figure 3b). This might be similar to the difference in the main target of action of quinolones against Gram-negative (DNA gyrase) and Gram-positive bacteria (Topoisomerase IV). Additionally, their antibacterial activities caused by the inhibition of DNA gyrase were weak overall, according to
Figure 2c and
Figure 3a.
Along with the increase in the LogP values, the antibacterial activities of the plant flavonoids against both Gram-negative and Gram-positive bacteria increased when the LogP values were higher than approximately 3.6 (
Figure 3a,b). However, their antibacterial activities against Gram-positive bacteria increased more quickly than against Gram-negative species. This might be attributed to the difference in the cell envelope structure and the quinone pool composition in the respiratory chain between Gram-negative and Gram-positive bacteria (explained in the following paragraph) [
22]. As mentioned in our previous research [
22], menaquinones (MKs) are the sole quinones for electron transfer in the respiratory chain of Gram-positive bacteria, while two types of quinones, namely, MKs and ubiquinones, have been discovered in that of Gram-negative bacteria. This may be an important reason why the plant flavonoids showed weaker antimicrobial activities against Gram-negative than against Gram-positive bacteria when the LogP values were more than approximately 3.6 (
Figure 3a).
Generally, when the LogP value is approximately between 0 and 3.0, the compound presents good cell membrane permeability to the lipid bilayer. As is well known, the cell envelope of Gram-negative bacteria contains a unique hydrophilic outer membrane and a lipophilic inner membrane which is similar to the cell membrane of Gram-positive bacteria. However, the cell structure of Gram-positive bacteria only contains a cell membrane, without an outer membrane. As shown in
Figure 4, when the LogP was less than zero, it was difficult for plant flavonoids to approach the inner membrane of Gram-negative bacteria and the cell membrane of Gram-positive species to reach the target site. This indicates that plant flavonoids exhibit little antibacterial activity against both Gram-negative and Gram-positive bacteria. When the LogP value was between 0 and 3.0, it was easy for the plant flavonoids to cross the inner membrane of Gram-negative bacteria and the cell membrane of Gram-positive species. Thus, plant flavonoids can inhibit the DNA gyrase of Gram-negative bacteria, exhibiting a certain antibacterial activity, while their antibacterial effects on Gram-positive bacteria mainly relies on the membrane action, showing less antibacterial activity. When the LogP value was more than 3.6, it was difficult for the plant flavonoids to cross the inner membrane of Gram-negative bacteria, and their effect gradually changed from the inhibition of DNA gyrase to cell membrane action, including damaging the membrane and inhibiting respiratory chain, thus having similar effects on Gram-positive species. Therefore, their antibacterial activities against both Gram-negative and Gram-positive bacteria show similar trends with the increase in the LogP values between 3.0 and 6.0. However, the antibacterial activity of plant flavonoids against Gram-positive bacteria with the increase in the LogP value increased more rapidly than that against Gram-negative species. This is likely due to the optional electron acceptors such as ubiquinones and MKs in the quinone pool of Gram-negative bacteria, while there are only MKs in that of Gram-positive species. When the logP value was more than 6.7, it was difficult for plant flavonoids to access the inner membrane (action site) due to the obstruction from the hydrophilic outer membrane, and their antibacterial activities against Gram-negative bacteria gradually decreased. Similarly, their antibacterial activities against Gram-positive species also decreased, when the logP value was more than 9.0, even though there is no outer membrane in Gram-positive bacteria. This is likely due to the difficulty that plant flavonoids have in approaching the polarity head of the cell membrane when the logP value is too high. The above statement provides a reasonable interpretation for the characteristics of both regression curves in
Figure 3, which, in turn, confirms the antibacterial mechanisms in this study and our previous report [
20,
21] of plant flavonoids against both Gram-negative and Gram-positive bacteria. Therefore, both the inhibition of DNA gyrase and action on the cell membrane, including membrane damage and the inhibition of the quinone pool of the respiratory chain, are two important mechanisms of plant flavonoids against Gram-negative bacteria. Additionally, their different antibacterial effects on Gram-negative and Gram-positive bacteria were attributed to these different antibacterial mechanisms and their cell envelope structure.
Although both regression curves (
Figure 2a,b) presented similar characters overall with the change in the LogP values, the MIC values of the tested plant flavonoids were obviously higher than those of the reported plant flavonoids. This was probably due to the different determination conditions, methods, details, and selected pathogenic strains. Among the tested and reported plant flavonoids, nine of them were the same. Thus their antibacterial activities against
E. coli ATCC 25922 were reorganized, and they are listed in
Table S1, allowing for a comparison of the differences in the tested and reported data [
25,
28,
29,
32,
34,
38]. The results indicate that only the MICs of compound
13 (Licoflavone C) presented more than a 10 times difference between the tested (3,026.36 μM) and reported (23.08 μM) values [
29]. Therefore, using the broth microdilution method in 96-well plates, the antibacterial activity of this compound against
E. coli ATCC 25922 was determined in six replicates, and all the MICs of compound
13 against
E. coli ATCC 25922 were 3,026.36 μM.
As mentioned in
Section 2.3, the antibacterial activities of plant flavonoids with logP values less than about 3.6 against Gram-negative bacteria were statistically stronger than those against Gram-positive species, and they also increased overall when the LogP value increased to more than about 3.6. However, the antibacterial activities of most plant flavonoids against Gram-negative bacteria were weak overall. Additionally, stronger antibacterial activities can be acquired through structural modification based on the antibacterial mechanism of plant flavonoids against Gram-negative bacteria, such as the inhibition of DNA gyrase and the quinone pool of the respiratory chain. However, due to the unique cell envelope structure of Gram-negative bacteria, it is challenging to discover plant flavonoid derivatives with very strong antibacterial activities against Gram-negative species. This is different from the case of Gram-positive bacteria, which does not have a hydrophilic outer membrane and only has the sole electron transfer acceptor menaquinone in the quinone pool [
22]. Furthermore, when attempting to discover flavonoids with remarkably inhibitory activity against Gram-negative bacteria, it should also be encouraged to pay an attention to the combinational use of plant flavonoids and clinical antibiotics for combating AMR [
7,
31,
34], as they have a good safety profile and can enhance susceptibility to many antibiotics.
4. Materials and Methods
4.1. Materials, Chemicals and Reagents
Thirty-seven plant flavonoids: icaritin (>98%), isoliquiritigenin (98%), formononetin (98%), isoliquiritigenin (98%), galangin (98%), baicalein (98%), diosmin (95%), hesperetin (97%), puerarin (98%), apigenin (≥95%), diosmetin (98%) and naringenin (97%) were purchased from Shanghai Macklin Biochemical Co., Ltd. (Shanghai, China); α-mangostin (>98.0%), licochalcone A (>98.0%), nobiletin (≥98.5%), tangeritin (≥98.5%), quercitrin (98%), sinensetin (98%), narirutin (98%), orientin (99%) and isoorientin (98%) were purchased from Chengdu Push Bio-technology Co., Ltd. (Chengdu, China); naringin (95%), neohesperidin (≥98%) and hesperidin (95%) were purchased from Shanghai Yuanye Bio-Technology Co., Ltd. (Shanghai, China); eriodictyol (≥98%), eriocitrin (≥98%), rhoifolin (≥98%) and licoflavone C (≥98%) were purchased from Wuhan ChemFaces Biochemical Co., Ltd. (Wuhan, China); glabridin (99.8%) and sophoraflavanone G (>98%) were purchased from Shanghai TopScience Co., Ltd. (Shanghai, China); quercetin (97%) was purchased from Meryer (Shanghai) Biochemical Technology Co., Ltd. (Shanghai, China); methyl-hesperidin (95%) was purchased from Shanghai Acmec Biochemical Co.,Ltd. (Shanghai, China); and didymin (≥98%), 5-demethylnobiletin (≥98%), 4′,5,7-trimethoxyflavone (≥98%), vitexin (≥98%) and isovitexin (≥98%) were purchased from Sichuan Weikeqi Biological Technology Co.,Ltd. (Sichuan, China). All the compounds were stored at −20°C. The stock solutions of the above plant flavonoids were prepared by dissolving them in a certain volume of dimethyl sulfoxide (DMSO) and diluting them with Mueller Hinton broth (MHB) to obtain a concentration of 4,096, 8,192, or 16,384 μg/mL. The stock solution was mixed well and then diluted to the desired concentration with MHB immediately before use. The concentrations of DMSO in all test systems were kept below 5.0%, and all those in the blank controls were 5.0%.
Casein hydrolysate (Qingdao Hope Bio-Technology Co., Ltd., Qingdao, China), starch soluble, beef extract, and agar powder (Sangon Biotech (Shanghai) Co., Ltd., Shanghai, China) were used to prepare the media. Mueller Hinton Agar (MHA) consisted of casein hydrolysate 17.5 g/L, starch soluble 1.5 g/L, beef extract 3.0 g/L, and agar powder 17.0 g/L, which were dissolved in purified water, and its pH value was maintained at 7.40 ± 0.20. MHB was prepared without agar powder according to the same composition and procedure as MHA. DMSO and 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) were purchased from Sangon Biotech (Shanghai) Co., Ltd. (Shanghai, China), and 96-well plates were purchased from Shanghai Excell Biological Technology Co., Ltd. (Shanghai, China). All reagents were of analytical or biochemical grade. All TopPette Pipettors (2~20 μL and 20~200 μL) were purchased from DLAB Scientific Co., Ltd., Beijing, China.
4.2. Bacterial Strains and Growth Conditions
E. coli ATCC 25922 and S. aureus ATCC 25923 were purchased from American Type Culture Collection, Manassas, VA, USA, and they were stored in MicrobankTM microbial storage (PRO-LAB diagnostics, Toronto, Canada) at −20°C. Prior to use, E. coli and S. aureus were cultured onto an MHA plate at 37°C, and then pure colonies from the culture plate were inoculated into MHB at 37°C for 24 h on a rotary shaker (160 rpm) of a CHA-2 gas bath constant-temperature oscillator (Jintan Yineng Experimental Instrument Factory, Changzhou, China). A 1:100 dilution of the overnight culture was made into fresh MHB, and then incubated at 37°C until the exponential phase for the following experiments. MHB was used for antimicrobial susceptibility tests.
4.3. Antimicrobial Susceptibility Assay and Prediction
According to the standard procedure recommended by the Clinical and Laboratory Standards Institute (CLSI) [
43], the exponential phase culture was diluted with MHB to achieve a bacterial concentration of approximately 1.0×10
6 CFU/mL, and then the susceptibility of the plant flavonoids against
E. coli ATCC 25922 or
S. aureus ATCC 25923 was determined using the broth microdilution method on the 96-well plates in triplicate [
7]. Based on the preliminary MIC values of the plant flavonoids, the initial concentration of each compound was set as 1,024, 2,048, or 4,096 μg/mL. After the 96-well plate was incubated at 35°C for 24 h, 20 μL of MTT (4.0 mg/mL) was added to each well, mixed well, and left for 30 min at ambient temperature. The MIC, defined as the lowest concentration of compounds that completely inhibited bacterial growth in the micro-wells, was assessed by observing no color change in the wells, while the blank wells had sufficient bacterial growth [
22].
According to the reported antimicrobial quantitative relationship [
20], the MICs of the plant flavonoids against Gram-positive bacteria
S. aureus were also calculated using the equation
y =−0.1285
x6 + 0.7944
x5 + 51.785
x4−947.64
x3 + 6,638.7
x2−21,273
x + 26,087. Here,
y is the antimicrobial activity (MIC) of the pant flavonoids against
S. aureus; and
x is the LogP value of the corresponding pant flavonoids, calculated using the software ACD/Labs 6.0 (Advanced Chemistry Development, Inc., Ontario, Canada). Generally, the LogP values of the equation range from around 2.3 to 9.0. However, an appropriate extension of this range is also workable. It was considered acceptable for the predicted MIC value to fall into the range of 1/4 to 4 of the tested MIC value. Leniently, those MIC values ranging from 1/8× to 8× of the tested MIC value were also considered acceptable events, especially when the LogP value was outside the predicted range while close to both boundary values. This equation was established from 92 data pairs consisting of the LogP values and the MICs of plant flavonoids against Gram-positive bacteria, with an
R2 value of 0.9413.
4.4. Structures and MICs of Reported Plant Flavonoids
As the antimicrobial activities of a certain compound against different pathogenic bacteria and strains vary, data on the chemical structures of plant flavonoids and their MIC against E. coli (a representative Gram-negative bacterium) were collected from the literature. Simultaneously, the tested strains of E. coli were sourced from American Type Culture Collection (ATCC) where possible. The structures, the MIC values against E. coli, and other related information of the plant flavonoids were unsystematically searched for on Google Scholar, and in several databases, such as SciFinder, Medline, Elsevier, ACS, ScienceDirect, Wiley Online Library, and Springer-Link, using the keywords “flavonoid” and “E. coli”, or and “antimicrobial”, or and “antibacterial”, and or and “Gram-negative bacteria”. Furthermore, the relevant references in the obtained literature were also tracked, and then the structures, the MIC values against E. coli, and other related information regarding the plant flavonoids were collected from the acquired literature. Finally, the structures of the collected compounds were drawn using the software ChemBioDraw Ultra 14.0.
4.5. Correlation and Regression Analyses
The lipophilic physicochemical parameters LogP (Log
10 of partition coefficient) and LogD
7.40 (Log
10 of distribution coefficient at pH 7.40) of the plant flavonoids measured in
Section 4.3 or previously reported were automatically calculated from their chemical structures using the software ACD/Labs 6.0 (Advanced Chemistry Development, Inc., Ontario, Canada). Subsequently, the MIC values (μM) and the LogP and LogD
7.40 values of the plant flavonoids measured or reported were listed in an Excel table. After this, correlation analyses between the calculated LogP or LogD
7.40 value (x) and the MIC (y) of the plant flavonoids in the table were performed using Microsoft Excel software (Microsoft Corporation, USA). Simultaneously, corresponding regression equations were established for further verification, and the value of the correlation coefficient (r) was automatically calculated. It is noteworthy that the compounds (
Table 2) without accurate MIC values (such as more than the measuring concentration) were not considered for the regression analyses, but their related information was used for the following discussion. Finally, the correlation between the antimicrobial activity MIC of the plant flavonoids against Gram-negative bacteria and their physicochemical parameter LogP or LogD
7.40 was validated using the r-test. The coefficient of determination (R
2) was also calculated to determine the fitting precisions according to equation (5) in previous work [
20].
4.6. Comparison of the Characteristics of Regression Curves
As referred to from our previous work [
20], the MIC was further transformed to the log
10(MIC), and then a regression analysis between the log
10(MIC) (y) and the LogP or LogD
7.40 (x) was performed. Then, the characteristics of the regression curves were compared on coordinate maps with the same vertical axis to explore the difference in the antibacterial effects and mechanisms of plant flavonoids against Gram-negative and Gram-positive bacteria, and to determine the possible reasons for the antibacterial effects and mechanisms changing with the lipophilic parameter LogP or LogD
7.40 value.
Author Contributions
Conceptualization, G.Y.; methodology, G.Y.; software, G.Y.; validation, G.Y., Y.Y. and X.X.; formal analysis, G.Y., Y.Y. and X.X.; investigation, Y.Y., X.X., A.F., L.Z., G.Y., F.L. and Y.W.; resources, G.Y.; data curation, Y.Y., X.X., L.Z., G.Y., F.L. and Y.W.; writing—original draft preparation, G.Y., Y.Y., A.F. and X.X.; writing—review and editing, G.Y., Y.Y., A.F., L.Z. and X.X.; visualization, G.Y. and Y.W.; supervision, G.Y.; project administration, G.Y.; funding acquisition, G.Y. All authors have read and agreed to the published version of the manuscript.