Next Article in Journal
Antioxidant Properties of Gluten-Free Pasta Enriched with Vegetable By-Products
Next Article in Special Issue
Impact of Thermal Processing on the Selected Biological Activities of Ginger Rhizome—A Review
Previous Article in Journal
Studies on the Selective Syntheses of Sodium Ditelluride and Dialkyl Ditellurides
Previous Article in Special Issue
Stimulation of Lignan Production in Schisandra rubriflora In Vitro Cultures by Elicitation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multidirectional Characterization of Phytochemical Profile and Health-Promoting Effects of Ziziphora bungeana Juz. Extracts

by
Karlygash Zhaparkulova
1,*,
Aigerim Karaubayeva
1,
Zuriyadda Sakipova
1,
Anna Biernasiuk
2,
Katarzyna Gaweł-Bęben
3,*,
Tomasz Laskowski
4,
Aliya Kusniyeva
1,
Azamat Omargali
5,
Tolkyn Bekezhanova
1,
Liliya Ibragimova
1,
Galiya Ibadullayeva
1,
Amangeldy Jakiyanov
1,
Karolina Czech
3,
Kuanysh Tastambek
1,6,7,8,
Kazimierz Głowniak
3,
Anna Malm
2 and
Wirginia Kukula-Koch
9
1
School of Pharmacy, S.D. Asfendiyarov Kazakh National Medical University, Tole-bi 94, Almaty 050012, Kazakhstan
2
Chair and Department of Pharmaceutical Microbiology, Faculty of Pharmacy, Medical University of Lublin, 1, Chodzki str., 20-093 Lublin, Poland
3
Department of Cosmetology, University of Information Technology and Management in Rzeszów, 2 Sucharskiego str., 35-225 Rzeszów, Poland
4
Department of Pharmaceutical Technology and Biochemistry and BioTechMed Centre, Faculty of Chemistry, Gdańsk University of Technology, Gabriela Narutowicza Str. 11/12, 80-233 Gdańsk, Poland
5
Edinburgh Dental Institute, The University of Edinburgh, Edinburgh EH8 9YL, UK
6
Department of Fundamental Medicine, Al-Farabi Kazakh National University, Almaty 050040, Kazakhstan
7
Department of Chemical and Biochemical Engineering, Geology and Oil-Gas Business Institute Named after K. Turyssov, Institute of Chemical and Biological Technologies, Satbayev University, Almaty 050043, Kazakhstan
8
Ecology Research Institute, Khoja Akhmet Yassawi International Kazakh-Turkish University, Turkistan 161200, Kazakhstan
9
Department of Pharmacognosy, Medical University of Lublin, 1, Chodźki str., 20-093 Lublin, Poland
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(24), 8994; https://doi.org/10.3390/molecules27248994
Submission received: 14 November 2022 / Revised: 9 December 2022 / Accepted: 11 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Plant Metabolites: Accumulation, Profiling and Bioactivity)

Abstract

:
Ziziphora species (Lamiaceae) have been used in traditional medicine as sedatives, antiseptics, carminatives, or expectorants. Despite their common applications in phytotherapy, there is still lack of evidence about the composition of their extracts and its impact on biological properties of the plants. The aim of this study was to evaluate the content of Ziziphora bungeana, a less studied species growing in Kazakhstan, using HPLC-ESI-QTOF-MS/MS instrumentation and to determine its antimicrobial, antioxidant, and cytotoxic activity together with inhibitory properties against tyrosinase and toxicity in erythrocyte lysis assay. Extracts from Z. bungeana were found to be sources of flavonoids, phenolic acids, organic acids, and terpenes that determined their antiradical activity. The minimum inhibitory concentrations of extracts were lower for Gram-positive bacteria (1.25–10 mg/mL) than for Gram-negative bacteria and fungi (5–20 mg/mL). The EC50 value calculated for antiradical activity ranged between 15.00 ± 1.06 µg/mL and 13.21 ± 3.24 µg/mL for ABTS and DPPH assays, respectively. Z. bungeana extracts were found to decrease the activity of tyrosinase by 50% (at 200 µg/mL) similarly to kojic acid and were slightly cytotoxic for human melanoma A375 cell line (at 200 µg/mL) with no effect on HaCaT keratinocytes. In the end, Z. bungeana did not reveal toxic effects in hemolytic assay as compared to the positive control Triton X-100. The performed tests show potential application of the plant in the treatment of infectious diseases, disorders caused by free radicals, and skin problems.

1. Introduction

Plant biodiversity is certainly an invaluable wealth of our planet. Herbal remedies have been used in the treatment of various diseases for centuries as plants are capable of synthesizing a wide range of metabolites of pharmacological significance. Traditionally plants were used as whole, in the form of the total extracts to enable the combined action of their ingredients. From a therapeutic point of view, the natural products, organic low-molecular structures called secondary metabolites, are beneficial to both the plants, as they protect them from harsh environmental factors or pathogens, and to humans for their medicinal or cosmetic applications. A multitude of compounds that are encountered in the world’s flora inspires the scientists to modify the original scaffoldings and to create more active, stable, or less toxic molecules to fight the diseases that are difficult to treat with conventional drugs. This work presents a compositional study and biological activity determinations that are performed on the Ziziphora bungeana species collected in Kazakhstan. This species will be tested for its composition and activity, the results of which will allow for determining of antimicrobial properties and potential application of the plant as an additive to drugs or foods, an ingredient of cosmetic products with anti-aging, anti-cancer, or whitening properties for the use in skin cancer therapy, and whether these extracts are safe. The undertaken direction of studies was inspired by former research data that proved a diversified composition of Ziziphora spp. and a marked biological potential of its extracts.
The genus Ziziphora L. (Lamiaceae) comprises about 30 species widespread all over Asia, Africa, and Europe that represent the prototypical example of the Lamiaceae family. Plants from this genus are known to produce monoterpenes, triterpenes and phenolic substances belonging mainly to the group of flavonoids [1,2,3]. Essential oil (EO) from Ziziphora clinopodioides and Z. tenuior was proven to be rich in pulegone, menthone, and limonene, whereas more polar extracts were characterized by a high phenolic content [4,5]. Rich composition in secondary metabolites allows for a multi-directional application of the plants. They are used in the form of infusions, decoctions and macerates as remedies for stomachache, common cold, inflammatory conditions, cough, migraine, fever, depression, diarrhea, and gastrointestinal diseases. Moreover, they are known from their sedative, expectorant, antiseptic, and carminative properties [2,4,6,7,8,9]. In Kazakh traditional medicine, Ziziphora species possess several medicinal uses. In particular Z. bungeana and Z. clinopodioides are administered in the treatment of cardiovascular system disorders or infections [10].
Despite the presence of polyphenols in the extracts that have impact on the antioxidant potential of the mentioned plant species, the constituents of less polar fractions may also influence the total antiradical potential of the plant, affect its antispasmodic, anti-inflammatory, anti-infective, and expectorant properties, which broadens their activity range [11].
Having in mind plentiful applications of different Ziziphora species, the aim of the study was to deliver information about another less known Ziziphora bungeana Juz., which is a synonym of Z. clinopodioides ssp. bungeana. This plant is distributed mainly in Kazakhstan, China, Central Asia, and Mongolia [12]. The scientific literature is still lacking sufficient information about its composition and beneficial actions, including the antimicrobial, antifungal, antioxidant, or skin whitening potential, whose determination will be performed in this study together with compositional analysis by HLPC-ESI-QTOF-MS/MS technique. Based on the scientific literature data on Ziziphora gender, this plant is expected to be a good natural antioxidant and antibacterial agent with high potential for its application in foods and cosmetics as a preservative and antiradical component. This study was designed to meet the diverse expectations towards plant extracts in the context of their potential use in the treatment of civilization diseases or in the skin care because of the harmful effects caused by the environment. The determination of the antioxidant, antimicrobial, and whitening properties of the extracts of different polarities will allow for the study on the potential use of Ziziphora bungeana in cosmetics. Moreover, the toxicity studies in relation to normal and cancer skin cells will bring evidence for a discussion about safety of its use. On the other hand, the qualitative analysis of different polarity extracts will develop a fingerprint responsible for the determined action of the plant.

2. Results and Discussion

2.1. Chemical Composition of Extracts by HPLC-ESI-QTOF-MS/MS

The applied chromatographic method was capable of separating metabolites present in the extracts from the aerial parts of Z. bungeana, whereas the applied mass spectrometer settings provided MS/MS spectra that helped in the identification of metabolites from different groups.
The HPLC-ESI-QTOF-MS/MS analysis confirmed the presence of phenolic acids, flavonoids, triterpenes, and monoterpene glucosides in the tested extracts. The list of 26 tentatively identified components is presented in the Table 1 and Figure 1, whereas the recorded mass chromatograms from positive and negative ionization modes are shown in the Supplementary File (Figures S1 and S2). It was proven that the components of Z. bungeana were previously described in other species from the same gender.
The applied LC-MS technology enabled the identification of components, e.g., thymol, carvacrol, or ziziphorosides, that were previously determined in Ziziphora spp. by GC-MS technique. Milder fragmentation conditions (fragmentation voltage of 110 V, capillary voltage of 3000 V, gas temperature of 275 °C and collision energy of 10 V) increased the chance to observe terpene compounds in the chromatogram in the liquid chromatography-based system. Ziziphoroside isomers, possibly ziziphoroside A, B, and C, were present in the mass spectra in the form of adducts with sodium ions. The remaining compounds were traced in the form of molecular ions with or without a proton.
Z. bungeana was proven to contain different types of metabolites whose presence was revealed in the HPLC-MS assessment of the extracts tested in this study. Among them flavonoids constituted the leading group of components, followed by phenolic acids and, interesting from the structural point of view, terpenes, e.g., ziziphorosides or schizonepetaside E. It is worthwhile to note that the scientific literature still lacks sufficient information about the composition of this plant species. For the moment, to the best of the authors’ knowledge, there is only one original manuscript that discusses the composition of extracts based on the HPLC-MS results. The researchers confirmed the presence of twelve flavonoids in Z. bungeana, that included: kaempferol-7-O-rutinoside, kaempferol-3-O-rutinoside; rutin; apigenin-7-O-rutinoside; 3′-hydroxyacacetin-7-O-rutinoside; acacetin-4′-O-rutinoside; pinocembrin-7-O-rutinoside; chrysin-7-O-rutinoside; linarin; 5,7,3′-trihydroxy-6,4′,5′-trimethoxyflavone, 5,4′-dihydroxy-6-methoxy-7,8-methylenedioxyflavone, and 5,7-dihydroxy-6-methoxyflavone. The above list of components presented in the Table 1 expands information on the composition of this species.
The analysis of previously published papers provided a more detailed list of components of Ziziphora genus that helped to enrich the list of the tentatively identified metabolites of Z. bungeana extracts. Previous investigation of the chemical profile of ethyl acetate, methanol, and water extracts from the aerial parts of Ziziphora taurica subsp. cleonioides showed that among the identified compounds, rosmarinic acid and chlorogenic acid were the most abundant components of the methanol extract with the calculated concentration of 3375.67 ± 38.02 and 3225.10 ± 16.44 µg/mL, respectively [32]. Both compounds were also determined in the studied species. Moreover, flavonoids constituted the major group of bioactive compounds present in Ziziphora clinopodioides Lam. [2] together with organic acids, alkaloids, and glycosides that were listed by other authors [33].
As mentioned above, Ziziphora species belong to the plants that synthesize secondary metabolites from different classes, which explains their various therapeutical applications. In the previous studies Ziziphora spp. were found to be rich sources of phenolic compounds. Interestingly, both extracts and the EO were sources of polyphenols. For example, the measured total phenolic content in Z. tenuior was equal to 49.0 ± 1.4 mg mg gallic acid/100 g of EO [15].
The metabolites that are mentioned above are important from the pharmacological point of view. Phenolic acids and flavonoids are known scavengers of free radicals that are efficiently inhibiting the progression of different inflammatory conditions and civilization diseases progressing with an important role of radicals [34]. Their presence in the final extracts from edible plants is certainly related to the type of the plant, but also to the extraction conditions [35].
Based on this information, the authors found it crucial to study the biological potential of Z. bungeana and focus on its antimicrobial, antioxidant, and anti-tyrosynase properties, as well as to evaluate safety.

2.2. Antimicrobial Activity Assessment

