Next Article in Journal
Association between Polymorphisms in Vitamin D Pathway-Related Genes, Vitamin D Status, Muscle Mass and Function: A Systematic Review
Next Article in Special Issue
Effect of Geumgwe-Sinkihwan on Renal Dysfunction in Ischemia/Reperfusion-Induced Acute Renal Failure Mice
Previous Article in Journal
WISER Survivor Trial: Combined Effect of Exercise and Weight Loss Interventions on Insulin and Insulin Resistance in Breast Cancer Survivors
Previous Article in Special Issue
Black Cumin (Nigella sativa L.): A Comprehensive Review on Phytochemistry, Health Benefits, Molecular Pharmacology, and Safety
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Plant-Derived and Dietary Hydroxybenzoic Acids—A Comprehensive Study of Structural, Anti-/Pro-Oxidant, Lipophilic, Antimicrobial, and Cytotoxic Activity in MDA-MB-231 and MCF-7 Cell Lines

by
Monika Kalinowska
1,*,
Ewelina Gołębiewska
1,
Grzegorz Świderski
1,
Sylwia Męczyńska-Wielgosz
2,
Hanna Lewandowska
2,
Anna Pietryczuk
3,
Adam Cudowski
3,
Aleksander Astel
4,
Renata Świsłocka
1,
Mariola Samsonowicz
1,
Anna Barbara Złowodzka
5,
Waldemar Priebe
6 and
Włodzimierz Lewandowski
7
1
Department of Chemistry, Biology and Biotechnology, Bialystok University of Technology, Wiejska 45E Street, 15-351 Bialystok, Poland
2
Institute of Nuclear Chemistry and Technology, 16 Dorodna Street, 03-195 Warsaw, Poland
3
Department of Water Ecology, Faculty of Biology, University of Bialystok, Ciolkowskiego 1J Street, 15-245 Bialystok, Poland
4
Environmental Chemistry Research Unit, Institute of Biology and Earth Sciences, Pomeranian University in Słupsk, Arciszewskiego 22a Street, 76-200 Słupsk, Poland
5
Faculty of Chemistry, Warsaw University of Technology, Noakowskiego 3 Street, 00-664 Warszawa, Poland
6
Department of Experimental Therapeutics, The University of Texas MD Anderson Cancer Center, 1901 East Rd., Houston, TX 77054, USA
7
Institute of Agricultural and Food Biotechnology—State Research Institute, Rakowiecka 36 Street, 02-532 Warsaw, Poland
*
Author to whom correspondence should be addressed.
Nutrients 2021, 13(9), 3107; https://doi.org/10.3390/nu13093107
Submission received: 7 August 2021 / Revised: 1 September 2021 / Accepted: 2 September 2021 / Published: 4 September 2021

Abstract

:
Seven derivatives of plant-derived hydroxybenzoic acid (HBA)—including 2,3-dihydroxybenzoic (2,3-DHB, pyrocatechuic), 2,4-dihydroxybenzoic (2,4-DHB, β-resorcylic), 2,5-dihydroxybenzoic (2,5-DHB, gentisic), 2,6-dihydroxybenzoic (2,6-DHB, γ-resorcylic acid), 3,4-dihydroxybenzoic (3,4-DHB, protocatechuic), 3,5-dihydroxybenzoic (3,5-DHB, α-resorcylic), and 3,4,5-trihydroxybenzoic (3,4,5-THB, gallic) acids—were studied for their structural and biological properties. Anti-/pro-oxidant properties were evaluated by using DPPH (2,2-diphenyl-1-picrylhydrazyl), ABTS•+ (2,2-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid), FRAP (ferric-reducing antioxidant power), CUPRAC (cupric-reducing antioxidant power), and Trolox oxidation assays. Lipophilicity was estimated by means of experimental (HPLC) and theoretical methods. The antimicrobial activity against Escherichia coli (E. coli), Pseudomonas aeruginosa (P. aeruginosa), Staphylococcus aureus (S. aureus), Bacillus subtilis (B. subtilis), Salmonella enteritidis (S. enteritidis), and Candida albicans (C. albicans) was studied. The cytotoxicity of HBAs in MCF-7 and MDA-MB-231 cell lines was estimated. Moreover, the structure of HBAs was studied by means of experimental (FTIR, 1H, and 13C NMR) and quantum chemical DFT methods (the NBO and CHelpG charges, electrostatic potential maps, and electronic parameters based on the energy of HOMO and LUMO orbitals). The aromaticity of HBA was studied based on the calculated geometric and magnetic aromaticity indices (HOMA, Aj, BAC, I6, NICS). The biological activity of hydroxybenzoic acids was discussed in relation to their geometry, the electronic charge distribution in their molecules, their lipophilicity, and their acidity. Principal component analysis (PCA) was used in the statistical analysis of the obtained data and the discussion of the dependency between the structure and activity (SAR: structure–activity relationship) of HBAs. This work provides valuable information on the potential application of hydroxybenzoic acids as bioactive components in dietary supplements, functional foods, or even drugs.

1. Introduction

Free radicals are molecules, atoms, or ions that contain an odd (unpaired) number of electrons in the valence orbital. Therefore, they are described as short-lived, highly unstable, and extremely reactive chemical species [1,2,3]. Naturally occurring metabolic processes—especially those related to oxygen respiration—cause the formation of certain amounts of oxygen free radicals, such as superoxide anions (O2•−), hydrogen peroxide (H2O2), or a particularly harmful hydroxyl radical (HO) [4,5]. Various lifestyle-related factors such as poor nutrition (e.g., dietary vitamin deficiencies or high-sugar diet), low physical activity, chronic stress, smoking, alcohol abuse, and air pollution can lead to excessive production of reactive oxygen species (ROS) in the human body [6]. A moderate ROS level in organisms has a positive effect on the immune system and cellular responses, while high concentrations can lead to cell dysfunction, changes in DNA structure, and even apoptosis [7]. Any imbalance between antioxidant defense and ROS production causes a phenomenon known as oxidative stress [8]; this process has been implicated in various degenerative diseases, such as cardiovascular diseases (e.g., coronary heart disease, atherosclerosis), autoimmune diseases (e.g., rheumatoid arthritis), neurodegenerative diseases (e.g., multiple sclerosis, Alzheimer’s disease, Parkinson’s disease), aging, and cancer (e.g., breast cancer) [9,10,11]. Providing antioxidants in food plays an important role in protecting the organism against oxidative stress [12].
Fruits, vegetables, and herbs are natural, generally cheap, and widely available sources of many valuable nutrients (e.g., minerals, fibers, vitamins) and phytochemicals (e.g., flavonoids, phenolic compounds, anthocyanins) [13]. Phenolic compounds are secondary plant metabolites synthesized during biotic (caused by insects and microorganisms) and abiotic (caused by extreme temperatures, excessive lighting, excess water, drought, heavy metals, etc.) environmental stress [14,15]. Phenolic acids (i.e., benzoic acids and cinnamic acids) are a large subclass of phenolic compounds isolated from plants that possess extremely important biological properties, including antioxidant, anti-inflammatory, antiviral, antibacterial, antifungal, and anticancer properties. Epidemiological studies have revealed positive effects of plant phenolic compounds on human health, including reduced risk of cancer, asthma, allergy, and cardiometabolic disorders [16]. Because the human body cannot produce polyphenols, they must be obtained from food [17]. For example, bioactive gallic acid (3,4,5-THB), which belongs to the group of hydroxybenzoic acids, can be found in common plants such as guava (Psidium guajava) [18], grapes (Vitis vinifera) [19], pomegranate (Punica granatum) [20], avocado (Persea americana) [21], blackcurrant (Ribes nigrum) [22], mango (Mangifera indica) [23], and others (Figure 1) [24,25]. Another phenolic compound—gentisic acid (2,5-DHB)—occurs in citrus fruits (Citrus spp.), sesame (Sesamum indicum), olive (Olea europaea), rose gum (Eucalyptus grandis), and Jerusalem artichoke (Helianthus tuberosus). Protocatechuic acid (3,4-DHB) is found in onion, garlic, and relatives (Allium spp.), as well as mulberry (Morus alba), Japanese pepper (Zanthoxylum piperitum), mandarin (Satsuma mandarin), sharp-leaf galangal (Alpinia oxyphylla), and buckwheat (Fagopyrum spp.) [26,27,28].
The significant biological, pharmacological, and medicinal properties of plant phenolic compounds have been reported in many papers. For example, in the study of Kalinowska et al. [29], isoferulic acid was shown to exhibit slight anticancer activity against HaCaT human immortalized keratinocyte cell line [29]. Ashidate et al. [30] reported that gentisic acid (2,5-DHB) inhibited LDL (low-density lipoprotein) oxidation and the formation of cholesterol ester hydroperoxides in human plasma, which may play a huge role in the treatment of atherogenesis [30,31]. The study conducted by Shi et al. [32] showed that protocatechuic acid (3,4-DHB) isolated from the dried fruits of Alpinia Oxyphylla Miq may be a potential neuroprotective agent; the obtained results proved that 3,4-DHB increased the activity of glutathione peroxidase and superoxide dismutase, reduced the content of lipid peroxide, prevented the oxidative damage induced by hydrogen peroxide, and improved the cognitive properties of tested elderly rats [32]. It is well known that oxidative DNA damage can be responsible for the development of cancer [33,34]. In the study conducted by Lin et al. [35], protocatechuic acid (3,4-DHB) showed an apoptotic effect on human gastric adenocarcinoma (AGS) cells [35]. Lodovici et al. [36] reported a potent inhibitory effect on iron-induced oxidative DNA damage by 2,3-DHB and 3,4-THB [36]. In another study, 2,3-DHB decreased the level of NF-κB induction in monocytes treated with hydrogen peroxide. Data show that 200 μM 2,3-dihydroxybenzoic acid increased cell survival to 50%. These anti-inflammatory properties were related to the fact that 2,3-DHB can act as an antioxidant and an iron chelator [37]. Sharma et al. [38] reported that extract of Hibiscus rosa-sinensis containing gentisic acid (2,5-DHB) caused inhibition of the rate-limiting enzyme in tumor promotion (ornithine decarboxylase; ODC) [38].
The antioxidant and antiradical efficiency of hydroxybenzoic acids depends on various factors, one of which is the number and position of the hydroxyl (-OH) group attached to the aromatic ring [39]. It was reported that the antioxidant properties of the phenolic compounds increase with an increase in the number of -OH groups in the phenolic ring. For example, in a study by Brand-Williams et al. [40], gallic acid with three hydroxyl groups showed higher antiradical activity via DPPH assay than protocatechuic acid with two –OH groups [40]. The position of the hydroxyl groups attached to the benzene skeleton also affects the properties of the phenolic acid. In the study by Velika and Kron [41], dihydroxybenzoic acids with the –OH groups substituted in the ortho and para positions to the carboxylic (-COOH) group (2,5-DHB) showed stronger activity against superoxide radicals than phenolic acids with the -OH groups in the meta position with respect to one another (2,4-DHB; 2,6-DHB; 3,5-DHB) [41]. Although the antioxidant properties of selected group(s) of hydroxybenzoic acid have already been described in several other papers, the comprehensive analysis of the biological (antioxidant, anticancer, antimicrobial, lipophilic) and structural (based on spectroscopic and theoretical electronic parameters) properties of the group of seven hydroxybenzoic acids has not previously been conducted. Understanding the dependency between the chemical structure and biological activity of molecules is of great importance in designing new biologically active compounds that can be applied as functional food ingredients, dietary supplements, food antioxidants, preservatives, or drugs.
The present study aimed to examine the biological activity of seven selected hydroxybenzoic acids constituting a logical series of benzoic acid derivatives. The selected compounds differ with respect to the number and position of the hydroxyl groups substituted to the aromatic ring. The research evaluates the anti- and pro-oxidative potential, as well as the antimicrobial, cytotoxic, and lipophilic properties of the selected di- (and one tri-) hydroxy derivatives of benzoic acid. The methods used in this work include (1) the pro-/antioxidant tests, with the use of DPPH radicals, ABTS•+ cation radicals, ferric-reducing antioxidant power (FRAP), cupric-reducing antioxidant capacity (CUPRAC), and Trolox oxidation assays; (2) a microbiological microdilution method with the estimation of the minimum inhibitory concentration (MIC); (3) cytotoxicity testing in the MCF-7 and MDA-MB-231 human breast cancer cell lines; and (4) lipophilicity measurement via the use of HPLC. The structural studies were carried out using spectroscopic methods (FTIR, 1H, and 13C NMR) and quantum chemical calculations in the Gaussian 09W program; the electronic parameters ionization potential, electron affinity, electronegativity, chemical hardness, chemical softness, and electrophilicity index, as well as NBO and ESP atomic charge distribution, were calculated. Principal component analysis (PCA) was used in the statistical analysis of the obtained experimental and theoretical data.

2. Materials and Methods

2.1. Materials

2,3-Dihydroxybenzoic acid (2,3-DHB, pyrocatechuic acid), 2,5-dihydroxybenzoic acid (2,5-DHB, gentisic acid), 2,6-dihydroxybenzoic acid (2,6-DHB, γ-resorcylic acid), 3,5-dihydroxybenzoic acid (3,5-DHB, α-resorcylic acid), 3,4,5-trihydroxybenzoic acid (3,4,5-THB, gallic acid), DPPH (2,2-diphenyl-1-picrylhydrazyl), ABTS (2,2-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid), potassium persulfate (K2S2O8), TPZT (2,4,6-tris(2-pyridil)-s-triazine), iron(III) chloride (FeCl3∙6H2O), sodium acetate (C2H3NaO2∙3H2O), copper(II) chloride (CuCl2), ammonium acetate (CH3COONH4), neocuproine (2,9-dimethyl-1,10-phenanthroline), iron(II) sulfate (FeSO4), Trolox (6-hydroxy-2,5,7,8-tetramethylchroman-2-carboxylic acid), hydrogen peroxide (H2O2), phosphate buffer pH 7, and horseradish peroxide (HRP) were purchased from Sigma-Aldrich Co. (St. Louis, MO, USA). 2,4-Dihydroxybenzoic acid (2,4-DHB, β-resorcylic acid) and 3,4-dihydroxybenzoic acid (3,4-DHB, protocatechuic acid) were purchased from Across Organics. Methanol was sourced from Merck (Darmstadt, Germany). All chemicals had an analytical purity and were used without further purification.
The following bacterial strains were used in this study: Escherichia coli (PCM 2857), Pseudomonas aeruginosa (PCM 2720), Bacillus subtilis (PCM 2850), Staphylococcus aureus (PCM 2267), Salmonella enteritidis (NCTC 4776), and Candida albicans (PCM 2566-FY). The microbial strains were bought from the Polish Collection of Microorganisms (Wroclaw, Poland) or the American Type Culture Collection. Bacteria (0.1 mL of the reconstituted suspension) were seeded onto sterile Mueller–Hinton agar plates (Biomaxima), to which previously suitable amounts of the tested compounds were added to give the desired concentration.
The human breast adenocarcinoma cell lines MDA-MB-231 and MCF-7 were purchased from the American Type Culture Collection (ATCC, Rockville, MD, USA) and maintained according to ATCC protocol. Briefly, MCF-7 cells were cultured in EMEM (Biological Industries, Israel) medium supplemented with 10% fetal calf serum (Gibco). MDA-MB-231 cells were cultured in Leibovitz‘s L15 (Biological Industries, Beit HaEmek, Israel) medium supplemented with 10% fetal calf serum (Gibco). The MCF-7 cells were incubated in a 5% CO2 atmosphere at 37 °C, whereas the MDA-231 cells were incubated in a 100% air atmosphere at 37 °C.

