UPLC‐ESI‐QTOF‐MS/MS Profiling, Antioxidant, and Cytotoxicity Potentials of Marrubium vulgare L. Extracts: Experimental Analysis and Computational Validation
Abstract
Plant extracts are emerging as valuable options for food additives and therapeutic treatments. This study evaluated the phytochemical profile, antioxidant activity, and cytotoxicity of aerial parts of Marrubium vulgare L. crude extract (MVCE) and its subfractions. The MVCE (80% ethanol) contained steroids, phenolic compounds, flavonoids, terpenes, and cardiac glycosides, with total phenolic content (TPC) and total flavonoid content (TFC) of 14.96 ± 0.12 mg GAE/g DW and 12.27 ± 0.63 mg RE/g DW, respectively. All MV extracts exhibited potent antioxidant activity against DPPH• (0.106–1.864 mg/mL) and ABTS+• (0.298–17.084 mg/mL). The MV residual aqueous fraction (MVRF) showed significant cytotoxicity against human cancer cell lines, including MCF7 (IC50 = 5.47 ± 1.32 µg/mL), HT29 (IC50 = 17.48 ± 1.47 µg/mL), and SW480 (IC50 = 7.51 ± 0.36 µg/mL). Ultra‐performance liquid chromatography‐mass spectrometry identified 26 bioactive compounds, with malic acid, caffeic acid, chlorogenic acid, kaempferol‐3‐glucuronide, and l‐tryptophan as the major ones. Molecular docking revealed strong binding affinities of the above compounds to breast (PDB ID: 6CHZ) and colorectal cancer (PDB ID: 1HVY) proteins. Pharmacokinetic and toxicological studies confirmed their safety and efficacy, supporting MVRF as a potential therapeutic agent. These findings highlight MV as a promising candidate for future anticancer research.
Article type: Research Article
Keywords: antioxidant, cytotoxicity, in silico study, UPLC‐ESI‐QTOF‐MS/MS
Affiliations: Chemistry Department, Faculty of Science Al‐Baha University Al‐Baha Kingdom of Saudi Arabia; Research Unit Advanced Materials, Applied Mechanics, Innovative Processes, and Environment, Higher Institute of Applied Sciences and Technology of Gabes University of Gabes Tunisia; Department of Chemistry and Industrial Chemistry, College of Applied and Industrial Sciences University of Bahri Sudan; Laboratory of Natural Substances National Institute of Research and Physico‐Chemical Analysis Ariana Tunisia; Department of Pharmacology and Toxicology, College of Pharmacy Umm Al‐Qura University Makkah Kingdom of Saudi Arabia; Department of Pharmacology and Toxicology Medicinal and Aromatic Plants Research Institute, National Center for Research Khartoum Sudan; Department of Chemistry, College of Science Jouf University Sakaka Kingdom of Saudi Arabia; Faculty of Science of Sfax, Department of Chemistry University of Sfax Sfax Tunisia; Postgraduate Program in Pharmaceutical Sciences (PPGCF) Institute of Health Sciences, Federal University of Pará Belém Brazil
License: © 2025 The Author(s). Chemistry & Biodiversity published by Wiley‐VHCA AG. CC BY 4.0 This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
Article links: DOI: 10.1002/cbdv.202500400 | PubMed: 40516110 | PMC: PMC12533807
Relevance: Moderate: mentioned 3+ times in text
Full text: PDF (661 KB)
Introduction
Cancer is a significant public health concern at the global level, identified as a noncommunicable and multifactorial burden [ref. 1, ref. 2]. As indicated by the Global Cancer Observatory (GCO) database, cancer represents one of the primary causes of mortality in Gulf Cooperation Council (GCC) countries [ref. 3]. In 2022, Saudi Arabia had the highest global accounting number of new cases (28 113) and the highest number of deaths (13 399) due to cancer. The estimated 5‐year prevalent cases were 46 288 in males and 48 663 in females [ref. 4]. Breast cancer (BC) and colorectal cancer (CRC) are the most commonly prevalent malignancies and the leading causes of cancer death for both sexes in Saudi Arabia, accounting for 3777 (13.4%)/1046 and 3750 (13.3%)/1883 new cases/deaths, respectively, of all cancer incidence [ref. 5]. An excessive and aberrant formation of reactive oxygen species (ROS) above the cellular tolerability threshold, followed by a weakening of the immune defense system, can interfere with cancer pathology [ref. 6, ref. 7, ref. 8].
The elimination of these cells through the enhancement of scavenging abilities represents a promising strategy for cancer therapies [ref. 6, ref. 7, ref. 8]. Conventional cancer therapies include surgical procedures, radiation, hormonal and targeted therapies, and chemotherapy, which often result in significant adverse side effects [ref. 9, ref. 10]. Doxorubicin (DOX) is an effective chemotherapeutic drug that is widely recognized as one of the most popular drugs used to treat metastatic human BC [ref. 11, ref. 12]. In contrast, contemporary traditional pharmacopoeias derived from plant‐derived biomolecules have emerged as promising strategies for safer and more effective therapeutic purposes, with the objective of preventing or curing cancer through the discovery of potential anticancer drugs [ref. 13, ref. 14, ref. 15, ref. 16].
Marrubium vulgare L. (MV), a member of the Lamiaceae family, is an annual indigenous plant commonly known as the white horehound. It has a naturalized distribution in North Africa, Europe, and Central Asia [ref. 17]. The plant has traditionally been reported to have a range of therapeutic effects, including the treatment of inflammation, pulmonary infections, rheumatoid arthritis, night blindness, and stomach disorders [ref. 18, ref. 19, ref. 20]. Phytochemical analysis revealed the abundance of diterpenes (labdanes), particularly marrubenol and marrubiin [ref. 18, ref. 19, ref. 20, ref. 21, ref. 22]. Extracts from this plant provide a valuable source of bioactive molecules that can be utilized in the development of novel drugs.
The Marrubium genus plants have been demonstrated to possess high phenolic contents, comprising the majority of classes, including phenylpropanoid and phenolic acids, esters, flavonoids, flavone derivatives (flavone aglycones and flavone glycosides), and tannins [ref. 23, ref. 24]. In addition, the presence of nitrogen‐containing compounds, polysaccharides, phytosterols, amino acids, organic acids, alkaloids, and minerals, in particular potassium salts in MV has been documented [ref. 25, ref. 26, ref. 27]. The high efficiency of these phytochemicals for various purposes and their known beneficial pharmacological attributes have led to this plant being designated a “plant drug” for use in cancer therapy and general health. Biological investigations have confirmed the analgesic, antioxidant, diuretic, hypoglycemic, hypolipidemic, antihypertensive, antidiabetic, diaphoretic, antinociceptive, antispasmodic, gastroprotective, anti‐inflammatory, antimicrobial, vasodilating, and antiproliferative effects of this plant [ref. 22, ref. 28, ref. 29]. As an herbal beverage, it was frequently incorporated into teas, syrups, and other preparations. It was also used as a bitter and choleretic, as well as for ovoid dyspeptic complaints and in food applications [ref. 30].
The present study was designed to elucidate the phytoconstituents and metabolite profiling of MV extract(s) through the use of an untargeted metabolomic approach, employing ultrahigh‐performance liquid chromatography with quadrupole time‐of‐flight and mass spectrometry (UHPLC‐QTOF‐MS/MS). The antioxidant and cytotoxic effects of the extracts were evaluated, while molecular docking and pharmacokinetic analysis of the selected major phytocompounds from MV residual aqueous fraction (MVRF), which underline the potent anticancer property, were performed.
Results and Discussion
Phytochemical Screening, Total Phenolic Content, and Total Flavonoid Content
The qualitative phytochemical screening of the MV crude extract (MVCE) revealed the presence of steroids, phenolic compounds, flavonoids, terpenes, and cardiac glycosides, but the absence of alkaloids and saponins (results not shown). The total phenolic content (TPC) and total flavonoid content (TFC) in the MVCE were found to be 14.96 ± 0.12 mg GAE/g dry weight (DW) and 12.27 ± 0.63 mg RE/g DW, respectively. Some recent studies have focused on the TPC and TFC of this species. Michalak et al. demonstrated a TPC of 55.72 mg GAE/mL and a TFC of 11.01 mg catechin equivalent/mL in a 50/50 ethanol/water (v/v) extract from aerial parts of MV collected from Poland using ultrasonication [ref. 31]. Emam et al. determined TPC by spectrophotometry and obtained 240.46 ± 12.19 mg of GAE/g DW of phenolic content in a 70% ethanol extract from macerated aboveground parts originating from Egypt for 3 weeks [ref. 32]. The hot infusion of aerial parts of MV from Tunisia exhibited a richness in phenolic content (73.619 mg GAE/g extract) [ref. 33]. Gavarić et al. reported that the TPC and TFC in MV extracts are influenced by the extraction method in the following descending order: microwave extraction > ultrasound‐assisted extraction > conventional solid‐liquid extraction. TPC ranged from 63.77 to 117.58 mg GAE/g DW, while TFC ranged from 42.51 to 65.80 mg catechin equivalent/g DW [ref. 34]. On the other hand, Khouchlaa et al. reported that the TPC and TFC of the methanolic extract of MV aerial parts from Morrocco did not exceed 0.0430 mg GAE/g DW and 1.644 mg quercetin/g DW, respectively, during three development stages [ref. 35]. Meanwhile, the methanolic extract of MV leaves collected from the Al‐Taif Governorate (KSA) showed a TPC of 36.8 mg GAE/g DW [ref. 36].
Overall, the above results reveal that the content of secondary metabolites is influenced by genetics, geographical and environmental factors, regional climatic conditions, the type of solid/liquid extraction (common/green), the kind of extract (organic/aqueous), and the polarity index of the solvent [ref. 34, ref. 37–ref. 39].
Antioxidant Activity
The MVCE and its subfractions exhibited concentration‐dependent DPPH• and ABTS+• radical scavenging activities expressed as IC50 values (mg/mL) summarized in Table 1. Our findings demonstrate that all extracts significantly reduce DPPH• radical to varying degrees (p < 0.05). MVPF (1.864 ± 0.069 mg/mL) is the least effective, followed by MVCF and MVRF, with no observable significant difference, followed by MVCE, significantly less effective than MVEF. The latter fraction showed the highest radical quenching potential among all extracts, but it was comparable to the positive control (Trolox, 0.012 ± 0.000 mg/mL) with values of 0.106 ± 0.008 mg/mL as the difference in IC50 values was not significant (p = 0.161). The effects of MVCE and subfractions on the ABTS+• radical were also studied. All extracts were able to quench the ABTS+• cationic radical except for MVPF (17.084 ± 1.101 mg/mL). The effectiveness of the crude extract and corresponding subfractions as ABTS+• quenchers were ranked in the following increasing order MVCF (3.375 ± 0.120 mg/mL) < MVRF (1.727 ± 0.018 mg/mL), MVCE (1.473 ± 0.075) (p = 0.980) ≤ MVEF (0.298 ± 0.027), Trolox (0.055 ± 0.000 mg/mL) (p = 0.984).
TABLE 1: DPPH (IC50, mg/mL) and ABTS+• (IC50, mg/mL) of Marrubium vulgare L. aerial parts crude extract and subfractions.
| Sample | Solvent | DPPH• | ABTS+• |
|---|---|---|---|
| MVCE | EtOH (80%) | 0.349 ± 0.072c | 1.473 ± 0.075cd |
| MVPF | Petroleum ether | 1.864 ± 0.069a | 17.084 ± 1.101a |
| MVCF | Chloroform | 0.656 ± 0.033b | 3.375 ± 0.120b |
| MVEF | Ethyl acetate | 0.106 ± 0.008d | 0.298 ± 0.027de |
| MVRF | Residual aqueous fraction | 0.544 ± 0.017b | 1.727 ± 0.018c |
| Trolox | — | 0.012 ± 0.000d | 0.055 ± 0.000e |
Note: Values are mean of three replicates ± standard deviation (n = 3). Different letters in the same column indicate a significant difference according to the Tukey test at p < 0.05.
Overall, we can conclude that only MVCE, MVEF, and MVRF worked excellently against both DPPH• and ABTS+• to reduce stable radicals among all studied extracts. The descending order of radical scavenging potential was as follows: MVEF > MVCE > MVRF, but they were still less effective than the positive control. This can be explained by the abundance of phytochemicals in these extracts serving as potential hydrogen donors that induce the reduction of free radicals [ref. 40]. DPPH• and ABTS+• assays are both hydrogen atom transfer processes working on an identical mechanism [ref. 41]. Therefore, it is evident that they revealed consistent results.
Recently, Duc et al. attributed the effectiveness of phenolic compounds as oxidation inhibitors to their chemical structure [ref. 42]. The presence and location of hydroxyl groups on the aromatic ring are important, as they contribute to the low bond dissociation energy and strength in both phenolic acids and flavonoids. Furthermore, the effectiveness of flavonoids in neutralizing free radicals is related to the presence of a double bond between C‐2 and C‐3 in conjugation with a carbonyl group at C‐4 and the absence of sugar units at C‐3. Previous studies reported the capacity of MV to quench DPPH• and ABTS+• radicals due to the presence of several potent antioxidants such as phenylethanoid glycosides, tannins, flavonoids, and phenolic acids [ref. 43].
Correlation Between Phenolic Compounds and Antioxidant Activity
In this study, a Pearson correlation analysis was conducted to explore the relationships between the TPC and TFC and antioxidant activity (DPPH• and ABTS+•) of the MVCE. The objective of this analysis was to ascertain the manner in which alterations in phenolic and flavonoid concentrations are associated with the plant’s antioxidant capabilities, thereby offering insight into the contribution of these compounds to the plant’s overall salutary effects. As illustrated in Table 2, there is a robust correlation between TPC and TFC (r = 0.887). The correlation coefficients between DPPH• and (TPC and TFC) were r = −0.852 (p > 0.05) and −0.998 (p < 0.05), respectively, indicating a significant negative correlation between phenolic acids and flavonoids and DPPH• inhibition. Furthermore, a moderate positive correlation was identified between ABTS+• and TPC (r = 0.439, p > 0.05), while a strong positive correlation was observed between ABTS+•and TFC (r = 0.805, p > 0.05).
TABLE 2: Pearson’s correlation test (r values) between TPC, TFC, and antioxidant capacity (DPPH• and ABTS+• assays) for MVCE.
| TPC | TFC | DPPH• | ABTS+• | |
|---|---|---|---|---|
| TPC | 1 | 0.887 | −0.852 | 0.439 |
| TFC | 1 | −0.998cbdv70100-tbl2-note-0001 | 0.805 | |
| DPPH• | 1 | −0.845 | ||
| ABTS+• | 1 |
* Significant correlation at p ≤ 0.05.
This contrasting antioxidant behavior may be explained by factors such as radical stability, reaction mechanism, and the structural properties of polyphenols, including basic structure, hydroxyl group number and position, and Bors criteria [ref. 44, ref. 45]. It is also important to consider other factors, such as pH, solvent, carbohydrate, and protein content, as these can influence the reported activity of natural extracts [ref. 46]. In conclusion, TPC and TFC exhibited a strong correlation with the DPPH• assay compared to the ABTS+• assay. This underscores the significant role of phenolics and flavonoids in antioxidant capacity, especially in terms of DPPH• radical scavenging activity.
Cytotoxicity on Normal and Cancerous Cells
The cytotoxicity of the crude extract of MV and its subfractions was evaluated after 72 h of incubation on three cancerous cell lines (MCF7, HT29, and SW480) and one non‐cancerous MRC5 cell line. The IC50 values (µg/mL) were determined through the use of the MTT assay. All extracts exhibited dose‐dependent cytotoxic activity against the three tested cancerous cell lines. As shown in Table 3, the crude extract and its five subfractions exhibited notable cytotoxicity against the aforementioned cell lines, with p values less than 0.05. The results demonstrated notable discrepancies among the cell lines and extracts, with some exhibiting comparatively diminished activity relative to DOX.
TABLE 3: Cytotoxic activity of Marrubium vulgare crude extract and subfractions against three cell lines and normal fibroblast (MTT 72 h, IC50 “µg/mL” ± SD, n = 3).
| Sample | Solvent | MCF7 | HT29 | SW480 | MRC5 |
|---|---|---|---|---|---|
| MVCE | EtOH (80%) | 13.60 ± 0.90d | 44.78 ± 2.70a | 4.92 ± 0.77d | 20.84 ± 0.34c |
| MVPF | Petroleum ether | 18.83 ± 0.46c | 15.85 ± 0.46c | 25.05 ± 1.93b | 21.52 ± 0.38c |
| MVCF | Chloroform | 22.06 ± 0.73b | 28.29 ± 2.07b | 35.35 ± 1.15a | 22.38 ± 0.37c |
| MVEF | Ethyl acetate | 24.68 ± 1.40a | 26.32 ± 2.31b | 14.81 ± 1.92c | 45.27 ± 1.90a |
| MVRF | Residual aqueous fraction | 5.47 ± 1.32 e | 17.48 ± 1.47c | 7.51 ± 0.36d | 35.77 ± 3.45b |
| Doxorubicin | — | 0.04 ± 0.00f | 0.05 ± 0.01d | 0.06 ± 0.01e | 2.48 ± 0.02d |
Note: Average cytotoxicity (IC50) of each extract against the three cancer cells. Values are mean of three replicates ± standard deviation (n = 3). Different letters in the same column indicate a significant difference according to the Tukey test (p < 0.05).
The IC50 values ranged from 5.47 ± 1.32 to 24.68 ± 1.40 µg/mL against MCF7. The order of the cytotoxic effects of the extracts is as follows: MVEF < MVCF < MVPF < MVCE < MVRF. With regard to HT29, the IC50 values exhibited a range of 15.85 ± 0.46–44.78 ± 2.70 µg/mL. The results demonstrated that MVCE was less potent than MVCF and MVEF (p = 0.756), which in turn were less potent than MVPF and MVRF (p = 0.866). Against SW480, the IC50 values ranged from 7.51 ± 0.36 to 35.35 ± 1.15 µg/mL. The order of the cytotoxic effects of the extracts was as follows: The order of the cytotoxic effects of the extracts was as follows: MVCF < MVPF < MVEF < MVCE, MVRF (p = 0.191). In contrast, the crude extract and its subfractions demonstrated a moderate to weak impact on the healthy MRC5 cell line, with IC50 values ranging from 13.29 ± 0.92 to 45.27 ± 1.90 µg/mL. The observed effect was found to diminish in the following descending order: MVCE, MVPF, and MVCF (p = 0.849) exhibited greater potency than MVRF and MVEF. As can be observed, MRC5 demonstrates reduced sensitivity to MV extract action, which is an anticipated outcome given that these normal cells are untransformed and of healthy connective tissue cell nature [ref. 47]. In general, human breast and colorectal cell lines demonstrate greater sensitivity to polar fractions than less polar fractions. MVRF demonstrated the greatest efficacy against all cancerous cell lines. It can thus be concluded that the adopted fractionation process resulted in the production of fractions with greater activity than the crude extract. These results can be attributed to the disparities in phytochemical composition, contingent upon the fraction’s polarity in relation to the solvent extraction type and polarity index.
UPLC‐ESI‐QTOF‐MS/MS Characterization of the Residual Aqueous Fraction of MV
The phytochemical constituents of the MVRF were tentatively identified by means of UHPLC‐QTOF‐MS/MS. In addition, the HMDB (http://www.hmdb.ca/) and Mass Bank (https://massbank.eu/MassBank/) databases were consulted for the purpose of identification, utilizing the precise masses of the molecules in question. In total, 26 peaks were identified, with 11 in the ESI negative mode, 12 in the ESI positive mode, and 3 compounds identified in both modes. The results are presented in Table 4, and the total ion chromatograms (TICs) for the ESI− and ESI+ modes are displayed in supplementary material Figure S1a,b, respectively. The identified secondary metabolites belong to various natural product classes, including three organic acids, two phenolic acids, six flavonoid glycosides, and other compounds, as detailed in Table 4. The MS/MS spectra and structures of some identified compounds are presented in Figures S2 and S3, respectively.
TABLE 4: UPLC‐ESI‐QTOF‐MS/MS metabolic analysis of MVRF.
| No. | Compound | RT (min) | Experimental (m/z) | Area | Error ppm | Ionization (ESI+/ESI−) | Theoretical (m/z) | Formula | MS2 ion products | References |
|---|---|---|---|---|---|---|---|---|---|---|
| Organic acids | ||||||||||
| 1 | Succinic acid | 1.008 | 117.0195 | 23 623.3 | −1.2 | [M−H]− | 117.01933 | C4H6O4 | 55.0201273.0300799.00807 | [ref. 48, ref. 49] |
| 2 | Malic acid | 1.046 | 133.0128 | 67 580.6 | −1.8 | [M−H]− | 133.0142 | C4H6O5 | 71.013489.0231115.0024 | [ref. 50, ref. 51, ref. 52] |
| 3 | Citraconic acid | 1.474 | 128.96 | 70 409.3 | −0.4 | [M−H]− | 129.01933 | C5H6O4 | 84.98897 | MassBank |
| Phenolic acids | ||||||||||
| 4 | Caffeic acid | 1.087 | 179.0569 | 33 793.7 | −0.7 | [M−H]− | 179.03499 | C9H8O4 | 71.0128475.0082289.031999.00691143.03537161.04589 | [ref. 49, ref. 53–ref. 57] |
| 5 | Chlorogenic acid | 6.087 | 355.1726 | 17 164.8 | 1.2 | [M+H]+ | 355.10236 | C16H18O9 | 163.0387 | [ref. 58, ref. 59, ref. 60, ref. 61] |
| Flavonoids | ||||||||||
| Flavonols | ||||||||||
| 6 | Quercitrin | 1.309 | 447.1172 | 6709.2 | 0.5 | [M−H]− | 447.09329 | C21H20O11 | 149.03605151.05052152.01283299.08168 | [ref. 62] |
| 7 | Kaempferol‐3‐glucuronide | 6.251 | 461.074 | 71 500.5 | 1.3 | [M−H]− | 461.07254 | C21H18O12 | 175.03933285.04057 | [ref. 62, ref. 63, ref. 64, ref. 65] |
| 9.756 | 463.0881 | 103 489.3 | 0.5 | [M+H]+ | 463.0871 | C21H18O12 | 287.05566 | [ref. 66] | ||
| 8 | Isorhamnetin‐3‐O‐rutinoside | 7.484 | 623.2025 | 27 863.4 | −0.4 | [M−H]− | 623.16174 | C28H32O16 | 161.02477315.10612461.16721623.19975 | [ref. 65, ref. 67] |
| Flavones | ||||||||||
| 9 | Vicenin‐2 | 5.785 | 593.1548 | 10 424.5 | 0.8 | [M−H]− | 593.15118 | C27H30O15 | 353.06696383.07783473.11029593.1507 | [ref. 68, ref. 69, ref. 70, ref. 71, ref. 72, ref. 73, ref. 74] |
| 6.117 | 595.1658 | 8988.8 | −0.4 | [M+H]+ | 595.16577 | C27H30O15 | 295.05884325.07021337.06948379.08152457.11172529.1370559.14537577.15467 | [ref. 75, ref. 76] | ||
| 10 | Baicalein‐7‐O‐glucuronide | 6.732 | 445.0813 | 10 455.9 | −0.7 | [M−H]− | 445.07764 | C21H18O11 | 175.02541269.04758 | [ref. 65, ref. 72, ref. 77] |
| Flavanones | ||||||||||
| 11 | Naringenin‐7‐O‐glucoside | 8.845 | 435.1279 | 8836.5 | 0.7 | [M+H]+ | 435.12857 | C21H22O10 | 273.07332 | [ref. 78] |
| Amino acids | ||||||||||
| 12 | l‐Tryptophan | 4.303 | 203.0833 | 17 746.0 | 0 | [M−H]− | 203.0826 | C11H12N2O2 | 72.0085274.02491116.05056142.06585159.09128186.0536 | [ref. 77] |
| 4.521 | 205.0973 | 17 706.8 | −1.6 | [M+H]+ | 205.09715 | C11H12N2O2 | 118.0647146.0593188.0700 | [ref. 79] | ||
| 13 | l‐Ornithine | 1.135 | 133.061 | 16 270.88 | −1 | [M+H]+ | 133.09715 | C5H12N2O2 | 74.0257 | HMDB |
| 14 | l‐Proline | 1.236 | 116.0704 | 14 577.7 | −1.4 | [M+H]+ | 116.0706 | C5H9NO2 | 70.0652 | [ref. 80] |
| 15 | Tyrosine | 1.261 | 182.0581 | 7844.0 | 0.2 | [M+H]+ | 182.08118 | C9H11NO3 | 91.0597 | [ref. 79] |
| 16 | Norvaline | 1.308 | 118.0866 | 18 930.7 | −2.9 | [M+H]+ | 118.08626 | C5H11NO2 | 55.057172.0847 | HMDB |
| Bile acids | ||||||||||
| 17 | Cholic acid | 14.050 | 409.1622 | 33 542.04 | −0.1 | [M+H]+ | 409.29486 | C24H40O5 | — | HMDB |
| Triterpene saponins | ||||||||||
| 18 | Glycyrrhizin | 23.218 | 821.5222 | 36 026.2 | 0.8 | [M−H]− | 821.39648 | C42H62O16 | — | [ref. 81, ref. 82] |
| Other compounds | ||||||||||
| 19 | Anserine | 11.195 | 239.0683 | 9384.6 | 0 | [M−H]− | 239.11496 | C10H16N4O3 | 176.03577207.04123 | HMDB |
| 20 | 5‐Aminoimidazole‐4‐Carboxamide‐1‐ribofuranosyl 5′‐monophosphate | 16.633 | 337.2076 | 12 561.2 | −0.4 | [M−H]− | 337.05548 | C9H15N4O8P | 96.96002 | HMDB |
| 21 | 1‐Myristoyl‐2‐hydroxy‐sn‐glycero‐3‐phosphate | 17.143 | 381.2335 | 7071.0 | −0.9 | [M−H]− | 381.20477 | C17H35O7P | 96.95987 | [83]HMDB |
| 22 | Nicotanamide | 1.646 | 122.9235 | 20 565.8 | 5.5 | [M+H]+ | 123.05529 | C6H6N2O | 80.0477 | [ref. 84] |
| 23 | S‐Lactoylglutathione | 1.024 | 380.0962 | 37 087.3 | 0.2 | [M+H]+ | 380.11221 | C13H21N3O8S | 260.0516 | [ref. 85, ref. 86] |
| 24 | l‐β‐Homoproline | 1.248 | 130.0505 | 8256.4 | −2.5 | [M+H]+ | 130.08626 | C6H11NO2 | 56.049484.0450 | [ref. 87] |
| 25 | Pipecolate | 1.309 | 130.0862 | 55 965.4 | −1.6 | [M+H]+ | 130.08626 | C6H11NO2 | 84.0845 | [88]MassBank |
| 26 | 4‐Methyl‐5‐thiazole ethanol | 1.647 | 144.1015 | 167 217.7 | 0.2 | [M+H]+ | 144.04776 | C6H9NOS | 58.064884.0802 | [ref. 89] |
Identification of Organic Acids and Derivatives
Succinic acid (1.008 min) was identified based on its [M─H]− ion at m/z 117.0199 corresponding molecular formula C4H6O4 and the resulting fragments ions at m/z 55.0201 [M─CO2─H2O]−, 73.0301 [M─CO2]− and 99.0081 [M─H2O]−. These findings are consistent with the data reported in the literature [ref. 48, ref. 49]. To the best of our knowledge, this is the first time that succinic acid has been reported in the Marrubium genus.
Malic acid (1.046 min) was identified based on its [M−H]− ion at m/z 133.0128 matching C4H6O5 Bursal et al. as well as the resulting product ions at m/z 115.0024, 89.0231, and 71.0134 corresponding to [M─H─H2O]−, [M─H─CO2]− and [M─H─H2O─CO2]−, respectively [ref. 90]. This is in excellent agreement with the previous reports [ref. 50, ref. 52]. Malic acid is found to be rich in Marrubium astracanicum subsp. macrodon [ref. 90]. Citraconic acid (RT = 1.474 min) was identified based on its deprotonated precursor [M−H]− at m/z 128.96 matching C5H6O4 and a base peak at m/z 84.98897 corresponding to [M─H─CO2]−. The fragmentation pattern matches the description of citraconic acid as documented in the MassBank database.
Identification of Phenolic Acids
In the present study, only two phenolic acids were identified in the residual aqueous fraction classified as hydroxycinnamic acids. Caffeic acid, also known as 3,4‐dihydroxycinnamic acid (RT = 1.087 min), was characterized in a negative mode based on its [M─H]− at m/z 179.0575 fitting C9H8O4 [ref. 49, ref. 55, ref. 57]. It was also confirmed with an MS/MS spectrum showing characteristic fragment ions resulting from continuous losses of water, CO, and CO2 moieties at m/z 161.0459 [M─H2O]−, 143.0354 [M─H─2H2O]−, 135 [M─H─CO2]−, 99.0069 [M─H─2H2O─CO2]−, and 89.0326 [M─H─CO2─CO─H2O]−. This fragmentation pattern aligned with previous studies [ref. 53, ref. 54, ref. 56, ref. 91]. Another fragment at m/z 71.0128 corresponds to the ion [M─H─C6H4O2]‐. Caffeic acid has previously been identified in MV [ref. 31, ref. 43, ref. 92, ref. 93] and other Marrubium species such as Marrubium trachyticum Boiss [ref. 94], Marrubium peregrinum L., and Marrubium friwaldskyanum Boiss [ref. 95].
Chlorogenic acid, also known as 5‐caffeoylquinic acid (RT = 6.087 min) was identified in a positive mode based on its protonated molecular ion [M+H]+ at m/z 355.1042 and fitting C9H8O4 [ref. 58, ref. 59, ref. 61]. Other characteristic peaks were detected at m/z 163.0403, 145.0298, 135.0466, and 117.0314 [ref. 58, ref. 59]. The fragment ion at m/z 163.0387 may be due to the loss of the C7H11O6 moiety resulting from the cleavage of the caffeoyl bond followed by the neutral loss of a water molecule [ref. 59, ref. 60]. The fragment ions at 145.0298 and 135.0466 are obtained from the neutral loss of H2O and CO from the ion at m/z 163.0387, respectively. Both ions yielded the fragment at m/z 117.0314 by neutral loss of H2O or CO, as suggested by Willems et al. [ref. 59]. Chlorogenic acid was identified in MV [ref. 20, ref. 96] and other Marrubium species [ref. 97] Marrubium lutescens and M. trachyticum Boiss [ref. 94, ref. 98]. Caffeic acid and chlorogenic acid are common hydroxy derivatives of cinnamic acid [ref. 99]. Both acids are recognized for their high therapeutic potential due to their antioxidant, anticancer, antimicrobial, and many other properties [ref. 85, ref. 100, ref. 101, ref. 102, ref. 103].
Identification of Flavonoids
In the current study, three subclasses of flavonoids in the residual aqueous fraction were identified. These subclasses include three flavonols, two flavones, and one flavanone, all of which are O‐ and C‐glycosides.
The compound eluted at 1.309 min showed a parent ion [M−H]− at m/z 447.1173. It was tentatively identified as quercitrin (quercetin 3‐rhamnoside) matching C21H20O11 [ref. 62]. The peaks at m/z 149 and 151 were obtained from the retro‐Diels Alder of the aglycone quercetin and corresponded to 1,3B− and 1,3A− ions, respectively. Many quercetin O‐glycosides were identified in some Marrubium species such as quercitrin‐O‐gallate, quercetin‐7‐O‐galloyl‐glucoside, and quercetin‐diglucuronide in Marrubium anisodon as well as quercetin‐3‐O‐glucoside‐6″‐acetate, quercetin‐3‐O‐pentosyl‐pentoside, quercetin‐3‐O‐glucuronide and quercetin‐3‐O‐glucoside6″‐acetate in Marrubium crassidens [ref. 104]. Ibrahim et al. [ref. 105] reviewed that quercetin 3‐O‐β‐d‐rutinoside, quercetin 3‐O‐β‐d‐glucoside, quercetin 3‐O‐α‐1‐rhamnosyl‐glucoside, and the aglycone quercetin were identified in MV but it seems that quercitrin was identified for the first time in this species. Chen et al. reviewed the wide spectrum of pharmacological activities exhibited by quercitrin, which include reducing oxidative stress and inflammation and inhibiting microorganisms [ref. 106]. In addition, quercitrin has been shown to protect vital organs and demonstrate therapeutic potential for various diseases such as gastric ulcers, immune diseases, and cancer.
Kaempferol‐3‐glucuronide (RT = 6.251 min) was identified due to its parent ion [M−H]− at m/z 461.0727 corresponding to C21H18O12 [ref. 62, ref. 64], which generated a predominant daughter ion at m/z 285.0406 indicating the loss of a glucuronic acid residue (176 Da) [ref. 65]. Another peak at m/z 175 was observed and explained by the loss of an aglycone portion [ref. 63]. In the positive mode, this compound was eluted at 9.756 min and showed a parent ion [M + H]+ at m/z 463.0873, which generated an intense daughter ion at m/z 287.0552 relative to the ion [M + H−176]+ [ref. 66]. Kaempferol‐3‐glucuronide was identified in M. crassidens [ref. 104]. It was reported that kaempferol‐3‐glucuronide has shown high scavenging activity and effectiveness in treating various types of cancers [ref. 51, ref. 107].
Isorhamnetin‐3‐O‐rutinoside (RT = 7.484 min) was identified as it yielded a negative precursor m/z = 623.1998 [M−H]− corresponding molecular formula C28H32O16 [ref. 72], which is confirmed by the detection of fragment ions at m/z values 461.1672 [M−H−162]− and 315.1061 [M−H−162−146]− corresponding to the aglycone part (isorhamnetin). The total loss of 308 Da is characteristic of the disaccharide moiety (rutinose) [ref. 65, ref. 67]. This compound was identified in Marrubium globosum [ref. 105] and M. friwaldskyanum [ref. 95]. As reported in previous research papers, it was found that isorhamnetin‐3‐O‐rutinoside exhibited high scavenging and antiproliferative effects [ref. 108, ref. 109].
The compound eluted at 5.785 min gave an [M−H]− ion at m/z 593.1507 in negative mode, and in positive mode, it eluted at 6.117 min and yielded an [M + H]+ ion at m/z 595.1673 indicating a molecular formula of C27H30O15, was identified as vicenin‐2 also known as apigenin 6,8‐di‐C‐glucoside. The fragmentation pattern aligns with the description of vicenin‐2 [ref. 70, ref. 73, ref. 74, ref. 110]. On the other hand, the fragmentation pattern of the protonated vicenin‐2 ion [M+H]+ is complex. Abundant water‐loss ions at m/z 577 ([M+H−H2O]+), m/z 559 [M + H−2H2O]+, and m/z 529 ([M + H−CH2O−2H2O]+) were observed [ref. 76, ref. 111]. In addition, other peaks derived from the loss of 120 and/or 150 Da from the protonated molecule, such as m/z 2950,1X+ 1 0,1X+ 2 [M+H−150−150]+; m/z 3250,1X+ 1 0,2X+ 2 [M+H−120−150]+, m/z 3370,2X+ 1 0,2X+ 2−H2O [M+H−120−120−18]+; m/z 3790,2X+ 1 0,4X+ 2−2H2O [M + H−120−96]+ and m/z 4570,2X+ 1−H2O [M + H−120−18]+. This fragmentation pattern is consistent with previous studies [ref. 69, ref. 75, ref. 76, ref. 112]. Fragmentation sites were displayed in the supporting information Figure S4. Apigenin and some derivatives have been identified in MV as O‐glycosides, O‐glucuronides, and C‐glycosides [ref. 113, ref. 114]. Vicenin‐2 has been noted to exhibit several pharmacological properties, including antioxidant liver‐protective, anti‐inflammatory, and anticancer effects [ref. 115].
Baicalein‐7‐O‐glucuronide, also known as baicalin (RT = 6.732 min) was tentatively identified based on the molecular ion [M−H]− at m/z = 445.0794 fitting C21H18O11 along with a predominant fragment ion at m/z value 269.0476 [M−H−176]− which corresponds to the loss of glucuronic acid [ref. 65, ref. 72]. Another peak at m/z 175.0254 was observed and explained by the loss of an aglycone portion [ref. 77]. This compound was previously found in the Marrubium genus and recently demonstrated in MV. In a recent review, Bao et al. emphasized the health benefits of baicalin, such as its anti‐inflammatory properties, ability to reduce oxidative stress, inhibit various viruses, and potentially suppress the growth of different tumor cells [ref. 116].
Naringenin‐7‐O glucoside, also known as prunin (RT = 8.845 min) was identified based on its protonated molecular ion [M+H]+ at m/z 435.1280 which corresponds C21H22O10 [ref. 78]. The peak detected at 273.0377 resulted from the loss of glycoside moiety (162 Da), which is an indicator of aglycon residue (naringenin). Recently, naringenin 7‐O‐glucoside was identified in Marrubium alysson L. and M. globosum [ref. 88, ref. 105]. In a recent study, naringenin‐7‐O‐glucoside was tested in vitro against a panel of cancer cell lines. It was demonstrated that this flavone glycoside has the potential to be a targeted therapeutic agent [ref. 117].
Implications of Major Compounds in Biological Activity
Among the 26 compounds identified in MVRF using UPLC‐ESI‐QTOF‐MS/MS, five major compounds were selected: malic acid, caffeic acid, chlorogenic acid, kaempferol‐3‐glucuronide, and l‐tryptophan. These compounds are known for their antioxidant activity and cytotoxic effects.
The antioxidant potential of malic acid [ref. 118, ref. 119], caffeic acid [ref. 102, ref. 120, ref. 121], chlorogenic acid [ref. 121, ref. 122], kaempferol‐3‐glucuronide (kaempferol linked to a glucuronic acid moiety via a glycosidic bond at the 3‐position of the flavonoid) [ref. 123], and l‐tryptophan [ref. 124] has been reported. In addition, the potential of malic acid to induce anticarcinogenic and genotoxic damage in human fibroblast cells (HDFa) and glioblastoma (U87‐MG) cell lines has been demonstrated [ref. 125]. The antitumor action of caffeic acid, including its inhibitory effect on cell migration and reducing metastases in tumor cells, has been investigated [ref. 126, ref. 127]. Also, the anticancer potential of caffeic acid in several human cancers is well documented [ref. 128, ref. 129, ref. 130]. Chlorogenic acid has been noted for its role in cancer prevention and therapy [ref. 131], with evidence suggesting its ability to induce apoptosis and cell‐cycle arrest [ref. 132]. The biological effects of kaempferol‐3‐glucuronide, similar to kaempferol, have been observed in several cancer cell lines. In addition, the role of kaempferol in inhibiting the invasion of human breast carcinoma has been explored [ref. 133].
Based on the above literature data, malic acid, caffeic acid, chlorogenic acid, kaempferol‐3‐glucuronide, and l‐tryptophan were selected as key compounds for subsequent molecular docking calculations and pharmacokinetic profiling.
Molecular Docking Analysis
Estrogen receptor alpha (ERα) is a nuclear hormone receptor activated by the hormone estrogen. It plays a crucial role in regulating the growth and development of tissues, particularly in the breast and reproductive organs. ERα is significantly involved in the pathogenesis of hormone‐driven cancers, especially BC, where its overexpression promotes tumor growth and survival [ref. 134, ref. 135]. In contrast, human thymidylate synthase (TS) is a critical enzyme in the folate pathway, responsible for catalyzing the conversion of deoxyuridylate (dUMP) to thymidylate (dTMP), a nucleotide essential for DNA synthesis and repair. The activity of TS is particularly vital for rapidly dividing cells, such as cancer cells, which have an increased demand for DNA replication [ref. 136]. The docking results for the phytocompounds against ERα (PDB ID: 6CHZ) and TS (PDB ID: 1HVY) are summarized in Table 5, with docking scores indicating the binding affinities of each compound. In the 6CHZ protein, malic acid exhibited a binding affinity of −6.254 kcal/mol, interacting with the active site of the receptor 6CHZ through key interactions (Figure 1A). A hydrogen bond was observed with residues Glu353 and Arg394, suggesting stabilization through electrostatic interactions. In addition, Arg394 was involved in an ionic interaction, which may enhance the binding strength and contribute to the overall stability of the complex. With a binding affinity of −6.992 kcal/mol, caffeic acid demonstrated stronger interactions, particularly through hydrogen bonding. Arg394 and Glu353 were key residues forming hydrogen bonds, indicating a firm and potentially selective binding mode (Figure 1B). These interactions could facilitate the proper positioning of the compound in the binding pocket, enhancing its potential biological activity.
TABLE 5: Docking score (in kcal/mol) of the phytocompounds in the active site of estrogen receptor alpha (PDB ID: 6CHZ) and human thymidylate synthase (PDB ID: 1HVY).
| Phytocompounds | 6CHZ | 1HVY |
|---|---|---|
| Malic acid | −6.254 | −6.848 |
| Caffeic acid | −6.992 | −6.228 |
| Chlorogenic acid | −6.164 | −5.228 |
| Kaempferol‐3‐glucuronide | −8.715 | −5.66 |
| l‐Tryptophan | −7.652 | −4.157 |

Chlorogenic acid displayed a moderate binding affinity of −6.164 kcal/mol, forming a critical interaction with Asp351 (Figure 2C). This single interaction may suggest a less robust binding compared to other compounds, but its affinity could still contribute to significant biological activity, depending on the overall dynamics of the binding site environment. Kaempferol‐3‐glucuronide emerged as the most potent among the compounds, with a binding affinity of −8.715 kcal/mol. The compound exhibited multiple interactions with key residues, including Leu387, Glu353, Asp351, and Leu346 (Figure 1D). This multitude of interactions suggests strong and stable binding, likely due to the complementarity between the molecule and the receptor pocket. The involvement of both hydrogen bonds and hydrophobic contact could explain the high binding affinity observed. l‐Tryptophan displayed a binding affinity of −7.652 kcal/mol and stood out due to its unique π–π interaction with Phe404, suggesting a key role for aromatic stacking in stabilizing the compound within the binding site (Figure 1E). In addition, ionic interactions with Arg394 and Glu353 further stabilized the complex, reinforcing the versatility of l‐tryptophan in engaging with multiple interaction types within the binding site. The co‐crystallized ligand (H3B‐9224) in this protein exhibited critical hydrogen bonding interactions with Glu353 and Arg394, which are key residues within the active site [ref. 137]. These interactions play a vital role in stabilizing the ligand and facilitating its optimal orientation for effective binding. A similar binding pattern was observed with the identified phytochemicals, including hydrogen bonding with these pivotal residues. This similarity suggests that the phytochemicals mimic the binding mechanism of the co‐crystallized ligand. The observed interactions highlight the structural and chemical compatibility of the phytochemicals with the binding site, reinforcing their potential as effective inhibitors of the 6CHZ protein.

The molecular docking study of selected phytocompounds against the protein target (PDB ID: 1HVY) revealed varying binding affinities and interactions with key residues. Among the tested compounds, malic acid demonstrated the highest binding affinity with a docking score of −6.848 kcal/mol. It exhibited multiple interactions, including hydrogen bonds with Tyr258, Ser216, and Arg50, and ionic interactions with Arg176, Arg215, and Arg50, suggesting a strong and stable binding to the active site (Figure 2A). Caffeic acid, with a docking score of −6.228 kcal/mol, showed hydrogen bonding with Asp218, Asn226, and Arg215, as well as ionic interactions with Arg50, Arg175, and Arg176. In addition, a π–π interaction was observed with Ala414, further stabilizing its binding within the active site (Figure 2B). Chlorogenic acid had a lower binding affinity (−5.228 kcal/mol), forming a single hydrogen bond with Lys107 (Figure 2C). Kaempferol‐3‐glucuronide, with a docking score of −5.66 kcal/mol, interacted through hydrogen bonding with Lys107 and Met309, indicating moderate binding strength (Figure 2D). l‐Tryptophan, which exhibited the lowest binding affinity (−4.157 kcal/mol), formed hydrogen bonds with Met309 and ionic interactions with Lys308 (Figure 2E). While the binding strength was relatively weaker, these interactions suggest their potential for stabilization within the active site. The co‐crystallized ligand raltitrexed, a well‐known inhibitor, served as a benchmark for evaluating the binding interactions of the tested phytochemicals. Raltitrexed exhibited critical interactions within the active site of the target protein, including ionic interactions with Lys308, which play a pivotal role in stabilizing its binding conformation [ref. 138]. Furthermore, it established a hydrogen bond with Lys77, contributing to its anchoring within the binding pocket and enhancing its specificity towards the target. In addition, a salt bridge interaction with Lys308 was observed, reinforcing the ligand’s overall stability and interaction strength. Overall, the results indicate that malic acid and caffeic acid have the strongest binding interactions with the target protein, primarily through hydrogen bonding and ionic interactions with critical residues, making them promising candidates for further investigation.
Pharmacokinetic and Drug‐Likeness Properties
As a preliminary step in drug discovery research and to forecast insights into their potential use as therapeutic agents, the selected phytoconstituents (Table 6) in MVRF were tested for their pharmacokinetic properties based on their absorption, distribution, metabolism, excretion, and toxicity (ADMET) profile [ref. 139, ref. 140, ref. 141, ref. 142, ref. 143]. Predicted values of Caco‐2 permeability indicate moderate to lower absorption. It was predicted that compounds 2, 3, and 5 would exhibit good intestinal absorption (in humans), with values exceeding 30%. With regard to the permeability of the selected phytoconstituents across the skin, all exhibited low permeability. The fraction unbound, which delineates the extravasate portion of the free drug in plasma, is within the range of 0.28 and 0.568, indicating that the drug exhibited a favorable distribution profile in the plasma and primarily interacted with the pharmacological target.
TABLE 6: Pharmacokinetic profile assessment of major MVRF phytoconstituents.
| Entry | 1 | 2 | 3 | 4 | 5 | Reference |
|---|---|---|---|---|---|---|
| Absorption | ||||||
| Water solubility | −1.381 | −2.33 | −2.449 | −2.866 | −2.891 | — |
| Caco‐2 permeability | −0.395 | 0.634 | −0.84 | −0.884 | 0.638 | > 0.9 |
| Intestinal absorption (human) | 13.831 | 69.407 | 36.377 | 25.165 | 77.224 | < 30% is poorly |
| Skin permeability (log Kp) | −2.735 | −2.722 | −2.735 | −2.735 | −2.735 | > −2.5 is low |
| P‐glycoprotein inhibitor | No | No | Yes | Yes | Yes | No |
| Distribution | ||||||
| VDss (human) | −0.998 | −1.098 | 0.581 | 1.295 | −0.081 | Low is < −0.15, high is > 0.45 |
| Fraction unbound (human) | 0.652 | 0.529 | 0.658 | 0.28 | 0.577 | — |
| BBB permeability | −0.788 | −0.647 | −1.407 | −1.441 | −0.495 | Poorly is < −1, high is > 0.3 |
| CNS permeability | −3.523 | −2.608 | −3.856 | −3.955 | −2.622 | Penetrate is > −2, unable is < −3 |
| Metabolism | ||||||
| CYP1A2 inhibitior | No | No | No | No | No | No |
| CYP2C19 inhibitior | No | No | No | No | No | No |
| CYP2C9 inhibitior | No | No | No | No | No | No |
| CYP2D6 inhibitior | No | No | No | No | No | No |
| CYP3A4 inhibitior | No | No | No | No | No | No |
| Excretion | ||||||
| Total clearance | 0.81 | 0.508 | 0.307 | 0.503 | 0.64 | — |
| Renal OCT2 substrate | No | No | No | No | No | — |
| Toxicity | ||||||
| AMES toxicity | No | No | No | No | No | No |
| Max. tolerated dose (human) | 1.212 | 1.145 | −0.134 | 0.46 | 0.793 | Low is ≤ 0.477, high is > 0.477 |
| hERG I inhibitor | No | No | No | No | No | No |
| hERG II inhibitor | No | No | No | No | No | No |
| Oral rat acute toxicity (LD50) | 1.818 | 2.383 | 1.973 | 2.513 | 2.357 | — |
| Oral rat chronic toxicity (LOAEL) | 3.104 | 2.092 | 2.982 | 4.641 | 1.231 | — |
| Hepatotoxicity | No | No | No | No | No | No |
| Skin sensitization | No | No | No | No | No | No |
| Tetrahymena pyriformis toxicity | 0.285 | 0.293 | 0.285 | 0.285 | 0.285 | — |
| Minnow toxicity | 3.348 | 2.246 | 5.741 | 6.898 | 1.835 | — |
Note: 1, malic acid; 2, caffeic acid; 3, chlorogenic acid; 4, kaempferol‐3‐glucuronide; 5: l‐tryptophan.
The volume of distribution (VDss) is a measure of the extent of drug distribution. The values of compounds 3 and 4 indicate that they will be distributed in tissue, whereas those of compounds 1 and 2 suggest that they will be distributed in plasma. With regard to blood–brain barrier (BBB) membrane permeability, compounds 3 and 4 exhibited a logBB value of less than −1, indicating that they were predicted to have difficulty crossing the BBB. In contrast, compounds 1, 2, and 5 demonstrated the ability to cross the BBB to a moderate extent. It is similarly postulated that compounds 2 and 5 will penetrate the central nervous system (CNS) with ease. The drug’s metabolic profile was predicted using cytochrome P450s, with the two main subtypes being CYP2D6 and CYP3A4 enzymes. The generated values indicated that none of the selected phytoconstituents exhibited inhibitory effects on any P450 isoform. The capacity of the body to eliminate the drug was determined by the total clearance, which encompasses both hepatic and renal processes. As illustrated in Table 6, our findings indicated that only compound 6 exhibited effective elimination. A toxicity assessment was also conducted for the selected phytoconstituents. The toxicity results were found to be promising, with none of the selected compounds exhibiting AMES, hepatotoxicity, or skin sensitization.
Target Prediction
The potential molecular targets for the major phytochemicals (chlorogenic acid, caffeic acid, malic acid, l‐tryptophan, and kaempferol‐3‐glucuronide) were predicted and the results revealed that they are categorized into several classes. As illustrated in Figure 3, the predominant target class for malic acid was enzyme (26.7%), for caffeic acid was lyase (66.7%), for chlorogenic acid was protease (33.3%), for kaempferol‐3‐glucuronide was enzyme (26.7%) and family AG protein‐coupled receptor (26.7%), and for l‐tryptophan was family AG protein‐coupled receptor (80%).

Conclusions
The extracts evaluated demonstrated promising potential as sources of antioxidant and anticancer compounds. The antioxidant activity, as demonstrated by the DPPH• and ABTS+• radical scavenging assays, exhibited significant cellular toxicity against the MCF7 (breast), HT29 (colon), and SW480 (colorectal) cell lines. This observation indicates a direct correlation between the observed biological effects and the presence of phenolic compounds. UPLC‐QTOF‐MS/MS analysis identified a variety of bioactive metabolites, including three organic acids, two phenolic acids, three flavonols, two flavones, one flavanone, five amino acids, one bile acid, and one triterpene saponin. In addition to these, other compounds were also identified and are considered coadjuvants in the observed biological activities. In silico analyses corroborated the in vitro findings, thereby reinforcing the therapeutic potential of the MVRF. However, for these results to be safely translated into clinical application, additional in vivo studies and clinical trials are needed to identify the bioactive compounds responsible for the observed effects and validate their efficacy and safety in humans.
Experimental Section
Chemical and Reagents
The HPLC‐grade methanol, formic acid (≥ 98%) and sodium hydroxide anhydrous pellets were sourced from Fisher Scientific (UK). Ammonium formate (mass spectrometer grade, 98% or greater) and HPLC‐grade acetonitrile were obtained from Sigma‐Aldrich (Germany). Ultrapure water for chromatographic purposes was prepared using a purification system (Millipore Milli‐Q, USA). The remaining chemicals and reagents were of analytical grade.
Plant Collection and Preparation of Extracts
Plant Collection
In April 2022, aerial parts of MV were collected (Figure 4). The plants were harvested from Buljurashi (Al‐Baha Province) in southern Saudi Arabia (coordinates 41.58′ N and 19.86′ E, altitude 2450 m) at the flowering stage. The plant material was authenticated by an expert taxonomist, Dr. Haidar A. Mohamed of the Department of Biology, Faculty of Science, Al‐Baha University. The available databases were consulted for the purpose of verifying the nomenclature [ref. 144, ref. 145]. The plant was archived with voucher specimen herbarium number BUH‐112 at the Pharmacognosy Department, Faculty of Clinical Pharmacy, Al‐Baha University, KSA.

Preparation of Extracts
The aerial parts of MV (500 g) were dried in a shaded environment, ground, and then successively extracted at room temperature with ethanol and water (8:2, 3 × 1.5 L), yielding 14% of the MVCE. The residue was subjected to fractionation, resulting in the isolation of five distinct fractions: petroleum ether (MVPF, 4.1 g), chloroform (MVCF, 12.9 g), ethyl acetate (MVEF, 1.3 g), and an aqueous residual fraction (MVRF, 10.9 g).
Phytochemical Screening, TPC, and TFC
A preliminary qualitative phytochemical screening was conducted to identify the phytochemical constituents of the crude extract of MV aerial parts using standard protocols [ref. 146]. The following phytochemicals were subjected to analysis: triterpenes, alkaloids, saponins, cardiac glycosides, steroids, flavonoids, tannins, and phenolic compounds.
The TPC of the crude extract (MVCE) was quantified in accordance with the methodology established by Pourmorad et al. [ref. 147] and expressed as mg GAE/g DW.
The TFC was also determined based on the method described by Pourmorad et al. [ref. 147], with some modifications, and was expressed in mg RE/g DW.
Antioxidant Activity
DPPH• and ABTS+• Radical Scavenging Activities
The method for determining antioxidant activity based on the free radical DPPH• (2,2‐diphenyl‐1‐picrylhydrazyl) was performed according to procedures previously described by [ref. 148]. The technique involves reducing the DPPH• radical in a methanolic solution in the presence of antioxidant compounds, which is indicated by a decrease in absorbance at 517 nm. Briefly, a DPPH• solution (2 mL, 0.1 mM) was prepared in methanol and mixed with 1 mL of various concentrations of MV extracts or the antioxidant standard compound (Trolox). The mixture was incubated at room temperature, in the dark, for 1 h. After this period, the absorbance was measured in a spectrophotometer. The percentage of inhibition of the DPPH• radical was calculated using the following Equation (1).
The results were expressed in terms of IC50 (concentration required to inhibit 50% of DPPH•), calculated by linear fitting of the tested concentrations.
The method for determining antioxidant activity based on the ABTS•+ radical cation discoloration assay was performed according to methods described in the literature [ref. 149]. Briefly, the ABTS•+ cation radical was generated by the reaction of an aqueous solution of ABTS (7 mM) with potassium persulfate (2.45 mM) and the mixture was kept in the dark at room temperature for 12–16 h before use. The formed ABTS•+ radical cation was diluted in ethanol to an absorbance of 0.700 ± 0.020 at 734 nm. The antioxidant activity was assessed by adding 10 µL of the sample or antioxidant standard to 1 mL of the ABTS•+ solution. The absorbance measurement was then performed after 30 min at 30°C. Appropriate solvents were used as controls in each assay. All determinations were carried out in triplicate. The percentage inhibition of the ABTS•+ radical was calculated using the following Equation (1). The results were expressed as Trolox‐equivalent antioxidant activity (TEAC), calculated in relation to a standard curve constructed with different concentrations of Trolox and IC50.
Cytotoxicity
Cell Lines
This study employed three cancer cell lines: MCF7 (human breast adenocarcinoma), HT29, and SW480 (human colorectal adenocarcinoma). Furthermore, MRC5 (normal human fetal lung fibroblast) cells were procured from the ATCC (USA). The three cancer cells were subcultured in RPMI‐1640 media in accordance with the methodology described by Abdalla et al. [ref. 150].
The Cytotoxicity Determination
The cytotoxicity of the six extracts, along with DOX, was evaluated using the MTT assay, as previously described [ref. 151].
Phytochemical Profile by UPLC‐ESI‐QTOF‐MS/MS Analysis
The UPLC‐ESI‐QTOF‐MS/MS analysis was conducted according to Mohammed et al. [ref. 72].
Sample Preparation
A stock solution of the extract was created using 50 mg of the lyophilized aqueous ethanolic extract, which was dissolved in 1000 µL of a solvent blend comprising water, methanol, and acetonitrile (H2O:MeOH:ACN) at a ratio of 2:1:1. The stock solution was rendered fully soluble by vortexing the sample and subjecting it to ultrasonication at 30 kHz for a period of 10 min. A 20 µL aliquot of the stock solution was diluted once more with 1000 µL of H2O:MeOH:ACN (2:1:1) and subjected to centrifugation at 10 000 rpm for 5 min. The resulting solution was then injected for analysis at a concentration of 1 µg/mL. Similarly, an LC‐MS analysis was conducted for the blank and quality control samples/internal standard (IS) to ensure the reliability of the experiment. The sample was introduced in both positive and negative modes.
Instruments and Acquisition Method
The isolation of small molecules was conducted using an ExionLC system (AB Sciex, Framingham, MA, USA), which was linked to an autosampler, an in‐line filter disk pre‐column (0.5 µm × 3.0 mm, Phenomenex, Torrance, CA, USA), and an Xbridge C18 column (3.5 µm, 2.1 × 50 mm) (Waters Corporation, Milford, MA, USA). The column was maintained at 40°C and operated at a flow rate of 300 µL/min. The mobile phase consisted of Solution A, which was 5 mM ammonium formate in 1% methanol, with the pH adjusted to 3.0 using formic acid, and Solution B, which was acetonitrile (100%) for the positive mode. The negative mode solution (C) consisted of 5 mM ammonium formate in 1% methanol, modified to pH 8 with sodium hydroxide. The gradient elution was conducted according to the following schedule: The initial phase of the elution process was conducted for 20 min at a concentration of 10% B, followed by a 5‐min period at 90% B. Subsequently, the concentration was reduced to 10% B for a further 5 min, after which the process was terminated, and the column was allowed to undergo a period of equilibration at 90% B.
The MS analysis was conducted using a Triple TOF 5600+ system, which is equipped with a Duo‐Spray source operating in ESI mode (AB SCIEX, Concord, Ontario, Canada). The sprayer capillary and declustering potential voltages were set to 4500 and 80 eV in positive mode and −4500 and −80 V in negative mode. The source temperature was set at 600°C, the curtain gas at 25 psi, and both Gas 1 and Gas 2 at 40 psi. The collision energy was set at 35 V (positive mode) and −35 V (negative mode), with a CE spread of 20 V and an ion tolerance of 10 ppm. The TripleTOF 5600+ was operated with an information‐dependent acquisition (IDA) protocol. Batches for the collection of MS and MS/MS data were generated with the use of Analyst‐TF 1.7.1. The IDA technique was employed to simultaneously gather full‐scan MS and MS/MS data. The methodology entailed the acquisition of high‐resolution survey spectra spanning a mass range of 50–1100 m/z, with the mass spectrometer employed in the aforementioned configuration [ref. 152].
LC‐MS Data Processing
The open‐source software MS‐DIAL 3.70 [ref. 153] was employed for a comprehensive analysis of small molecules in the absence of any targeted sampling. The ReSpect positive database (2737 records) and the ReSpect negative database (1573 records) were selected as reference databases based on the acquisition method. The search settings were established with MS1 and MS2 mass tolerance at 0.01 and 0.05 Da for data acquisition, with a minimum peak height of 100 amplitude, a mass slice width of 0.05 Da, a smoothing level of two scans, and a minimum peak width of six scans. For identification, the MS1 and MS2 tolerances were set as follows: The alignment tolerance was set at 0.2 Da, while the retention time tolerance was 0.05 min, and the MS1 tolerance was set at 0.25 Da. The MS‐DIAL results were employed to reprocess the data in PeakView 2.2 with the MasterView 1.1 package (AB SCIEX) for the purpose of confirming features (peaks) from the TIC in accordance with the following criteria: aligned features with a signal‐to‐noise ratio exceeding 5 and sample intensities relative to the blank exceeding 5.
In Silico Molecular Docking
The structures of the phytocompounds were obtained from PubChem and processed using Schrödinger’s LigPrep, which optimizes ligand geometry and generates low‐energy 3D structures with accurate chirality. Molecular docking [ref. 154] of the phytocompounds was performed using Schrödinger’s Glide into the active sites of ERα (PDB ID: 6CHZ) and human TS (PDB ID: 1HVY). The protein crystal structures were preprocessed as follows [ref. 155, ref. 156]. Finally, docking was conducted in Glide using the standard precision (SP) method [ref. 157].
Pharmacokinetics Study
In silico molecular studies of the selected phytoconstituents were investigated to predict their pharmacokinetic (https://biosig.lab.uq.edu.au/pkcsm/) and target prediction (http://www.swisstargetprediction.ch/) characteristics.
Statistical Analysis
All assays were conducted in triplicate, and the data was expressed as mean ± standard deviation (SD). Differences between the data were analyzed using Tukey’s one‐way analysis of variance (ANOVA) with SPSS software (version 22.0).
Author Contributions
Ines El Mannoubi: led the conception and design of the study, supervised the execution of the research, and coordinated the writing of the manuscript. Mozaniel Santana de Oliveira: contributed significantly to the methodological structuring and data analysis. Nuha M. Alghamdi: Contributed to data collection and analysis. Seham H. Bashir: contributed to data collection and analysis. Suada Alsaied Mohamed: contributed to data collection and analysis. Hedia Chaabane: participated in the interpretation of results and critical review of the manuscript. Ashraf N. Abdalla: participated in the interpretation of results and critical review of the manuscript. Majdi Abid: participated in the interpretation of results and critical review of the manuscript. Adel Kadri: coordinated the writing of the manuscript, provided technical support and academic supervision.
Conflicts of Interest
The authors declare no conflicts of interest.
References
- In Vitro and In Silico Evaluation of Antiproliferative Activity of New Isoxazolidine Derivatives Targeting EGFR: Design, Synthesis, Cell Cycle Analysis, and Apoptotic Inducers,”. Pharmaceuticals, 2023. [PubMed]
- The Rise of Gastrointestinal Cancers as a Global Phenomenon: Unhealthy Behavior or Progress?,”. International Journal of Environmental Research and Public Health, 2023. [PubMed]
- The Current and Future Cancer Burden in the Gulf Cooperation Council (GCC) Countries,”. Cancer Med, 2024. [DOI | PubMed]
- Drug Repurposing for Cancer Treatment: A Comprehensive Review,”. International Journal of Molecular Sciences, 2024. [PubMed]
- Demographic and Clinical Characteristics Associated With Advanced Stage Colorectal Cancer: A Registry‐Based Cohort Study in Saudi Arabia,”. BMC Cancer, 2024. [PubMed]
- Reactive Oxygen Species‐Upregulating Nanomedicines Towards Enhanced Cancer Therapy,”. Biomaterials Science, 2023. [PubMed]
- Promoting Reactive Oxygen Species Accumulation to Overcome Tyrosine Kinase Inhibitor Resistance in Cancer,”. Cancer Cell International, 2024. [PubMed]
- Oxidative Cell Death in Cancer: Mechanisms and Therapeutic Opportunities,”. Cell Death & Disease, 2024. [PubMed]
- Cancer Treatment and Survivorship Statistics, 2022,”. CA: A Cancer Journal for Clinicians, 2022. [PubMed]
- Cancer Chemotherapy and Beyond: Current Status, Drug Candidates, Associated Risks and Progress in Targeted Therapeutics,”. Genes & Diseases, 2023. [PubMed]
- A Review on the Efficacy and Toxicity of Different Doxorubicin Nanoparticles for Targeted Therapy in Metastatic Breast Cancer,”. Biomedicine & Pharmacotherapy, 2017. [PubMed]
- Synergy Between PEDF and Doxorubicin in Breast Cancer Cells: Effects on Metastatic and Metabolic Pathways,”. International Journal of Molecular Sciences, 2024. [PubMed]
- The Chemotherapeutic Potentials of Compounds Isolated From the Plant, Marine, Fungus, and Microorganism: Their Mechanism of Action and Prospects,”. Journal of Tropical Medicine, 2022. [PubMed]
- Phytochemical Profiling, Antimicrobial, Antibiofilm, Insecticidal, and Anti‐Leishmanial Properties of Aqueous Extract From Juglans regia L. Root Bark: In Vitro and In Silico Approaches,”. International Journal of Food Properties, 2023
- Anticancer Drug Discovery Based on Natural Products: From Computational Approaches to Clinical Studies,”. Biomedicines, 2024. [PubMed]
- Calocedrus formosana Essential Oils Induce ROS‐Mediated Autophagy and Apoptosis by Targeting SIRT1 in Colon Cancer Cells,”. Antioxidants, 2024. [PubMed]
- The In‐Vitro Evaluation of Antibacterial, Antifungal and Cytotoxic Properties of Marrubium vulgare L. essential Oil Grown in Tunisia,”. Lipids in Health and Disease, 2011. [PubMed]
- Chemical Composition and Antioxidant Activity of Marrubium vulgare L. Essential Oil From Tunisia,”. African Journal of Biotechnology, 2011
- Essential Oil of Marrubium vulgare: Chemical Composition and Biological Activities. A Review,”. Natural Product Sciences, 2019
- The Bioactive Ingredients and Therapeutic Effects of Marrubium vulgare—A Review,”. International Journal of Biological and Pharmaceutical Sciences Archive, 2021
- The Vasorelaxant Activity of Marrubenol and Marrubiin From Marrubium vulgare ,”. Planta Medica, 2003. [PubMed]
- Marrubium vulgare L.: A Phytochemical and Pharmacological Overview,”. Molecules, 2020. [PubMed]
- A Review of the Chemical and Pharmacological Aspects of the Genus Marrubium ,”. Current Pharmaceutical Design, 2010. [PubMed]
- Marrubium vulgare L. Leave Extract: Phytochemical Composition, Antioxidant and Wound Healing Properties,”. Molecules, 2017. [PubMed]
- Comparative Phytochemical Investigation of the Composition and Content of Biologically Active Substances in Marrubium vulgare and M. alternidens ,”. Chemistry of Natural Compounds, 2003
- Lead Finding From Whole Plant of Marrubium vulgare L. With Hepatoprotective Potentials Through In Silico Methods,”. Asian Pacific Journal of Tropical Biomedicine, 2012
- Evaluation of Antioxidant Activity of Alkaloid Extracts Obtained From Marrubium vulgare L.,”. Current Developments in Nutrition, 2022
- Antiproliferative Activity of Marrubium persicum Extract in the MCF‐7 Human Breast Cancer Cell Line,”. Asian Pacific Journal of Cancer Prevention, 2012. [PubMed]
- Marrubium vulgare L.: A Review on Phytochemical and Pharmacological Aspects,”. Journal of Intercultural Ethnopharmacology, 2017
- A Comprehensive Assessment of the Chemical Composition, Antioxidant, Genoprotective and Antigenotoxic Activities of Lamiaceae Species Using Different Experimental Models In Vitro,”. Food & Function, 2021. [PubMed]
- Biological Activity of Horehound (Marrubium vulgare L.) Herb Grown in Poland and Its Phytochemical Composition,”. Pharmaceuticals, 2024. [PubMed]
- Anti‐Alzheimer’s Potency of Rich Phenylethanoid Glycosides Extract From Marrubium vulgare L.: In Vitro and In Silico Studies,”. Pharmaceuticals, 2024. [PubMed]
- Tunisian Horehound (Marrubium vulgare) Aqueous Extract Improves Treatment of Bleomycin‐Induced Lung Fibrosis in Rat,”. Dose‐Response, 2022. [DOI]
- The Application of Ultrasonic Waves and Microwaves to Improve Antihyperglycaemic and Antimicrobial Activities of Marrubium vulgare Extracts,”. Antibiotics, 2022. [PubMed]
- Phytochemical Screening, Total Phenolic Content and Antioxidant Activities of Marrubium vulgare L. Collected in Three Development Stages,”. Notulae Scientia Biologicae, 2023
- HPLC‐Analysis, Biological Activities and Characterization of Action Mode of Saudi Marrubium vulgare Against Foodborne Diseases Bacteria,”. Molecules, 2021. [PubMed]
- Phytochemical Screening, Polyphenols, Flavonoids and Tannin Content, Antioxidant Activities and FTIR Characterization of Marrubium vulgare L. From 2 Different Localities of Northeast of Morocco,”. Heliyon, 2020. [PubMed]
- Chemical Composition and In Vitro Antioxidant and Antimicrobial Activities of Marrubium vulgare L. ,”. Scientific World Journal, 2021. [PubMed]
- Variation of the Chemical Composition of Essential Oils and Total Phenols Content in Natural Populations of Marrubium vulgare L. ,”. Plants, 2022. [PubMed]
- Phytochemical Composition, Antioxidant and Antibacterial Activities of Crude Extracts From the Species Euphorbia atlantica Coss. ,”. Pharmaceutical Chemistry Journal, 2019
- The Antimicrobial, Antioxidative, Anti‐Inflammatory Activity and Cytotoxicity of Different Fractions of Four South African Bauhinia Species Used Traditionally to Treat Diarrhoea,”. Journal of Ethnopharmacology, 2012. [PubMed]
- Isolation and Evaluation of the Antioxidant Capacity of Compounds From Ehretia asperula Zoll. & Moritzi,”. Indonesian Journal of Chemistry, 2024
- Phytochemical Profiling, In Vitro Evaluation of Total Phenolic Contents and Antioxidant Properties of Marrubium vulgare (Horehound) Leaves of Plants Growing in Algeria,”. Industrial Crops and Products, 2014
- High‐Throughput Screening and Characterization of Phenolic Compounds in Stone Fruits Waste by LC‐ESI‐QTOF‐MS/MS and Their Potential Antioxidant Activities,”. Antioxidants, 2021. [PubMed]
- Radical Scavenging Mechanisms of Phenolic Compounds: A Quantitative Structure‐Property Relationship (QSPR) Study,”. Frontiers in Nutrition, 2022. [DOI | PubMed]
- Common Trends and Differences in Antioxidant Activity Analysis of Phenolic Substances Using Single Electron Transfer Based Assays,”. Molecules, 2021. [PubMed]
- The Effects of Resveratrol‐Rich Extracts of Vitis vinifera Pruning Waste on HeLa, MCF‐7 and MRC‐5 Cells: Apoptosis, Autophagia and Necrosis Interplay,”. Pharmaceutics, 2022. [PubMed]
- Qualitative HPLC‐DAD/ESI‐TOF‐MS Analysis, Cytotoxic, and Apoptotic Effects of Croatian Endemic Centaurea ragusina L. Aqueous Extracts,”. Chemistry and Biodiversity, 2017. [DOI]
- Metabolic Profiling of Vitex pubescens Vahl Bark via UPLC‐ESI‐QTOF/MS/MS Analysis and Evaluation of Its Antioxidant and Acetylcholinesterase Inhibitory Activities,”. BMC Complementary Medicine and Therapies, 2024. [PubMed]
- Chemical Profiling and Quantitation of Bioactive Compounds in Platycladi Cacumen by UPLC‐Q‐TOF‐MS/MS and UPLC‐DAD,”. Journal of Pharmaceutical and Biomedical Analysis, 2018. [PubMed]
- Extraction and Characterization of Flavonoid Kaempferol‐3‐Glucuronide from Mango for Treatment of Different Types of Cancers,”. London Journal of Research in Science: Natural, 2024
- Assessment of the Antioxidant and Enzyme Inhibition Activities of Cousinia Iconica With Focus on Phytochemical Investigation by LC‐MS/MS,”. Plant Biosystems, 2021
- Qualitative Analysis and Componential Differences of Chemical Constituents in Taxilli herba From Different Hosts by UFLC‐Triple TOF‐MS/MS,”. Molecules, 2021. [PubMed]
- In Vitro Bioactivities of Commonly Consumed Cereal, Vegetable, and Legume Seeds as Related to Their Bioactive Components: An Untargeted Metabolomics Approach Using UHPLC–QTOF‐MS2,”. Antioxidants, 2023. [PubMed]
- Chemical Characterization and Quantitative Determination of Flavonoids and Phenolic Acids in Yerba Santa (Eriodictyon spp.) Using UHPLC/DAD/Q‐ToF,”. Journal of Pharmaceutical and Biomedical Analysis, 2023. [PubMed]
- LC‐ESI‐QTOF‐MS2 Characterization of Phenolic Compounds in Different Lentil (Lens culinaris M.) Samples and Their Antioxidant Capacity,”. Front Biosci, 2023. [DOI]
- UPLC‐ESI‐MS/MS Profiling of Secondary Metabolites From Methanol Extracts of In Vivo and In Vitro Tissues of Daucus capillifolius Gilli (A Comparative Study),”. Molecules, 2024. [PubMed]
- Main Iridoid Glycosides and HPLC/DAD‐Q‐TOF‐MS/MS Profile of Glycosides From the Antioxidant Extract of Eucommia Ulmoides Oliver Seeds,”. Industrial Crops and Products, 2016
- Analysis of a Series of Chlorogenic Acid Isomers Using Differential Ion Mobility and Tandem Mass Spectrometry,”. Analytica Chimica Acta, 2016. [PubMed]
- Quantitative and Qualitative Analysis of Flavonoids and Phenolic Acids in Snow Chrysanthemum (Coreopsis tinctoria Nutt.) by HPLC‐DAD and UPLC‐ESI‐QTOF‐MS,”. Molecules, 2016. [PubMed]
- Quality Evaluation of Ilex asprella Based on Simultaneous Determination of Five Bioactive Components, Chlorogenic Acid, Luteoloside, Quercitrin, Quercetin, and Kaempferol, Using UPLC–Q‐TOF MS Study,”. Journal of AOAC International, 2019. [PubMed]
- TLC‐MALDI‐TOF‐MS‐Based Identification of Flavonoid Compounds Using an Inorganic Matrix,”. Phytochemical Analysis, 2019. [PubMed]
- Identification of Isomeric Flavonoid Glucuronides in Urine and Plasma by Metal Complexation and LC‐ESI‐MS/MS,”. Journal of Mass Spectrometry, 2006. [PubMed]
- LC‐DAD‐ESI‐MS/MS Analysis and Cytotoxic and Antiproliferative Effects of Chlorogenic Acid Derivative Rich Extract From Nerium oleander L. pink Flowers,”. Food & Function, 2021. [PubMed]
- LC/MS/MS and GC/MS/MS Metabolic Profiling of Leontodon hispidulus, In Vitro and in Silico Anticancer Activity Evaluation Targeting Hexokinase 2 Enzyme,”. Scientific Reports, 2024. [PubMed]
- Phenolic Characterization of a Purple Maize (Zea mays cv. “Moragro”) by HPLC–QTOF‐MS and Study of Its Bioaccessibility Using a Simulated In Vitro Digestion/Caco‐2 Culture Model,”. Journal of Agricultural and Food Chemistry, 2024. [PubMed]
- A Red‐Berry Mixture as a Nutraceutical: Detailed Composition and Neuronal Protective Effect,”. Molecules, 2021. [PubMed]
- Comprehensive Characterization of C‐Glycosyl Flavones in Wheat (Triticum aestivum L.) Germ Using UPLC‐PDA‐ESI/HRMSn and Mass Defect Filtering,”. Journal of Mass Spectrometry, 2016. [PubMed]
- Detection of Flavone C‐Glycosides in the Extracts From the Bark of Prunus avium L. and Prunus cerasus L. ,”. European Journal of Mass Spectrometry, 2020. [PubMed]
- Comparative Study of the Phytochemical Profiles of the Rhizomes of Cultivated and Wild‐Grown Polygonatum sibiricum ,”. Separations, 2022
- Characterization of Potentially Health‐Promoting Constituents in Sea Fennel (Crithmum maritimum) Cultivated in the Conero Natural Park (Marche Region, Central Italy),”. Pharmaceutical Biology, 2023. [PubMed]
- Phytochemical Profiling, In Vitro and In Silico Anti‐Microbial and Anti‐Cancer Activity Evaluations and Staph GyraseB and h‐TOP‐IIβ Receptor‐Docking Studies of Major Constituents of Zygophyllum coccineum L. Aqueous‐Ethanolic Extract and Its Subsequent Fractions: An Approach to Validate Traditional Phytomedicinal Knowledge,”. Molecules, 2021. [PubMed]
- Qualitative and Quantitative Analysis of Phytochemicals in Sayeok‐Tang via UPLC‐Q‐Orbitrap‐MS and UPLC‐TQ‐MS/MS,”. Pharmaceuticals, 2024. [PubMed]
- LC‐MS/MS‐Based Chemical Profiling of Water Extracts of Moringa Oleifera Leaves and Pharmacokinetics of Their Major Constituents in Rat Plasma,”. Food Chemistry: X, 2024. [PubMed]
- Identification and Mass Spectrometric Characterization of Glycosylated Flavonoids in Triticum Durum Plants by High‐Performance Liquid Chromatography With Tandem Mass Spectrometry,”. Rapid Communications in Mass Spectrometry, 2005. [PubMed]
- Approach to the Study of Flavone Di‐C‐Glycosides by High Performance Liquid Chromatography‐Tandem Ion Trap Mass Spectrometry and Its Application to Characterization of Flavonoid Composition in Viola yedoensis ,”. Journal of Mass Spectrometry, 2014. [PubMed]
- Phytochemical Elucidation and Effect of Maesa indica (Roxb.) Sweet on Alleviation of Potassium Dichromate‐Induced Pulmonary Damage in Rats,”. Plants, 2024. [PubMed]
- UPLC‐ESI/MS/MS Profiling and Anti‐Inflammatory Activity of Gleditsia caspica,”. Archives of Pharmaceutical Sciences Ain Shams University, 2020
- Metabolomic Workflow for the Accurate and High‐Throughput Exploration of the Pathways of Tryptophan, Tyrosine, Phenylalanine, and Branched‐Chain Amino Acids in Human Biofluids,”. Journal of Proteome Research, 2022. [PubMed]
- Linking MS1 and MS2 Signals in Positive and Negative Modes of LC‐HRMS in Untargeted Metabolomics Using the ROIMCR Approach,”. Analytical and Bioanalytical Chemistry, 2023. [PubMed]
- Chemometric Analysis and Chemical Characterization for the Botanical Identification of Glycyrrhiza Species (G. glabra, G. uralensis, G. inflata, G. echinata and G. lepidota) Using Liquid Chromatography‐Quadrupole Time of Flight Mass Spectrometry (LC‐QToF),”. Journal of Food Composition and Analysis, 2022
- LC‐MS/MS Standardization and Validation of Glycyrrhizin From the Roots of Taverniera cuneifolia: A Potential Alternative Source of Glycyrrhiza glabra ,”. Heliyon, 2022. [PubMed]
- Phospholipid Synthesis inside Phospholipid Membrane Vesicles,”. Communications Biology, 2022. [PubMed]
- Nicotinamide Deficiency in Primary Open‐Angle Glaucoma,”. Investigative Ophthalmology & Visual Science, 2019. [PubMed]
- Association of Serum Levels of Deoxyribose 1–Phosphate and S–lactoylglutathione With Neoadjuvant Chemotherapy Sensitivity in Patients With Gastric Cancer: A Metabolomics Study,”. Oncology letters, 2020. [DOI | PubMed]
- Metabolic Disturbance in Hippocampus and Liver of Mice: A Primary Response to Imidacloprid Exposure,”. Scientific Reports, 2020. [PubMed]
- Phytochemical Analysis of Ononis arvensis L. by Liquid Chromatography Coupled With Mass Spectrometry,”. Journal of Mass Spectrometry, 2019. [PubMed]
- Marrubium alysson L. Ameliorated Methotrexate‐Induced Testicular Damage in Mice Through Regulation of Apoptosis and miRNA‐29a Expression: LC‐MS/MS Metabolic Profiling,”. Plants, 2022. [PubMed]
- Metabolomic Profiling of Burkholderia pseudomallei Using UHPLC‐ESI‐Q‐TOF‐MS Reveals Specific Biomarkers Including 4‐Methyl‐5‐Thiazoleethanol and Unique Thiamine Degradation Pathway,”. Cell & Bioscience, 2015. [PubMed]
- Evaluation of Antioxidant Capacity of Endemic Plant Marrubium astracanicum Subsp. macrodon: Identification of Its Phenolic Contents by Using HPLC‐MS/MS,”. Natural Product Research, 2019. [PubMed]
- Identification of the Free Phenolic Profile of Adlay Bran by UPLC‐QTOF‐MS/MS and Inhibitory Mechanisms of Phenolic Acids Against Xanthine Oxidase,”. Food Chemistry, 2018. [PubMed]
- Isolation, Identification and Activity of Natural Antioxidants From Horehound (Marrubium vulgare L.) Cultivated in Lithuania,”. Food Chemistry, 2012
- LC‐TOFMS Analysis, In Vitro and In Silico Antioxidant Activity on NADPH Oxidase, and Toxicity Assessment of an Extract Mixture Based on Marrubium vulgare L. and Dittrichia viscosa L.,”. Journal Biology Biomedical Research, 2024
- Phytochemical Composition, Biological Activity and Molecular Docking Studies of the Endemic Marrubium trachyticum Boiss. ,”. Türkiye Tarımsal Araştırmalar Dergisi, 2024
- Metabolome Profiling of Marrubium peregrinum L. and Marrubium friwaldskyanum Boiss Reveals Their Potential as Sources of Plant‐Based Pharmaceuticals,”. International Journal of Molecular Sciences, 2023. [PubMed]
- Marrubium Vulgare Ethanolic Extract Induces Proliferation Block, Apoptosis, and Cytoprotective Autophagy in Cancer Cells In Vitro,”. Cell. Mol. Biol., 2016. [PubMed]
- The Phytochemical and Antimicrobial Examination of Phenolic Acids Contained in Extracts From Four Marrubium Species,”. Acta Pol. Pharm ‐ Drug Res, 2018
- Phytochemical Analysis and In Vitro Anti‐Inflammatory, Anticancer Activities of Marrubium lutescens on Melanoma Cancer Cell Line and Molecular Docking Studies,”. Journal of Herbal Medicine, 2024
- Antioxidant Properties of Hydroxycinnamic Acids: A Review of Structure‐Activity Relationships,”. Current Medicinal Chemistry, 2013. [PubMed]
- The Study of Anti‐/Pro‐Oxidant, Lipophilic, Microbial and Spectroscopic Properties of New Alkali Metal Salts of 5‐O‐Caffeoylquinic Acid,”. International Journal of Molecular Sciences, 2018. [PubMed]
- Caffeic Acid on Metabolic Syndrome: A Review,”. Molecules, 2021. [PubMed]
- Antioxidant Activity of Caffeic Acid: Thermodynamic and Kinetic Aspects on the Oxidative Degradation Pathway,”. Free Radical Research, 2022. [PubMed]
- Therapeutic Implications of Caffeic Acid in Cancer and Neurological Diseases,”. Front. Oncol., 2022. [DOI]
- Identification of Flavonoids from Marrubium and Ballota Species (Lamiaceae) and Determination of Chemotaxonomic Markers Using High Performance Liquid Chromatography Mass Spectrometer,”. {Journal of Sciences, Islamic Republic of Iran, 2022
- Chemical Review on Genus Marrubium,”. Records of Pharmaceutical and Biomedical Sciences, 2023
- Chemistry, Pharmacokinetics, Pharmacological Activities, and Toxicity of Quercitrin,”. Phytotherapy Research, 2022. [PubMed]
- Superoxide Scavenging by Polyphenols: Effect of Conjugation and Dimerization,”. Redox Repor, 2002
- Biological Activities of Flavonoids From Nitraria retusa (Forssk.) Asch. and Their Acylated Derivatives,”. Food Chemistry, 2011
- Structures, Sources, Identification/Quantification Methods, Health Benefits, Bioaccessibility, and Products of Isorhamnetin Glycosides as Phytonutrients,”. Nutrients, 2023. [PubMed]
- UPLC‐ESI‐MS polyphenolic profiling and the antimicrobial potential of Ailanthus altissima (Mill.) Swingle leaves growing in Egypt”. Journal of Research in Pharmacy, 2024
- Modulatory Effect of Polyphenolic Compounds From the Mangrove Tree Rhizophora mangle L. on Non‐Alcoholic Fatty Liver Disease and Insulin Resistance in High‐Fat Diet Obese Mice,”. Molecules, 2018. [PubMed]
- New Features on the Fragmentation and Differentiation of C‐Glycosidic Flavone Isomers by Positive Electrospray Ionization and Triple Quadrupole Mass Spectrometry,”. Rapid Communications in Mass Spectrometry, 2008. [PubMed]
- Compositional Analysis and In Vivo Anti‐Diabetic Activity of Wild Algerian Marrubium vulgare L. Infusion,”. Fitoterapia, 2012. [PubMed]
- Antioxidant and Antifungal Activities of Marrubiin, Extracts and Essential Oil From Marrubium vulgare L. Against Pathogenic Dermatophyte Strains,”. Journal de Mycologie Médicale, 2020. [PubMed]
- Vicenin‐2 Inhibits Wnt/β‐catenin Signaling and Induces Apoptosis in HT‐29 Human Colon Cancer Cell Line,”. Drug Design, Development and Therapy, 2018. [PubMed]
- Research Progress on Pharmacological Effects and New Dosage Forms of Baicalin,”. Veterinary Medicine and Science, 2022. [PubMed]
- In‐Vitro Cytotoxicity, Apoptotic Property, and Gene Expression Changes Induced by Naringenin‐7‐O‐Glucoside in Triple‐Negative Breast Cancer,”. Cureus, 2024. [DOI | PubMed]
- Effects of Dietary Malic Acid Supplementation on Growth Performance, Antioxidant and Immunological Parameters, and Intestinal Gene Expressions in Rainbow Trout, Oncorhynchus mykiss ,”. Aquaculture, 2023
- Bio‐Fermented Malic Acid Facilitates the Production of High‐Quality Chicken via Enhancing Muscle Antioxidant Capacity of Broilers,”. Antioxidants, 2022. [PubMed]
- Caffeic Acid: An Antioxidant With Novel Antisickling Properties,”. FEBS Open Bio, 2021
- In Vitro and In Vivo Antioxidant Properties of Chlorogenic Acid and Caffeic Acid,”. International Journal of Pharmaceutics, 2011. [PubMed]
- Antioxidant Potential of Chlorogenic Acid in Age‐Related Eye Diseases,”. Pharmacology Research & Perspectives, 2024. [DOI | PubMed]
- Kaempferol‐3‐O‐Glucuronide Ameliorates Non‐Alcoholic Steatohepatitis in High‐Cholesterol‐Diet‐Induced Larval Zebrafish and HepG2 Cell Models via Regulating Oxidation Stress,”. Life, 2021. [PubMed]
- Evaluation of the Antioxidant Properties of Tryptophan and Its Metabolites in In Vitro Assay,”. Journal of Complementary and Integrative Medicine, 2016. [PubMed]
- Anticarcinogenic Properties of Malic Acid on Glioblastoma Cell Line Through Necrotic Cell Death Mechanism,”. MANAS Journal of Engineering, 2021
- Chemopreventive Effects of Dietary Phytochemicals Against Cancer Invasion and Metastasis: Phenolic Acids, Monophenol, Polyphenol, and Their Derivatives,”. Cancer Treatment Reviews, 2012. [PubMed]
- 127 M. Alam , G. M. Ashraf , K. Sheikh , et al., "Potential Therapeutic Implications of Caffeic Acid in Cancer Signaling: Past, Present, and Future", Frontiers in Pharmacology 13 (2022): 845871, 10.3389/fphar.2022.845871.35355732
- The Effects of Caffeic, Coumaric and Ferulic Acids on Proliferation, Superoxide Production, Adhesion and Migration of Human Tumor Cells In Vitro,”. European Journal of Pharmacology, 2015. [PubMed]
- Caffeic Acid Expands Anti‐Tumor Effect of Metformin in Human Metastatic Cervical Carcinoma HTB‐34 Cells: Implications of AMPK Activation and Impairment of Fatty Acids De Novo Biosynthesis,”. International Journal of Molecular Sciences, 2017. [PubMed]
- Inhibitory Effect of Caffeic Acid on Cancer Cell Proliferation by Oxidative Mechanism in Human HT‐1080 Fibrosarcoma Cell Line,”. Molecular and Cellular Biochemistry, 2011. [PubMed]
- Chlorogenic Acid for Cancer Prevention and Therapy: Current Status on Efficacy and Mechanisms of Action,”. Pharmacological Research, 2022. [PubMed]
- Chlorogenic Acid Induces Apoptosis and Cell‐Cycle Arrest in Colorectal Cancer Cells,”. Molecular Biology Reports, 2023. [PubMed]
- Inhibitory Effects of Kaempferol on the Invasion of Human Breast Carcinoma Cells by Downregulating the Expression and Activity of Matrix Metalloproteinase‐9,”. Biochemistry and Cell Biology, 2015. [PubMed]
- Estrogen Receptors Alpha (ERα) and Beta (ERβ): Subtype‐selective Ligands and Clinical Potential,”. Steroids, 2014. [PubMed]
- 135 N. Fuentes and P. Silveyra , “Estrogen Receptor Signaling Mechanisms,” in Advances in Protein Chemistry and Structural Biology, edited by R. Donev , 116 (2019): 135–170. Elsevier.
- Structural Analyses of human Thymidylate Synthase Reveal a Site That May Control Conformational Switching Between Active and Inactive States,”. Journal of Biological Chemistry, 2017. [PubMed]
- Discovery of Selective Estrogen Receptor Covalent Antagonists for the Treatment of ERαWT and ERαMUT Breast Cancer,”. Cancer Discovery, 2018. [PubMed]
- Human Thymidylate Synthase Is in the Closed Conformation When Complexed With dUMP and Raltitrexed, an Antifolate Drug,”. Biochemistry, 2001. [PubMed]
- Design, Synthesis, Molecular Properties and In Vitro Antioxidant and Antibacterial Potential of Novel Enantiopure Isoxazolidine Derivatives,”. Arabian Journal of Chemistry, 2020
- Novel Enantiopure Isoxazolidine and C‐Alkyl Imine Oxide Derivatives as Potential Hypoglycemic Agents: Design, Synthesis, Dual Inhibitors of α‐Amylase and α‐Glucosidase, ADMET and Molecular Docking Study,”. Bioorganic Chemistry, 2020. [PubMed]
- Synthesis, Structure‐Activity Relationship and in Silico Studies of Novel Pyrazolothiazole and Thiazolopyridine Derivatives as Prospective Antimicrobial and Anticancer Agents,”. ChemistrySelect, 2021
- Design, Synthesis, Biological Evaluation and In Silico Studies of Novel 1,2,3‐Triazole Linked Benzoxazine‐2,4‐Dione Conjugates as Potent Antimicrobial, Antioxidant and Anti‐Inflammatory Agents,”. Arabian Journal of Chemistry, 2022
- Identification of Dual‐Target Isoxazolidine‐Isatin Hybrids With Antidiabetic Potential: Design, Synthesis, In Vitro and Multiscale Molecular Modeling Approaches,”. Heliyon, 2024. [PubMed]
- “R Bot Gard Kew, Harvard Univ Herb Libr Aust Natl Herb,” Marrubium vulgare L. in GBIF Secretariat (. 2023. [DOI]
- World Flora Online,” Published on the Internet, (. 2024
- Qualitative Tests for Preliminary Phytochemical Screening: An Overview,”. International Journal of Chemical Studies, 2020
- Antioxidant activity, phenol and flavonoid contents ofsome selected Iranian medicinal plants,”. African Journal of Biotechnology, 2006
- Use of a Free Radical Method to Evaluate Antioxidant Activity,”. LWT—Food Science and Technology, 1995
- 149 R. Re , N. Pellegrini , A. Proteggente , et al. Free Radical Biology and Medicine 26 (1999): 1231–1237.10381194
- Synergistic Anti Leukemia Effect of a Novel Hsp90 and a Pan Cyclin Dependent Kinase Inhibitors,”. Molecules, 2020. [PubMed]
- Co‐Inhibition of P‐gp and Hsp90 by an Isatin‐Derived Compound Contributes to the Increase of the Chemosensitivity of MCF7/ADR‐Resistant Cells to Doxorubicin,”. Molecules, 2021. [PubMed]
- Comparative Metabolite Profiling of Four Citrus Peel Cultivars via Ultra‐Performance Liquid Chromatography Coupled With Quadrupole‐Time‐of‐Flight‐Mass Spectrometry and Multivariate Data Analyses,”. Journal of Chromatographic Science, 2019. [PubMed]
- MS‐DIAL: Data‐Independent MS/MS Deconvolution for Comprehensive Metabolome Analysis,”. Nature Methods, 2015. [PubMed]
- Lapachol and (α/β)‐Lapachone as Inhibitors of SARS‐CoV‐2 Main Protease (Mpro) and hACE‐2: ADME Properties, Docking and Dynamic Simulation Approaches,”. Pharmacognosy Magazine, 2022
- Nobiletin as a Neuroprotectant Against NMDA Receptors: An In Silico Approach,”. Pharmaceutics, 2022. [PubMed]
- Identification of TBK1 inhibitors against breast cancer using a computational approach supported by machine learning,”. Frontiers in Pharmacology, 2024. [DOI | PubMed]
- 157 A. J. Siddiqui , R. Badraoui , S. Jahan , et al., "Targeting NMDA receptor in Alzheimer’s disease: identifying novel inhibitors using computational approaches", Frontiers in Pharmacology 14 (2023): 1208968, 10.3389/fphar.2023.1208968.37416066