The data presented in Table 2 and Table 3 indicate that the extracts from Ziziphora bungeana showed some potential antimicrobial activity. They were more effective against reference Gram-positive bacteria than towards Gram-negative bacteria and yeasts. The lowest concentration of extracts which inhibited the growth of these microorganisms or killed them ranged from 1.25 mg/mL to 20 mg/mL and from 2.5 mg/mL to 20 mg/mL, respectively.
As shown in Table 2, in the case of Gram-positive bacteria, the MIC values of the extracts were in the range of 1.25–10 mg/mL. Their activity was the same towards staphylococci, both Staphylococcus aureus ATCC 43,300 (MRSA—Methicillin Resistant S. aureus), S. aureus ATCC 29,213 (MSSA—Methicillin Susceptible S. aureus) and Staphylococcus epidermidis ATCC 12,228 with MIC = 2.5 mg/mL and MBC = 2.5–5 mg/mL (except S. aureus ATCC 29213; MIC = 5 mg/mL for Z2). In turn, Micrococcus luteus ATCC 10,240 was the most (MIC = 1.25 mg/mL and MBC = 5 mg/mL) and the least sensitive (MIC = 5 mg/mL and MBC = 10 mg/mL) to Z2 and Z3, respectively. The antibacterial effect against two reference Bacillus spp. strains was lower. MIC values were mainly 5 mg/mL. However, MIC = 2.5 mg/mL was shown for the Z3 against B. subtillis ATCC 6633 and MIC = 10 mg/mL for Z1 towards B. cereus ATCC 10876. MBCs of all extracts for bacilli were the same—10 mg/mL.
The activity of extracts towards Gram-negative rods-shaped, was slightly weaker with MIC = 5–20 mg/mL and MBC = 10–20 mg/mL. Among them, Z2 showed the highest effect towards Bordetella bronchiseptica ATCC 4617 (MIC = 5 mg/mL and MBC = 20 mg/mL). For Z1, MIC = 10 mg/mL against B. bronchiseptica and Pseudomonas aeruginosa ATCC 9027 was shown. In the case of Z2, the same MIC value towards reference Klebsiella pneumoniae strain was indicated. The growth of remaining Gram-negative bacteria was inhibited by these extracts at a concentration of 20 mg/mL.
Taking into account the MBC/MIC and MFC/MIC ratios, as presented at Figure 2, it was shown that extracts from Z. bungeana had a beneficial bactericidal or fungicidal effect towards reference microorganisms. The values of MBC/MIC or MFC/MIC were in the range 1–4. MICs equal to MBC or MFC (MBC/MIC = 1 and MFC/MIC = 1) were shown for most, i.e., 12 (75%) and 11 (68.75%) strains in the case of Z1 and Z3, respectively. For Z2, these values were different. However, these ratios were mainly 1 (for six (37.5%) strains) and 2 (for seven (43.75%) strains). The value of 4 was found the least frequently (only for Z2 and Z3). The bacteriostatic effect (MBC/MIC > 4 or MFC/MIC > 4) of the tested extracts was not demonstrated.
In the next stage of this study, the total antimicrobial activity (TAA) was assessed. The total antibacterial activity or total antifungal activity of the studied extracts Z1–Z2 was shown in Table 4. The MIC and TAA, both total antibacterial activity and total antifungal activity values are important pharmacological tools. They are useful in determining the activity of extracts in mg/mL (potency) of plants extracts for isolating bioactive compounds and total activity on mL/g (efficacy) for the selection of appropriate plant species [36]. Generally, their TAA values were the highest against Gram-positive bacteria: Staphylococcus spp., Micrococcus luteus, followed by Bacillus spp. (3.35 ± 0.0–20.20 ± 0.0 mL/g) and the lowest against Gram-negative bacteria and fungi belonging to Candida spp. (1.09 ± 0.0–4.21 ± 1.46 mL/g). Z1 had higher TAA towards reference strains of S. aureus, S. epidermidis and M. luteus (13.39 ± 0.0 mL/g). TAA values of Z2 varied slightly and were in the range 4.21 ± 1.46 to 20.20 ± 0.0 mL/g against these bacteria. In turn, TAA values of Z3 were slightly lower (4.38 ± 0.0–8.75 ± 0.0 mL/g).
As shown in Figure 3, Z1, Z2, and Z3 extracts had the mean total antibacterial activities of 6.49, 6.01 and 4.41 mL/g, respectively. In turn, the mean total antifungal activity of extracts was lower in the range 1.24–2. mL/g. In conclusion, Z1 and Z2 had a similar and better efficacy against both bacteria and fungi compared to Z3. The higher the TAA value, the more efficacious the plant extract.
The results indicated that Z1–Z3 extracts from Z. bungeana showed some antimicrobial activity with bactericidal or fungicidal effect. Among all studied reference microorganisms, Gram-positive bacteria were the most sensitive to them. The lowest concentrations of Z1–Z3 extracts which inhibited the growth of the tested microorganisms or killed them ranged from 1.25 mg/mL to 20 mg/mL and from 2.5 mg/mL to 20 mg/mL, respectively. Overall, the Gram-positive bacteria were more sensitive to the extracts than the Gram-negative bacteria and yeasts from Candida. The difference in the sensitivity between these microorganisms may be due to the variation in their cell wall structure. The Gram-positive bacterial cell wall consists of 70–100 layers of peptidoglycans. Peptidoglycan is comprised of two polysaccharides, N-acetyl-glucosamine and N-acetyl-muramic acid cross-linked by peptide side chains and cross bridges [37]. It is possible that active compounds from extracts can easier break important bonds in cell wall structure in these bacteria. On the other hand, the cell wall of Gram-negative bacteria is far more complex, and it is among other things the reason they are more resistant for biologically active compounds.
There is little information in the literature on the biological activity of Z. bungeana extracts. However, there are reports on other Ziziphora species. Some authors showed antimicrobial effect of different extracts, EO or selected compounds derived from Ziziphora gender. The results and findings described herein are in accordance with some other studies.
The antibacterial activity of EO and its two main components (pulegone and 1,8-cineole) obtained from the aerial flowering parts of Ziziphora clinopodioides subsp. bungeana (Juz.) Rech. f. was analyzed by Sonboli A. et al. [38] against seven bacteria. It was found that the EO exhibited interesting activity against S. epidermidis, S. aureus, E. coli, and B. subtilis with MIC values of 3.75 mg/mL. These results were similar to ours for Gram-positive bacteria.
In turn, the inhibitory effect of methanol extract and EO from Ziziphora persica was tested against 98 strains belonging to 51 bacteria species by standard dilution methods. The results showed that both extract and EO had antibacterial activity against many tested bacteria. The lowest MIC values (7.81 µg/mL) of EO were obtained against Bacillus dipsauri, Corynebacterium cystitidis, and Corynebacterium flavescens [39].
Z. clinopodioides was studied by subsequent researchers. The LC-MS/MS results of Özkan E.E. et al. [3] indicated that quinic acid, malic acid and rhoifolin are the abundant compounds in aerial and root ethanol extracts of Z. clinopodioides. Both extracts exhibited moderate antifungal activity with MIC = 39.06 μg/mL against Candida tropicalis. Moreover, these extracts showed some better or the same antibacterial effect against reference S. aureus, S. epidermidis, and E. faecalis strains (MIC = 0.312–1.25 mg/mL) as our extracts of Z. bungeana. In the case of Gram-negative bacteria (P. aeruginosa, E. coli, K. pneumoniae and P. mirabilis) and C. albicans, no activity was showed.
Moreover, the studies of Anzabi Y. et al. [40] showed that the Z. clinopodioides EOs was effective on many tested bacteria and can be used as natural antimicrobial drug against microorganisms causing urogenital tract infections in women. The aerial parts of Z. clinopodioides were also screened by other authors [41] for their possible antimicrobial activities. Methanol extract was found to have moderate antimicrobial activity against some microorganisms tested. Acinetobacter lwoffii and Candida krusei were the most sensitive for this extract. The antimicrobial properties were also found in Z. clinopodioides EOs collected from provinces in western Iran. The studied EO inhibited the growth of Listeria monocytogenes, S. typhimurium, E. coli O157:H7, B. subtilis, B. cereus, and S. aureus at MIC values between 0.03% and 0.04%. The Gram-positive bacteria were the most susceptible to it, while Gram-negative bacteria were resistant [1]. The interesting antibacterial activity against seven Gram-positive or Gram-negative bacteria exhibited also EO and methanol extract of Z. clinopodioides subsp. rigida (BOISS.) RECH. f. from Iran. The obtained results indicated that B. subtilis was the most sensitive microorganism to this EO, with the lowest MIC = 3.8 mg/mL. The growth inhibition of S. epidermidis and S. aureus was observed at similar MIC = 7.5 mg/mL. The inhibitory activity of EO against E. faecalis, K. pneumoniae, and E. coli was also determined with MIC values equal to or greater than 15 mg/mL. No activity was observed against P. aeruginosa [42].
The subsequent results of Hazrati et al. [2] showed 17 and 21 different compounds (comprising 99.7% of total EO) in Z. clinopodioides and Z. tenuior, respectively. The major identified compounds in EO analysis reported as pulegone and menthone for Z. clinopodioides, or pulegone and limonene for Z. tenuior. Both Ziziphora species were also rich in phenolic compounds. These authors investigated the antibacterial activity of EOs against important foodborne pathogenic bacteria and showed that they could be considered as good sources of natural antibacterial material as well as food preservative [2,15].
Additionally, Celiket al. [43] evaluated the antimicrobial and anti-biofilm properties of Z. tenuior. EO against multi-drug resistant Acinetobacter baumannii with MIC = 0.6–1.25 μL/mL and MBC = 2.5–5.0 μL/mL. Furthermore, minimal biofilm inhibition concentration (MBIC) values of 0.3–1.25 μL/mL and minimal biofilm eradication concentration (MBEC) values of 5–10 μL/mL were observed.
Considering all above information, it can be concluded that Ziziphora plants may deliver extracts that are important from a pharmacological point of view, as they may help combat the occurring bacterial and fungal infections.

2.3. Antioxidant Activity Assessment

Antioxidant activity of Ziziphora bungeana extracts was compared using DPPH and ABTS radical scavenging assays (Figure 4A) and determination of superoxide dismutase (SOD) activity (Figure 4B). In respect of the radical scavenging potential, extract Z3 showed the most significant activity with EC50 values of 15.00 ± 1.06 µg/mL and 13.21 ± 3.24 µg/mL for ABTS and DPPH assays, respectively. All tested extracts also showed significant SOD activity, dependent on the extract concentration. The most effective was extract Z1, showing > 90% SOD activity in all three tested concentrations. Extract Z2 was the least effective. At the concentration of 50 µg/mL, the mean SOD activity detected for this extract was 64.4 ± 0.55%.
The antioxidant activity of 50% (v/v) ethanolic extract from Z. bungeana was recently compared with other plants from Lamiaceae family by measuring its influence on the level of lipid peroxidation in the liver microsome and the membrane-stabilizing properties [44]. The antioxidant potential of the extracts was significant in both assays but moderate in comparison with other Lamiaceae plants.
However in the publication of Gursoy and co-investigators [41], the aerial parts of Z. clinopodioides were found to be the strongest radical scavengers among the tested species from Lamiaceae family, namely: Z. clinopodioides, Cyclotrichium niveum, and Mentha longifolia subsp. typhoides var. typhoides in the DPPH and beta-carotene/linoleic acid assays. The calculated IC50 values for the tested extracts were 37.73 ± 1.18 µg/mg for DPPH and 83.56 ± 1.19% in the inhibition capacity of the linoleic acid. Moreover, the total phenolic content of its methanolic extract was the highest among the tested species and was equal to 129.55 +/− 2.26 µg/mg.
Recently, the antioxidant activity of the aqueous, ethyl acetate, and methanolic extracts from other Ziziphora species, Ziziphora taurica subsp. taurica, were compared using DPPH and ABTS scavenging assays. In both assays, the methanolic ectract from Z. taurica was the most effective with IC50 values of 5.74 ± 0.08 mg/mL and 2.74 ± 0.10 mg/mL for DPPH and ABTS scavenging, respectively. The EC50 values obtained in our study suggest that the antioxidant potential of Z. bungeana extracts is higher than that of Z. taurica [45].

2.4. Tyrosinase Activity Assay