2.2. Methods

DPPH assay was measured according to the method described by Rice-Evans [42]. Appropriate volumes of methanolic solutions of the tested hydroxybenzoic acids (C = 30 μM) were mixed with methanol in test tubes to prepare a series of solutions of various concentrations (final volume 1 mL). Then, 2 mL of methanolic solution of DPPH (C = 60 μM) was added to each tube; the control contained methanol instead of the sample. The tubes with reaction mixtures were capped, vortexed, and left to stand for 60 min at 23 °C in the dark. The absorbance of the solutions was measured by the use of a Cary 5000 UV–Vis spectrophotometer (Santa Clara, CA, USA) at λ = 516 nm against pure methanol. The ability to scavenge DPPH radicals (%I) was calculated according to the equation:
%I = [(Acontrol − Asample)/Acontrol] × 100%
where %I = the % inhibition of DPPH radicals; Acontrol = the absorbance of the control; and Asample = the absorbance of the tested sample. The concentration of the tested hydroxybenzoic acids was plotted against the %I, and the IC50 values (antioxidant concentration that inhibited 50% of DPPH radicals) were calculated from the obtained scavenging curves.
ABTS assay was conducted according to the procedure of Re et al. [43]. The reagent was prepared by mixing aqueous solutions of ABTS (C = 7 mM) and K2S2O8 (C = 2.45 mM) at a volumetric ratio of 1:1, and leaving them for 12–16 h at 23 °C to generate ABTS•+ cation radicals. Then, 1 mL of the obtained ABTS•+ solution was mixed with 60 mL of methanol. Samples containing 1 mL of ABTS•+ solution and 1 mL of methanolic solution of the tested compound (C = 50 μM) were mixed and incubated for 7 min at 23 °C; the control contained methanol instead of the sample. The absorbance of the solutions was measured with a UV–Vis spectrophotometer at 734 nm against methanol (the blank). The ability to scavenge ABTS•+ cation radicals was expressed as the %I of ABTS•+ and calculated according to the aforementioned equation.
FRAP assay was carried out as described in [42]. The working FRAP reagent was prepared by mixing 250 mL of acetate buffer (C = 300 mM; pH 3.6), 25 mL of TPZT in 40 mM HCl (C = 10 mM), and 25 mL of FeCl3 aqueous solution (C = 20 mM) (10:1:1 v/v/v). Then, 3 mL of the freshly prepared FRAP reagent was added to 0.4 mL of methanolic solution of the tested compound (C = 50 μM) and mixed thoroughly. After 8 min of incubation at 23 °C, the absorbance at 595 nm was measured against the blank (containing 0.4 mL of methanol instead of the sample). Ferric-reducing antioxidant activity was expressed as Fe(II) ion equivalents (μM) using the calibration curve obtained for FeSO4 (the standard curve equation: y = 3.4303x − 0.1058; R2 = 0.9984).
CUPRAC assay was carried out according to the method described in [44]. First, 250 mL of CuCl2 (C = 10 mM; in water) was mixed with 250 mL of ammonium acetate (pH 7; in water) and neocuproine (C = 75 mM; in ethanol). Then, 3 mL of the CUPRAC mixture was added to 0.5 mL of methanolic solution of the tested compound (C = 50 μM) and 0.6 mL of distilled water. The samples were mixed and incubated for 60 min at room temperature (23 °C). The absorbance of the samples was measured with a UV–Vis spectrophotometer at λ = 450 nm against a blank (containing 0.5 mL of methanol instead of the sample). Cupric-reducing antioxidant activity was expressed as Trolox equivalents (μM) by using the calibration curve prepared over the range of 10–200 μM of Trolox (y = 4575.8x + 0.0271; R2 = 0.9919).
The pro-oxidant activity was expressed as the rate of Trolox oxidation according to the protocol described in [45]. First, a mixture of 0.5 mL of Trolox (C = 0.4 mM; in water) with 0.5 mL of H2O2 (C = 0.2 mM; in water), 0.5 mL of HRP (C = 0.04 μM; in phosphate buffer; pH = 7), 5 μL of the tested compounds (C = 50 μM; in methanol), and 495 μM of distilled water was prepared as the tested sample. The absorbance was measured at λ = 272 nm every 10 min for 60 min. The control contained methanol instead of the tested substance.
All anti-/pro-oxidant assays were conducted in five repetitions for three independent experiments for each tested compound. An Agilent Cary 5000 spectrophotometer (Agilent, Santa Clara, CA, USA) was used to measure the absorbance. The results were expressed as the means of the values obtained for the replications. Average, standard deviation calculation, and graphs were completed using Microsoft Excel 2019.
The MIC (minimal inhibitory concentration) was determined by serial dilution of the tested compound in an agar medium, to which the appropriate inoculum of microorganisms was then added and incubated. Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, Bacillus subtilis, Salmonella enteritidis, and Candida albicans were grown overnight and then resuspended in physiological saline to an optical density at 600 nm (OD = 600) of 0.60, corresponding to 5.0 × 108 CFU/mL. Bacteria (0.1 mL of the reconstituted suspension) were seeded onto sterile Mueller–Hinton agar plates, to which previously suitable amounts of the tested compounds were added to reach the desired concentration. Tested chemicals were dissolved in DMSO solution. Negative controls were agar plates, to which DMSO was added, while positive controls were plates with gentamycin (in the case of bacteria) or flucanozole (in the case of fungi). The plates were incubated at 37 °C for 24 h. The lowest concentration without visible bacterial growth was determined as the MIC.
Cytotoxicity evaluation was conducted to assess the impact of hydroxybenzoic acids on the proliferation of MCF-7 and MDA-MB-231 cells, which was measured via neutral red (NR) assay. In brief, cells were seeded in 96-well microplates (TPP, Switzerland) at a density of 2 × 104 cells/well in 100 μL of culture medium. At least three independent experiments in six replicate wells were conducted per concentration. Twenty-four hours after cell seeding, cells were treated with increasing concentrations of hydroxybenzoic acids (0.02–5 mM) for forty-eight hours. After treatment, cells were incubated for 3 h at 37 °C with 100 µL of neutral red solution (stock solution of NR) (Sigma-Aldrich), and then 5 mg/mL of PBS was diluted 1:100 in cell culture medium, incubated for 12 h at 37 °C, and centrifuged to remove any undissolved NR powder. Next, NR solution was aspirated, cells were washed with 150 µL of PBS, and 200 µL of an acid–ethanol solution (49% water, 50% ethanol, and 1% acetic acid) was added to each well. After 15 min of gentle shaking, optical density was read at 540 nm with an Infinite M200 plate reader spectrophotometer (Tecan, Austria). Based on the survival curve, the IC50 (the concentration of a compound that is required for the 50% inhibition of free radical viability) was calculated by fitting data to the logistic equation.
The lipophilicity of the studied compounds was determined via RP-HPLC analysis using a Waters Alliance 2695 HPLC separation module (Milford, MA, USA) and a Waters 2996 photodiode array detector (Milford, MA, USA) (λ = 254 nm). The analyses were carried out on the following chromatographic columns: RP-18e Purospher STAR (C18), 150 mm × 4.0 mm, 5 µm, (Merck, Darmstadt, Germany); RP-8e Purospher STAR (C8), 150 mm × 4.6 mm, 5 µm, (Merck, Darmstadt, Germany); Zorbax Eclipse XDB-CN (CN), 150 mm × 4.6 mm, 5 µm, (Agilent Technologies, Santa Clara, CA, USA); Nucleosil Phenyl (PHE), 250 mm × 4.6 mm, 7 µm, (Supelco, Bellefonte, PA, USA); and Regis Rexchrom IAM.PC.DD.2 (IAM), 100 mm × 4.6 mm, 10 µm, (Regis Technologies Inc., Morton Grave, IL, USA). The experimental conditions and the calculations were as described in [46].

2.3. Calculations in ACD/Labs and Gausssian09 Programs

The values of the partition coefficients in the octan-1-ol/water solvent system and pKa were calculated using the ACD/Labs program [47].
The molecules of dihydroxybenzoic acids and trihydroxybenzoic (gallic) acid were optimized using the B3LYP/6-311++G(d,p) method in the gas phase, in water and methanol. Based on the obtained structural data (bond lengths), the aromaticity indices Aj [48], BAC [49], HOMA, GEO, EN [50], and Bird’s I6 index [51] were calculated on the basis of the equations described in the literature. The index values were normalized in such a way that for a pure aromatic system (benzene) they assumed a maximum value of 1 (Aj, BAC, HOMA index) [48,49,50] or 100 (Bird’s I6 index) [51], and for a non-aromatic system (cyclohexatriene) took values of 0. Using the GIAO method in B3LYP/6-311++G(d,p), the NICS magnetic indices [52] were calculated. The NICS index describes the value of the chemical shift of a given chemical compound calculated in the center of the aromatic ring [52]. The distribution of electronic charges on the atoms of the analyzed molecules was calculated using the NBO [53] and CHelpG [54] methods. The electrostatic potential (ESP) distribution maps were calculated using the CHelpG method. The electronic parameters ionization potential, electron affinity, electronegativity, chemical hardness, chemical softness, and electrophilicity index were calculated on the basis of the values of the energy of HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital). Theoretical calculations were performed using the Gaussian09 program [55].

2.4. FTIR, 1H, and 13C NMR Study

The FTIR spectra were recorded with an Alfa (Bruker) spectrometer (Bremen, Germany) within the range of 400–4000 cm−1. Samples in the solid state were measured in KBr matrix pellets and using the ATR technique. The 1H and 13CNMR spectra were measured on a Bruker Avance II 400 spectrometer (Bremen, Germany) in DMSO solution at 25 °C. TMS was used as an internal reference.

2.5. Statistical Analysis

For overall exploration of the datasets including lipophilic (LogPC18, LogPGalas, LogPexp) and electronic parameters (pKa1, ΔE(LUMO-HOMO), IP), aromaticity indices (BAC, HOMA, I6, NICS), FTIR (νC-(OH), β(CH) 18b, γ(CH) 11), and NMR (δC2, δC3) parameters as features designed to build the model—as well as the anti-/pro-oxidant and antimicrobial potential of the selected dihydroxybenzoates as accompanied features—principal component analysis (PCA) was used. PCA produces a low-dimensional representation of a dataset, finds a sequence of linear combinations of the variables that have maximal variance, and serves as an efficient tool for data visualization. Prior to its use, the impact of the implementation of the PCA on a dataset was checked using analysis of Pearson’s correlation matrix. The PCA itself was calculated in its default mode.

3. Results

3.1. DPPH Assay

To determine the antioxidant potential of the analyzed hydroxybenzoic acids to neutralize free radicals, one of the most popular and most frequently employed tests was used: the DPPH (2,2-diphenyl-1-picrylhydrazyl) assay. DPPH is a stable, synthetic free radical that has a free electron on the valence shell of the nitrogen atom that forms a nitrogen bridge [56]. Due to the delocalization of the unpaired electron, the DPPH molecule does not dimerize [57]. When the DPPH radical reacts with a compound capable of donating a hydrogen atom, it is reduced, and the solution turns from purple to yellow, as measured spectrophotometrically [58]. The antioxidant activity data for studied compounds, expressed as IC50 parameters (the concentration of antioxidant necessary to inactivate 50% of free radicals), are displayed in Table 1. The 2,4-DHB and 3,5-DHB showed the lowest antiradical properties (or no antiradical activity), with IC50 values higher than 120,000 and 1000 μM, respectively. Among the tested hydroxybenzoic acids, 3,4,5-THB—with three hydroxyl groups attached to the aromatic ring—possessed the highest DPPH radical scavenging activity (IC50 = 2.42 ± 0.08 μM). The inhibition of DPPH radicals increased with the increase in the concentration of the tested compound, as shown in Figure 2. The 2,3-DHB, 2,5-DHB, and 3,4-DHB showed a slightly lower ability to capture the radicals than 3,4,5-THB. The activity of the tested compounds can be ranked according to their increasing antioxidant potential, as follows: 2,4-DHB < 3,5-DHB < 2,6-DHB < 2,3-DHB < 3,4-DHB < 2,5-DHB < 3,4,5-THB.
The ability of the analyzed phenolic acids to capture the DPPH radicals depends on many factors, one of which is the number of hydroxyl groups in the structure and their position on the aromatic ring [59]. The compounds with the lowest antioxidant potential in this study were 2,4-, 3,5-, and 2,6-DHB. Due to the presence of two hydroxyl groups in the meta- position to one another, they showed negligible activity at the analyzed concentrations. The most effective polyphenols include those with a hydroxyl group in para- or ortho- positions towards one another [60]. Therefore, the fact that 3,4,5-THB shows the greatest activity towards the generated DPPH radicals is influenced by the presence of three hydroxyl groups in its structure, and by two of them being in the ortho- position in relation to the third one. Mishra et al. [61] also reported the high antioxidant activity of 3,4,5-THB in the DPPH assay. The IC50 parameter after 30 min of incubation with a free radical was 4.2 μM (which gave 0.07 μmol of antioxidant/μmol of DPPH versus 0.04 μmol of antioxidant/μmol of DPPH reported herein). In the above work, the activity of 3,4-DHB and 2,3-DHB was similar to the results obtained for 3,4,5-THB, due to the presence of hydroxyl groups in their structures in the para- and ortho- positions, respectively; this was also related to them being in the ortho- position to one another, as well as the presence of an internal hydrogen bond between the -OH substituents. This facilitates the transfer of hydrogen from the -OH group of the phenols to the radical [62].
A review of the available literature indicates that there are some papers describing the antioxidant activity of hydroxybenzoic acids (Table 2), but due to the different measurement procedures used (e.g., concentrations of reagents and samples, type of solvents, reaction time, temperature) it is difficult to compare the data obtained by different researchers.

3.2. ABTS Assay

