Antibacterial activity and mechanism of plant flavonoids to gram-positive bacteria predicted from their lipophilicities

Antimicrobial resistance seriously threatened human health, and new antimicrobial agents are desperately needed. As one of the largest classes of plant secondary metabolite, flavonoids can be widely found in various parts of the plant, and their antibacterial activities have been increasingly paid attention to. Based on the physicochemical parameters and antibacterial activities of sixty-six flavonoids reported, two regression equations between their ACD/LogP or LogD7.40 and their minimum inhibitory concentrations (MICs) to gram-positive bacteria were established with the correlation coefficients above 0.93, and then were verified by another sixty-eight flavonoids reported. From these two equations, the MICs of most flavonoids against gram-positive bacteria could be roughly calculated from their ACD/LogP or LogD7.40, and the minimum MIC was predicted as approximately 10.2 or 4.8 μM, more likely falls into the range from 2.6 to 10.2 μM, or from 1.2 to 4.8 μM. Simultaneously, both tendentiously concave regression curves indicated that the lipophilicity is a key factor for flavonoids against gram-positive bacteria. Combined with the literature analyses, the results also suggested that the cell membrane is the main site of flavonoids acting on gram-positive bacteria, and which likely involves the damage of phospholipid bilayers, the inhibition of the respiratory chain or the ATP synthesis, or some others.