Tyrosinase (EC 1.14.18.1) is a cooper containing metalloenzyme catalyzing the first two, rate-limiting steps of mammalian melanogenesis. Neither increased nor decreased activation of tyrosinase is desirable as it may lead to hyper- or hypopigmentation disorders, respectively. Natural extracts and compounds with tyrosinase inhibitory activity are particularly desired by the cosmetic industry as they serve as active ingredients in skin lightening cosmetics and rituals [46]. On the other hand, the compounds increasing the activity of tyrosinase might be considered as topical treatment for vitiligo [47].
Investigating the influence of novel extracts and compounds on tyrosinase activity is commonly performed using the assay utilizing commercially available mushroom tyrosinase, incubated with its substrate L-3,4-dihydroxyphenylalanine (L-DOPA), in the presence or absence of tested compound. Despite the low costs, simplicity, and high throughput of the procedure, the assay has several limitations, resulting from substantial differences between mushroom and mammalian tyrosinase [48]. Therefore, the influence of plant-derived extracts and compounds on the activity of mushroom and mammalian tyrosinases may vary significantly [49,50].
As shown in Figure 5, none of the analyzed Ziziphora extracts significantly inhibited mushroom tyrosinase up to the concentration of 200 µg/mL (Figure 5B). Extract Z2 slightly increased the activity of this enzyme at 25 and 50 µg/mL. In respect of the murine tyrosinase all Ziziphora extracts showed significant, dose-dependent inhibitory potential (Figure 5A). The most potent inhibitor of murine tyrosinase was extracts Z2, decreasing the activity of tyrosinase by 50% at 200 µg/mL which was comparable with the inhibitory activity of kojic acid (KA), a tyrosinase inhibitor widely used in skin lightening cosmetics [49]. Extracts Z3 showed the lowest activity, significantly decreasing the activity of tyrosinase only at the highest analyzed concentration (200 µg/mL). To our knowledge, this is the first study investigating the effect of Ziziphora spp. extract on the activity of mammalian tyrosinase.
Several phytochemicals identified in Ziziphora extracts were described in scientific literature as effective mushroom tyrosinase inhibitors, including acetophenone [51] identified in Z. tenuior [52] and cuminyl aldehyde (syn. cumaldehyde) [53] from Z. clinopodioides subsp. rigida [54]. The last compound was also shown to suppress melanin synthesis in B16F10 murine melanoma cells [55].
Z. clinopodioides extracts were previously found to exhibit weak tyrosinase inhibitory potential. The extracts from the overground parts of the plant exhibited weak inhibitory potential against the enzyme at the concentration of 200 μg/mL with 8.60 ± 0.87% inhibition compared to kojic acid (KA) for whom the inhibition percentage was calculated as 95.26 ± 0.23% [3].
Mushroom tyrosinase activity was also analyzed by Tomczyk and co-workers in respect of Z. taurica extracts [45]. The IC50 values for aqueous, ethyl acetate and methanolic extracts were 2.29 ± 0.13, 1.37 ± 0.07 and 1.46 ± 0.06 mg/mL, respectively. These values suggest that Z. bungeana extracts Z1, Z2, and Z3 might be effective against mushroom tyrosinase, but at higher concentrations than tested in this study.

2.5. Cytotoxic Activity

Ziziphora spp. were shown to contain several compounds with broad range of anticancer activities, including pulegone, menthol, menthone, cineole, piperitone, isomenthol, and curcumin [10]. However, only a few recent studies described the anti-cancer potential of whole Ziziphora extracts and EOs [56,57,58].
In this study, the authors focused on the assessment of Z. bungeana cytotoxic effect on human and murine melanoma cells (Figure 6B–D) in comparison with a known chemotherapeutic agent, 5′fluorouracil (5′FU). Human keratinocytes HaCaT served as noncancerous control cells (Figure 5A).
Extract Z1 at 200 µg/mL was slightly cytotoxic for human melanoma A375 cell line, reducing the number of viable cells by ca. 20%. It was not cytotoxic for HaCaT keratinocytes, B16F10, and SKMEL-3 melanoma cells. Extract Z2 at 200 µg/mL significantly reduced the number of viable A375 and SK-MEL3 melanoma cells by ca 28% and 23%. However, it showed comparable cytotoxicity towards HaCaT keratinocytes. Extract Z3 was cytotoxic only for B16F10 murine melanoma cells, reducing their viability by 15% at 200 µg/mL.
The cytotoxicity of Ziziphora spp. Extracts towards murine and human melanoma cells as well as human noncancerous skin cells has not been described in the scientific literature to date. Several compounds found in Ziziphora preparations, such as extracts and EO, including linalool and α-terpineol (Z. clinopodioides), carvacrol (Z. tenuior, Z. clinopodioides), thymol (Z. tenuior), and terpinen-4-ol (Z. clinopodioides), are known to induce apoptosis in melanoma cell lines [11]. In the study of Yousefbeyk and co-investigators [59] Z. clinopoides n-hexane extract that was found rich in pulegone, menthone and menthol exhibited strong cytotoxic activity against K-562 and T-47D cell lines with EC50 values of 80 ± 2.56 μg/mL and 77.41 ± 12.89, respectively. Interestingly, more polar fractions did not show cytotoxic effects.
Available scientific data on the in vitro cytotoxic effect of Ziziphora spp. preparations were obtained using EOs. Azimi and co-workers showed that the EO from Z. tenuior induces apoptosis in human colorectal cancer cells HT-29 in a concentration range of 50–200 µg/mL. The apoptotic effect was mediated by increased caspase 3 and 9 expression at mRNA and protein levels and decreased levels of Bcl-2 [56]. Ghavan et al. showed that Z. clinopodioides subsp. rigida EO is cytotoxic for human ovarian cancer cells (OVCAR-3) [60].

2.6. The Hemolytic Activity Assay (Toxicity towards Erythrocytes)

In the present studies, the toxicity of Z1–Z3 extracts from Z. bungeana towards red blood cells was calculated in vitro hemolytic assay. The erythrocyte model (erythrocyte lysis assay; ELA) was used to analyzed their effect on cell membrane [60]. The results revealed that studied extracts exhibit negligible toxicity as compared to the positive control Triton X-100 (100% erythrocyte lysis). As presented in Figure 7, hemolytic activity of each extract was related to their concentration. The highest concentrations of the studied extracts (20 mg/mL) showed some hemolytic activity in the range 6.1–30%. In turn, their concentrations that did not exert any hemolytic effect ranged from 1.25 to 2.5 mg/mL and the percentage of lysed red blood cells of 0–4.5 was within the permissible limit of 5% hemolysis [61]. The Z1 and Z3 extracts exhibited lower hemolytic activity (0–18.9%) than Z2 (2.5–30%) and did not affect the stability of the erythrocyte membrane. Data obtained using ELA confirm that antimicrobial effect, especially against Gram-positive bacteria (MIC = 1.25–5 mg/mL), was observed at non-cytotoxic concentrations of extracts from Z. bungeana.
The erythrocyte model presents general indication of membrane toxicity. The red blood cell membrane shows similarity to other cell membranes. Hemolysis is due to erythrocyte cells destruction resulting from lysis of the membrane lipid bilayer [60,61]. The obtained data using ELA confirm antibacterial activity of Z1–Z3 extracts at their non-cytotoxic concentrations (MIC = 1.25–5 mg/mL) against staphylococci, micrococci, and some bacilli. Therefore, it seems practical to use these extracts in the future in the prevention and treatment of infections caused by selected microorganisms, mainly Gram-positive bacteria.

2.7. Chemometric Assessment

Principal component analysis (PCA) was conducted separately for the Z1–Z3 relative compositions (Crel) and the Z1–Z3 relative activities (Arel) (See Tables S2 and S3 in the Supplementary File), whereas the extracts were treated as vectors defined on Crel or Arel values. The resulting spaces were two-dimensional, since in both cases the first two principal components extracted almost 100% of the information (expressed as variance) of studied ensembles. In the case of the Crel system (Figure 8A), the relative composition of the Z3 and Z1 extracts were reversely correlated, constituting the first dimension of the studied dataset, whereas the Z2 extract exhibited rather unique proportions of the selected eleven analyzed compounds (named in the Table 1 as 1, 2, 4, 6, 9, 10, 11, 13, 15, 23, 26—see Table S2 in the Supplementary File), defining the second dimension of the vector space. The obtained conclusions are logical, as dichloromethane is characterized by a much lower polarity than ethanol and water, and that is why the extracts obtained using dichloromethane can show a different fingerprint from alcoholic or water ones. In the case of Arel dataset (Figure 8B), the relative biological activities of Z2 and Z3 extracts were reversely correlated (Dimension 1), while the Z1 extract exhibited different properties, possibly thanks to other metabolites, e.g., peptides, sugars, or proteins whose identity was not analyzed in this study, solely explaining the second dimension of the studied space.
The linear maps of the studied compounds (Figure 9A) and the activity tests (Figure 9B) were presented in the same spaces as the respective linear maps of the Z1–Z3 vectors in Figure 8. In the case of Crel ensemble, the relative amounts of the 1, 2, 6, 9, 11, and 13 compounds were very similar for all three extracts. Z3 clearly excelled in the relative amounts of 23 and 26, whereas it was very low on 10′s concentration. On the contrary, Z1 contained impressive amounts of 10 while lacking 23 and 26. Z2 was quite low in 4 and 10 yet exhibited higher-than-average amounts of 15 and 26. While taking into account the relative values of activity tests, Z3 excelled with the B16F10 cell line (I and II), on the contrary to Z1. Moreover, Z2 was quite good at IV (A375 cell line) and V (SKMEL-3 cell line), yet toxic to HaCaT cells (III), whereas Z1 exhibited poor activity at VIII (SOD assay).
While comparing the above with the distribution of the relative compositions of the eleven analyzed compounds (named in Table 1 as 1, 2, 4, 6, 9, 10, 11, 13, 15, 23, 26) within the Z1–Z3 extracts, one might conclude that eventual toxic effects towards HaCat cell line (III) at the highest dose, exhibited solely by Z2, could result from the presence of very high, relative amounts of the compounds 15 and 26. These findings may be due to the fact, that the analyzed concentration used in the calculations was high and exceeded safe doses for both diosmin and ursolic acid, respectively. In the meantime, Z2 utterly failed at the tests I and II (B16F10 cell line), similarly to Z1. Since Z3 succeeded at I and II, while it was also rich with the compounds 23 and 26, the good result at I and II could be directly associated with high relative amounts of 23 (acacetin). Finally, Z1 extract did relatively well at the tests IV (A375 cell line) and VII (mushroom tyrosinase assay), while it exhibited high amounts of a ziziphoroside isomer 2 (10). Possibly, its presence influences the total activity of the extract. Previously, other species of Ziziphora were proven to inhibit tyrosinase [10]. On the basis of the resulting images, no other compounds could be related to the biological activities of Z1–Z3 extracts in a straightforward manner.

3. Material and Methods

3.1. Materials

3.1.1. Plant Material

The aerial parts of Ziziphora bungeana Lam. were collected in the summer of 2021 in the flowering stage in the Turkestan region of the Republic of Kazakhstan and identified by the Institute of Botany and Phytointroduction, Science Committee, Ministry of Education and Science of the Republic of Kazakhstan. A voucher sample (№01-05/337 from 5 October 2021) has been deposited in the herbarium of the Institute of Botany and Phytointroduction, Almaty, Republic of Kazakhstan.

3.1.2. Microorganisms

The reference strains of microorganisms from American Type Culture Collection (ATCC) (Manassas, VA, USA) were used in the study. The representative Gram-positive bacteria were: Staphylococcus aureus ATCC 29,213 (Methicillin Susceptible Staphylococcus aureus—MSSA), Staphylococcus aureus ATCC 43,300 (Methicillin Resistant Staphylococcus aureus—MRSA), Staphylococcus epidermidis ATCC 12228, Enterococcus faecalis ATCC 29212, Micrococcus luteus ATCC 10240, Bacillus subtilis ATCC 6633 and Bacillus cereus ATCC 10876), while those of Gram-negative bacteria: Escherichia coli ATCC 25922, Klebsiella pneumoniae ATCC 13883, Salmonella typhimurium ATCC 14028, Pseudomonas aeruginosa ATCC 9027 and Bordetella bronchiseptica ATCC 4617. Moreover, the fungi belonging to yeasts: Candida albicans ATCC 10231, Candida albicans ATCC 2091, Candida parapsilosis ATCC 22019, Candida glabrata ATCC 90030, and Candida krusei ATCC 14,243 were used.