During the ABTS assay measurement, the freshly prepared 2,2-azino-bis(3-ethylbenzothiazoline-6-sulfonic acid (ABTS•+) cation radical was reduced by the antioxidant, which caused the discoloration of the blue-green solution and a reduction in the absorbance measured at 734 nm. The obtained results indicated that at the tested concentration (50 μM) the hydroxybenzoic acids showed low-to-significant scavenging properties against ABTS•+, as their average percentage of inhibition (%I) was over 60% (Figure 3). The highest antiradical activity was observed for 2,3-DHB, 2,5 -DHB, 3,4,5-THB, and 3,4-DHB, with %I values of 86.40, 80.11, 79.50, and 74.51%, respectively. These compounds also showed the best DPPH radical scavenging ability. 3,5-DHB was found to be less effective against ABTS•+ radicals, with %I = 60.39%, but much more active than 2,4-DHB and 2,6-DHB (16.17 and 8.12%, respectively). In contrast to the ABTS test, 3,5-DHB showed a negligible ability to capture the synthetic DPPH radicals.

3.3. FRAP Assay

In the FRAP assay, the Fe(III)-2,4,6-tris(2-pyridil)-s-triazine (TPTZ) complex was reduced by the HBA to Fe(II)-TPTZ at an acidic pH (3.6), and the color of the solution changed to intense blue. The increase in the absorbance at 595 nm was measured. The obtained results are summarized in Table 3. Among all tested compounds, 2,5-DHB possessed the strongest FRAP antioxidant capacity (236.00 µM Fe2+). For 2,3-DHB and 3,4,5-THB, lower reduction potentials were reported (173.79 and 158.10 µM Fe2+, respectively). The FRAP values obtained for compounds at a concentration of 50 μM increased in the following order: 44.22 (3,4-DHB) < 158.10 (3,4,5-THB) < 173.79 (2,3-DHB) < 236.00 (2,5-DHB). The above list did not include other phenolic acids, due to the negative values of the FRAP units obtained.

3.4. CUPRAC Assay

In this reducing assay, the light-blue Cu(II)-neocuproine complex was reduced by the antioxidants to yellow-orange Cu(I)-neocuproine at a close-to-physiological pH (7.4). The reaction was monitored at 450 nm. This method is based on the reduction of Cu(II) to Cu(I) by DHB and the formation of the yellow complex of Cu(I) with neocuproine (2,9-dimethyl-1,10-phenanthroline), which shows an absorption maximum at 450 nm [75]. The results of the above assay demonstrated that all tested compounds showed cupric-reducing antioxidant activity. The CUPRAC antioxidant capacity of the tested samples ranged from 20.51 to 73.85 μM of Trolox equivalents (Figure 4). The highest antioxidant capacity in the CUPRAC method was observed for 3,4,5-THB (73.85 μM of Trolox equivalents). Comparably high values were also obtained for 2,5-DHB, 2,3-DHB, and 3,4-DHB (68.77, 60.83, and 60.53 μM of Trolox equivalents, respectively). As in the DPPH, ABTS, and FRAP assays, these compounds were characterized by much higher antioxidant activity compared to the other tested compounds. 2,4-DHB, 2,6-DHB, and 3,5-DHB showed the lowest cupric-reducing antioxidant ability. According to these results, the total antioxidant capacity determined in CUPRAC assay follows the order 2,4-DHB < 2,6-DHB < 3,5-DHB< 3,4-DHB ~2,3-DHB < 2,5-DHB < 3,4,5-THB.
Summarizing the results of the four antioxidant assays based on different mechanisms of action (DPPH, ABTS•+, FRAP, and CUPRAC assays), 2,3-DHB, 2,5-DHB, 3,4-DHB, and 3,4,5-THB showed the strongest antioxidant activity (Table 2 and Table 3, Figure 3 and Figure 4). The higher antioxidant activity of these DHBs can be explained by the resonance structures of their radicals [76]. In the case of DHBs with the -OH groups located in the ortho- and para- positions to one another, both oxygen atoms from the -OH groups participate in the delocalization of the electronic charge in the radical molecules. Meanwhile, the -OH substituents being in the meta- position of in the ring limits the delocalization of the charge with the participation of the second -OH [77]. 3,4,5-THB showed the highest antioxidant activity in the DPPH and CUPRAC assays; it contains the greatest number of hydroxyl (-OH) groups (three) bonded to the benzoic ring in an ortho- position to one another. It was previously shown that the antioxidant activity of phenolic compounds increases with the increasing number of hydroxyl groups substituted in the aromatic ring [39]. The same relationship was observed in the study of Sroka and Cisowski [69], where 3,4,5-THB showed a higher DPPH radical scavenging ability than 3,4-DHB and, 2,3-DHB which each contain two hydroxyl groups in their structures [69]. It should be stressed that, in our study, 3,4,5-DHB showed the highest antioxidant activity in the DPPH and CUPRAC assays, but not in the ABTS or FRAP assays. This may be related to differences in the mechanisms of action and reaction environments of the individual methods. The tests conducted in the above study can be divided into two types of reaction mechanisms: in the DPPH and ABTS•+ assays antioxidants can react with radicals by a combination of SET (single electron transfer) and HAT (hydrogen atom transfer) mechanisms (i.e., SPLET mechanism), whereas the FRAP and CUPRAC assays are based only on the SET reaction mechanism [57]. Moreover, the results of reactions based on the SET mechanism are related to the environment, e.g., the type of solvent and the pH value of the solution [78]. Other factors, such as the acidity of the studied molecules, may also affect the reaction pathway. The ABTS and FRAP assays are conducted predominantly in water solutions, where the hydroxybenzoic acids may be more ionized than in methanolic or ethanolic solutions. Phenoxide anions oxidize much more easily and faster than corresponding neutral molecules of phenols. This may accelerate the rate of electron transfer reaction in the SPLET mechanism. Moreover, the final results of ABTS assay may be affected by the reduction of the ABTS•+ by phenolic compounds to the parent substrate ABTS, the degradation of ABTS•+ to nonreactive products, or the formation of additional adducts between phenolics and ABTS•+ [79,80].
In conclusion, the experimental results show the following: (1) the antioxidant activity of the studied phenolic acids is directly related to the number and position of hydroxyl –OH groups bonded to the benzene ring; (2) hydroxybenzoic acids with hydroxyl groups substituted in the ortho- (2,3-DHB and 3,4-DHB) and para- (2,5-DHB) positions show higher antioxidant activity; (3) the hydroxyl groups substituted in the meta- (2,4-DHB and 3,5-DHB) position reduce the antioxidant and antiradical properties of hydroxybenzoic acids; and (4) the reaction environment and differences in the mechanisms of action of individual methods affect the activity of the compounds and the test results.

3.5. Pro-Oxidant Activity

HBAs were tested for their pro-oxidative effects in the Trolox oxidation assay. The radicals of HBAs were produced in their reaction with H2O2 catalyzed by horseradish peroxide. Then, the phenoxyl radicals react with Trolox, which undergoes oxidation to Trolox radicals, and then Trolox quinones, whereas the phenoxyl radicals are transformed to phenolic compounds. The maximum absorption for Trolox quinone is 272 nm. The pro-oxidant activity of the tested benzoic acid derivatives was measured at two concentrations (0.15 and 0.35 μM) and expressed as the rate of Trolox oxidation; the results are presented in Figure 5. 2,6-DHB, 3,4-DHB, and 3,5-DHB showed the highest pro-oxidant activity (degree of Trolox oxidation after 120 min of incubation for compounds at the concentrations of 35 μM: 9.52, 8.44, and 8.03%, respectively). The pro-oxidant activity of 2,4-DHB, 2,5-DHB, 2,6-DHB, and 3,4-DHB increased with an increase in their concentration; in the case of the other studied compounds, this cannot be clearly stated. The phenolic acids that possessed the strongest activity in the ABTS•+, FRAP, and CUPRAC assays (i.e., 2,3-DHB, 2,5-DHB, and 3,4,5-THB) showed the lowest pro-oxidant properties. Moreover, a low value of Trolox oxidation was obtained for 2,4-DHB.

3.6. Antimicrobial Activity

The antibacterial and antifungal activities of hydroxybenzoic acids were assayed using the microdilution method of minimum inhibitory concentration (MIC). The MIC values of the tested phenols for all of the tested strains were between 2 and 6 mg/mL (Table 4). The studies revealed that 2,4-DHB and 3,4-DHB possessed the strongest antimicrobial properties towards E. coli, P. aeruginosa, S. aureus, B. subtilis, S. enteritidis, and C. albicans at the concentration of 2 mg/mL. Similarly good results were obtained for 3,5-DHB, with only one higher MIC (3 mg/mL) for P. aeruginosa. 2,5-DHB, 2,6-DHB, and 3,4,5-THB presented the same antimicrobial activity against E. coli, B. subtilis and S. enteritidis at the concentration of 3 mg/mL. The lowest antimicrobial properties against the tested microorganisms were reported for 2,3-DHB (MIC = 5 mg/mL).
The antimicrobial activities of hydroxybenzoic acids have been widely described in previous research. In a study performed by Borges et al. [81], 3,4,5-THB (gallic acid) showed antimicrobial activity (MIC = 0.5–2.0 mg/mL) against E. coli, P. aeruginosa, S. aureus, and L. monocytogenes [81]. Sanchez-Maldonado et al. [82] obtained the MIC values for 3,4,5-THB and 3,4-DHB against L. plantarum, L. hammesii, E. coli, and B. subtilis, which ranged from 0.49 to 4.56 mg/mL for gallic acid and from 0.31 to 3.87 mg/mL for protocatechuic acid; these authors found that the increase in the number of -OH groups increased the growth-inhibitory activity of hydroxybenzoic acid, but this is not reflected in the results of our work [82]. In another study, 3,4-DHB and 3,4,5-THB isolated from T. javanicum was found to have significant activity against Gram-negative bacteria (E. coli, P. aeruginosa, V. cholera, S. typhi, and S. typhimurium) and some Gram-positive bacteria (M. luteus and S. aureus) [66]. Stojkovic et al. [83] reported the antibacterial activity of protocatechuic acid against S. aureus, B. cereus, M. flavus, P. aeruginosa, E. coli, and E. cloacae. The MIC value for all mentioned strains was 3.0 mg/mL [83]. 2,4-DHB (β-resorcylic acid) has been found to exhibit significant antimicrobial properties against various pathogens, including Salmonella [84], E. coli [85], and L. monocytogenes [86]. In the study of Jayesh et al. [87], 2,3-DHB electrospun into a nanofiber blend of poly(D,L-lactide) and poly(ethylene oxide) inhibited biofilm formation by P. aeruginosa [87]. Moreover, 2,3-DHB isolated from the fruit extract of Flacourtia inermis Roxb exhibited antifungal activity; at a concentration of 30 mg/mL it inhibited the growth of A. fumigatus, A. flavus, A. niger, and tested Chrysosporium species (zone of inhibition: 30, 24, 22, and 25 mm, respectively) [88]. Lattanzio et al. [89] reported that 2,5-DHB (gentisic acid) showed antifungal activity; at a concentration of 5 mM, it completely inhibited spore germination and mycelial growth of Botrytis cinerea and Rhizopus stolonifer [89]. In the work of Merkl et al. [30], 2,5-DHB inhibited the growth of E. coli, B. aureus, and L. monocytogenes (MIC = 2.5 mM) [90]. According to Campos et al. [91], the mechanism of the antimicrobial activity of phenolic acids may be related to their partially lipophilic character. It can be concluded that the protein denaturation is caused by an undissociated form of phenol that passes through the cell membrane via passive diffusion, disrupting its structure and possibly acidifying the cytoplasm [91].

3.7. Cytotoxicity Studies in MCF-7 and MDA-MB-231 Cell Lines

The effects of seven hydroxybenzoic acids on two types of human breast cancer cell lines (MCF-7 and MDA-MB-231) were investigated; the results are shown in Table 5 and Figure 6. Among the tested compounds, only 3,4,5-THB (gallic acid) turned out to be toxic to the MCF-7 cell line. The IC50 parameter (concentration of the sample that inhibited cell viability by 50% compared to control cells) for gallic acid was 0.44 mM. The other tested hydroxybenzoic acids were found to be nontoxic to MCF-7 cells. 3,4,5-THB was found to be the most toxic among the tested compounds (IC50 = 0.36 mM) towards both MDA-MB-231 (estrogen-independent) and MCF-7 human breast cancer cells. At a concentration of 5 mM, 3,4,5-THB inhibited the growth of MCF-7 and MDA-MB-231 cells almost completely. Compared to the MCF-7 cell line, the MDA-MB-231 cell line turned out to be more sensitive to the tested hydroxybenzoic acids. Except for 3,4,5-THB, the studied compounds were only slightly toxic to the cancer cell lines. For example, the IC50 parameter for 2,4-DHB was 4.77 mM. 3,4-DHB and 3,5-DHB turned out to be nontoxic to both the MDA-MB-231 and MCF-7 cell lines.
Numerous studies also indicate the antitumor potential of gallic acid. For example, in the study of Hsu et al. [92], 3,4,5-THB at a low concentration (5 μg/mL) significantly reduced the cell proliferation of MCF-7 to 42.73% of control (incubation time: 24 h) [92]. In another study [93], 3,4,5-THB was found to inhibit 50% of MCF-7 cancer cells at a concentration of 7.5 μg/mL (incubation time: 24 h) [93]. Maurya et al. [94] tested the effect of gallic acid on the A549 human lung adenocarcinoma cell line. Exposure to gallic acid for 24 h inhibited the A549 cell growth and decreased cell viability [94]. Moghtaderi et al. [95] reported that 3,4,5-THB inhibits the growth of MDA-MB-231 cells in a dose-dependent manner (IC50 ≈ 150 μM in 24 h); the combination of gallic acid with curcumin enhanced the apoptotic effect [95]. Another study [96] showed that 3,4,5-THB has anticancer activity against A2780 cancer cells (IC50 = 45.25 μg/mL; incubation time: two days), and this antiproliferative effect was greater for the synthesized derivatives of 3,4,5-THB (IC50 < 39 μg/mL) [96]. This fact may be related to the presence of alkyl chains, which increase the lipophilicity of the compound (easier transfer through the cell membrane) [97].
Other hydroxybenzoic acids that exhibited anticancer effects on the MDA-MB-231 cell line in our work were 2,3-DHB (IC50 = 4.10 mM) and 2,5-DHB (IC50 = 4.39 mM) (after 48 h of incubation); these compounds are often referred to as aspirin metabolites. Sankaranarayanan et al. [98] also reported their chemopreventive effects against that breast cancer line. 2,3-DHB was found to significantly inhibit the formation of colon cancer cells at a concentration of ~500 µM in MDA-MB-231 cells, whereas 2,5-DHB inhibited the colony at a concentration of 100 µM (incubated for 48 h) [98]. The effect of 2,5-DHB on MCF-7 cell proliferation was also examined in our previous study [62].
In our study, after 48 h of incubation, 3,4-DHB and 3,5-DHB did not show cytotoxic activity in the MCF-7 and MDA-MB-231 cell lines, but in the antimicrobial assay they possessed the strongest properties. In previous literature, protocatechuic acid has been referred to as a potent anticancer agent. In the study of Yin et al. [28], 3,4-DHB decreased the viability of MCF-7 human breast cancer cells, A549 lung cancer cells, HepG2 cells, HeLa cervical cancer cells, and LNCaP prostate cancer cells [28]. Yip et al. [99] reported that 3,4-DHB at a concentration of 100 µM inhibited HepG2 hepatocellular carcinoma cells’ viability to 40% (incubation time: 4 days) [99]. The results obtained by Reis et al. [100] showed that Leccinum vulpinum Watling extract—rich in protocatechuic acid (77 µg/g extract)—decreased cellular proliferation and induced apoptosis in MCF-7 human breast cancer cells [100].

3.8. Lipophilicity of Hydroxybenzoic Acids

The chromatographic lipophilicity parameters (logkw) obtained under different chromatographic conditions are shown in Table 6. The alkyl (C18 and C8)-modified silica are stationary phases commonly applied for the estimation of lipophilicity by the use of HPLC, where the separation mechanism is mainly based on hydrophobic (van der Waals) interactions. In the case of the phenyl-modified silica phase (PHE), the π⟶π interactions are assumed. In the case of a more polar cyano-bonded phase (CN), the additional hydrogen bonds between cyano groups and a hydrogen atom from the hydroxyl group exist. The IAM (immobilized artificial membrane) stationary phase is used to mimic the distribution of compounds in the phospholipid membrane. The most distinct differences between particular values of logkw were obtained for the IAM column, whereas the lowest differences were observed in the case of the CN and PHE columns. The highest linear correlation was obtained for the logkw obtained by the use of the C18 and PHE columns (R2 = 0.772). The experimental logP values for HBAs (Table S1 and Table 7) highly correlate with the logkw obtained on C18 and the theoretical logP classic, logP galas, pKa1, and pKa2 parameters. This means that both theoretical and chromatographic logP parameters can be successfully used to predict the hydrophilic–lipophilic properties of the hydroxybenzoic acid series. The obtained (experimental and theoretical) logP values for HBAs are in the range 0 < logP < 3; this means that these compounds possess from slightly hydrophilic (logP: 0–1) to moderately lipophilic properties (logP: 1–3) [101]. Generally, the disubstituted hydroxybenzoic acids were more lipophilic than the trisubstituted one (gallic acid has hydrophilic properties, and dissolves well in water and other polar solvents). According to the experimental logP values, the compounds can be ordered by increasing lipophilicity as follows: 3,4,5-THB→3,4-DHB~3,5-DHB→ 2,3-DHB→2,4-DHB →2,5-DHB →2,6-DHB. Dissociation and association phenomena affect the value of the partition coefficient. Therefore, we observed a relationship between the acid dissociation constants (expressed as pKa) and the values of partition coefficients. Each functional group that can be a hydrogen bond donor or acceptor increases the hydrophilic nature of the compound. Hydroxyl and carboxylic groups present in the structures of the compounds can form hydrogen bonds with water molecules in the aqueous environment, which affects their solubility in water. Therefore, in the case of 3,4,5-trihydroxybenzoic acid, an increase in the hydrophilic properties of the acid is noticeable. The more hydrogen bonds can be formed between a molecule and water molecules, the greater its solubility in water. In the case of ortho- substituted benzoates, the -OH group in the ortho- position is mostly engaged in hydrogen bond formation with the -COOH group, facilitating dissociation of the H+ carboxylic group and increasing the acidity of 2,6-DHB (pKa = 1.30) and, to a lesser extent, 2,3-DHB (pKa = 2.91), 2,5-DHB (pKa = 2.97), and 2,4-DHB (pKa = 3.11). In the case of 3,5-DHB (pKa = 4.04), 3,4-DHB (pKa = 4.26), and 3,4,5-THB (pKa = 4.40), the -OH groups are too far from -COOH to interact with this functional group, but may form hydrogen bonds with the solvent. As a result, the acidity of these acids decreases (higher values of pKa) compared to the ortho-substituted benzoates. The dependency between acidity (pKa exp.) and lipophilicity (logPexp.) in the series of studied hydroxybenzoic acids is shown in Figure 7. With the increase in the acidity of the compounds, the lipophilicity increases as well. Therefore, the hydroxybenzoic acids can be divided into three groups characterized by: (a) lower acidity and lipophilicity (3,4,5-THB; 3,4-DHB; 3,5-DHB), (b) moderate acidity and lipophilicity (2,3-DHB; 2,4-DHB; 2,5-DHB), and (c) higher acidity and lipophilicity (2,6-DHB). There is no strict dependency between lipophilicity or acidity of hydroxybenzoic acids and their biological properties (antioxidant, antimicrobial, and cytotoxic), but general conclusions can be drawn. Dihydroxybenzoic acids with the highest antimicrobial activity are characterized by low-to-moderate lipophilicity (logPexp. in the range 0.86–1.63), and acidity at the level of pKa 3.11–4.26 (Table 7). The strongest cytotoxic activity was observed in trihydroxybenzoic acid, with the lowest lipophilicity and acidity in the group of studied compounds.

3.9. DFT Study

The energy levels of the highest occupied molecular orbitals (EHOMOs) and the lowest unoccupied molecular orbitals (ELUMOs) for the optimized structures of the seven selected hydroxybenzoic acids were determined using the B3LYP/6-311++G(d,p) calculation method in the Gaussian 09W software. Using these data, other energy parameters such as HOMO–LUMO energy gap (ΔE), ionization potential (IP), electron affinity (A), electronegativity (χ), electrochemical potential (μ), electrophilicity index (ω), and chemical hardness (η) and softness (σ) were calculated, and the results are summarized in Table S2. The energy value of the HOMO orbital is a measure of the ability of one molecule to donate electrons to another, while the energy value of the LUMO orbital determines the tendency of the molecule to accept electrons. The lower the ionization potential of a compound (the minimum energy that is needed to detach an electron from the molecule) and the higher the energy of HOMO, the better the antioxidant (electron donor). The chemical activity of the substance also depends on other parameters, e.g., the value of the HOMO–LUMO energy gap [111,112]. The lower the value of ΔE, the easier the electron enters the excited state, and the better the antioxidant. Many physicochemical descriptors, such as those listed above, are often used to evaluate the biological and antioxidant activity of phenolic compounds [113]. According to the calculations (Table S2), we can see that the EHOMO values for three phases did not show significant differences; the values of EHOMO in the gas phase ranged from −9.2813 to −8.7735 eV, in aqueous solution from −9.2753 to −8.7555 eV, and in methanolic solution from −9.2755 to −8.7550 eV. Moreover, taking into account the energy of the HOMO orbitals, we can classify the studied compounds into two groups: group I—better electron-donor molecules (phenolic acids that presented the highest EHOMO values, above −9 eV), i.e., 3,4,5-THB, 2,5-DHB, 2,3-DHB, and 3,4-DHB; and group II—better electron-acceptor molecules (phenolic acids that presented the lowest EHOMO values, below −9 eV), i.e., 2,4-DHB, 3,5-DHB, and 2,6-DHB. In this work, the determination of the antioxidant activity of phenolic compounds was carried out in two environments: methanol (DPPH) and water (ABTS•+, FRAP, CUPRAC). Hydroxybenzoic acids from group I showed better antioxidant activity in the above-mentioned assays than compounds from group II. According to the calculated parameters, 3,4,5-THB—with the highest energy of HOMO (-8.7550 and -8.7555 eV) and electronic chemical potential (−7.0637 and −7.0638 eV), and with the lowest values of ionization potential (8.7550 and 8.7555 eV) and electronegativity (−7.0637 and 7.0638 eV)—should be the best antioxidant among the tested phenolic acids (data for methanolic and aqueous solutions, respectively). In our study, gallic acid was found to have the highest antioxidant activity only in the DPPH and CUPRAC assays. Moreover, the relationship between the increase in antioxidant properties and the increase in the compound’s ability to donate electrons was observed only in the CUPRAC test. The energy values of the HOMO orbitals of hydroxybenzoic acids in the aqueous phase increased in the following order: 2,4-DHB (−9.2753 eV) < 2,6-DHB (−9.2475 eV) < 3,5-DHB (−9.2290 eV) < 3,4-DHB (−8.9340 eV) < 2,3-DHB (−8.9153 eV) < 2,5-DHB (−8.7749 eV) < 3,4,5-THB (−8.7555 eV).

3.10. Aromaticity

According to the values of the calculated geometric indices of hydroxybenzoic acids in which the hydroxyl groups were substituted in position 2,n-(where n = 3, 4, 5, 6), aromaticity increased in the sequence 2,3-DHB < 2,4-DHB < 2,5-DHB < 2,6-DHB (Table S3). This is indicated by the values of all indices calculated both for molecules modeled in the gas phase and in solvents (water, methanol). The values of the calculated aromaticity indices in water and methanol were very similar to one another. The aromaticity of the studied molecules increased with the change in the substitution site of the second carboxyl group (moving away) relative to the carboxyl group substituted in position 2. In the case of 3,4-DHB and 3,5-DHB, it was also observed that the molecules in which the hydroxyl groups were further apart had higher aromaticity. According to the calculated values of the aromaticity indices, among all tested systems, 3,5-DHB acid showed the highest aromaticity. High values of aromaticity indices were also noted for 3,4,5-trihydroxybenzoic acid (slightly lower than for 3,5-DHB and 2,6-DHB). When comparing the aromaticity of all studied molecules, they can be ordered by increasing aromaticity as follows: 2,3-DHB < 2,4-DHB < 2,5-DHB < 3,4-DHB < 2,6-DHB < 3,4,5-THB < 3,5-DHB. The NICS magnetic aromaticity index is negative for aromatic compounds; the lower its value, the higher the aromaticity of a given molecule. The highest aromaticity according to the NICS calculations was demonstrated for 3,4,5-trihydroxybenzoic acid, and the lowest for 2,4-DHB. The values were similar for calculations carried out in the gas phase and in solvents (water and methanol).

3.11. The Distribution of Electronic Charges

The electron charge distribution in the analyzed molecules was calculated using the CHelpG and NBO methods. The charge values are presented in Table S4. The electronic charge values correlate well with the values of chemical shifts of 13C NMR (presented in the next part of the manuscript). The positive electronic charges (NBO and CHelpG) were gathered around the carbon atom from the ring substituted by the -OH group (Figure S3). The highest electronic density was reported for the C1 carbon atom in 2,4- and 2,6-dihydroxybenzoic acids, dropped slightly for 2,3- and 2,5-dihydroxybenzoic acids, and was much lower for 3,4- and 3,5-dihydroxybenzoic acids and 3,4,5-trihydroxybenzoic acid. The electronic density around the second carbon atom became positive due to the pulling out of electrons by the highly electronegative oxygen atom of the hydroxyl group in 2,n-dihydroxybenzoic acids, while in the remaining acids these values were negative. The electronic density around the sixth atom was negative in almost all acids, except for 2,6-dihydroxybenzoic acid. In the case of the carbon of the carboxyl group (C7), the electronic density was positive in all analyzed acids, while for 2,6-dihydroxybenzoic acid it was much lower than in the other structures. This acid is characterized by high values of aromaticity indices, similar to 3,5-dihydroxybenzoic acid, in which the value of electronic density around the carbon atom of the carboxylic group is also lower than in other molecules.
The symmetric position of the hydroxyl groups in the aromatic ring with respect to the carboxylic group significantly influences the decrease in the electronic density of the carbon atom of the carboxylic group, as well as the increase in symmetrization of the electronic charge distribution in the aromatic ring and, consequently, the increase in the aromaticity of a given system. Moreover, the calculated NBO and CHelpG total electronic charges of the ring ( ē ring) increase in the order 2,6- < 3,5- < 2,4- < 2,3- < 3,4- < 2,5-DHB < 3,4,5-THB. For all dihydroxybenzoates, the value of the ē ring was negative, whereas for 3,4,5-THB the total electronic charge of the ring had a positive value. Surprisingly, the increase in the ē ring correlates well with the increase in the antioxidant activity of the studied molecules (specifically, it is in line with the results of the DPPH and CUPRAC assays). For 2,6-DHB (the molecule with the lowest antioxidant activity), the lowest electronic charge in the ring and the highest electronic charge in the COO were reported. The NBO charges calculated for the oxygen of the -OH group for 3,4,5-THB and 2,3-, 2,5-, and 3,4-DHB were more negative than for the other DHBs, facilitating the donation of protons in the reaction with radicals (the ionization potential also points to easier electron donation by 3,4,5-THB and 2,3-, 2,5-, and 3,4-DHB; Table S2).

3.12. ESP Maps

The electrostatic potential map shows the areas of a molecule related to its electrophilic (red) and nucleophilic (blue) reactivity (Figure 8). In 2,3-, 2,4-, 2,5-, and 2,6-dihydroxybenzoic acids, it was observed that the oxygen atom of the carbonyl group was highly susceptible to electrophilic attack. As the second hydroxyl group moves away from the substitution position of the first hydroxyl group (the second position on the aromatic ring), the charge density in the region of the carboxylic group increases. Moreover, the area of the oxygen atom of the hydroxyl group (from the second position) is characterized by high electrophilic reactivity. The second position of the hydroxyl group in these acids is characterized by a significant susceptibility to nucleophilic attack, and as this hydroxyl group is moved away from the position of the first hydroxyl group, its nucleophilic reactivity increases. In 3,4-DHB, 3,5-DHB, and 3,4,5-THB, the electrophilic reactivity of the carbonyl group was lower than in the remaining acids. Moreover, the electrophilic reactivity of the hydroxyl group located closer to the carbonyl group was lower than in the other acids. In molecules with a symmetrical distribution of hydroxyl groups in the aromatic ring with respect to the carboxyl group—i.e., in 2,6-DHB and 3,5-DHB—a symmetric distribution of electrostatic potential was observed. These molecules are characterized by the highest aromaticity among the tested compounds. The distribution of the electrostatic potential for the molecules optimized in water and methanol was analogous to that in the molecules tested in the gas phase.

3.13. FTIR, 1H, and 13C NMR Study

The wavenumbers, intensities, and assignments of the bands occurring in the FTIR spectra of dihydroxybenzoic acids and 3,4,5-trihydroxybenzoic acid are presented in Table S5 [114]. The bands were numbered along with the notation used by Varsányi [115]. In the spectra of hydroxybenzoic acids (Figure S1), the bands assigned to the stretching vibrations of the carbonyl group νC=O were located in the range of 1639–1704 cm−1. These were very intense and characteristic bands, whose location in the spectra was affected by the number and position of hydroxyl substituents in the ring. In the range of 3214–3375 cm−1, the bands derived from -OH stretching occurred. The narrow band at 3374 cm−1 corresponded to the free, unassociated -OH group present only in the spectra of 2,4-DHB. The hydrogen atom of the -OH group can easily form intramolecular hydrogen bonds with the oxygen atoms of other functional groups, which affects the shape and position of the absorption bands assigned to the -OH stretching vibrations; this causes widening of the bands. The formation of hydrogen bonds leads to the formation of dimers and polyassociates. The spectra of carboxylic acids were characterized by very wide bands (with 3–5 maxima) assigned to the stretching of the -OH group (related to the presence of intermolecular hydrogen bonds), the maxima of which were in the range of 2500–2878 cm−1. In the spectra of all studied hydroxybenzoic acids, strong bands assigned to the stretching vibrations of νC-OH occurred in the range of 1194–1263 cm−1, and deformations of the -OH group at 1385–1446, 1208–1241 (in-plane β(OH)), and 576–633 (out-of-plane ν(OH)). Based on the FTIR spectra, the stabilization/destabilization of the aromatic system of ligands can be discussed. The destabilization of the aromatic system occurs via a decrease in the wavenumbers as well as the number and intensities of the bands assigned to the vibrations of the aromatic ring [116,117]. In the spectra of the 2,n-dihydroxybenzoic acids (n = 3,4,5,6), the increase in the wavenumbers of selected bands derived from the vibrations of the aromatic ring (no. 8a, 19a, 14, 18b, 7b, 9b) was observed in the following sequence: 2,3-DHB < 2,4-DHB < 2,5-DHB < 2,6-DHB. This may show the increase in the aromaticity of the acids in the same sequence (the calculation of the aromaticity indices confirmed the observation). The analysis of the bands assigned to the aromatic ring vibrations in the spectra of 3,n-dihydroxybenzoic acids (n = 4), revealed that 3,5-DHB showed higher aromaticity than 3,4-DHB. This can be concluded based on the higher wavenumber values of the aromatic bands and the occurrence of some additional aromatic bands in the spectrum of 3,5-DHB compared with the spectrum of 3,4-DHB. The highest wavenumber values of some aromatic bands were observed in the spectrum of 2,6-DHB, compared with other hydroxybenzoic acids, which may suggest the highest stability of the aromatic system in the case of this acid.
The change in the positions of signals from the 1H and 13C NMR spectra of hydroxybenzoic acids (Figure S2) showed that the change of the substituted position in the ring affects the electronic charge distribution in the molecule and the aromaticity of the system (Table S6). Strongly electronegative oxygen atoms of the -OH groups decrease the π-electron density at the nucleus, deshielding the nucleus and causing a larger chemical shift in the 13C NMR. The values of the chemical shifts of the carbon atoms substituted by the -OH group are highlighted in bold in Table S6. The highest decrease in the electronic charge was observed in the case of the second and fourth carbon atoms in 2,4-DHB, the second and fifth in 2,5-DHB, and the third and fourth in 3,4-DHB. These acids may reveal a greater ability to donate electrons or protons to radical molecules and, therefore, may show higher antioxidant activity than the other hydroxybenzoic acids.
The number and position of the -OH substituents in the ring affect the -COOH substitution capability by metal cations, as well as the strength of the hydroxybenzoic acids. 3,4- and 3,5-DHB, as well as 3,4,5-THB, are weaker acids than the other studied compounds.
The electronic density around the C7 atom in 3,4- and 3,5-DHB, as well as 3,4,5-THB, is higher than in the other acids. Significant differences in the electronic density values around the C1 carbon atom were observed between the 2,n-dihydroxybenzoic acids and the 3,n-dihydroxybenzoic acids. In the case of the former, the electronic density was much higher (lower values of the chemical shifts in the 13C NMR spectra) than in the latter group of acids, as a result of the substitution site of the hydroxyl groups. As the position of the hydroxyl groups moves away from the carboxyl group, the electronic density around the carbon atom of the carboxyl group decreases (the value of 13C chemical shifts increases).
The comparison of the values of chemical shifts on the 13C NMR spectra (differences between the maximum value of the chemical shift and the lowest value of the chemical shift within a given molecule) of the tested acids shows that the most stable distribution of electronic charge occurs in the aromatic ring of 2,5-dihydroxybenzoic acid. The least stable electron charge distribution of the aromatic ring was observed in 2,6- and 3,5-dihydroxybenzoic acids. The values of chemical shifts were compared with the electronic charge distribution calculated by the NBO method. For electron densities determined theoretically by the NBO method and the experimental values of 13C on the C1, C2, C5, and C6 atoms, a high correlation R > 0.9 was noted.
The electronic stability of the aromatic ring (aromaticity of the π-electron system under study) can be assessed by analyzing the values of chemical shifts of aromatic protons in the 1H NMR spectra. Along with the change in the position of the hydroxyl groups, a change in the value of the chemical shifts of the protons attached to the aromatic ring was observed. The increase in the value of chemical shifts of aromatic protons in the 1H NMR spectrum indicates an increase in the aromaticity of the benzene ring. In the case of the tested acids, due to the different sites of substitution of the hydroxyl groups in the ring, the assessment of aromaticity of these compounds based on the criterion of the value of protic chemical shifts is ambiguous.
Changes in the values of chemical shifts of the proton of the carboxyl group and the protons of the hydroxyl groups in individual acids were also observed. The highest values of chemical shifts of hydroxyl groups were observed in the 1H NMR spectrum of 2,4-dihydroxybenzoic acid, and the lowest for 2,5-dihydroxybenzoic acid. The electronic density around the protons of the hydroxyl groups according to the values of the proton chemical shifts observed in the 1H NMR spectra increases as follows: 2,6-DHB < 2,4-DHB < 3,5-DHB < 3,4-DHB < 2,3-DHB < 2,5-DHB < 3,4,5-THB. The highest electronic density was observed around the protons of the hydroxyl groups in 3,4,5-THB and 2,5-DHB; these acids showed the highest antioxidant potential in the DPPH, CUPRAC, and FRAP tests.
In the case of the proton of the carboxyl group, higher values of chemical shifts were observed in the 1H NMR spectra of 3,4 and 3,5-dihydroxybenzoic acids than in the case of other acids.

3.14. Statistical Analysis

The relationships between the lipophilic and electronic parameters, aromaticity indices, FTIR and NMR parameters, and anti-/pro-oxidant and antimicrobial potential of the selected dihydroxybenzoates were evaluated by principal component analysis. The list of variables correlated with one another is summarized in Table S7. Predominant highly significant correlation coefficients between variables proved the application of PCA to be justified and promising; however, lipophilic and electronic parameters, aromaticity indices, and FTIR and NMR parameters (those that significantly change position in the spectra depending on the number and location of the substituents in the ring) were used as input variables to the PCA, while anti-/pro-oxidant and antimicrobial potential were used as accompanying variables. According to the PCA and the Kaiser criteria, three principal components with eigenvalues > 1 were obtained, accounting for 89.9% of the total variance. PC1, with an eigenvalue equal to 8.16 explaining 51.0% of the total variance, was contributed to by logPGalas, logPexp., δC2, logPC18, and ΔE(LUMO-HOMO) with negative correlation, and by pKa1, γ(CH) 11, δC3, νC-(OH), and β(CH) 18b with positive correlation (Table S8). Based on the characteristics contributing to the first principal component, this reflects the hydrophilic–hydrophobic and spectral characteristics of dihydroxybenzoates. PC2, with an eigenvalue equal 3.61, explains around 22.6% of the total variance, and shows directly proportional correlation between I6, HOMA, BAC, IP, and ΔE(LUMO-HOMO). This could be accepted as a component indicating some aromaticity due to the energy difference between the frontier orbitals. PC3, with an eigenvalue 2.61 explaining 16.3% of the total variance, was contributed to by positive correlation between NICS, νC-(OH), and IP, as well as negative correlation between BAC and NBO, mainly reflecting the atomic charge and electric characteristics. Biplots of factor loadings and sample scores of all combination of PCs with corresponding anti-/pro-oxidant and antimicrobial potential are presented on Figure 9a–d, respectively. The antioxidant activity of the studied compounds described by the CUPRAC and ABTS parameters correlates with PC1, DPPH parameters (IC50) correlate with PC3, whereas pro-oxidant activity correlates with PC2. The antioxidant parameters obtained via the ABTS and CUPRAC assays negatively correlate with the DPPH assay’s parameter. Compounds with high scavenging activity against DPPH (lower IC50 values) possess high antioxidant activity as measured by the ABTS and CUPRAC assays. Moreover, the pro-oxidant activity of HBAs is negatively correlated with antioxidant properties (compounds with pro-oxidant activity—2,6-, 3,5-, and 3,4-DHB—are weaker antioxidants compared with the other HBAs). The antimicrobial properties of HBAs correlate with PC2 and PC3. 2,3- and 2,5-DHB possess the lowest antimicrobial activity (compared with other HBAs, as expressed by their generally having the highest MIC values), and are characterized by lower IP, ΔE(LUMO-HOMO), and BAC values than the other HBAs. PCA analysis showed differences between the HBAs that can be divided into groups taking into account the parameters used for calculations of the PCs: (a) 3,4,5-THB, 3,4-DHB, and 3,5-DHB (higher values of spectral parameters, lower values of logP, higher pKa), (b) 2,5-DHB, 2,3-DHB, and 2,4-DHB.

4. Conclusions

The number and position of hydroxyl substituents in the benzoic acid moiety strongly affect the electronic charge distribution over the molecules, as well as their biological properties. PCA analysis was successfully applied to describe the structure–activity relationship in the group of HBAs. Three principal components described 89.9% of the total variance. The applied variables describing the structure of HBAs included lipophilic and electronic parameters, aromaticity indices, NBO charges, and spectral (FTIR and NMR) parameters. The antioxidant activity—described by parameters obtained via the DPPH, ABTS, FRAP, and CUPRAC assays—strongly correlates with the calculated electronic parameters (Table S2) and NBO and CHelpG charges. With the decrease in the total electronic charge of the aromatic ring and the increase in electronic charge gathered on the -OH group, the electron-/hydrogen-donating properties of molecules also increase. Because the antioxidant assays were based on different mechanisms and factors, such as type of solvent, pH may affect the results. The PCA analysis showed that the CUPRAC and ABTS parameters correlate with PC1, the DPPH parameter (IC50) correlates with PC3, and pro-oxidant activity correlates with PC2. This suggests that different structural parameters can play a significant role in assessing antioxidant activity, depending on the test used. Therefore, it is not possible to rank the HBAs according to antioxidant activity, but taking into account the obtained parameters it can be stated that 3,4,5-THB and 3,4-, 2,3-, and 2,5-DHB showed the highest antioxidant activity. The position of the hydroxyl groups in the meta- position in relation to one another reduces the antioxidant and antiradical properties of hydroxybenzoic acids. 2,4-, 3,4-, and 3,5-DHB showed the highest antimicrobial activity against E. coli, P. aeruginosa, S. aureus, B. subtilis, S. enteritidis, and C. albicans. These compounds are characterized by low-to-moderate lipophilicity (logPexp. in the range 0.86–1.63) and acidity at the level of pKa 3.11–4.26. Two of them—3,4- and 3,5-DHBs—showed pro-oxidant activity in the Trolox oxidation test, but were non-cytotoxic to the MCF-7 and MDA-MB-231 cell lines. The highest cytotoxicity was observed in 3,4,5-THB (IC50 = 0.44 ± 1.80 and IC50 = 0.36 ± 0.50 respectively for MCF-7, and MDA-MB-231), which also possessed the highest antioxidant activity in the DPPH and CUPRAC assays, and no pro-oxidant activity in the Trolox oxidation test. 3,4,5-THB had the lowest lipophilicity and acidity among the studied compounds. Taking into account these two parameters, the HBAs can be divided into three groups that are characterized by: (a) lower acidity and lipophilicity (3,4,5-THB; 3,4-DHB; 3,5-DHB), (b) moderate acidity and lipophilicity (2,3-DHB; 2,4-DHB; 2,5-DHB), and (c) higher acidity and lipophilicity (2,6-DHB).
The HBAs commonly present in plant-based foods have significant biological properties. Some of them can be applied as strong antioxidants, antimicrobials, or even as chemopreventive or anticancer agents. Their mechanism of action is still not elucidated, and requires further study. It is worth undertaking research on the mechanism of action of HBAs at the cellular level, because some of them have the potential to be used as dietary supplements, functional food ingredients, or in the prevention and treatment of diseases.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/nu13093107/s1: Figure S1: The FTIR spectra of hydroxybenzoic acids; Figure S2: 1H and 13C NMR spectra of hydroxybenzoic acids; Figure S3. The atom numbering scheme; Table S1: Pearson’s correlation matrix (p < 0.05) of different theoretical and experimental logP and pKa parameters; Table S2: Calculated parameters of chemical reactivity of tested hydroxybenzoic acids in gas, aqueous, and methanolic phases; Table S3: The values of aromaticity indexes of the tested acids (Aj, BAC, HOMA, GEO, EN, I6, NICS) calculated for structures optimized in the gas and aqueous phases using the B3LYP/6-311++G(d,p) method; Table S4: The distribution of electronic charges (calculated via ChelpG and NBO) for dihydroxybenzoic acids and trihydroxybenzoic (gallic) acid; Table S5: Wavenumbers (cm−1), intensities, and assignments of bands occurring in the experimental FTIR spectra of hydroxybenzoic acids; Table S6: Chemical shifts δ (ppm) in the 1H and 13C NMR spectra of hydroxybenzoic acids; Table S7: Pearson’s correlation coefficient matrix between lipophilic and electronic parameters, aromaticity indices, and FTIR and NMR parameters of the selected dihydroxybenzoates; Table S8: PCA factor loadings, eigenvalues, and explained percentage of variance.

Author Contributions

Conceptualization, M.K. and W.L.; methodology, M.K., S.M.-W., A.P., A.C., and A.A.; software, G.Ś., A.A. and E.G.; formal analysis, M.K., G.Ś., E.G., S.M.-W., H.L., A.A. and R.Ś.; investigation, A.B.Z., A.P., A.C., M.K., A.A. and E.G.; resources, M.K. and W.L.; data curation, M.K., S.M.-W. and A.P.; writing—original draft preparation, E.G., M.K., G.Ś. and R.Ś.; writing—review and editing, M.K., H.L., M.S., W.P. and W.L.; visualization, E.G., A.A., G.Ś. and S.M.-W.; supervision, M.K.; project administration, M.K.; funding acquisition, W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Centre, Poland, grant no. UMO-2018/31/B/NZ7/03083.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Martínez-Cayuela, M. Oxygen free radicals and human disease. Biochimie 1995, 77, 147–161. [Google Scholar] [CrossRef]
  2. Phaniendra, A.; Jestadi, D.B.; Periyasamy, L. Free Radicals: Properties, Sources, Targets, and Their Implication in Various Diseases. Indian J. Clin. Biochem. 2015, 30, 11–26. [Google Scholar] [CrossRef] [Green Version]
  3. Di Meo, S.; Venditti, P. Evolution of the Knowledge of Free Radicals and Other Oxidants. Oxid. Med. Cell. Longev. 2020, 2020, 9829176. [Google Scholar] [CrossRef]
  4. Jurgowiak, M.; Oliński, R. Wolne rodniki a starzenie się. Kosmos 1995, 44, 71–88. [Google Scholar]
  5. Pizzino, G.; Irrera, N.; Cucinotta, M.; Pallio, G.; Mannino, F.; Arcoraci, V.; Squadrito, F.; Altavilla, D.; Bitto, A. Oxidative Stress: Harms and Benefits for Human Health. Oxid. Med. Cell. Longev. 2017, 2017, 8416763. [Google Scholar] [CrossRef]
  6. Niemann, B.; Rohrbach, S.; Miller, M.R.; Newby, D.E.; Fuster, V.; Kovacic, J.C. Oxidative Stress and Cardiovascular Risk: Obesity, Diabetes, Smoking, and Pollution: Part 3 of a 3-Part Series. J. Am. Coll. Cardiol. 2017, 70, 230–251. [Google Scholar] [CrossRef] [PubMed]
  7. Pham-Huy, L.A.; He, H.; Pham-Huy, C. Free radicals, antioxidants in disease and health. Int. J. Biomed. Sci. 2008, 4, 89–96. [Google Scholar]
  8. Umeno, A.; Biju, V.; Yoshida, Y. In vivo ROS production and use of oxidative stress-derived biomarkers to detect the onset of diseases such as Alzheimer’s disease, Parkinson’s disease, and diabetes. Free Radic. Res. 2017, 51, 413–427. [Google Scholar] [CrossRef] [PubMed]
  9. Valko, M.; Leibfritz, D.; Moncol, J.; Cronin, M.T.D.; Mazur, M.; Telser, J. Free radicals and antioxidants in normal physiological functions and human disease. Int. J. Biochem. Cell Biol. 2007, 39, 44–84. [Google Scholar] [CrossRef]
  10. Kang, D.H. Oxidative stress, DNA damage, and breast cancer. AACN Adv. Crit. Care 2002, 13, 540–549. [Google Scholar] [CrossRef]
  11. Sugamura, K.; Keaney, J.F. Reactive oxygen species in cardiovascular disease. Free Radic. Biol. Med. 2011, 51, 978–992. [Google Scholar] [CrossRef] [Green Version]
  12. Kurutas, E.B. The importance of antioxidants which play the role in cellular response against oxidative/nitrosative stress: Current state. Nutr. J. 2016, 15, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Septembre-Malaterre, A.; Remize, F.; Poucheret, P. Fruits and vegetables, as a source of nutritional compounds and phytochemicals: Changes in bioactive compounds during lactic fermentation. Food Res. Int. 2018, 104, 86–99. [Google Scholar] [CrossRef]
  14. Guardado, E.; Molina, E.; Joo, M.; Uriarte, E. Antioxidant and Pro-Oxidant Effects of Polyphenolic Compounds and Structure-Activity Relationship Evidence. Nutr. Well-Being Health 2012, 2, 23–48. [Google Scholar]
  15. Kowalczewski, P.Ł.; Radzikowska, D.; Ivanišová, E.; Szwengiel, A.; Kačániová, M.; Sawinska, Z. Influence of abiotic stress factors on the antioxidant properties and polyphenols profile composition of Green Barley (Hordeum vulgare L.). Int. J. Mol. Sci. 2020, 21, 397. [Google Scholar] [CrossRef] [Green Version]
  16. Ruskovska, T.; Budić-Leto, I.; Corral-Jara, K.F.; Ajdžanović, V.; Arola-Arnal, A.; Bravo, F.I.; Deligiannidou, G.-E.; Havlik, J.; Janeva, M.; Kistanova, E.; et al. Systematic Bioinformatic Analyses of Nutrigenomic Modifications by Polyphenols Associated with Cardiometabolic Health in Humans—Evidence from Targeted Nutrigenomic Studies. Nutrients 2021, 13, 2326. [Google Scholar] [CrossRef] [PubMed]
  17. Taamalli, A.; Contreras, M.D.M.; Abu-Reidah, I.M.; Trabelsi, N.; Ben Youssef, N. Quality of phenolic compounds: Occurrence, health benefits, and applications in food industry. J. Food Qual. 2019, 2019, 9594646. [Google Scholar] [CrossRef]
  18. Dos Santos, W.N.L.; da Silva Sauthier, M.C.; dos Santos, A.M.P.; de Andrade Santana, D.; Almeida Azevedo, R.S.; da Cruz Caldas, J. Simultaneous determination of 13 phenolic bioactive compounds in guava (Psidium guajava L.) by HPLC-PAD with evaluation using PCA and Neural Network Analysis (NNA). Microchem. J. 2017, 133, 583–592. [Google Scholar] [CrossRef]
  19. Yilmaz, Y.; Toledo, R.T. Major Flavonoids in Grape Seeds and Skins: Antioxidant Capacity of Catechin, Epicatechin, and Gallic Acid. J. Agric. Food Chem. 2004, 52, 255–260. [Google Scholar] [CrossRef]
  20. Li, R.; Chen, X.G.; Jia, K.; Liu, Z.P.; Peng, H.Y. A systematic determination of polyphenols constituents and cytotoxic ability in fruit parts of pomegranates derived from five Chinese cultivars. SpringerPlus 2016, 5, 914. [Google Scholar] [CrossRef] [Green Version]
  21. Di Stefano, V.; Avellone, G.; Bongiorno, D.; Indelicato, S.; Massenti, R.; Lo Bianco, R. Quantitative evaluation of the phenolic profile in fruits of six avocado (Persea americana) cultivars by ultra-high-performance liquid chromatography-heated electrospray-mass spectrometry. Int. J. Food Prop. 2017, 20, 1302–1312. [Google Scholar] [CrossRef] [Green Version]
  22. Nour, V.; Trandafir, I.; Cosmulescu, S. Antioxidant capacity, phenolic compounds and minerals content of blackcurrant (ribes nigrum L.) leaves as influenced by harvesting date and extraction method. Ind. Crops Prod. 2014, 53, 133–139. [Google Scholar] [CrossRef]
  23. Palafox-Carlos, H.; Yahia, E.M.; González-Aguilar, G.A. Identification and quantification of major phenolic compounds from mango (Mangifera indica, cv. Ataulfo) fruit by HPLC-DAD-MS/MS-ESI and their individual contribution to the antioxidant activity during ripening. Food Chem. 2012, 135, 105–111. [Google Scholar] [CrossRef]
  24. Hallmann, E.; Ponder, A.; Aninowski, M.; Narangerel, T.; Leszczyńska, J. The interaction between antioxidants content and allergenic potency of different raspberry cultivars. Antioxidants 2020, 9, 256. [Google Scholar] [CrossRef] [Green Version]
  25. Butkhup, L.; Samappito, W.; Samappito, S. Phenolic composition and antioxidant activity of white mulberry (Morus alba L.) fruits. Int. J. Food Sci. Technol. 2013, 48, 934–940. [Google Scholar] [CrossRef]
  26. Khadem, S.; Marles, R.J. Monocyclic phenolic acids; hydroxy- and polyhydroxybenzoic acids: Occurrence and recent bioactivity studies. Molecules 2010, 15, 7985–8005. [Google Scholar] [CrossRef] [PubMed]
  27. Dludla, P.V.; Nkambule, B.B.; Jack, B.; Mkandla, Z.; Mutize, T.; Silvestri, S.; Orlando, P.; Tiano, L.; Louw, J.; Mazibuko-Mbeje, S.E. Inflammation and oxidative stress in an obese state and the protective effects of gallic acid. Nutrients 2019, 11, 23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Yin, M.C.; Lin, C.C.; Wu, H.C.; Tsao, S.M.; Hsu, C.K. Apoptotic effects of protocatechuic acid in human breast, lung, liver, cervix, and prostate cancer cells: Potential mechanisms of action. J. Agric. Food Chem. 2009, 57, 6468–6473. [Google Scholar] [CrossRef] [PubMed]
  29. Kalinowska, M.; Gołębiewska, E.; Mazur, L.; Lewandowska, H.; Pruszyński, M.; Świderski, G.; Wyrwas, M.; Pawluczuk, N.; Lewandowski, W. Crystal Structure, Spectroscopic Characterization, Antioxidant and Cytotoxic Activity of New Mg(II) and Mn(II)/Na(I) Complexes of Isoferulic Acid. Materials 2021, 14, 3236. [Google Scholar] [CrossRef]
  30. Ashidate, K.; Kawamura, M.; Mimura, D.; Tohda, H.; Miyazaki, S.; Teramoto, T.; Yamamoto, Y.; Hirata, Y. Gentisic acid, an aspirin metabolite, inhibits oxidation of low-density lipoprotein and the formation of cholesterol ester hydroperoxides in human plasma. Eur. J. Pharmacol. 2005, 513, 173–179. [Google Scholar] [CrossRef]
  31. Abedi, F.; Razavi, B.M.; Hosseinzadeh, H. A review on gentisic acid as a plant derived phenolic acid and metabolite of aspirin: Comprehensive pharmacology, toxicology, and some pharmaceutical aspects. Phyther. Res. 2020, 34, 729–741. [Google Scholar] [CrossRef] [PubMed]
  32. Shi, G.F.; An, L.J.; Jiang, B.; Guan, S.; Bao, Y.M. Alpinia protocatechuic acid protects against oxidative damage in vitro and reduces oxidative stress in vivo. Neurosci. Lett. 2006, 403, 206–210. [Google Scholar] [CrossRef] [PubMed]
  33. Valko, M.; Izakovic, M.; Mazur, M.; Rhodes, C.J.; Telser, J. Role of oxygen radicals in DNA damage and cancer incidence. Mol. Cell. Biochem. 2004, 266, 37–56. [Google Scholar] [CrossRef]
  34. Lewandowska, H.; Kalinowska, M. New polyphenol-containing LDL nano-preparations in oxidative stress and DNA damage: A potential route for cell-targeted PP delivery. Materials 2020, 13, 5106. [Google Scholar] [CrossRef] [PubMed]
  35. Lin, H.-H.; Chen, J.-H.; Huang, C.-C.; Wang, C.-J. Apoptotic effect of 3,4-dihydroxybenzoic acid on human gastric carcinoma cells involving JNK/p38 MAPK signaling activation. Int. J. Cancer 2007, 120, 2306–2316. [Google Scholar] [CrossRef]
  36. Lodovici, M.; Guglielmi, F.; Meoni, M.; Dolara, P. Effect of natural phenolic acids on DNA oxidation in vitro. Food Chem. Toxicol. 2001, 39, 1205–1210. [Google Scholar] [CrossRef]
  37. Sappey, C.; Boelaert, J.R.; Legrandpoels, S.; Grady, R.W.; Piette, J. NF-κB transcription factor activation by hydrogen peroxide can be decreased by 2, 3-Dihydroxybenzoic acid and its ethyl ester derivative. Arch. Biochem. Biophys. 1995, 321, 263–270. [Google Scholar] [CrossRef] [PubMed]
  38. Sharma, S.; Khan, N.; Sultana, S. Study on prevention of two-stage skin carcinogenesis by Hibiscus rosa sinensis extract and the role of its chemical constituent, gentisic acid, in the inhibition of tumour promotion response and oxidative stress in mice. Eur. J. Cancer Prev. 2004, 13, 53–63. [Google Scholar] [CrossRef]
  39. Rafat Husain, S.; Cillard, J.; Cillard, P. Hydroxyl radical scavenging activity of flavonoids. Phytochemistry 1987, 26, 2489–2491. [Google Scholar] [CrossRef]
  40. Brand-Williams, W.; Cuvelier, M.E.; Berset, C. Use of a free radical method to evaluate antioxidant activity. LWT Food Sci. Technol. 1995, 28, 25–30. [Google Scholar] [CrossRef]
  41. Velika, B.; Kron, I. Antioxidant properties of benzoic acid derivatives against Superoxide radical. Free Rad. Antiox. 2012, 2, 62–67. [Google Scholar] [CrossRef] [Green Version]
  42. Rice-Evans, C.A.; Diplock, A.T.; Symons, M.C.R. Techniques in Free Radical Research; Elsevier: New York, NY, USA, 1991; pp. 1–291. [Google Scholar]
  43. Re, R.; Pellegrini, N.; Proteggente, A.; Pannala, A.; Yang, M.; Rice-Evans, C. Antioxidant activity applying an improved ABTS radical cation decolorization assay. Free Radic. Biol. Med. 1999, 26, 1231–1237. [Google Scholar] [CrossRef]
  44. Apak, R.; Güçlü, K.; Özyürek, M.; Karademir, S.E. Novel total antioxidant capacity index for dietary polyphenols and vitamins C and E, using their cupric ion reducing capability in the presence of neocuproine: CUPRAC method. J. Agric. Food Chem. 2004, 52, 7970–7981. [Google Scholar] [CrossRef]
  45. Zeraik, M.L.; Petrônio, M.S.; Coelho, D.; Regasini, L.O.; Silva, D.H.S.; da Fonseca, L.M.; Machado, S.A.S.; Bolzani, V.S.; Ximenes, V.F. Improvement of Pro-Oxidant Capacity of Protocatechuic Acid by Esterification. PLoS ONE 2014, 9, e110277. [Google Scholar] [CrossRef] [Green Version]
  46. Kalinowska, M.; Bajko, E.; Matejczyk, M.; Kaczyński, P.; Łozowicka, B.; Lewandowski, W. The study of anti-/pro-oxidant, lipophilic, microbial and spectroscopic properties of new alkali metal salts of 5-o-caffeoylquinic acid. Int. J. Mol. Sci. 2018, 19, 463. [Google Scholar] [CrossRef] [Green Version]
  47. PhysChem, ADME & Toxicity Calculations with Percepta Software. Available online: https://www.acdlabs.com/products/percepta/ (accessed on 7 July 2021).
  48. Julg, A.; Francoics, P. Structural studies of some non-alternant hydrocarbons; a new definition of aromaticity. Theor. Chim. Acta 1967, 7, 249–261. [Google Scholar] [CrossRef]
  49. Krygowski, T.M.; Ciesielski, A.; Bird, C.W.; Kotschy, A. Aromatic Character of the Benzene Ring Present in Various Topological Environments in Benzenoid Hydrocarbons. Nonequivalence of Indices of Aromaticity. J. Chem. Inf. Comput. Sci. 1995, 35, 205. [Google Scholar] [CrossRef]
  50. Cyrañski, M.K.; Krygowski, T.M.; Katritzky, A.R.; Schleyer, P.V.R. To what extent can aromaticity be defined uniquely? J. Org. Chem. 2002, 67, 1333–1338. [Google Scholar] [CrossRef] [PubMed]
  51. Bird, C.W. A new aromaticity index and its application to five-membered ring heterocycles. Tetrahedron 1985, 41, 1409–1414. [Google Scholar] [CrossRef]
  52. Schleyer, P.V.R.; Jiao, H. What is aromaticity? Pure Appl. Chem. 1996, 68, 209–218. [Google Scholar] [CrossRef]
  53. Weinhold, F.; Landis, C.R. Natural bond orbitals and extensions of localized bonding concepts. Chem. Educ. Res. Pract. 2001, 2, 91–104. [Google Scholar] [CrossRef]
  54. Breneman, C.M.; Wiberg, K.B. Determining atom-centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis. J. Comput. Chem. 1990, 11, 361–373. [Google Scholar] [CrossRef]
  55. Frish, M.J.; Trucks, G.W.; Schlegel, H.; Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; et al. Gaussian 09, Revision A.02; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  56. Sharma, O.P.; Bhat, T.K. DPPH antioxidant assay revisited. Food Chem. 2009, 113, 1202–1205. [Google Scholar] [CrossRef]
  57. Francenia Santos-Sánchez, N.; Salas-Coronado, R.; Villanueva-Cañongo, C.; Hernández-Carlos, B. Antioxidant Compounds and Their Antioxidant Mechanism. In Antioxidants; IntechOpen: London, UK, 2019; pp. 1–28. [Google Scholar]
  58. Holtz, R.W. In Vitro Methods to Screen Materials for Anti-aging Effects. In Skin Aging Handbook; Elsevier Inc.: Amsterdam, The Netherlands, 2009; pp. 329–362. [Google Scholar]
  59. Wojdyło, A.; Oszmiański, J.; Czemerys, R. Antioxidant activity and phenolic compounds in 32 selected herbs. Food Chem. 2007, 105, 940–949. [Google Scholar] [CrossRef]
  60. Cieślik, E.; Grȩda, A.; Adamus, W. Contents of polyphenols in fruit and vegetables. Food Chem. 2006, 94, 135–142. [Google Scholar] [CrossRef]
  61. Mishra, K.; Ojha, H.; Chaudhury, N.K. Estimation of antiradical properties of antioxidants using DPPH- assay: A critical review and results. Food Chem. 2012, 130, 1036–1043. [Google Scholar] [CrossRef]
  62. Kalinowska, M.; Mazur, L.; Jabłońska-Trypuć, A.; Lewandowski, W. A new calcium 2,5-dihydroxybenzoate: Synthesis, characterization and antioxidant studies and stress mediated cytotoxity in MCF-7 cells. J. Saudi Chem. Soc. 2018, 22, 742–756. [Google Scholar] [CrossRef]
  63. Prihantini, A.I.; Tachibana, S.; Itoh, K. Antioxidant Active Compounds from Elaeocarpussylvestris and their Relationship between Structure and Activity. Procedia Environ. Sci. 2015, 28, 758–768. [Google Scholar] [CrossRef] [Green Version]
  64. Santoso, S.P.; Ismadji, S.; Angkawijaya, A.E.; Soetaredjo, F.E.; Go, A.W.; Ju, Y.H. Complexes of 2,6-dihydroxybenzoic acid with divalent metal ions: Synthesis, crystal structure, spectral studies, and biological activity enhancement. J. Mol. Liq. 2016, 221, 617–623. [Google Scholar] [CrossRef] [Green Version]
  65. Mardani-Ghahfarokhi, A.; Farhoosh, R. Antioxidant activity and mechanism of inhibitory action of gentisic and α-resorcylic acids. Sci. Rep. 2020, 10, 19487. [Google Scholar] [CrossRef]
  66. Nofrizal; Putra, D.P.; Arbain, D. Antioxidant and antibacterial constituents from two Sumatran ferns, Trichomanes javanicum and Oleandra pistillaris. Nat. Prod. Commun. 2017, 12, 1263–1264. [Google Scholar]
  67. Karamać, M.; Kosińska, A. Comparison of radical-scavenging activities for selected phenolic acids. Pol. J. Food Nutr. Sci. 2005, 55, 165–170. [Google Scholar]
  68. Czerniak, K.; Walkiewicz, F. Synthesis and antioxidant properties of dicationic ionic liquids. New J. Chem. 2017, 41, 530–539. [Google Scholar] [CrossRef]
  69. Sroka, Z.; Cisowski, W. Hydrogen peroxide scavenging, antioxidant and anti-radical activity of some phenolic acids. Food Chem. Toxicol. 2003, 41, 753–758. [Google Scholar] [CrossRef]
  70. Tung, Y.T.; Wu, J.H.; Kuo, Y.H.; Chang, S.T. Antioxidant activities of natural phenolic compounds from Acacia confusa bark. Bioresour. Technol. 2007, 98, 1120–1123. [Google Scholar] [CrossRef]
  71. Mathew, S.; Abraham, T.E.; Zakaria, Z.A. Reactivity of phenolic compounds towards free radicals under in vitro conditions. J. Food Sci. Technol. 2015, 52, 5790–5798. [Google Scholar] [CrossRef] [Green Version]
  72. Yen, G.C.; Der Duh, P.; Tsai, H.L. Antioxidant and pro-oxidant properties of ascorbic acid and gallic acid. Food Chem. 2002, 79, 307–313. [Google Scholar] [CrossRef]
  73. Biskup, I.; Golonka, I.; Gamian, A.; Sroka, Z. Antioxidant activity of selected phenols estimated by ABTS and FRAP methods. Adv. Hyg. Exp. Med. 2013, 67, 958–963. [Google Scholar] [CrossRef]
  74. Li, X.; Wang, X.; Chen, D.; Chen, S. Antioxidant activity and mechanism of protocatechuic acid in Vitro. Funct. Foods Health Dis. 2011, 1, 232–244. [Google Scholar] [CrossRef]
  75. Kukoc-Modun, L.; Tsikas, D.; Biocic, M.; Radić, N. Flow injection analysis of N-acetyl-l-cysteine based on the reduction of copper(ii)-neocuproine reagent. Anal. Lett. 2016, 49, 607–617. [Google Scholar] [CrossRef]
  76. Spiegel, M.; Kapusta, K.; Kołodziejczyk, W.; Saloni, J.; Żbikowska, B.; Hill, G.A.; Sroka, Z. Antioxidant Activity of Selected Phenolic Acids–Ferric Reducing Antioxidant Power Assay and QSAR Analysis of the Structural Features. Molecules 2020, 25, 3088. [Google Scholar] [CrossRef]
  77. Rice-Evans, C.A.; Miller, N.J.; Paganga, G. Structure-antioxidant activity relationships of flavonoids and phenolic acids. Free Radic. Biol. Med. 1996, 20, 933–956. [Google Scholar] [CrossRef]
  78. Owusu Apenten, R. A FRAP Assay at pH 7 unveils Extra Antioxidant Activity from Green, Black, White and Rooibos Tea but not Apple Tea. Food Nutr. Rep. 2015, 1, 16–23. [Google Scholar] [CrossRef]
  79. Litwinienko, G.; Ingold, K.U. Abnormal solvent effects on hydrogen atom abstraction. 2. Resolution of the curcumin antioxidant controversy. The role of sequential proton loss electron transfer. J. Org. Chem. 2004, 69, 5888–5896. [Google Scholar] [CrossRef]
  80. Foti, M.C.; Daquino, C.; Geraci, C. Electron-transfer reaction of cinnamic acids and their methyl esters with the DPPH(*) radical in alcoholic solutions. J. Org. Chem. 2004, 69, 2309–2314. [Google Scholar] [CrossRef] [PubMed]
  81. Borges, A.; Ferreira, C.; Saavedra, M.J.; Simões, M. Antibacterial activity and mode of action of ferulic and gallic acids against pathogenic bacteria. Microb. Drug Resist. 2013, 19, 256–265. [Google Scholar] [CrossRef] [PubMed]
  82. Sánchez-Maldonado, A.F.; Schieber, A.; Gänzle, M.G. Structure-function relationships of the antibacterial activity of phenolic acids and their metabolism by lactic acid bacteria. J. Appl. Microbiol. 2011, 111, 1176–1184. [Google Scholar] [CrossRef] [PubMed]
  83. Stojković, D.S.; Živković, J.; Soković, M.; Glamočlija, J.; Ferreira, I.C.F.R.; Janković, T.; Maksimović, Z. Antibacterial activity of Veronica montana L. extract and of protocatechuic acid incorporated in a food system. Food Chem. Toxicol. 2013, 55, 209–213. [Google Scholar] [CrossRef]
  84. Mattson, T.E.; Johny, A.K.; Amalaradjou, M.A.R.; More, K.; Schreiber, D.T.; Patel, J.; Venkitanarayanan, K. Inactivation of Salmonella spp. on tomatoes by plant molecules. Int. J. Food Microbiol. 2011, 144, 464–468. [Google Scholar] [CrossRef]
  85. Baskaran, S.A.; Upadhyay, A.; Kollanoor-Johny, A.; Upadhyaya, I.; Mooyottu, S.; Roshni Amalaradjou, M.A.; Schreiber, D.; Venkitanarayanan, K. Efficacy of plant-derived antimicrobials as antimicrobial wash treatments for reducing enterohemorrhagic escherichia coli O157: H7 on apples. J. Food Sci. 2013, 78, M1399–M1404. [Google Scholar] [CrossRef] [PubMed]
  86. Upadhyay, A.; Upadhyaya, I.; Kollanoor-Johny, A.; Ananda Baskaran, S.; Mooyottu, S.; Karumathil, D.; Venkitanarayanan, K. Inactivation of Listeria monocytogenes on frankfurters by plant-derived antimicrobials alone or in combination with hydrogen peroxide. Int. J. Food Microbiol. 2013, 163, 114–118. [Google Scholar] [CrossRef] [PubMed]
  87. Ahire, J.J.; Dicks, L.M.T. 2,3-dihydroxybenzoic acid-containing nanofiber wound dressings inhibit biofilm formation by Pseudomonas aeruginosa. Antimicrob. Agents Chemother. 2014, 58, 2098–2104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Pj, B.; Shibumon, G.; Sunny, K.; Cincy, G. 2,3-Dihydroxybenzoic Acid: An Effective Antifungal Agent Isolated from Flacourtia inermis Fruit. Int. J. Pharm. Clin. Res. 2010, 2, 101–105. [Google Scholar]
  89. Lattanzio, V.; Di Venere, D.; Linsalata, V.; Lima, G.; Ippolito, A.; Salerno, M. Antifungal activity of 2,5-dimethoxybenzoic acid on postharvest pathogens of strawberry fruits. Postharvest Biol. Technol. 1996, 9, 325–334. [Google Scholar] [CrossRef]
  90. Merkl, R.; Hrádková, I.; Filip, V.; Šmidrkal, J. Antimicrobial and antioxidant properties of phenolic acids alkyl esters. Czech J. Food Sci. 2010, 28, 275–279. [Google Scholar] [CrossRef] [Green Version]
  91. Campos, F.M.; Couto, J.A.; Figueiredo, A.R.; Tóth, I.V.; Rangel, A.O.S.S.; Hogg, T.A. Cell membrane damage induced by phenolic acids on wine lactic acid bacteria. Int. J. Food Microbiol. 2009, 135, 144–151. [Google Scholar] [CrossRef]
  92. Hsu, J.D.; Kao, S.H.; Ou, T.T.; Chen, Y.J.; Li, Y.J.; Wang, C.J. Gallic acid induces G2/M phase arrest of breast cancer cell MCF-7 through stabilization of p27Kip1 attributed to disruption of p27 Kip1/Skp2 complex. J. Agric. Food Chem. 2011, 59, 1996–2003. [Google Scholar] [CrossRef] [PubMed]
  93. Marcelia, M.; Arsianti, A.; Chan, J.O.T.; Stevano, J.W.; Fadilah, F.; Putrianingsih, R.; Azizah, N.N.; Bahtiar, A.; Tanimoto, H. Kiyomikakiuchi In vitro cytotoxicity of the synthesized gallic acid derivatives (N-Alkyl Gallamide) against breast MCF-7 cancer cells. Orient. J. Chem. 2018, 34, 2268–2272. [Google Scholar] [CrossRef]
  94. Maurya, D.K.; Nandakumar, N.; Devasagayam, T.P.A. Anticancer property of gallic acid in A549, a human lung adenocarcinoma cell line, and possible mechanisms. J. Clin. Biochem. Nutr. 2011, 48, 85–90. [Google Scholar] [CrossRef] [Green Version]
  95. Moghtaderi, H.; Sepehri, H.; Delphi, L.; Attari, F. Gallic acid and curcumin induce cytotoxicity and apoptosis in human breast cancer cell MDA-MB-231. BioImpacts 2018, 8, 185–194. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Sherin, L.; Sohail, A.; Shujaat, S. Time-dependent AI-Modeling of the anticancer efficacy of synthesized gallic acid analogues. Comput. Biol. Chem. 2019, 79, 137–146. [Google Scholar] [CrossRef]
  97. Roleira, F.M.F.; Tavares-Da-Silva, E.J.; Varela, C.L.; Costa, S.C.; Silva, T.; Garrido, J.; Borges, F. Plant derived and dietary phenolic antioxidants: Anticancer properties. Food Chem. 2015, 183, 235–258. [Google Scholar] [CrossRef] [PubMed]
  98. Sankaranarayanan, R.; Valiveti, C.K.; Dachineni, R.; Kumar, D.R.; Lick, T.; Bhat, G.J. Aspirin metabolites 2,3-DHBA and 2,5-DHBA inhibit cancer cell growth: Implications in colorectal cancer prevention. Mol. Med. Rep. 2020, 21, 20–34. [Google Scholar] [CrossRef] [Green Version]
  99. Yip, E.C.H.; Chan, A.S.L.; Pang, H.; Tam, Y.K.; Wong, Y.H. Protocatechuic acid induces cell death in HepG2 hepatocellular carcinoma cells through a c-Jun N-terminal kinase-dependent mechanism. Cell Biol. Toxicol. 2006, 22, 293–302. [Google Scholar] [CrossRef] [Green Version]
  100. Reis, F.S.; Sousa, D.; Barros, L.; Martins, A.; Morales, P.; Ferreira, I.C.F.R.; Vasconcelos, M.H. Leccinum vulpinum Watling induces DNA damage, decreases cell proliferation and induces apoptosis on the human MCF-7 breast cancer cell line. Food Chem. Toxicol. 2016, 90, 45–54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Chmiel, T.; Mieszkowska, A.; Kempińska-Kupczyk, D.; Kot-Wasik, A.; Namieśnik, J.; Mazerska, Z. The impact of lipophilicity on environmental processes, drug delivery and bioavailability of food components. Microchem. J. 2019, 146, 393–406. [Google Scholar] [CrossRef]
  102. Hansch, C.; Leo, A.; Hoekman, D.; Livingstone, D.J. Exploring QSAR: Hydrophobic, Electronic, and Steric Constants; American Chemical Society: Washington, DC, USA, 1995; Volume 1, pp. 155–156. [Google Scholar]
  103. Zevatskii, Y.E.; Lysova, S.S. Empirical Procedure for the Calculation of Ionization Constants of Organic Compounds in Water from Their Molecular Volume. Russ. J. Org. Chem. 2009, 45, 825–834. [Google Scholar] [CrossRef]
  104. Chemical Book. Available online: https://www.chemicalbook.com/ProductIndex_EN.aspx (accessed on 8 July 2021).
  105. Human Metabolome Database. Available online: https://hmdb.ca/ (accessed on 8 July 2021).
  106. Kortum, G.; Vogel, W.; Andrussow, K. Dissociation Constants of Organic Acids in Aqueous Solution; Butterworths: London, UK, 1961. [Google Scholar]
  107. Sangster, J.M. Octanol-Water Partition Coefficients: Fundamentals and Physical Chemistry; Wiley: New York, NY, USA, 1997; Volume 2. [Google Scholar]
  108. Taylor, T.; Seitz, E.; Fox, P.; Fuls, J. Antimicrobial Compositions Containing an Aromatic Carboxylic Acid and a Hydric Solvent. U.S. Patent 7,592,300, 22 September 2009. [Google Scholar]
  109. Serjeant, E.P.; Bempsey, B. Ionisation constants of organic acids in aqueous solution. IUPAC Chem. Data Ser. 1979, 23, 169–190. [Google Scholar]
  110. Patty, F.A. Industrial Hygiene and Toxicology; John Wiley & Sons Inc.: New York, NY, USA, 1963. [Google Scholar]
  111. Honório, K.M.; Da Silva, A.B.F. An AM1 study on the electron-donating and electron-accepting character of biomolecules. Int. J. Quantum Chem. 2003, 95, 126–132. [Google Scholar] [CrossRef]
  112. Marković, Z.; Crossed, D.; Signorović, J.; Dimitrić Marković, J.M.; Živić, M.; Amić, D. Investigation of the radical scavenging potency of hydroxybenzoic acids and their carboxylate anions. Mon. Chem. Chem. Mon. 2014, 145, 953–962. [Google Scholar] [CrossRef]
  113. Chen, Y.; Xiao, H.; Zheng, J.; Liang, G. Structure-thermodynamics-antioxidant activity relationships of selected natural phenolic acids and derivatives: An experimental and theoretical evaluation. PLoS ONE 2015, 10, e0121276. [Google Scholar] [CrossRef] [PubMed]
  114. National Institute of Advanced Industrial Science and Technology. Available online: https://sdbs.db.aist.go.jp/ (accessed on 10 July 2021).
  115. Varsanyi, G.; Kovner, M.A.; Lang, L. Assignments for Vibrational Spectra of 700 Benzene Derivatives; Akademiai Kiado: Budapest, Hungary, 1973. [Google Scholar]
  116. Lewandowski, W.; Kalinowska, M.; Lewandowska, H. The influence of metals on the electronic system of biologically important ligands. Spectroscopic study of benzoates, salicylates, nicotinates and isoorotates. Review. J. Inorg. Biochem. 2005, 99, 1407–1423. [Google Scholar] [CrossRef]
  117. Kalinowska, M.; Siemieniuk, E.; Kostro, A.; Lewandowski, W. The application of Aj, BAC, I6, HOMA indexes for quantitative determination of aromaticity of metal complexes with benzoic, salicylic, nicotinic acids and benzene derivatives. J. Mol. Struct. THEOCHEM 2006, 761, 129–141. [Google Scholar] [CrossRef]
Figure 1. Some plants containing gallic acid (on the basis of data from [18,19,20,21,22,23,24,25,26,27,28]).
Figure 1. Some plants containing gallic acid (on the basis of data from [18,19,20,21,22,23,24,25,26,27,28]).
Nutrients 13 03107 g001
Figure 2. Scavenging of DPPH radicals by 2,3-DHB, 2,5-DHB, 2,6-DHB, 3,4-DHB, and 3,4,5-THB. Mean values from three independent experiments ± SD are shown.
Figure 2. Scavenging of DPPH radicals by 2,3-DHB, 2,5-DHB, 2,6-DHB, 3,4-DHB, and 3,4,5-THB. Mean values from three independent experiments ± SD are shown.
Nutrients 13 03107 g002aNutrients 13 03107 g002b
Figure 3. Comparison of the percentage of inhibition of ABTS•+ cation radicals measured after 7 min of reaction time for the tested hydroxybenzoic acids at a concentration of 50 μM (methanolic solutions). Mean values from three independent experiments ± SD are shown.
Figure 3. Comparison of the percentage of inhibition of ABTS•+ cation radicals measured after 7 min of reaction time for the tested hydroxybenzoic acids at a concentration of 50 μM (methanolic solutions). Mean values from three independent experiments ± SD are shown.
Nutrients 13 03107 g003
Figure 4. Comparison of the antioxidant activity of the tested hydroxybenzoic acids at a concentration of 50 μM (methanolic solutions), as measured by the CUPRAC assay. Mean values from three independent experiments ± SD are shown.
Figure 4. Comparison of the antioxidant activity of the tested hydroxybenzoic acids at a concentration of 50 μM (methanolic solutions), as measured by the CUPRAC assay. Mean values from three independent experiments ± SD are shown.
Nutrients 13 03107 g004
Figure 5. The degree of Trolox oxidation (%) by tested hydroxybenzoic acids at concentrations of 0.15 and 0.35 μM (methanolic solutions), measured within 90 or 120 min. Mean values from three independent experiments ± SD are shown.
Figure 5. The degree of Trolox oxidation (%) by tested hydroxybenzoic acids at concentrations of 0.15 and 0.35 μM (methanolic solutions), measured within 90 or 120 min. Mean values from three independent experiments ± SD are shown.
Nutrients 13 03107 g005
Figure 6. Survival curves of MDA-MB-231 and MCF-7 cells in the presence of different concentrations of the tested hydroxybenzoic acids.
Figure 6. Survival curves of MDA-MB-231 and MCF-7 cells in the presence of different concentrations of the tested hydroxybenzoic acids.
Nutrients 13 03107 g006aNutrients 13 03107 g006b
Figure 7. The dependency between the lipophilicity (logPexp.) and acidity (pKa exp.) in the series of studied hydroxybenzoic acids (yellow-marked: hydroxybenzoates with the highest antioxidant activity in the group; circled with a red square: compounds with the highest antimicrobial properties against selected microorganisms; circled with a blue square: compounds with the highest cytotoxic activity against selected cancer cell lines).
Figure 7. The dependency between the lipophilicity (logPexp.) and acidity (pKa exp.) in the series of studied hydroxybenzoic acids (yellow-marked: hydroxybenzoates with the highest antioxidant activity in the group; circled with a red square: compounds with the highest antimicrobial properties against selected microorganisms; circled with a blue square: compounds with the highest cytotoxic activity against selected cancer cell lines).
Nutrients 13 03107 g007
Figure 8. Electrostatic potential maps (calculated via the CHelpG method in B3LYP/6-311++G(d,p) for dihydroxybenzoic acids and 3,4,5-trihydroxybenzoic acid (in the gas phase and water solution).
Figure 8. Electrostatic potential maps (calculated via the CHelpG method in B3LYP/6-311++G(d,p) for dihydroxybenzoic acids and 3,4,5-trihydroxybenzoic acid (in the gas phase and water solution).
Nutrients 13 03107 g008
Figure 9. Principal component analysis (PCA) of lipophilic and electronic parameters, aromaticity indices, FTIR and NMR parameters, and the anti-/pro-oxidant (a,b) and antimicrobial potential (c,d) of the studied compounds.
Figure 9. Principal component analysis (PCA) of lipophilic and electronic parameters, aromaticity indices, FTIR and NMR parameters, and the anti-/pro-oxidant (a,b) and antimicrobial potential (c,d) of the studied compounds.
Nutrients 13 03107 g009aNutrients 13 03107 g009bNutrients 13 03107 g009c
Table 1. DPPH antioxidant activity of the tested hydroxybenzoic acids (methanolic solutions) measured after 60 min of reaction time. Mean values from three independent experiments ± SD are shown.
Table 1. DPPH antioxidant activity of the tested hydroxybenzoic acids (methanolic solutions) measured after 60 min of reaction time. Mean values from three independent experiments ± SD are shown.
Compound2,3-DHB2,4-DHB2,5-DHB2,6-DHB3,4-DHB3,5-DHB3,4,5-THB
DPPH/IC50 (µM)5.86 ± 0.64>120,0004.99 ± 0.79209.02 ± 0.985.24 ± 0.95>10002.42 ± 0.08
Table 2. Antioxidant activity of selected hydroxybenzoic acids on the basis of literature data.
Table 2. Antioxidant activity of selected hydroxybenzoic acids on the basis of literature data.
Method2,3-DHB2,4-DHB2,5-DHB2,6-DHB3,4-DHB3,5-DHB3,4,5-THB
DPPH,
IC50 (μM)
3.52 * (CDPPH = 1 mM,
t = 30 min) [63]
4.95 (CDPPH = 60 μM,
t = 60 min) [62]
231.51 (CDPPH = 500 μM,
t = 30 min) [64]
3.91 * (CDPPH = 1 mM,
t = 30 min) [63]
384.6 (CDPPH = 90 μM,
t = 30 min) [65]
3.78 * (CDPPH = 1 mM,
t = 30 min) [63]
5.18 * (CDPPH = 1 mM,
t = 30 min) [63]
15.38 * (CDPPH = 1 mM,
t = 30 min) [63]
3.83 (CDPPH = 50 μM,
t = 30 min) [66]
0.0237 (CDPPH = 1 mM,
t = 20 min) [67]
37.04 (CDPPH = 90 μM,
t = 30 min) [65]
0.0574 (CDPPH = 1 mM,
t = 20 min) [67]
0.0292 (CDPPH = 1 mM,
t = 20 min) [67]
15.83 (CDPPH = 60 μM,
t = 30 min) [68]
6.30 (CDPPH = 60 μM,
t = 30 min) [68]
DPPH radical scavenging
activity (%)
46 (C2,3-DHB = 0.11 mg/mL, CDPPH = 100 μM, t = 1 min) [69]0.11 (C2,4-DHB = 0.11 mg/mL, CDPPH = 100 μM, t = 1 min) [69]30.5 (C2,5-DHB = 0.11 mg/mL, CDPPH = 100 μM, t = 1 min) [69] 41.2 (C3,4-DHB = 0.11 mg/mL, CDPPH = 100 μM, t = 1 min) [69]0.6 (C3,5-DHB = 0.11 mg/mL, CDPPH = 100 μM, t = 1 min) [69]75 (C3,4,5-THB = 0.11 mg/mL, CDPPH = 100 μM, t = 1 min) [69]
≈ 70 (C3,4-DHB = 10 μM, CDPPH = 100 μM, t = 30 min) [70] ≈ 70 (C3,4,5-THB = 10 μM, CDPPH = 100 μM, t = 30 min) [70]
≈ 88.5 (C3,4-DHB = 15 μM, CDPPH = 60 μM, t = 30 min) [71] 40 (C3,4,5-THB≈ 0.7 mM, CDPPH = 200 μM, t = 30 min) [72]
88.5 (C3,4,5-THB = 15 μM, CDPPH = 60 μM, t = 30 min) [71]
ABTS,
IC50 (μM)
10.73 (CABTS = 7 mM,
t = 6 min) [73]
169.44 (CABTS = 7 mM,
t = 6 min) [73]
15.56 (CABTS = 7 mM,
t = 6 min) [73]
20.67 (CABTS = 7 mM,
t = 6 min) [73]
41.85 (CABTS = 7 mM,
t = 6 min) [73]
4.33 (CABTS = 7 mM,
t = 6 min) [73]
≈ 6.50 (CABTS = 7 mM,
t = 6 min) [68]
0.89 * (CABTS = n.d.,
t = 6 min) [74]
10.07 (CABTS = 7 mM,
t = 6 min) [68]
ABTS radical scavenging
activity (%)
≈ 38 (C3,4,5-THB = 10 μM, CABTS = 7 mM, t = 7 min) [71] ≈ 80 (C3,4,5-THB = 10 μM, CABTS = 7 mM, t = 7 min) [71]
FRAP (μM Fe2+) ≈ 0.012 (C = 0,001 mg/mL), λ = 593 nm, t = 30 min)
[73]
≈ 0.98 (C = 0,001 mg/mL),
λ = 593 nm, t = 30 min)
[73]
≈ 0.65 (C = 0,001 mg/mL),
λ = 593 nm, t = 30 min)
[73]
≈ 0.03 (C = 0,001 mg/mL),
λ = 593 nm, t = 30 min)
[73]
≈ 1.03 (C = 0,001 mg/mL),
λ = 593 nm, t = 30 min)
[73]
75.39 (C = 5 μM, λ = 594 nm, t = n.d.) [62]
236.86 (C = 50 μM, λ = 594 nm, t = n.d.) [62]
* Results expressed as μg/L; ≈ the exact value was not given; n.d.: no data.
Table 3. FRAP antioxidant activity of tested hydroxybenzoic acids at a concentration of 50 μM (methanolic solutions). Mean values from three independent experiments ± SD are shown.
Table 3. FRAP antioxidant activity of tested hydroxybenzoic acids at a concentration of 50 μM (methanolic solutions). Mean values from three independent experiments ± SD are shown.
Compound2,3-DHB2,4-DHB2,5-DHB2,6-DHB3,4-DHB3,5-DHB3,4,5-THB
FRAP (µM Fe2+)173.79 ± 39.06(-) *236.00 ± 6.31(-) *44.22 ± 9.38(-) *158.10 ± 14.81
* Negative values were obtained for these acids.
Table 4. The minimum inhibitory concentration (MIC) values (mg/mL) for the tested hydroxybenzoic acids against selected microorganisms.
Table 4. The minimum inhibitory concentration (MIC) values (mg/mL) for the tested hydroxybenzoic acids against selected microorganisms.
MicroorganismsCompound
2,3-DHB2,4-DHB2,5-DHB2,6-DHB3,4-DHB3,5-DHB3,4,5-THB
MIC (mg/mL)
Escherichia coli5233223
Pseudomonas aeruginosa5254234
Staphylococcus aureus5255226
Bacillus subtilis5233223
Salmonella enteritidis5233223
Candida albicans5253223
Table 5. IC50 values for the tested hydroxybenzoic acids on breast cancer cells. Mean values from three independent experiments ± SD are shown.
Table 5. IC50 values for the tested hydroxybenzoic acids on breast cancer cells. Mean values from three independent experiments ± SD are shown.
CompoundCytotoxicity (mM) in Cell Line
MDA-MB-231MCF-7
2,3-DHBIC50 = 4.10 ± 0.76NT *
2,4-DHBIC50 = 4.77 ± 1.20NT *
2,5-DHBIC50 = 4.39 ± 1.50NT *
2,6-DHBIC50 = 3.97 ± 2.50NT *
3,4-DHBNT *NT *
3,5-DHBNT *NT *
3,4,5-THBIC50 = 0.36 ± 0.50IC50 = 0.44 ± 1.80
* NT: nontoxic.
Table 6. Lipophilicity parameters determined by chromatographic methods (the logarithm of the retention factor, logkw) and octanol/water partition coefficients logP obtained for hydroxybenzoic acids.
Table 6. Lipophilicity parameters determined by chromatographic methods (the logarithm of the retention factor, logkw) and octanol/water partition coefficients logP obtained for hydroxybenzoic acids.
CompoundC18C8CNIAMPHE
log kwϕ0R2log kwϕ0R2log kwϕ0R2log kwϕ0R2log kwϕ0R2
2,3-DHB2.470.600.9501.700.870.9941.591.160.9972.171.080.9961.250.930.972
2,4-DHB2.660.570.9831.811.380.9812.061.150.9991.141.000.9951.391.050.996
2,5-DHB2.390.550.9851.370.900.9931.641.150.9971.601.110.9971.171.010.986
2,6-DHB2.990.680.9261.171.050.9271.491.080.9751.311.140.9890.871.050.998
3,4-DHB2.320.450.9791.000.830.9881.581.080.9960.441.130.9931.070.500.996
3,5-DHB2.560.450.9830.900.830.9781.781.070.9940.511.160.9931.030.550.985
3,4,5-THB1.910.370.9770.500.910.9791.241.060.9770.151.440.9890.830.330.961
Table 7. Theoretical (calculated in ACD/Labs or Percepta: logP classic, logP galas, pKa classic) and experimental logP and pKa values.
Table 7. Theoretical (calculated in ACD/Labs or Percepta: logP classic, logP galas, pKa classic) and experimental logP and pKa values.
CompoundLogP ClassicLogP GalasLogPexp.ACD/pKa ClassicpKaexp.
2,3-DHB1.80 ± 0.261.421.20 [102]pKa1 = 3.0 ± 0.1
pKa2 = 10.1 ± 0.1
pKa3 = 15.6 ± 0.1
2.91 (at 25 °C) [103]
2.98 [104]
2.56 [105]
2,4-DHB1.60 ± 0.261.821.63 [102]pKa1 = 3.3 ± 0.1
pKa2 = 9.1 ± 0.2
pKa3 = 13.8 ± 0.1
3.11 (at 25 °C) [106]
3.11 [104]
3.1 [105]
2,5-DHB1.56 ± 0.261.551.74 [107]pKa1 = 3.0 ± 0.1
pKa2 = 11.0 ± 0.2
pKa3 = 14.7 ± 0.2
2.97 (at 25 °C) [104]
2.84 [105]
2.97 [108]
2,6-DHB2.24 ± 0.272.322.20 [102]pKa1 = 1.3 ± 0.1
pKa2 = 12.5 ± 0.1
pKa3 = 14.1 ± 0.1
1.30 [104]
1.51 [105]
3,4-DHB1.16 ± 0.240.970.86 [102]pKa1 = 4.5 ± 0.1
pKa2 = 9.1 ± 0.2
pKa3 = 12.8 ± 0.1
4.26 (at 25 °C) [109]
4.48 [104]
4.16 [105]
3,5-DHB1.12 ± 0.240.990.86 [102]pKa1 = 4.0 ± 0.1
pKa2 = 9.8 ± 0.1
pKa3 = 11.3 ± 0.1
4.04 [104]
3.61 [105]
4.04 [108]
3,4,5-THB0.91 ± 0.330.610.70 [102]pKa1 = 4.3 ± 0.1
pKa2 = 9.0 ± 0.2
pKa3 = 12.3 ± 0.1
pKa4 = 14.4 ± 0.2
4.40 [110]
4.41 [104]
3.94 [105]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kalinowska, M.; Gołębiewska, E.; Świderski, G.; Męczyńska-Wielgosz, S.; Lewandowska, H.; Pietryczuk, A.; Cudowski, A.; Astel, A.; Świsłocka, R.; Samsonowicz, M.; et al. Plant-Derived and Dietary Hydroxybenzoic Acids—A Comprehensive Study of Structural, Anti-/Pro-Oxidant, Lipophilic, Antimicrobial, and Cytotoxic Activity in MDA-MB-231 and MCF-7 Cell Lines. Nutrients 2021, 13, 3107. https://doi.org/10.3390/nu13093107