Antimicrobial resistance (AMR) has been seriously threatened human public health and global economic development, and new antimicrobial agents are desperately needed 1,2 . Antibiotics, as the secondary metabolites produced by many bacteria, actinomycetes and fungi, showed remarkably antimicrobial activities, while they also bring some toxic side effects to human body, and are unavoidable to lead to the resistance 3 . Many plant ingredients present weaker antimicrobial activities, while some of them can reverse the resistance of antimicrobial agents 4 . Simultaneously, most of them are considered nontoxic to human body because of their ubiquity in all sorts of plant derived foods and beverages.
As one of the largest classes of plant secondary metabolite, flavonoids can be widely found in various parts of the plants, such as fruit, vegetables, nuts and tea 4 . These compounds have a wide range of pharmacological activities including antibiosis, antioxidation, and coronary heart disease prevention, etc. It is worth noting that some flavonoids can enhance the sensitivity of bacteria to antibiotics, and even reverse the AMR 4,5 . Thereout, the antibacterial activities of flavonoids have been paid more and more attention to. Recently, several investigations were performed for the antimicrobial activities of flavonoids, and the probable relationships between their chemical structures and antimicrobial activities were also summarized [4][5][6] . However, the regularity conclusions on the structure-activity relationships of flavonoids against bacteria still need to be further explored.
During our researches on antimicrobial agents [7][8][9] , it is vaguely found that the antimicrobial activities of flavonoids are not related to their special structure, while may be related to their polarities or lipid-water partition coefficients. Many data of plant flavonoids, involving their chemical structures and antibacterial activities reported in previous papers, were searched and analyzed for proving it. The inhibitory activities of plant flavonoids against gram-positive bacteria especially Staphylococcus aureus can be widely searched, while those against gram-negative ones and fungi were reported too few to carry out statistical analyses 4,6 . Thereby, the former was our focus in this research. As the inhibitory activities of a compound against different pathogenic bacteria are Data analysis and correlation establishment. The regression analyses for the physicochemical parameters CLogP, ACD/LogP, or LogD 7.40 and the antimicrobial activities (MIC or MIC 90 ) of these flavonoids to a certain pathogenic bacterium were respectively performed, and their regression curves were showed on Fig. S1 to S6 in Supplementary Information. From these figures, nearly all regression curves indicate that the antibacterial activities of these flavonoids present similar change characteristics along with the increase of their LogP or LogD 7. 40 . First, the antibacterial activities will dramatically increase when the LogP or LogD 7. 40     www.nature.com/scientificreports/    www.nature.com/scientificreports/ to a specific value. Along with the further increase of LogP or LogD 7.40 , the antibacterial activities will first increase tendentiously and then decrease. Simultaneously, their regression equations between the physicochemical parameter (x) and the MIC (y), together with the correlation coefficients (r), were respectively presented on Fig. S1 to S6, and summarily listed in Table 7. Most correlation coefficients (r) were more than 0.90 (Table 7). This indicated that there is a good correlation between the physicochemical parameter CLogP, ACD/LogP, or LogD 7.40 and the antimicrobial activities (MIC), of these flavonoids to a certain pathogenic bacterium.  www.nature.com/scientificreports/ As we pointed out above, the antimicrobial activities of a compound against different pathogenic bacteria were varied, and even against the same one in different determination conditions. Thereby, the regression analyses were respectively performed for these flavonoids reported in different papers. Considering that the pathogenic bacteria used for antibacterial experiments mainly involved S. aureus, S. epidermidis, and B. subtilis, the same compound should present similar inhibitory activities and identical antibacterial mechanism to these gram-positive bacteria. Thereby, we put the physicochemical parameters and the average MICs to S. aureus, S. epidermidis, or/and B. subtilis (Tables 1, 2, 3, 4, 5 and 6), of these flavonoids together for further regression analyses. The results indicated that the correlation between CLogP and antibacterial activities (MICs) is weak with a correlation coefficient of 0.8412, while that between ACD/LogP or LogD 7.40 (x) and MICs (y) is more reliable (Fig. 7). The regression equations were respectively expressed as y = − 1.6745x 5 + 56.143x 4 − 741.93x 3 + 4831.8x 2 − 15531x + 19,805 and y = − 1.1474x 5 + 38.802x 4 − 515.39x 3 + 3361.9x 2 − 10789x + 13,706, with the correlation coefficients of 0.9349 and 0.9309, respectively. These further proved, by a larger sample, that the inhibitory activities of these flavonoids to gram-positive bacteria will nonlinearly increase as the ACD/LogP or LogD 7.40 increase to approximately 7.0, and then decrease along with the further increase of ACD/LogP or LogD 7.40 .
Verification. To verify the above correlations, other sixty-eight flavonoids (Fig. 8) including flavone, isoflavone, flavonol, flavanonol, dihydroflavone, dihydroisoflavone, flavane, and chalcone subclasses etc., reported in seven papers 4,16-21 , were selected for the comparison of theoretical and reported MICs. Using above two regres-  7.40 , and y is the antimicrobial activities (MICs)), the theoretical MICs of these flavonoids can be calculated. Considering that many factors, such as determination method, concentration of bacterial suspension, and test medium used, may influence on the determination of MIC 5 , the results reported would fluctuate within a reasonable range of the actual values. Thereout, the predicted MICs ranged from 1/4 × to 4 × the determined one were acceptable (marked as A), especially those ranged from 1/2 × to 2 × the determined one, were considered as complete coincidence (marked as C) since the MICs were generally determined by double dilution method 22 . Simultaneously, those more than or equal to the minimum value when the determined MICs were no upper limit were also regarded as complete coincidence (marked as C). Otherwise, those were unacceptable (marked as U). The results (Table 8) indicated that the predicted MICs were in acceptable or complete coincidence with the measured ones for approximate 85.3% flavonoids. Although the antibacterial activities of ten flavonoids (14.7%) are unsatisfactorily predicted, there are six compounds with the predicted MICs falling into the range of 1/8 × to 8 × determined ones. This together indicated that the MICs of most flavonoids against gram-positive bacteria can be roughly calculated from their ACD/LogP or LogD 7.40 although the predicted values are not in accordance with their tested ones for a few flavonoids. At least, these indicated that the ACD/LogP or LogD 7.40 is a key factor for the inhibitory activities of plant flavonoids against gram-positive bacteria.

Discussion and conclusion
Flavonoids can be widely found in various parts of the plant, and their antibacterial activities have been paid more and more attention to, especially after some of them were discovered to have the potency to enhance the susceptibility of some antibiotics to bacteria 4,5 . Based on the related data of plant flavonoids reported, many related physicochemical parameters were calculated, using software ChemBioDraw Ultra 12.0 and ACD/Labs 6.0, for the discovery of the correlations between the physicochemical parameters and the MICs of flavonoids against gram-positive bacteria. Two regression equations between the ACD/LogP or LogD 7.40 (x) and the antimicrobial activities (MICs) (y) were established as y = − 1.6745x 5  After all, the antibacterial activities of a compound to different pathogens are varied, and so these two regression equations, mainly valuable for Staphylococcus and Bacillus, may not always be suitable for flavonoids to other gram-positive bacteria. However, the acceptable range from 1/4 × to 4 × the determined MICs will increase the applicability of these two equations used for the prediction of plant flavonoids to other gram-positive bacteria. To say the least, if necessary, similar regression equations can be also established from the physicochemical parameters and the MICs to other gram-positive bacteria, of flavonoids. Thereby, we concluded that the MICs of most flavonoids against gram-positive bacteria can be roughly calculated from their physicochemical parameters ACD/LogP or LogD7.40.
Lipophilicity is a very important descriptor indicating membrane permeation 23 , and generally expressed as LogP which is valid only for a single electrical species. For ionizable drugs, LogD that refers to a pH-dependent mixture of all electrical species presented at any given pH was regarded as a better descriptor reflecting the actual partitioning and lipophilicity 24,25 . Generally, most flavonoids contain two or more phenolic hydroxyl groups 4-6 , and present similar weak acidity with the pKa of 7.0 to 10.0. Thereby, their LogD will correspondingly decrease along with the increase of environmental pH from about 5.0. Considering the pH in human blood or in the media of MIC determination was approximately 7.40, their LogD at pH 7.40 were selected. These together above indicate that the lipophilicity of plant flavonoids is a key factor for their inhibitory activities to gram-positive bacteria. As the lipophilicity is closely related to membrane permeability 26 , the tendentiously concave regression curves between the antibacterial activity and the LogP or LogD 7.40 also indicate that the cell membrane is probably an important site of flavonoids acting on gram-positive bacteria.
Different antibacterial mechanisms of plant flavonoids were reported 4-6 , such as causing cell-membrane damage, inhibition on various synthase involving the nucleic acid synthesis, the bacterial respiratory chain, or the cell envelope synthesis. However, the results above suggested that the antibacterial activities of these plant flavonoids had no obvious relationship with the specific fragments of their structures, while presented great relationship with their lipophilicities. Simultaneously, the antibacterial activities of plant flavonoids will dramatically increase as the LogP or LogD increases from 2.5 to 4.0 which range the membrane permeability remarkably decrease while the affinity to lipid bilayer greatly increase [27][28][29] . According to this, plant flavonoids may not target specific synthases, but more likely to nonspecifically act on the cell-membrane bilayer or the respiratory chain to kill www.nature.com/scientificreports/ bacteria. This deduction was indirectly supported by many researches which were reviewed in three paper 4-6 , such as follows: (1)   www.nature.com/scientificreports/ Although many other antibacterial mechanisms acting on various synthase for the nucleic acid or cell envelope syntheses were mentioned in these reviews 4,6 , two facts found from the researches of the cited literature are worth further discussing. First, most flavonoids used for mechanism exploration have the cLogP ranged from about 2.0 to 4.0, and are easy to infiltrate into the bacterial cell, while they present very weak antibacterial activities with the MICs more than 250 μg/mL. Second, most experiments were achieved by the determination of enzyme activities in vitro 35,36 , the molecular docking of flavonoids with various synthases 37 , the proteomics technology without the combination of related experiments and the consideration of first the chicken or the egg 38 . Another thing should be considered is whether some molecules can pass through the cell membrane and infiltrate into the bacterial cell or not. Moreover, previous works indicated the antibacterial activity to gram-positive bacteria was observed only four of fourteen flavonoids, while only four of seven flavonoids with DNA gyrase inhibition showed weak inhibitory activity to gram-positive bacteria 20 . Simultaneously, the authors pointed out that mechanisms other than DNA gyrase inhibition may also play a role in the antibacterial activity. Thereby, the conclusion that some of these flavonoids studied are potent inhibitors of DNA gyrase is worth reconsidering 20 . In fact, this work just right indicated that the inhibitory activity of flavonoids against gram-positive bacteria did not correlate with their in vitro DNA gyrase inhibition to a large extent. This was also supported by previous publication 39 . These together further confirmed that the cell-membrane should be the main region of plant flavonoids acting on Gram-positive bacteria, and which likely involving the disruption or damage of phospholipid bilayers, the inhibition of the respiratory chain or ATP synthesis, or some others.
According to the regression equations and above conclusions, the inhibitory activities of flavonoids to grampositive bacteria will increase when the alkyl especially isopentyl were introduced into the structures of flavonoids no matter carbon position it is introduced into. This can be interpreted that the introduction of alkyl would increase the lipophilicity of flavonoids or the LogP, and thereout increase their interactions with phospholipids of cell membrane. However, the introduction of too many alkyls will overmuch increase the LogP of these flavonoids, and which will lead their lipophilicities too large to pass through the hydrophilic region of phospholipid bilayers. This was proved by previous similar work 26,32,40 . On the contrary, the inhibitory activities of flavonoids to gram-positive bacteria will decrease when polar groups, such as hydroxyl and glycosyl, were introduced into their structures. This can be interpreted as that the excessive hydrophilicity of flavonoids will hinder its infiltration into phospholipid bilayers and interaction with hydrophobic region of cell membrane.
Based on the physicochemical parameters and MICs of various flavonoids, the regression equations and above conclusions were achieved. For a certain subclass of flavonoids, the regression equations with larger correlation coefficient can be established for their more accurate MIC predictions, and then can be further used for the structural design and optimization to obtain more efficient antibacterial activity.
As the inhibitory activities of plant flavonoids against gram-negative bacteria were reported less, it is difficult to draw a statistical conclusion. Considering that the cell envelope of gram-negative bacteria was different from that of gram-positive ones, it is worth further exploring whether the above regression equations and above conclusions are suitable for plant flavonoids against gram-negative bacteria. However, these can provide good references for their related researches. Referring to the above conclusions, the anti-MRSA activities of trimethylhydroquinone, vitamin K 3 and carnosic acid were successfully predicted and verified by our laboratory 9,41 .
In conclusion, the MICs of most flavonoids against gram-positive bacteria can be roughly calculated from their physicochemical parameters ACD/LogP or LogD 7.40 , and the lipophilicity is a key factor of plant flavonoids against gram-positive bacteria. Combined with the analyses of previous publications, the results also suggest that the cell membrane may be the main site of plant flavonoids acting on gram-positive bacteria, and which likely involves the damage of phospholipid bilayers, the inhibition of the respiratory chain or ATP synthesis, or some others. Base on this, the inhibitory activities and mechanisms of plant flavonoids to gram-positive bacteria were diagrammatically presented as Fig. 9.

Methods
Information and data. The structures, antimicrobial activities and other related information of plant flavonoids were unsystematically searched from Google academic search engine, and several databases SciFinder, Medline, Elsevier, ACS, ScienceDirect, Wiley Online Library, Springer-Link, and RSC, using keywords flavonoid and antimicrobial, or and antibacterial, and or and anti-MRSA. Furthermore, the relevant references in the obtained literature were also tracked. The structures, antibacterial activities, and other related information of flavonoids were collected from the obtained literature that can provide more than five or more flavonoids. As the antimicrobial activities of a certain compound against different pathogenic strains were varied, compounds reported in different papers were independently collected for the following analyses. Finally, the structures of selected compounds were drawn using software ChemBioDraw Ultra 12.0.
Simulation calculation of physicochemical parameters. The physicochemical parameters Gibbs energy, LogP, CLogP, MR, CMR and tPSA were calculated using software ChemBioDraw Ultra 12.0. Moreover, another software ACD/Labs 6.0 was also used for the calculations of physicochemical parameters LogP, LogD 7.40 and solubility (SolDB).

Data analysis and correlation establishment. The physicochemical parameters and antibacterial
activities of flavonoids reported in the same paper were respectively listed in a table, even those of the same compound. The regression analyses between the calculated values of each parameter and the antimicrobial activities (expressed as MICs) of all compounds in a table were respectively performed using Microsoft Excel software. It is noting that compounds without related antimicrobial information were not considered for the regression analyses, while they can be used for the following discussion. The physicochemical parameters significantly correlating with the antimicrobial activities were selected for the further analyses of correlations between the physicochemical parameters and antimicrobial activities of flavonoids.
Verification. Some other flavonoids were searched from above several databases, and the chemical structures of various flavonoids presented in previous publications were also drawn using software ChemBioDraw Ultra 12.0. The physicochemical parameters LogP and LogD 7.40 of these flavonoids were respectively calculated by software ACD/Labs 6.0, and then their antimicrobial activities (MICs) were respectively predicted using the above regression equations. Comparing with the predicted MICs with the determined one, the regression equations can be verified.  www.nature.com/scientificreports/ Correspondence and requests for materials should be addressed to G.Y.
Reprints and permissions information is available at www.nature.com/reprints.
Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.