3.2. Methods

3.2.1. Extraction Procedure

First, the aerial parts of the plant were powdered using an electric mill (type WZ-1, ZBPP, Poland). Next, 10 g portions of the aerial parts of the plant were divided into three parts to provide three extracts using the following extracting solvents—water (Z1), dichloromethane (Z2) and 96% ethanol (Z3). After adding 50 mL of solvents, the extraction was performed three times, 30 min each, at room temperature using an ultrasonic bath with no heating. Then, the extracts were centrifuged at 3500 rpm for 10 min, filtered through a nylon syringe filter (pore diameter 0.22 µm), and evaporated in the weighted vials using the Eppendorf Concentrator Plus (Hamburg, Germany) at the temperature of 45 °C. Weighted samples were kept in the freezer before chromatographic studies and bioactivity evaluations.

3.2.2. The HPLC-ESI-QTOF-MS/MS Analysis

Compositional studies of Z. bungeana extracts were performed using an HPLC-MS platform produced by Agilent Technologies (Santa Clara, CA, USA) which was composed of an HPLC chromatograph equipped in a binary pump (G1312C), a degasser (G1322A), an autosampler (G1329B), a photodiode array detector (DAD) (G1315D), and a QTOF-MS/MS mass spectrometer (G6530B).
The extracts’ constituents were separated in gradient method composed of 0.1% formic acid (solvent A) and acetonitrile with the addition of 0.1% formic acid (solvent B) in the following program: 0 min: 10% B, 10 min: 20% B, 15 min: 40% B, 17–22 min: 95% B, 22.10 min: 10% B. The run lasted 30 min, the flow rate was set at 0.200 mL/min and the injection volume was set at 5 µL, and the concentration of the extracts was 10 mg/mL. Chromatographic separation was performed on the RP-18 chromatographic column (dimensions: 150 mm × 2.1 mm; dp = 3.5 µm) (Zorbax Eclipse Plus by Agilent Technologies, Santa Clara, CA, USA).
The detection on the mass spectrometer was achieved in the following settings, using both negative and positive ionization mode: m/z range of 100–1700 Da, capillary voltage of 3000 V, gas and sheath gas temperatures of 275 and 325 °C, gas flows of 12 L/min, respectively, fragmentation voltage of 110 V, skimmer voltage of 65 V, and collision energies of 10 and 20 V. In the used method, the MS/MS spectra were recorded for the two most intense peaks per scan The structure determination was based on the fragmentation spectra, literature data, retention times, and open databases (Metlin).

3.2.3. In Vitro Antimicrobial Activity Assay

The three extracts Z1–Z3 from Ziziphora bungeana were investigated in vitro for antibacterial and antifungal activities. In these studies, the broth microdilution was used. The tests were performed in accordance with the guidelines of the European Committee on Antimicrobial Susceptibility Testing (EUCAST) [14,62,63] and Clinical and Laboratory Standards [64,65]. The used microbial cultures were first subcultured and on nutrient agar (for bacteria) or Sabouraud agar (for fungi) (BioMaxima S.A., Lublin, Poland) and incubated at 35 °C for 18–24 h. Microbial suspensions were prepared in sterile saline (0.85% NaCl) with an optical density of 0.5 McFarland standard scale (1.5 × 108 CFU/mL (CFU—Colony Forming Units/mL) for bacteria and 5 × 106 CFU/mL for yeasts). Samples containing examined extracts were first dissolved in dimethyl sulfoxide (DMSO) to the concentration of 200 mg/mL. The minimal inhibitory concentration (MIC) of these extracts was evaluated by the microdilution broth method in 96-well polystyrene plates. In this study, two-fold dilutions of the extracts in selective broth, Mueller-Hinton (MH) (BioMaxima S.A., Lublin, Poland) for bacteria and RPMI (Roswell Park Memorial Institute) 1640 with MOPS (3-(N-Morpholino)propanesulfonic acid) (Sigma-Aldrich Chemicals, St. Louis, MO, USA), were performed. The final concentrations of extracts (diluted in broth) ranged from 20 to 0.156 mg/mL.
Next, the bacterial or fungal suspensions were introduced into each well of the microplate to obtain final density of 1.5 × 106 CFU/mL for bacteria and 5 × 104 CFU/mL for yeasts. After 18–24 h incubation at 35 °C, the MIC value was assessed in the BioTek spectrophotometer (Biokom, Janki, Poland) as the minimal concentration of the samples that showed complete microbial growth inhibition. The inhibition of bacterial and fungal growth was assessed by comparison with control cultures in media without any sample tested. Standard drugs: ciprofloxacin (antibacterial chemotherapeutic) and nystatin (antifungal antibiotic) (Sigma-Aldrich Chemicals, St. Louis, MO, USA) were used as reference substances. Appropriate DMSO, sterile, and growth controls were prepared. The media with and without tested extracts/DMSO were used as controls [13,62,63,64,65,66,67].
Subsequently, minimal bactericidal concentration (MBC) or minimal fungicidal concentration (MFC) values of extracts were determined by transferring the cultures from each MIC determination well to the appropriate solid medium. After incubation, the lowest concentrations of extracts with no visible bacterial or fungal growth were evaluated as MBC or MFC. All the experiments were repeated three times as independent assays, and representative data are presented [13,62,63,64,65,66]. The MBC/MIC or MFC/MIC ratios were calculated in order to determine bactericidal/fungicidal (MBC/MIC ≤ 4, MFC/MIC ≤ 4) or bacteriostatic/fungistatic (MBC/MIC > 4, MFC/MIC > 4) effect of the tested extracts [65].

3.2.4. Total Antimicrobial Activity (TAA) Assay

The total antibacterial activity or total antifungal activity (TAA) of each of Z1, Z2, and Z3 extracts from Z. bungeana was obtained by dividing the quantity extracted from one gram of each plant extract by the MIC value. TAA was calculated following a standard formula:
T A A   m L / g = M a s s   o f   e x t r a c t   f r o m   1   g r a m   o f   p o w d e r   m g   p e r   g r a m           M I C     m g   p e r   m L
The total antibacterial activity (TAA) is a function of the extraction yield in milligram per 1 g of plant material and the minimal inhibitory concentration (MIC), expressed in milliliter per gram (mL/g). TAA indicates the volume of water or solvent, when added to 1 g of the extract, that will still inhibit the growth of the pathogen [68,69,70].

3.2.5. Antioxidant Activity

DPPH Scavenging Assay

The DPPH radical scavenging assay was performed as described by Matejic et al. [71]. Briefly, 100 µL of Z1, Z2 or Z3 diluted extracts (0.48–1000 µg/mL) was mixed with equal volume DPPH working solution (25 mM DPPH in 99.9% methanol; A540 ≈ 1). 100 µL of the solvent mixed with 100 µL DPPH served as a control sample. After 20 min incubation at RT in darkness, the absorbance of the samples was measured at λ = 540 nm using a FilterMax F5 microplate reader (Molecular Devices, San Jose, CA, USA). Obtained values of measurements were corrected by the absorbance values of the samples without DPPH. The percentage of DPPH radical scavenging was calculated based on the equation:
% of DPPH˙ scavenging = [1 − (Abs(S)/Abs(C))] × 100%
where Abs(S) is the absorbance of the sample and Abs(C) is the absorbance of the control sample (DPPH + solvent).
Obtained results were used to calculated EC50 values defined as the concentration of dried extract/fraction that is required to scavenge 50% of the DPPH radical activity.

ABTS Scavenging Assay

ABTS radical scavenging assay was performed according to Re and co-workers [72] with some modifications. Briefly, 135 µL of ABTS working solution (7 mM ABTS in 2.45 mM K2S2O8 diluted in distilled H2O up to A405 ≈ 1) was mixed with 15 µL of Z1, Z2 or Z3 diluted extract (0.48–1000 µg/mL) or solvent control. Following 15 min incubation at RT in darkness, the absorbance of the samples was measured at λ = 405 nm using a microplate reader (FilterMax F5 Molecular Devices, USA). The obtained values were corrected by the absorbance value of the sample without ABTS. The percentage of ABTS radical scavenging was calculated based on the equation:
% of ABTS scavenging = [1 − (Abs(S)/Abs(C))] × 100%
where Abs(S) is the absorbance of the extract and Abs(C) is the absorbance of the control sample (ABTS + solvent).
Obtained results were used to calculated EC50 values defined as the concentration of dried extract/fraction that is required to scavenge 50% of the ABTS radical activity.

SOD Inhibitory Assay

The influence of Z1, Z2, and Z3 extracts on the activity of superoxide dismutase (SOD) was measured using SOD Determination Kit (cat. No. 19160, Sigma Aldrich, Merck, Darmstadt, Germany), according to manufacturer’s instructions.

3.2.6. Tyrosinase Inhibitory Activity

The tyrosinase inhibitory activity of Z1, Z2, Z3 extracts was compared using commercially available mushroom tyrosinase (Sigma Aldrich) and murine tyrosinase contained in the lysate of B16F10 murine melanoma cells (ATCC CRL-6475; LGC Standards, Łomianki, Poland), prepared as previously described [73].
Mushroom tyrosinase activity assay was performed according to the protocol described by Uchida and co-workers [74] For this analysis, 120 μL phosphate buffer (100 mM, pH = 6.8) was mixed with 20 μL of diluted extracts (final concentrations 10–200 µg/mL) and 20 μL of mushroom tyrosinase (500 U/mL) and pre-incubated at room temperature for 10 min. Following the addition of 40 μL 4 mM L-DOPA, the samples were incubated for another 20 min at RT.
The activity of murine tyrosinase was assessed by Incubating the volume of B16F10 cell lysate containing 20 µg protein with 20 µL of diluted extracts (final concentrations 10–200 µg/mL), 40 µL 4 mM l-DOPA and 100 mM phosphate buffer pH 6.8 (up to 200 µL). The reaction was carried out for 4 h at 37 °C. Control samples (100% tyrosinase activity) for both assays contained an appropriate volume of the solvent instead of the extract. In both assays, the dopachrome formation was measured spectrophotometrically at λ = 450 nm using FilterMax F5 microplate reader (Molecular Devices, USA). The obtained values were corrected by the absorbance value of the extracts without mushroom or murine tyrosinase and l-DOPA. Each sample was analyzed in 3 independent repetitions. Kojic acid was used as a known tyrosinase inhibitor control.

3.2.7. In Vitro Cytotoxicity Assay

The cytotoxicity of Z1, Z2 and Z3 extracts was established by Neutral Red Uptake Test, as described by Repetto et al. [75] using human immortalized keratinocytes HaCaT (CLS Cell Lines Service GmbH, Eppelheim, Germany) [76], murine melanoma B16F10 (ATCC CRL-6475) and human melanoma A375 (ATCC CRL-1619) and SK-MEL3 (ATCC HTB-69) (LGC Standards, Łomianki, Poland). All cell lines were maintained in Dulbecco’s modified Eagle’s medium (DMEM)/high glucose supplemented with 10% fetal bovine serum (FBS, Pan Biotech, Aidenbach, Germany) at 37 °C in a humidified atmosphere with 5% CO2. For the experiments 3000 cells were plated per well onto a 96-well plate and grown overnight. Then, the cells were treated with various concentrations of Z1, Z2, or Z3 extracts (12.5–200 µg/mL) or an equal volume of the solvent control. Following 48 h of culture, the cells were incubated for 3 h in DMEM containing 1% FBS and 33 µg/mL neutral red, following by washing in PBS and lysis using acidified ethanol solution (50% v/v ethanol, 1% v/v acetic acid). The absorbance of the released neutral red was measured using FilterMax F5 microplate reader (Molecular Devices, San Jose, CA, USA) at λ= 540 nm. The mean measurement value for the lysate from control cells was set as 100% cellular viability and used to calculate the percentage of viable cells following extracts treatment.

3.2.8. Toxicity to Erythrocyte Assay

The erythrocyte lysis assay (ELA) was performed to study the toxicity of the extracts Z1, Z2 and Z3 from Ziziphora bungeana on red blood cells. In the first, erythrocytes were harvested from 5.0 mL fresh sheep blood (BioMaxima S.A., Poland) by centrifugation for 10 min at 1000× g and washed with 0.85% NaCl. Subsequently, 2% erythrocyte suspension was prepared in sterile phosphate buffer saline and in a volume of 100 μL was added to each well of a 96-well microtiter plate. The serial dilutions of these extracts ranging from 0.01 to 20 mg/mL were performed. To estimate the relative hemolytic potential of Z1, Z2, and Z3, the appropriate controls, i.e., 100% erythrocyte lysis using 4% Triton X-100 (Pol-Aura, Różnowo, Poland) and 0% lysis in saline solution, were used. Plates with samples were incubated for 1 h at 37 °C, then centrifuged for 10 min at 1000× g to separate the unlysed erythrocytes, and subsequently, the supernatant was transferred to a new plate. The absorbance was measured spectrophotometrically at 450 nm. The ELA represents an advantageous bioassay, because the lytic response can be measured photometrically by the amount of released hemoglobin. The hemolysis percentage was calculated according to the equation: % hemolysis = [(A450 of tested extract treated sample-A450 of buffer treated sample)/(A450 of 4% Triton X-100 treated samples-A450 of buffer treated sample)] × 100 [60,77,78,79].

3.3. Chemometric Analysis

All the chemometric analyses and visualizations were performed using R v4.2.0 [80] programming language in RStudio [81] software with pracma [82], factoextra [83], matlib [84], and corrplot [85] packages installed. After the standard, formal decomposition of the covariance matrices was calculated for the Crel (relative compositions) and Arel (relative activities) autoscaled datasets, and two principal components (PCs) were considered relevant in both cases. After the selection of the relevant PCs, their vectors were rotated in space in order to maximize the values of correlation coefficients between the original variables and the two orthogonal factors using the VARIMAX algorithm. In every case, compound/activity test scores in the space of the resulting varivectors (dimensions) were calculated by multiplying the matrix of the autoscaled Crel/Arel dataset by the matrix of the original variables’ loadings in the space of the resulting varivectors.

4. Conclusions

The presented results show the significance of Ziziphora bungeana extracts in terms of their composition and bioactivity. Twenty-six secondary metabolites were identified in the prepared extracts from Z. bungeana in the HPLC-ESI-QTOF-MS/MS analysis, that belonged to flavonoids, phenolic acids, terpenes, and organic acids. The results of antimicrobial studies indicated that extracts Z1, Z2, and Z3 showed potential activity with bactericidal or fungicidal effects. Among reference microorganisms, Gram-positive bacteria strains Staphylococcus spp., Micrococcus luteus, followed by Bacillus spp. were the most susceptible to the tested extracts (3.347–20.202 mL/g) in comparison with Gram-negative bacteria and fungi. Spectrophotometric assays proved the strongest antiradical properties of Z3 (EC50 values of 15.00 ± 1.06 µg/mL and 13.21 ± 3.24 µg/mL for ABTS and DPPH assays, respectively) and a marked SOD stimulatory action (>90% SOD activity) for Z1. In the murine tyrosinase assay all Ziziphora extracts showed significant, dose-dependent whitening properties. The most potent inhibitor of murine tyrosinase was extract Z2, decreasing the activity of tyrosinase by 50% at 200 µg/mL which was comparable with the inhibitory activity of kojic acid. All extracts were slightly cytotoxic for melanoma cells. However, Z2 at the concentration of 200 µg/mL showed a comparable cytotoxicity towards HaCaT keratinocytes. Moreover, our data suggest that the extracts Z1 and Z3 are not toxic for HaCaT cell lines or for erythrocyte membranes at the tested concentrations, which gives hope for its potential internal and external administration. The chemometric analysis performed to deliver the connections between the composition and biological properties of the extracts confirmed a different identity of all three extracts. According to the obtained results, the presence of the ziziphoroside isomer could induce anti-tyrosinase properties to the highest extent, whereas the presence of a higher quantity of acacetin could increase the anticancer potential of an extract.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27248994/s1. Figure S1: Fingerprints of the analyzed extracts in the negative ionization mode; Figure S2. Fingerprints of the analyzed extracts in the positive ionization modes; Table S1. MS/MS spectra of the identified compounds. Table S2: The relative compositions of the Z1–Z3 extracts. For the each studied compound (126), the HPLC peak areas representing the absolute amounts of a given compound within Z1, Z2 and Z3 extracts were rescaled, in order to sum to 1 in every row of the table below, expressing the relative compositions of the Z1–Z3 extracts; Table S3. The relative biological activities of the Z1–Z3 extracts. Biological activities represented by ‘survivability’ of cells/enzymes (S) were taken from Figure 6 for the Z1–Z3 extracts of concentrations equal to 200 ug/mL and used to calculate biological activity (A) using formula A = 100% − S. For each biological activity test (I–VIII), the values of A for the Z1, Z2 and Z3 extracts were rescaled in order to sum to 1 in every row of the table below.

Author Contributions

Conceptualization, K.Z., W.K.-K., K.G.-B., A.B., K.G., T.L. and Z.S.; methodology, K.Z., A.K. (Aigerim Karaubayeva)., A.K. (Aliya Kusniyeva), L.I., T.L., W.K.-K., K.G.-B., K.C., A.B. and A.M.; formal analysis, W.K.-K., A.B., K.G.-B., K.C., A.K. (Aigerim Karaubayeva), K.Z., T.B., K.T. and A.J.; investigation, K.Z., Z.S., A.K. (Aigerim Karaubayeva), W.K.-K., A.O., T.L., A.B. and K.G.-B.; resources, K.Z., L.I., Z.S., K.G.-B., K.G., W.K.-K., A.B., A.M. and G.I.; writing—original draft preparation, all authors; writing—review and editing, K.Z., W.K.-K., A.B. and K.G.-B.; visualization, K.C., W.K.-K., A.B., K.G.-B., T.L. and K.T.; supervision, K.Z.; project administration, K.Z.; funding acquisition, K.Z., K.G.-B., K.G., Z.S., W.K.-K. and A.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant IRN № AP09259196 from the Science Committee of the Ministry of Education and Science of the Republic of Kazakhstan, by the statutory project of the University of Information Technology and Management in Rzeszow, Poland (DS 503-07-01-38), and by the statutory project of the Medical University of Lublin (DS 24).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The obtained data are presented in the manuscript and supplementary File.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Sample Availability

Samples of the plant are available from the authors.

References

  1. Shahbazi, Y. Chemical compositions, antioxidant and antimicrobial properties of Ziziphora clinopodioides Lam. essential oils collected from different parts of Iran. J. Food Sci. Technol. 2017, 54, 3491–3503. [Google Scholar] [CrossRef] [PubMed]
  2. Hazrati, S.; Govahi, M.; Sedaghat, M.; Beyraghdar Kashkooli, A. A comparative study of essential oil profile, antibacterial and antioxidant activities of two cultivated Ziziphora species (Z. clinopodioides and Z. tenuior). Ind. Crops Prod. 2020, 157, 112942. [Google Scholar] [CrossRef]
  3. Eroglu Ozkan, E.; Boga, M.; Yilmaz, M.; Mataraci, E.; Yeşil, Y. LC-MS/MS analyses of Ziziphora clinopodioides Lam. from Turkey: Antioxidant, anticholinesterase, antimicrobial and, anticancer activities. İstanbul J. Pharm. 2020, 50, 33–41. [Google Scholar] [CrossRef]
  4. Ding, W.; Yang, T.; Liu, F.; Tian, S. Effect of different growth stages of Ziziphora clinopodioides Lam. on its chemical composition. Pharm. Mag. 2014, 10, S1–S5. [Google Scholar] [CrossRef] [Green Version]
  5. Hulya, D. Compositional diversity in essential oil of Ziziphora tenuior L. ecotypes. Genetika 2021, 53, 717–727. [Google Scholar] [CrossRef]
  6. Masrournia, M. Elemental Determination and Essential Oil Composition of Ziziphora clinopodioides and Consideration of its Antibacterial Effects. Asian J. Chem. 2013, 25, 6553. [Google Scholar] [CrossRef]
  7. Ozturk, Y.; Aydm, S.; Tecik, B.; Baser, K.H.C. Effects of essential oils from certain Ziziphora species on swimming performance in mice. Phytother. Res. 1995, 9, 225–227. [Google Scholar] [CrossRef]
  8. Tarakci, Z.; Coşkun, H.; Tuncturk, Y. Some Properties of Fresh and Ripened Herby Cheese, a Traditional Variety Produced in Turkey. Food Technol. Biotechnol 2004, 42, 47–50. [Google Scholar]
  9. Ozturk, S. Antibacterial activity and chemical constitutions of Ziziphora Clinopodioides. Food Control. 2007, 18, 535–540. [Google Scholar] [CrossRef]
  10. Šmejkal, K.; Malaník, M.; Zhaparkulova, K.; Sakipova, Z.; Ibragimova, L.; Ibadullaeva, G.; Žemlička, M. Kazakh Ziziphora Species as Sources of Bioactive Substances. Molecules 2016, 21, 826. [Google Scholar] [CrossRef] [Green Version]
  11. Hosseinzadeh, M.H.; Ebrahimzadeh, M.A. Antioxidant Potential of Ziziphora Clinopodioides Lam: A Narrative Review. mazu tbsrj 2020, 2, 1–7. [Google Scholar] [CrossRef]
  12. Srivedavyasasri, R.; Zhaparkulova, K.A.; Sakipova, Z.B.; Ibragimova, L.; Ross, S.A. Phytochemical and Biological Studies on Ziziphora bungeana. Chem. Nat. Compd. 2018, 54, 195–197. [Google Scholar] [CrossRef] [PubMed]
  13. Abdykerimova, S.; Sakipova, Z.; Nakonieczna, S.; Koch, W.; Biernasiuk, A.; Grabarska, A.; Malm, A.; Kozhanova, K.; Kukula-Koch, W. Superior Antioxidant Capacity of Berberis iliensis—HPLC-Q-TOF-MS Based Phytochemical Studies and Spectrophotometric Determinations. Antioxidants 2020, 9, 504. [Google Scholar] [CrossRef] [PubMed]
  14. Ahmadi, A.; Gandomi, H.; Derakhshandeh, A.; Misaghi, A.; Noori, N. Phytochemical composition and in vitro safety evaluation of Ziziphora clinopodioides Lam. ethanolic extract: Cytotoxicity, genotoxicity and mutagenicity assessment. J. Ethnopharmacol. 2021, 266, 113428. [Google Scholar] [CrossRef]
  15. Awwad, A.; Poucheret, P.; Idres, Y.; Tshibangu, D.; Servent, A.; Karine, F.; Lazennec, F.; Bidel, L.; Cazals, G.; Tousch, D. In Vitro Tests for a Rapid Evaluation of Antidiabetic Potential of Plant Species Containing Caffeic Acid Derivatives: A Validation by Two Well-Known Antidiabetic Plants, Ocimum gratissimum L. Leaf and Musanga cecropioides R. Br. ex Tedlie (Mu) Stem Bark. Mol. Cells 2021, 26, 5566. [Google Scholar] [CrossRef]
  16. Liu, C.; Wahefu, A.; Lu, X.; Abdulla, R.; Dou, J.; Zhao, H.; Aisa, H.A.; Xin, X.; Liu, Y. Chemical Profiling of Kaliziri Injection and Quantification of Six Caffeoyl Quinic Acids in Beagle Plasma by LC-MS/MS. Pharmaceuticals 2022, 15, 663. [Google Scholar] [CrossRef]
  17. de Oliveira, J.K.; Ronik, D.F.V.; Ascari, J.; Mainardes, R.M.; Khalil, N.M. A stability-indicating high performance liquid chromatography method to determine apocynin in nanoparticles. J. Pharm. Anal. 2017, 7, 129–133. [Google Scholar] [CrossRef]
  18. Zengin, G.; Mahomoodally, F.; Ibrahime, S.; Ak, G.; Etienne, O.; Jugreet, S.; Brunetti, L.; Leone, S.; Cristina, S.; Di Simone, S.; et al. Chemical Composition and Biological Properties of Two Jatropha Species: Different Parts and Different Extraction Methods. Antioxidants 2021, 10, 792. [Google Scholar] [CrossRef]
  19. Wang, X.; Qian, Y.; Li, X.; Jia, X.; Yan, Z.; Han, M.; Qiao, M.; Ma, X.; Chu, Y.; Zhou, S.; et al. Rapid determination of rosmarinic acid and its two bioactive metabolites in the plasma of rats by LC-MS/MS and application to a pharmacokinetics study. Biomed. Chromatogr 2021, 35, e4984. [Google Scholar] [CrossRef]
  20. Micucci, M.; Protti, M.; Aldini, R.; Frosini, M.; Corazza, I.; Marzetti, C.; Mattioli, L.; Tocci, G.; Chiarini, A.; Mercolini, L.; et al. Thymus vulgaris L. Essential Oil Solid Formulation: Chemical Profile and Spasmolytic and Antimicrobial Effects. Biomolecules 2020, 10, 860. [Google Scholar] [CrossRef]
  21. Huwait, E.; Mobashir, M. Potential and Therapeutic Roles of Diosmin in Human Diseases. Biomedicines 2022, 10, 1076. [Google Scholar] [CrossRef] [PubMed]
  22. Li, Y.; Guang, C.; Zhao, N.; Feng, X.; Qiu, F. LC-MS/MS Method for Simultaneous Determination of Linarin and Its Metabolites in Rat Plasma and Liver Tissue Samples: Application to Pharmacokinetic and Liver Tissue Distribution Study After Oral Administration of Linarin. Molecules 2019, 24, 3342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Wianowska, D.; Dawidowicz, A.L.; Bernacik, K.; Typek, R. Determining the true content of quercetin and its derivatives in plants employing SSDM and LC–MS analysis. Eur. Food Res. Technol. 2016, 243, 27–40. [Google Scholar] [CrossRef] [Green Version]
  24. Shen, J.; Jia, Q.; Huang, X.; Yao, G.; Ma, W.; Zhang, H.; Ouyang, H.; He, J. Development of a HPLC-MS/MS Method to Determine the 13 Elements of Semen cuscutae and Application to a Pharmacokinetic Study in Rats. Evid. Based Complement. Altern. Med. 2019, 2019, 6518528. [Google Scholar] [CrossRef] [Green Version]
  25. Cheruvu, H.S.; Yadav, N.K.; Valicherla, G.R.; Arya, R.K.; Hussain, Z.; Sharma, C.; Arya, K.R.; Singh, R.K.; Datta, D.; Gayen, J.R. LC-MS/MS method for the simultaneous quantification of luteolin, wedelolactone and apigenin in mice plasma using hansen solubility parameters for liquid-liquid extraction: Application to pharmacokinetics of Eclipta alba chloroform fraction. J. Chromatogr. B 2018, 1081–1082, 76–86. [Google Scholar] [CrossRef]
  26. Kim, S.; Kim, J.; Kim, N.; Lee, D.; Lee, H.; Lee, D.Y.; Kim, K.H. Metabolomic Elucidation of the Effect of Sucrose on the Secondary Metabolite Profiles in Melissa officinalis by Ultraperformance Liquid Chromatography-Mass Spectrometry. ACS Omega 2020, 5, 33186–33195. [Google Scholar] [CrossRef]
  27. Silvestro, L.; Tarcomnicu, I.; Dulea, C.; Attili, N.R.; Ciuca, V.; Peru, D.; Rizea Savu, S. Confirmation of diosmetin 3-O-glucuronide as major metabolite of diosmin in humans, using micro-liquid-chromatography-mass spectrometry and ion mobility mass spectrometry. Anal. Bioanal Chem. 2013, 405, 8295–8310. [Google Scholar] [CrossRef] [Green Version]
  28. Knez Hrnčič, M.; Cör, D.; Simonovska, J.; Knez, Ž.; Kavrakovski, Z.; Rafajlovska, V. Extraction Techniques and Analytical Methods for Characterization of Active Compounds in Origanum Species. Molecules 2020, 25, 4735. [Google Scholar] [CrossRef]
  29. Stebounova, L.; Ebert, S.M.; Murry, L.T.; Adams, C.M.; Murry, D.J. Rapid and Sensitive Quantification of Ursolic Acid and Oleanolic Acid in Human Plasma Using Ultra-performance Liquid Chromatography–Mass Spectrometry. J. Chromatogr. Sci. 2018, 56, 644–649. [Google Scholar] [CrossRef] [Green Version]
  30. Castellano, J.M.; Ramos-Romero, S.; Perona, J.S. Oleanolic Acid: Extraction, Characterization and Biological Activity. Nutrients 2022, 14, 623. [Google Scholar] [CrossRef]
  31. Sánchez-González, M.; Lozano-Mena, G.; Juan, M.E.; García-Granados, A.; Planas, J.M. Liquid chromatography-mass spectrometry determination in plasma of maslinic acid, a bioactive compound from Olea europaea L. Food Chem 2013, 141, 4375–4381. [Google Scholar] [CrossRef] [PubMed]
  32. Sarikurkcu, C.; Kakouri, E.; Sarikurkcu, R.T.; Tarantilis, P.A. Study on the Chemical Composition, Enzyme Inhibition and Antioxidant Activity of Ziziphora taurica subsp. cleonioides. Appl. Sci. 2019, 9, 5515. [Google Scholar] [CrossRef] [Green Version]
  33. Zhu, Y.; Xiong, Y.; Wang, H.; Li, P. Pharmacognostical and phytochemical studies on Ziziphora clinopodioides Lam.—A Kazakh and Uygur ethnomedicinal plant. J. Pharm. Pharmacogn. Res. 2017, 5, 345–353. [Google Scholar]
  34. Koch, W. Dietary Polyphenols—Important Non-Nutrients in the Prevention of Chronic Noncommunicable Diseases. A Systematic Review. Nutrients 2019, 11, 1039. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Koch, W.; Baj, T.; Kukula-Koch, W.; Marzec, Z. Dietary intake of specific phenolic compounds and their effect on the antioxidant activity of daily food rations. Open Chem. 2015, 13, 869–876. [Google Scholar] [CrossRef]
  36. Eloff, J.N. Quantification the bioactivity of plant extracts during screening and bioassay guided fractionation. Phytomedicine 2004, 11, 370–371. [Google Scholar] [CrossRef]
  37. Henley-Smith, C.J.; Steffens, F.E.; Botha, F.S.; Lall, N. Predicting the influence of multiple components on microbial inhibition using a logistic response model—A novel approach. BMC Complement. Altern. Med. 2014, 14, 190. [Google Scholar] [CrossRef] [Green Version]
  38. Sonboli, A.; Mirjalili, M.H.; Hadian, J.; Ebrahimi, S.N.; Yousefzadi, M. Antibacterial activity and composition of the essential oil of Ziziphora clinopodioides subsp. bungeana (Juz.) Rech. f. from Iran. Z. Nat. C J. Biosci. 2006, 61, 677–680. [Google Scholar] [CrossRef]
  39. Ozturk, S.; Ercisli, S. The chemical composition of essential oil and in vitro antibacterial activities of essential oil and methanol extract of Ziziphora persica Bunge. J. Ethnopharmacol. 2006, 106, 372–376. [Google Scholar] [CrossRef]
  40. Anzabi, Y.; Khaki, A.; Rasoli, A.; Ebrahimpour, S.; Fallah Rostami, F. Antibacterial properties of essential oils and methanol extracts of Ziziphora tenuior Lam. (a native plant) in pre-flowering stage against isolated bacteria from urogenital tract infections. Bulg. Chem. Commun. 2016, 48, 120–125. [Google Scholar]
  41. Gursoy, N.; Sihoglu-Tepe, A.; Tepe, B. Determination of in vitro antioxidative and antimicrobial properties and total phenolic contents of Ziziphora clinopodioides, Cyclotrichium niveum, and Mentha longifolia ssp. typhoides var. typhoides. J. Med. Food 2009, 12, 684–689. [Google Scholar] [CrossRef] [PubMed]
  42. Salehi, P.; Sonboli, A.; Eftekhar, F.; Nejad-Ebrahimi, S.; Yousefzadi, M. Essential oil composition, antibacterial and antioxidant activity of the oil and various extracts of Ziziphora clinopodioides subsp. rigida (BOISS.) RECH. f. from Iran. Biol. Pharm. Bull. 2005, 28, 1892–1896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Celik, C.; Tutar, U.; Karaman, I.; Hepokur, C.; Atas, M. Evaluation of the Antibiofilm and Antimicrobial Properties of Ziziphora tenuior L. Essential Oil Against Multidrug-resistant Acinetobacter baumannii. Int. J. Pharmacol. 2016, 12, 28–35. [Google Scholar] [CrossRef]
  44. Ydyrys, A.; Zhaparkulova, N.; Aralbaeva, A.; Mamataeva, A.; Seilkhan, A.; Syraiyl, S.; Murzakhmetova, M. Systematic Analysis of Combined Antioxidant and Membrane-Stabilizing Properties of Several Lamiaceae Family Kazakhstani Plants for Potential Production of Tea Beverages. Plants 2021, 10, 666. [Google Scholar] [CrossRef] [PubMed]
  45. Tomczyk, M.; Ceylan, O.; Locatelli, M.; Tartaglia, A.; Ferrone, V.; Sarikurkcu, C. Ziziphora taurica subsp. taurica: Analytical Characterization and Biological Activities. Biomolecules 2019, 9, 367. [Google Scholar] [CrossRef] [Green Version]
  46. Zolghadri, S.; Bahrami, A.; Hassan Khan, M.T.; Munoz-Munoz, J.; Garcia-Molina, F.; Garcia-Canovas, F.; Saboury, A.A. A comprehensive review on tyrosinase inhibitors. J. Enzym. Inhib. Med. Chem. 2019, 34, 279–309. [Google Scholar] [CrossRef] [Green Version]
  47. Niu, C.; Aisa, H.A. Upregulation of Melanogenesis and Tyrosinase Activity: Potential Agents for Vitiligo. Molecules 2017, 22, 1303. [Google Scholar] [CrossRef] [Green Version]
  48. Pillaiyar, T.; Manickam, M.; Namasivayam, V. Skin whitening agents: Medicinal chemistry perspective of tyrosinase inhibitors. J. Enzym. Inhib. Med. Chem. 2017, 32, 403–425. [Google Scholar] [CrossRef] [Green Version]
  49. Yoshimori, A.; Oyama, T.; Takahashi, S.; Abe, H.; Kamiya, T.; Abe, T.; Tanuma, S. Structure-activity relationships of the thujaplicins for inhibition of human tyrosinase. Bioorg. Med. Chem. 2014, 22, 6193–6200. [Google Scholar] [CrossRef]
  50. Strzępek-Gomółka, M.; Gaweł-Bęben, K.; Kukula-Koch, W. Achillea Species as Sources of Active Phytochemicals for Dermatological and Cosmetic Applications. Oxid. Med. Cell Longev. 2021, 2021, 6643827. [Google Scholar] [CrossRef]
  51. Zhu, H.; Ghoufi, A.; Szymczyk, A.; Balannec, B.; Morineau, D. Zhu et al. Reply. Phys. Rev. Lett. 2013, 111, 089802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Pakniyat, A.; Caines, P.E. Time Optimal Hybrid Minimum Principle and the Gear Changing Problem for Electric Vehicles. IFAC Pap. 2015, 48, 187–192. [Google Scholar] [CrossRef]
  53. Nath, R.; Beyer, D.; Friedland, J.; Grimm, P.; Nag, S. Regarding, Kubo, Coursey, Hanson et al., IJROBP 1998; 40:697–702. Int. J. Radiat. Oncol. Biol. Phys. 1999, 44, 469–470. [Google Scholar] [PubMed]
  54. Salehi, T.; Mahzounieh, M.; Saeedzadeh, A. Detection of InvA Gene in Isolated Salmonella from Broilers by PCR Method. Int. J. Poult. Sci. 2005, 4, 557–559. [Google Scholar] [CrossRef] [Green Version]
  55. Nitoda, T.; Fan, M.D.; Kubo, I. Effects of cuminaldehyde on melanoma cells. Phytother. Res. 2008, 22, 809–813. [Google Scholar] [CrossRef]
  56. Azimi, M.; Mehrzad, J.; Ahmadi, A.; Ahmadi, E.; Ghorbani Ranjbary, A. Apoptosis Induced by Ziziphora tenuior Essential Oil in Human Colorectal Cancer Cells. BioMed. Res. Int. 2021, 2021, 5522964. [Google Scholar] [CrossRef]
  57. Taghizadeh, M.S.; Niazi, A.; Moghadam, A.; Afsharifar, A.R. The potential application of the protein hydrolysates of three medicinal plants: Cytotoxicity and functional properties. J. Food Sci. 2020, 85, 3160–3167. [Google Scholar] [CrossRef]
  58. Ghavam, M. A GC-MC analysis of chemical compounds and identification of the antibacterial characteristics of the essential oil of two species exclusive to Iranian habitats: New chemotypes. PLoS ONE 2022, 17, e0273987. [Google Scholar] [CrossRef]
  59. Yousefbeyk, F.; Ostad, S.; Sourmaghi, M.; Amin, G. Investigation of chemical composition and cytotoxic activity of aerial parts of Ziziphora clinopodioides Lam. Res. J. Pharmacogn. 2016, 3, 41–57. [Google Scholar]
  60. Zohra, M.; Fawzia, A. Hemolytic activity of different herbal extracts used in Algeria. Int. J. Pharm. Sci. Res. 2014, 5, 495–500. [Google Scholar]
  61. Luna-Vázquez-Gómez, R.; Arellano-García, M.E.; García-Ramos, J.C.; Radilla-Chávez, P.; Salas-Vargas, D.S.; Casillas-Figueroa, F.; Ruiz-Ruiz, B.; Bogdanchikova, N.; Pestryakov, A. Hemolysis of Human Erythrocytes by Argovit™ AgNPs from Healthy and Diabetic Donors: An In Vitro Study. Materials 2021, 14, 2792. [Google Scholar] [CrossRef] [PubMed]
  62. European Committee for Antimicrobial Susceptibility Testing (EUCAST) of the European Society of Clinical Microbiology and Infectious Diseases (ESCMID). EUCAST Definitive Document E.DEF 3.1, June 2000: Determination of minimum inhibitory concentrations (MICs) of antibacterial agents by agar dilution. Clin. Microbiol. Infect. 2000, 6, 509–515. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Wiegand, I.; Hilpert, K.; Hancock, R.E. Agar and broth dilution methods to determine the minimal inhibitory concentration (MIC) of antimicrobial substances. Nat. Protoc. 2008, 3, 163–175. [Google Scholar] [CrossRef] [PubMed]
  64. O’Donnell, F.; Smyth, T.J.; Ramachandran, V.N.; Smyth, W.F. A study of the antimicrobial activity of selected synthetic and naturally occurring quinolines. Int. J. Antimicrob. Agents 2010, 35, 30–38. [Google Scholar] [CrossRef] [Green Version]
  65. Wald-Dickler, N.; Holtom, P.; Spellberg, B. Busting the Myth of “Static vs. Cidal”: A Systemic Literature Review. Clin. Infect. Dis. 2018, 66, 1470–1474. [Google Scholar] [CrossRef] [PubMed]
  66. Wayne, P. Reference method for broth dilution antifungal susceptibility testing of yeasts. Clin. Lab. Stand. Inst. 2008, 3, M27-A2. [Google Scholar]
  67. Abdrahimov, R. Annual river runoff of the ile-balkash basin and prospects of its assessment due to climatic changes and water economy activities. Int. J. Geomate 2020, 18, 230–239. [Google Scholar] [CrossRef]
  68. Eloff, J.; McGaw, L. Using African Plant Biodiversity to Combat Microbial Infections. In Novel Plant Bioresources: Applications in Food, Medicine and Cosmetics; Wiley: Hoboken, NJ, USA, 2014; pp. 163–173. [Google Scholar] [CrossRef]
  69. Elisha, I.L.; Botha, F.S.; McGaw, L.J.; Eloff, J.N. The antibacterial activity of extracts of nine plant species with good activity against Escherichia coli against five other bacteria and cytotoxicity of extracts. BMC Complement. Altern. Med. 2017, 17, 133. [Google Scholar] [CrossRef] [Green Version]
  70. Famuyide, I.M.; Aro, A.O.; Fasina, F.O.; Eloff, J.N.; McGaw, L.J. Antibacterial and antibiofilm activity of acetone leaf extracts of nine under-investigated south African Eugenia and Syzygium (Myrtaceae) species and their selectivity indices. BMC Complement. Altern. Med. 2019, 19, 141. [Google Scholar] [CrossRef] [Green Version]
  71. Matejić, J.; Dzamic, A.; Mihajilov-Krstev, T.; Randjelovic, V.; Krivošej, Z.D.; Marin, P. Total phenolic and flavonoid content, antioxidant and antimicrobial activity of extracts from Tordylium maximum. J. Appl. Pharm. Sci. 2013, 3, 55–59. [Google Scholar]
  72. 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] [PubMed]
  73. Strzępek-Gomółka, M.; Gaweł-Bęben, K.; Angelis, A.; Antosiewicz, B.; Sakipova, Z.; Kozhanova, K.; Głowniak, K.; Kukula-Koch, W. Identification of Mushroom and Murine Tyrosinase Inhibitors from Achillea biebersteinii Afan. Extract. Molecules 2021, 26, 964. [Google Scholar] [CrossRef] [PubMed]
  74. Uchida, R.; Ishikawa, S.; Tomoda, H. Inhibition of tyrosinase activity and melanine pigmentation by 2-hydroxytyrosol. Acta Pharm. Sin. B 2014, 4, 141–145. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Repetto, G.; del Peso, A.; Zurita, J.L. Neutral red uptake assay for the estimation of cell viability/cytotoxicity. Nat. Protoc. 2008, 3, 1125–1131. [Google Scholar] [CrossRef] [PubMed]
  76. Boukamp, P.; Petrussevska, R.T.; Breitkreutz, D.; Hornung, J.; Markham, A.; Fusenig, N.E. Normal keratinization in a spontaneously immortalized aneuploid human keratinocyte cell line. J. Cell Biol. 1988, 106, 761–771. [Google Scholar] [CrossRef] [Green Version]
  77. Biernasiuk, A.; Berecka-Rycerz, A.; Gumieniczek, A.; Malm, M.; Łączkowski, K.Z.; Szymańska, J.; Malm, A. The newly synthesized thiazole derivatives as potential antifungal compounds against Candida albicans. Appl. Microbiol. Biotechnol. 2021, 105, 6355–6367. [Google Scholar] [CrossRef]
  78. Turecka, K.; Chylewska, A.; Kawiak, A.; Waleron, K.F. Antifungal Activity and Mechanism of Action of the Co(III) Coordination Complexes With Diamine Chelate Ligands Against Reference and Clinical Strains of Candida spp. Front. Microbiol. 2018, 9, 1594. [Google Scholar] [CrossRef] [Green Version]
  79. Silva, S.; Rodrigues, C.F.; Araújo, D.; Rodrigues, M.E.; Henriques, M. Candida Species Biofilms’ Antifungal Resistance. J. Fungi 2017, 3, 8. [Google Scholar] [CrossRef]
  80. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2020. [Google Scholar]
  81. RStudio Team. RStudio: Integrated Development Environment for R.; RStudio, Inc.: Boston, MA, USA, 2019. [Google Scholar]
  82. Borchers, H.W. Pracma: Practical Numerical Math Functions 2022, Version 2.3.8 from R-Forge; R Package Vignette: Madison, WI, USA, 2022. [Google Scholar]
  83. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses, R Package Version 1.0.7; R Foundation for Statistical Computing: Vienna, Austria, 2020; pp. 337–354. [Google Scholar]
  84. Friendly, M.; Fox, J.; Chalmers, P. Matlib: Matrix Functions for Teaching and Learning Linear Algebra and Multivariate Statistics; R Package Vignette: Madison, WI, USA, 2021. [Google Scholar]
  85. Wei, T.; Simko, V. R Package “Corrplot”: Visualization of a Correlation Matrix; R Package Vignette: Madison, WI, USA, 2021. [Google Scholar]
Figure 1. Structure of the components identified in Ziziphora bungeana extracts.
Figure 1. Structure of the components identified in Ziziphora bungeana extracts.
Molecules 27 08994 g001aMolecules 27 08994 g001b
Figure 2. The frequency of occurrence of particular MBC/MIC or MFC/MIC ratios of Z. bungeana extracts against the reference strains of bacteria and fungi used in the study.
Figure 2. The frequency of occurrence of particular MBC/MIC or MFC/MIC ratios of Z. bungeana extracts against the reference strains of bacteria and fungi used in the study.
Molecules 27 08994 g002
Figure 3. The efficacy (mean TAA values, mL/g) of Z. bungeana extracts against all the studied reference bacteria and fungi.
Figure 3. The efficacy (mean TAA values, mL/g) of Z. bungeana extracts against all the studied reference bacteria and fungi.
Molecules 27 08994 g003
Figure 4. Antioxidant activity of Z1, Z2 and Z3 extracts from Z. bungeana: (A) neutralization of DPPH and ABTS free radicals, displayed as EC50; (B) relative activity of superoxide dismutase (SOD) of Z1, Z2 and Z3 extracts; histograms show mean values ±SD, n = 3.
Figure 4. Antioxidant activity of Z1, Z2 and Z3 extracts from Z. bungeana: (A) neutralization of DPPH and ABTS free radicals, displayed as EC50; (B) relative activity of superoxide dismutase (SOD) of Z1, Z2 and Z3 extracts; histograms show mean values ±SD, n = 3.
Molecules 27 08994 g004
Figure 5. The influence of Z1, Z2 and Z3 extracts from Z. bugeana on the activity of murine (A) and mushroom (B) tyrosinase; histograms show mean tyrosinase activity ±SD, * p < 0.05, ** p < 0.01, *** p < 0.001; KA—kojic acid.
Figure 5. The influence of Z1, Z2 and Z3 extracts from Z. bugeana on the activity of murine (A) and mushroom (B) tyrosinase; histograms show mean tyrosinase activity ±SD, * p < 0.05, ** p < 0.01, *** p < 0.001; KA—kojic acid.
Molecules 27 08994 g005
Figure 6. In vitro cytotoxicity of Z1, Z2 and Z3 Z. bungeana extracts on human keratinocytes HaCaT (A), murine melanoma B16F10 (B) and human melanoma cell lines A375 (C) and SK-MEL3 (D) following 48 h culture; graphs show mean viability of the cells ± SD in comparison with appropriate solvent controls; * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 6. In vitro cytotoxicity of Z1, Z2 and Z3 Z. bungeana extracts on human keratinocytes HaCaT (A), murine melanoma B16F10 (B) and human melanoma cell lines A375 (C) and SK-MEL3 (D) following 48 h culture; graphs show mean viability of the cells ± SD in comparison with appropriate solvent controls; * p < 0.05, ** p < 0.01, *** p < 0.001.
Molecules 27 08994 g006
Figure 7. Hemolytic effect (% of hemolysis) of the studied extracts from Z. bungeana.
Figure 7. Hemolytic effect (% of hemolysis) of the studied extracts from Z. bungeana.
Molecules 27 08994 g007
Figure 8. Relations of the Z1Z3 vectors in the space of the first two principal components (subjected to VARIMAX rotation), regarding their relative compositions (a) and relative biological activities (b).
Figure 8. Relations of the Z1Z3 vectors in the space of the first two principal components (subjected to VARIMAX rotation), regarding their relative compositions (a) and relative biological activities (b).
Molecules 27 08994 g008
Figure 9. Linear maps of the selected compounds in the space of the first two principal components (subjected to VARIMAX rotation), regarding their relative compositions (a) and relative biological activities (b).
Figure 9. Linear maps of the selected compounds in the space of the first two principal components (subjected to VARIMAX rotation), regarding their relative compositions (a) and relative biological activities (b).
Molecules 27 08994 g009
Table 1. The list of tentatively identified compounds present in the analyzed samples that were obtained from the HPLC-ESI-QTOF-MS/MS analysis in positive and negative ionization modes (DBE—double bond equivalent, error—error of measurement in ppm, Ion.—ionisation type, Rt—retention time).
Table 1. The list of tentatively identified compounds present in the analyzed samples that were obtained from the HPLC-ESI-QTOF-MS/MS analysis in positive and negative ionization modes (DBE—double bond equivalent, error—error of measurement in ppm, Ion.—ionisation type, Rt—retention time).
NoIon.+/−Rt
[min]
Molecular Formulam/z
Theoretical
m/z
Experimental
Error DBEMS/MS SpectrumProposed CompoundDistribution References
13.03C16H18O9353.0878353.0884−1.688191, 179, 173, 154Chlorogenic acidZ1, Z2, Z3[13]
23.9C16H18O9353.0878353.0880−0.558179, 191Neochlorogenic acidZ1, Z2, Z3[13]
3+3.9C16H28O8349.1857349.18511.713281, 163Schizonepetaside E Z1, Z2, Z3
44.4C16H18O9353.0878353.0884−1.688191, 179, 173, 135, 155(Z)-Chlorogenic acidZ1, Z2, Z3[13]
55.2C14H18O7297.098297.1011−10.486NDPiceinTraces Z3[14]
66.5C9H8O4179.0350179.0357−3.996135, 117, 107Caffeic acidZ1, Z2, Z3[15]
77.7C11H12O4207.0663207.0664−0.576192, 179, 174, 163, 135Ethyl ester of caffeic acidZ1, Z3[16]
89.5C9H10O3165.0557165.054010.345NDApocynin Traces: Z1, Z2, Z3[17]
9+10.9C16H26O7331.1751
353.1571 (+Na)
331.1775
353.1607 (+Na)
−7.18
−10.3
4201Ziziphoroside isomer 1Z1, Z2, Z3[14]
10+12.9C16H26O7331.1751
353.1571 (+Na)
331.1761
353.1607 (+Na)
−2.94
−10.3
4201, 147, 119Ziziphoroside isomer 2Z1, Z2, Z3[14]
11+14.3C16H26O7331.1751
353.1571 (+Na)
331.1757
353.1609 (+Na)
−1.73
−10.86
4201, 165, 147Ziziphoroside isomer 3Z1, Z2, Z3[14]
1218.5C21H18O11445.0776445.07672.113269, 175, 1135,7,21-trihydroxyflavone-21-
O-glucopyranoside
Z1, Z3[18]
1319.1C18H16O8359.0772 359.0780−2.1111197, 179, 161, 135Rosmarinic acidZ1, Z2, Z3[19]
14+19.8C10H14O151.1117151.1129−7.724133, 123, 109, 105Thymol Z1, Z2, Z3[20]
1520.1C28H32O15607.1668607.1671−0.4213561, 253Diosmin Z1, Z2, Z3[21]
1620.4C28H32O14591.1719591.1727−1.313NDLinarin Traces: Z1, Z2, Z3 [22]
1720.9C15H10O7301.0354301.0363−3.0611NDQuercetin Z1, Z2, Z3 (traces)[23]
1821.1C15H10O6285.0405285.03963.0111241, 151, 133Luteolin Z1, Z2, Z3[24]
1921.6 C15H10O5269.0455269.04622.0311225, 151Apigenin Z1, Z2, Z3[25]
20 21.8C18H16O8359.0772359.07720.1111344, 329Thymonin Z1, Z2, Z3[26]
2122.3C16H12O6299.0561299.05552.0411284, 256, 165, 135Diosmetin Z1, Z2, Z3[27]
22+22.6C10H14O151.1117151.1135−11.714136, 123, 117, 105Carvacrol Z1, Z2, Z3[28]
2322.7C16H12O5283.0612283.0620−2.8311268, 240Acacetin Z1, Z2, Z3[22]
2423.0C30H48O3455.3531455.3538−1.67455Oleanolic acid Z1, Z3[22,29]
2523.3C30H48O4471.3480471.34790.187337Maslinic acidZ1, Z2, Z3[30]
2624.0C30H48O3455.3531455.35280.597455Ursolic acidZ1, Z3[31]
Table 2. The activity data of Z. bungeana extracts expressed as MIC (Minimum Inhibitory Concentration), MBC (Minimum Bactericidal Concentration) [mg/mL] and MBC/MIC value against the reference strains of microogranisms. (CIP—ciprofloxacin (MIC and MBC) [µg/mL]).
Table 2. The activity data of Z. bungeana extracts expressed as MIC (Minimum Inhibitory Concentration), MBC (Minimum Bactericidal Concentration) [mg/mL] and MBC/MIC value against the reference strains of microogranisms. (CIP—ciprofloxacin (MIC and MBC) [µg/mL]).
Species of MicroorganismZ1Z2Z3CIP
MICMBCMBC
/MIC
MICMBCMBC
/MIC
MICMBCMBC
/MIC
MICMBCMBC
/MIC
Gram-positiveStaphylococcus aureus
ATCC 29213
2.52.515512.52.510.240.241
Staphylococcus aureus
ATCC 43300
2.52.512.5522.5520.240.241
Staphylococcus epidermidis
ATCC 12228
2.52.512.5522.5520.120.121
Micrococcus luteus
ATCC 10240
2.5521.255451020.981.962
Bacillus subtilis
ATCC 6633
510251022.51040.030.031
Bacillus cereus
ATCC 10876
10101510251020.060.122
Gram-negativeBordetella bronchiseptica
ATCC 4617
101015204202010.980.981
Klebsiella pneumoniae
ATCC 13883
2020110202202010.120.121
Salmonella typhimurium
ATCC 14028
2020120201202010.060.061
Escherichia coli
ATCC 25922
2020120201202010.0040.0041
Pseudomonas aeruginosa
ATCC 9027
1020220201202010.480.982
The representative (modal) data are presented. The sensitivity of the fungi belonging to Candida spp. to the tested extracts Z1–Z3 was similar to that of Gram-negative bacteria (MIC = 5–20 mg/mL and MFC = 20 mg/mL). Candida parapsilosis ATCC 22,019 was the most susceptible to Z2 and Z1 at MIC = 5 mg/mL and 10 mg/mL, respectively. Z2 showed also activity towards other Candida spp. with MIC = 10 mg/mL, except Candida glabrata ATCC 90,030 (MIC = 20 mg/mL). Moreover, the minimal concentrations of these extracts, which inhibited growth or killed these microorganisms were 20 mg/mL (Table 3).
Table 3. The activity data of Z. bungeana extracts expressed as MIC (Minimum Inhibitory Concentration), MFC (Minimum Fungicidal Concentration) [mg/mL] and MFC/MIC value against the reference strains of fungi (NYS —nystatin (MIC and MFC) [µg/mL]).
Table 3. The activity data of Z. bungeana extracts expressed as MIC (Minimum Inhibitory Concentration), MFC (Minimum Fungicidal Concentration) [mg/mL] and MFC/MIC value against the reference strains of fungi (NYS —nystatin (MIC and MFC) [µg/mL]).
Species of MicroorganismZ1Z2Z3NYS
MICMFCMFC
/MIC
MICMFCMFC
/MIC
MICMFCMFC
/MIC
MICMFCMFC
/MIC
Candida albicans
ATCC 10231
2020110202202010.480.481
Candida albicans
ATCC 2091
2020110202202010.240.241
Candida parapsilosis
ATCC 22019
102015204202010.240.482
Candida glabrata
ATCC 90030
2020120201202010.240.482
Candida krusei
ATCC 14243
2020110202202010.240.241
The representative (modal) data are presented. As shown our results (Table 2 and Table 3), the most common MIC value of 20 mg/mL was found for Z1 (7 (43.75%) strains) and Z3 extracts (10 (62.5%) strains). In the case of Z2 extract, the values of MIC = 20 mg/mL and MIC = 10 mg/mL, occurred with the same frequency (25% each) against reference strains of microorganisms. The same frequency of MIC = 10 mg/mL was shown for Z1. MIC values of 5 mg/mL were shown for 5 (31.25%), 2 (12.5%) and 1 (6.25%) strains in the case of Z2, Z3 and Z1 extracts, respectively. The Z1 and Z3 inhibited the growth of microorganisms at the minimum concentration of 2.5 mg/mL (4 strains (25%) each). Additionally, Z2 inhibited the growth of 1 (6.25%) and 2 (12.5%) strains with MIC = 1.25 mg/mL and 2.5 mg/mL, respectively.
Table 4. The activity data of Z. bungeana extracts expressed as TAA (Total Antibacterial Activity or Total Antifungal Activity) [mL/g] against the reference strains of bacteria and fungi.
Table 4. The activity data of Z. bungeana extracts expressed as TAA (Total Antibacterial Activity or Total Antifungal Activity) [mL/g] against the reference strains of bacteria and fungi.
Species of MicroorganismTAA (mL/g)
Z1Z2Z3
Gram-positive bacteriaStaphylococcus aureus
ATCC 29213
13.39 ± 0.06.74 ± 2.928.75 ± 0.0
Staphylococcus aureus
ATCC 43300
11.16 ± 3.868.42 ± 2.927.30 ± 2.53
Staphylococcus epidermidis
ATCC 12228
13.39 ± 0.010.10 ± 0.08.75 ± 0.0
Micrococcus luteus
ATCC 10240
13.39 ± 0.020.20 ± 0.05.84 ± 2.53
Bacillus subtilis
ATCC 6633
5.58 ± 1.935.05 ± 0.07.30 ± 2.53
Bacillus cereus
ATCC 10876
3.35 ± 0.04.21 ± 1.464.38 ± 0.0
Gram-negative bacteriaBordetella bronchiseptica
ATCC 4617
2.79 ± 0.974.21 ± 1.461.46 ± 0.63
Klebsiella pneumoniae
ATCC 13883
1.67 ± 0.02.53 ± 0.01.09 ± 0.0
Salmonella typhimurium
ATCC 14028
1.67 ± 0.01.68 ± 0.731.09 ± 0.0
Escherichia coli
ATCC 25922
2.23 ± 0.971.68 ± 0.731.46 ± 0.63
Pseudomonas aeruginosa
ATCC 9027
2.79 ± 0.971.26 ± 0.01.09 ± 0.0
FungiCandida albicans
ATCC 10231
2.23 ± 0.972.53 ± 0.01.46 ± 0.63
Candida albicans
ATCC 2091
1.67 ± 0.02.53 ± 0.01.09 ± 0.0
Candida parapsilosis
ATCC 22019
3.35 ± 0.04.21 ± 1.461.46 ± 0.63
Candida glabrata
ATCC 90030
2.23 ± 0.971.68 ± 0.731.09 ± 0.0
Candida krusei
ATCC 14243
2.23 ± 0.972.10 ± 0.891.09 ± 0.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhaparkulova, K.; Karaubayeva, A.; Sakipova, Z.; Biernasiuk, A.; Gaweł-Bęben, K.; Laskowski, T.; Kusniyeva, A.; Omargali, A.; Bekezhanova, T.; Ibragimova, L.; et al. Multidirectional Characterization of Phytochemical Profile and Health-Promoting Effects of Ziziphora bungeana Juz. Extracts. Molecules 2022, 27, 8994. https://doi.org/10.3390/molecules27248994

AMA Style

Zhaparkulova K, Karaubayeva A, Sakipova Z, Biernasiuk A, Gaweł-Bęben K, Laskowski T, Kusniyeva A, Omargali A, Bekezhanova T, Ibragimova L, et al. Multidirectional Characterization of Phytochemical Profile and Health-Promoting Effects of Ziziphora bungeana Juz. Extracts. Molecules. 2022; 27(24):8994. https://doi.org/10.3390/molecules27248994

Chicago/Turabian Style

Zhaparkulova, Karlygash, Aigerim Karaubayeva, Zuriyadda Sakipova, Anna Biernasiuk, Katarzyna Gaweł-Bęben, Tomasz Laskowski, Aliya Kusniyeva, Azamat Omargali, Tolkyn Bekezhanova, Liliya Ibragimova, and et al. 2022. "Multidirectional Characterization of Phytochemical Profile and Health-Promoting Effects of Ziziphora bungeana Juz. Extracts" Molecules 27, no. 24: 8994. https://doi.org/10.3390/molecules27248994

APA Style

Zhaparkulova, K., Karaubayeva, A., Sakipova, Z., Biernasiuk, A., Gaweł-Bęben, K., Laskowski, T., Kusniyeva, A., Omargali, A., Bekezhanova, T., Ibragimova, L., Ibadullayeva, G., Jakiyanov, A., Czech, K., Tastambek, K., Głowniak, K., Malm, A., & Kukula-Koch, W. (2022). Multidirectional Characterization of Phytochemical Profile and Health-Promoting Effects of Ziziphora bungeana Juz. Extracts. Molecules, 27(24), 8994. https://doi.org/10.3390/molecules27248994

Article Metrics

Back to TopTop