AMA Style

Kalinowska M, Gołębiewska E, Świderski G, Męczyńska-Wielgosz S, Lewandowska H, Pietryczuk A, Cudowski A, Astel A, Świsłocka R, Samsonowicz M, et al. Plant-Derived and Dietary Hydroxybenzoic Acids—A Comprehensive Study of Structural, Anti-/Pro-Oxidant, Lipophilic, Antimicrobial, and Cytotoxic Activity in MDA-MB-231 and MCF-7 Cell Lines. Nutrients. 2021; 13(9):3107. https://doi.org/10.3390/nu13093107

Chicago/Turabian Style

Kalinowska, Monika, Ewelina Gołębiewska, Grzegorz Świderski, Sylwia Męczyńska-Wielgosz, Hanna Lewandowska, Anna Pietryczuk, Adam Cudowski, Aleksander Astel, Renata Świsłocka, Mariola Samsonowicz, and et al. 2021. "Plant-Derived and Dietary Hydroxybenzoic Acids—A Comprehensive Study of Structural, Anti-/Pro-Oxidant, Lipophilic, Antimicrobial, and Cytotoxic Activity in MDA-MB-231 and MCF-7 Cell Lines" Nutrients 13, no. 9: 3107. https://doi.org/10.3390/nu13093107

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop