Selection of Mexican Medicinal Plants by Identification of Potential Phytochemicals with Anti-Aging, Anti-Inflammatory, and Anti-Oxidant Properties through Network Analysis and Chemoinformatic Screening

Many natural products have been acquired from plants for their helpful properties. Medicinal plants are used for treating a variety of pathologies or symptoms. The axes of many pathological processes are inflammation, oxidative stress, and senescence. This work is focused on identifying Mexican medicinal plants with potential anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects through network analysis and chemoinformatic screening of their phytochemicals. We used computational methods to analyze drug-like phytochemicals in Mexican medicinal plants, multi-target compounds, and signaling pathways related to anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence mechanisms. A total of 1373 phytochemicals are found in 1025 Mexican medicinal plants, and 148 compounds showed no harmful functionalities. These compounds displayed comparable structures with reference molecules. Based on their capacity to interact with pharmacological targets, three clusters of Mexican medicinal plants have been established. Curatella americana, Ximenia americana, Malvastrum coromandelianum, and Manilkara zapota all have anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects. Plumeria rubra, Lonchocarpus yucatanensis, and Salvia polystachya contained phytochemicals with anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence reported activity. Lonchocarpus guatemalensis, Vallesia glabra, Erythrina oaxacana, and Erythrina sousae have drug-like phytochemicals with potential anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects. Between the drug-like phytochemicals, lonchocarpin, vallesine, and erysotrine exhibit potential anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects. For the first time, we conducted an initial virtual screening of selected Mexican medicinal plants, which was subsequently confirmed in vivo, evaluating the anti-inflammatory activity of Lonchocarpus guatemalensis Benth in mice.

Inflammation is a manifestation of the innate immune system's defensive and adaptive response against injuries and harmful agents, such as bacteria, viruses, and toxins, in order to re-establish homeostasis.A dysregulated inflammatory response activates phytological processes, which can cause sepsis and organ failure.Chronic inflammation can occur Biomolecules 2023, 13, 1673 2 of 29 because of pro-inflammatory cytokine secretion stimulated by senescent cells [5].The phenomenon of low-grade chronic inflammation is characteristic of human aging and is termed "inflammaging" [2,6].The para-senescence process, which is triggered by secretory phenotypes linked to senescence (SASPs), leads to the generation of persistent inflammation [7].As a result, the inflammatory cascade causes an increase in extracellular ROS concentrations and oxidative stress.Elevated oxidative stress can cause the breakdown of the extracellular matrix (ECM) and activation of cell necrosis and apoptosis [5].
Aging and chronic degenerative pathologies exhibit common features, including elevated bioavailability of reactive oxygen species (ROS) and oxidative stress, chronic/persistent inflammation, glycation, and mitochondrial abnormalities.The production of excessive ROS leads to the destruction of nucleic acids and proteins, resulting in alterations to the cellular structure and functional outcome [5].
Oxidative stress arises from diminished anti-oxidant levels and the disturbance of the dynamic redox circuitry system, leading to the accumulation of free radicals within the body.The occurrence of oxidative stress initiates a negative sequence of events, leading to the modification of the chemical structure of cells, deterioration of the cell membrane, obstruction of crucial enzyme activities and energy production, and hindrance of vital cellular processes necessary for the body's regular functioning [8].Anti-oxidants have a vital function in the delay or prevention of cellular damage and the provision of essential protection against oxidative stress.Among the noteworthy enzymatic anti-oxidants are superoxide dismutase (SOD), glutathione peroxidase (GPX), catalase (CAT), and thioredoxin (Trx).The enzyme SOD facilitates the transformation of O 2 •− into the less reactive H 2 O 2 .GPX catalyzes the decomposition of H 2 O 2 and lipid hydroperoxide (LOOH), while CAT catalyzes the oxidation of H 2 O 2 to H 2 O and molecular oxygen (O 2 ), and Trx catalyzes the conversion of H 2 O 2 to produce H 2 O.To counteract or neutralize the effects of free radicals, the human body generates anti-oxidants.The stability of free radicals and anti-oxidants is vital for the appropriate physiological function of the body [5].
Natural Products have been used since ancient times in drug development.Between 1981 and 2019, 36.3% of the small molecules were developed for the treatment of antimicrobial, antiparasitic, and anticancer conditions based on natural products or their derivatives [9][10][11].Less than 20% of it is sourced from plants, as stated in reference [12].Cognitive and cardiovascular health benefits can come from consuming flavonoids, terpenoids, saponins, and polysaccharides like astragaloside, ginkgolide, ginsenoside, and gipenoside.The U.S. FDA has designated the term "Botanical Medicines" to refer to plant components that comprise complex mixtures.These mixtures require both chemical quality control and a manufacturing validation process [9].
For centuries, traditional medicine has been a fundamental part of the primary healthcare of populations worldwide [13,14].According to anthropological and archaeological data, medicinal plants were used in traditional Mexican medicine for over 5000 years.Using medicinal plants has been passed down through oral tradition, and many of these plants continue to be used as their historical usage [15].Mexico has 28,906 types of vascular plants that are grouped in over 2900 genera, 319 families, and 76 orders, according to Enci-cloVida [16].Its number of endemic species (approximately 50%) is second only surpassed by South Africa [17].Today, the number of plant species used in traditional medicine in Mexico exceeds 3000; the work "Atlas de las Plantas de la Medicina Traditional Mexicana" records 3103 [18], while an independent study recorded 3352 species [19].Enciclovida records 1025 species of native medicinal plants [20] consulted in June 2023.
Medicinal plant species, either native or introduced, are used by over 52 different ethnic groups.Herbal remedies are preferred by over 90% of Mexicans for their effectiveness, affordability, and fewer side effects compared to allopathic medicine [21].
Approximately 85% of health professionals are aware of herbal medicines, and approximately 75% recommend their use [22,23].Traditional herbal medicine remains valuable in rural and indigenous communities despite modern healthcare.Herbal medicines are important in traditional medicine for culturally appropriate and accessible healthcare.
Because of a lack of insurance and cultural beliefs, Mexican immigrants in the U.S. use alternative medicines.In Texas-Mexico border cities, 68.3% of pharmacists find patients using complementary and alternative medicines [21].
In addition, some states maintain botanical and medicinal plant gardens.Specialized stores sell traditional medicinal plants and medicines without proper quality control and regulation.In the 1990s, the Ministry of Health introduced a classification system to evaluate the safety and efficacy of herbal remedies used in traditional Mexican medicine.The NOM-072-SSA1-2012 [24] standardized and defined legal uses and packaging of herbal remedies since 2013.However, this regulatory law is not in force throughout Mexico [25].
Computational tools and novel algorithms have accelerated and optimized the drug discovery process.It is estimated at 20 years and saves around USD 1.3 billion [24].Several methods have been shown to reduce drug development costs by up to 50%.Methods include chemoinformatics, QSAR, docking, molecular similarity, network pharmacology, and computational design [21,26].Bioinformatic development has improved computer and web tools for designing new drugs.The analysis of compound databases using chemoinformatics has led to the discovery of new drugs.
The aim is to identify Mexican medicinal plants (MMPs) with one or more potential anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence (AOX-AINF-AAG-ASEN) effects.Identification of drug-like phytochemicals in MMPs eases plant selection with one or more potential AOX-AINF-AAG-ASEN effects.To achieve the objectives, (I) we will create an MMP database through the literature examination to be analyzed for reported AOX-AINF-AAG-ASEN activity; (II) we will search for AOX-AINF-AAG-ASEN phytochemicals candidates in MMPs using chemoinformatics and network pharmacology; and (III) the selection process for superior plants involves the thorough assessment of their multiple activities and compounds while also considering their synergistic properties that interact with different action mechanisms, including AOX-AINF-AAG-ASEN actions.

Identification of MMPs
To delimit the size of the population of MMPs, a further review was performed up to June 2023.The list of MMPs was bought from diverse, credible sources, including the Enci-cloVida National Biodiversity Platform [20], Atlas de las Plantas de la Medicina Traditional Mexicana, and the Index and Synonymic of Mexican Medicinal Plants (IMEPLAN) [18].

Search of AOX-AINF-AAG-ASEN Activity from MMPs
We reviewed bibliographies published in PubMed and Scopus up to July 2023 to find MMPs with or without scientific evidence related to AOX-AINF-AAG-ASEN activities.

Bibliographic Review of the Phytochemicals Reported to MMPs
To obtain a complete dataset of MMPS-reported molecules, we searched PubMed, chemical abstracts, and Scopus until July in the year 2023.

Creating MMP Dataset
We have created a dataset of MMPs, categorizing them according to family, genus, species, study type, synonymies, AOX-AINF-AAG-ASEN activities, other activities, diseases, and phytochemicals.

Identifying Plants with One or More AOX-AINF-AAG-ASEN Effects and Their Corresponding Phytochemicals
We analyzed the MMP dataset using scientific evidence to create a network that classified them by family, genus, and species and identified their AOX-AINF-AAG-ASEN capabilities.Using Cytoscape 3.8 and the Cytohubba plugin [27,28], the highly interconnected plants and their activity were identified.

MMPs Are Identified through Network Analysis Based on Taxonomic Features and AOX-AINF-AAG-ASEN Activities
We constructed a network derived from the MMP database based on family, genus, and species to locate plants without AOX-AINF-AAG-ASEN activity.These plants contain compounds with these reported activities.The network analysis was executed using the Cytoscape software version 3.8 [27].The Cytohubba plugin [28] was employed to find the most connected plants and their reported activity in the network.

Collection of Drugs with AOX-AINF-AAG-ASEN Activity
A review of the approved drugs by the FDA in the Drug Bank (https://go.drugbank.com/) [29] accessed on 15 August 2022 concerning compounds reported as anti-oxidants and anti-inflammatories was performed.To create the reference dataset, we use senolytic and anti-aging compounds obtained from previous work [30,31].

Identifying Drug-like Phytochemicals from MMPs
Estimation of the Drug-likeness Index of Phytochemicals in the MMP Dataset Based on the Quantitative Estimate of Drug-Likeness (QED) We were interested in detecting which phytochemicals from the literature had druglike properties.To identify phytochemicals with favorable characteristics to oral administration.The accomplishment of the aim was facilitated by the application of DruLito's druglikeness estimator, a Quantitative Estimate of Drug-Likeness (QED) (http://www.niper.gov.in/pi_dev_tools/DruLiToWeb/DruLiTo_index.htmlaccessed on 15 August 2022) [36].The DruLito software (version 1.0.0)used SDF files from PubChem to estimate the drug likeness of NMMPs.We analyzed only those molecules available on this server.
We used QED to evaluate the drug likeness of the compounds.QED is a technique for evaluating drug likeness by concurrently considering the primary molecular characteristics.QED score is determined by molecular descriptors, such as molecular weight, Lap, HBA, HBD, rotatable bonds, TPSA, and aromatic bond count.QED goes from zero to one.
Zero shows an unsuitable molecule, while one corresponds to molecules with favorable characteristics.We selected the phytochemicals with high QED scores (over 0.5).
2.6.Drug-like Phytochemicals Are Excluded Because of Toxicity and Lack of Targeting for AOX-AINF-AAG-ASEN Drug-like phytochemicals were sieved according to the following two criteria: (1) exclusion of undesirable effects (namely, tumorigenic, reproductive, and irritant functions) and ( 2) preservation of AOX-AINF-AAG-ASEN compounds with target information accessible in the PubChem database (retrieved February 2023).The signaling pathways used to construct this network are cited as the inflammatory response pathway, the relationship between inflammation, COX and EGFR, interleukin-1 family signaling, TP53 regulates transcription of cell death genes, the intrinsic pathway of apoptosis, TNF-related weak inducer of apoptosis (TWEAK) signaling pathway, cell cycle, an overview of the interferons-mediated signaling pathway, DNA damage/telomere stress-induced senescence, oncogene-induced senescence, DNA damage response, oxidative stress-induced senescence, senescence-associated secretory phenotype (SASP), and apoptosis.Check the Supplementary Materials for more information on plant data, targets, compounds, and network genes.

Identifying Drug-like Phytochemicals via Fingerprint Analysis with Reference Compounds
To determine their fingerprints and perform a comparison, we used ChemmineR and rcd in R-Studio version 3.4 [37].We used the extended value with a default length of 1024 (number of bits), taking rings and atomic properties.We used AOX-AINF-ASEN compounds as references and phytochemicals as tested datasets.Then, we performed a cluster analysis using the Tanimoto coefficient to compare them.We used Ward's clustering method to classify the molecules into three groups.To confirm the clustering results, the Dunn index and the silhouette coefficient were used [37].We performed different clustering methods; the reference compound cluster was analyzed to obtain molecules with one or more potential AOX-AINF-AAG-ASEN activities.

Structural Network Analysis Identifies MMPs with One or More Potential AOX-AINF-AAG-ASEN Activities
A structural network was employed, which was constructed using drug-like phytochemicals having similarity to references and target genes related to inflammation, oxidative stress, and senescence pathways.This network was also built based on signaling pathways found in Homo sapiens, which are associated with inflammation, oxidative stress, senescence, and their related processes.The signaling pathways involved in these processes included the inflammatory pathway, the relationship between inflammation, COX and EGFR, interleukin-1 family signaling, TP53 regulates transcription of cell death genes, the intrinsic pathway of apoptosis, TNF-related weak inducer of apoptosis (TWEAK) signaling pathway, cell cycle, overview of interferons-mediated signaling pathway, DNA damage/telomere stress-induced senescence, oncogene-induced senescence, DNA damage response, oxidative stress-induced senescence, senescence-associated secretory phenotype (SASP), and apoptosis.

Enrichment Analysis of the Pharmacological Targets from the Most Relevant MMPs
To determine the most important pharmacological target, drug-like phytochemical pathways involved in the structural network, an analysis was performed.They were built using Cytoscape software v 3.8.The JEPETTO plugin of Cytoscape [38] was employed to identify the major signaling pathways involved in MMPs and that interact with their phytochemicals.The significantly enriched pathways (p < 0.05, p-values) were revised using the XD-score, agreeing with previous reports [39,40].The KEGG database was used for this enrichment [41] (Supplementary File contains the results from the enrichment analysis).
Highly connected genes were identified in the network using the Cytohubba plugin score.

Animal Care
The mice were obtained from the medicine school of UNAM with the support of project # 008-CIC-2022.Their care was in line with Mexican standards (NOM-062-ZOO-1999; NOM-087-SEMARNAT-SSA1-2002) [42,43], and international official guidelines (Guide for the Care and Use of Laboratory Animals, revised in 2011) [44].
2.9.1.Experimental Testing L. guatemalensis Benth (Leguminosae) was collected in Santiago Pinotepa Nacional Km 8 road Pinotepa Nacional-Jamiltepecsws, Oaxaca, Mexico.A voucher on the plant is registered in the Mexican National Herb (IBUNAM: MEXU: 1139179).For more details, see Supplementary Materials.The dried root material was converted into chips and extracted with hexane for 72 h [45].Male BALB/c 8-week-old mice (26.67 g) were separated randomly into six experimental groups of six mice each (n = 6).Each group received one of these pretreatments: indomethacin (Sigma-Aldrich Cat.I7378, St. Louis, MO, USA) 10 mg/kg, as previously reported in inflammation tests [46,47].L. guatemalensis extract dissolved in corn oil and administered at 10, 31, and 100 mg/kg independently in each group, corn oil alone, and distilled water.L. guatemalensis doses were based on a chemotaxonomic study of a published article where the analgesic activity of another member of the Lonchocarpus genus (Lonchocarpus sericeus) was evaluated in mice [48].The LD50 of the extract was used as a reference for the establishment of the rest of the doses [48].The pretreatments were executed 30 min before the administration of carrageenan.In each group, inflammation was induced by administering a sub-plantar injection of 50 µL of 1% carrageenan (Sigma-Aldrich Cat.22049, St. Louis, MO, USA) into the right paw, as earlier described by Winter et al. [49][50][51].Simultaneously, in the left paw, distilled water was being administered.In the paw edema model, indomethacin, a potent non-steroidal anti-inflammatory drug (NSAID), is considered an excellent positive control due to its potent inhibition of vascular events [52].The edema of the paw was determined using a digital micrometer produced by Mitutoyo.Measurements were conducted before carrageenan was introduced and at various time points after inoculation [53].The edema rate was determined as per the previously specified method, employing the following equation: Edema (%) = (M f − M i /M i ) (100).Where M f is the paw measurement at different times, and M i is the paw measurement before inflammation was induced, both obtained from the same right paw [54].

Statistical Analysis
The data are presented as the mean ± standard error of the mean (s.e.m) of the percentage of edema induced by the administration of carrageenan.Multiple group comparisons were performed using one-way analysis of variance (ANOVA) post hoc Dunnett's test to detect inter-group differences.The difference was considered statistically significant when p < 0.05.The Graph Pad Prism 8.4.3 (GraphPad Software Inc., San Diego, CA, USA) software was used to perform the statistical analysis.
Figure 1 summarizes the methodology followed in this work.
The data are presented as the mean ± standard error of the mean (s.e.m) of the percentage of edema induced by the administration of carrageenan.Multiple group comparisons were performed using one-way analysis of variance (ANOVA) post hoc Dunnett's test to detect inter-group differences.The difference was considered statistically significant when p < 0.05.The Graph Pad Prism 8.4.3 (GraphPad Software Inc., San Diego, CA, USA) software was used to perform the statistical analysis.
Figure 1 summarizes the methodology followed in this work.

Identifying MMPs with and without Evidence of One or More AOX-AINF-AAG-ASEN Activities
Out of the 1025 MMPs that were consulted in EncicloVida [19], Atlas de las Plantas de la Medicina Traditional Mexicana, and Index and Synonymic of Mexican Medicinal Plants (IMEPLAN) [27], 364 demonstrated properties of AOX-AINF-AAG-ASEN, while 661 NMMPs lacked any such background.We list the results in Table 1.The clustering and relationship between families, genera, and species of all MMPs and their activities of anti-oxidation, anti-inflammation, anti-aging, and anti-senescence

Identifying MMPs with and without Evidence of One or More AOX-AINF-AAG-ASEN Activities
Out of the 1025 MMPs that were consulted in EncicloVida [19], Atlas de las Plantas de la Medicina Traditional Mexicana, and Index and Synonymic of Mexican Medicinal Plants (IMEPLAN) [27], 364 demonstrated properties of AOX-AINF-AAG-ASEN, while 661 NMMPs lacked any such background.We list the results in Table 1.The clustering and relationship between families, genera, and species of all MMPs and their activities of anti-oxidation, anti-inflammation, anti-aging, and anti-senescence are presented in Figure 2A.The analysis by the Cytohubba plugin showed that Asteraceae, Fabaceae, Chromolaena, Malvaceae, Piperaceae, and some genera like Granulosa, Peperomia, Chromolaena, Mimosa, and Malvaviscus (Figure 2B) have been extensively researched.Also, a large part of the MMPs lack any background related to AOX-AINF-AAG-ASEN activity.
are presented in Figure 2A.The analysis by the Cytohubba plugin showed that Asteraceae, Fabaceae, Chromolaena, Malvaceae, Piperaceae, and some genera like Granulosa, Peperomia, Chromolaena, Mimosa, and Malvaviscus (Figure 2B) have been extensively researched.Also, a large part of the MMPs lack any background related to AOX-AINF-AAG-ASEN activity.

MMPs with AOX-AINF-AAG-ASEN Multi-Effect and Their Corresponding Phytochemicals
The Figure 2A network was employed with Cytohubba to identify successful multieffect MMPS.The results are shown in Table 2.The most relevant families with multieffects included Dilleniaceae, Lauraceae, Olacaceae, Fabaceae, Arecaceae, Meliaceae, Asteraceae, Malvaceae, Solanaceae, and Sapotaceae.Ximenia americana, Curatella americana, Malvastrum coromandelianum, and Manilkara zapota exhibited the most significant multi-effect activity because they fulfilled all the activities considered in this research.The plants under investigation display diverse chemical structures that are classified in multiple clusters, which were previously noted for their AOX-AINF-AAG-ASEN properties.The flavonoids present in these plants are predominantly active compounds.3A, categorized by family, genus, and species without AOX-AINF-AAG-ASEN.The phytochemical profile of the phytochemicals is included.
Figure 3B shows a network analysis considering those MMPs with a greater number of phytochemicals that showed activity.In Table 3 are listed the most relevant families of MMPs (Apocynaceae, Fabaceae, and Lamiaceae).Rutin, stigmasterol, apigenin, and oleanolic acid are the phytochemicals that display the highest level of interconnectedness and are also emphasized.    3 shows the top 10 of the most interconnected nodes.

Identification of Reference Compounds with Some Reported AOX-AINF-AAG-ASEN Activity through Hierarchical Analysis
Figure 4 shows a heatmap of the main molecular similarities found among compounds with some AOX-AINF-AAG-ASEN properties.We analyzed the chemical space using the k-means algorithm and hierarchical grouping.Three clusters of reference compounds with individual AOX-AINF-AAG-ASEN activity were identified through the organization of the molecular descriptors' matrix.The top cluster in the heatmap (marked in a red rectangle) revealed the most similar compounds, which included homologous nucleoside antibiotics, macrolides, anti-oxidants, aminoglycosides, and cardiac glycosides.The analysis conducted allowed us to use these molecules as reference compounds for structural analysis, which was performed using fingerprints.Table 4 summarizes the chemical structure and pharmacological activity of some compounds from cluster one. of MMPs with potential AOX-AINF-AAG-ASEN activity and phytochemicals with reported studied activities.Table 3 shows the top 10 of the most interconnected nodes.

Identification of Reference Compounds with Some Reported AOX-AINF-AAG-ASEN Activity through Hierarchical Analysis
Figure 4 shows a heatmap of the main molecular similarities found among compounds with some AOX-AINF-AAG-ASEN properties.We analyzed the chemical space using the k-means algorithm and hierarchical grouping.Three clusters of reference compounds with individual AOX-AINF-AAG-ASEN activity were identified through the organization of the molecular descriptors' matrix.The top cluster in the heatmap (marked in a red rectangle) revealed the most similar compounds, which included homologous nucleoside antibiotics, macrolides, anti-oxidants, aminoglycosides, and cardiac glycosides.The analysis conducted allowed us to use these molecules as reference compounds for structural analysis, which was performed using fingerprints.Table 4 summarizes the chemical structure and pharmacological activity of some compounds from cluster one. of MMPs with potential AOX-AINF-AAG-ASEN activity and phytochemicals with reported studied activities.Table 3 shows the top 10 of the most interconnected nodes.

Identification of Reference Compounds with Some Reported AOX-AINF-AAG-ASEN Activity through Hierarchical Analysis
Figure 4 shows a heatmap of the main molecular similarities found among compounds with some AOX-AINF-AAG-ASEN properties.We analyzed the chemical space using the k-means algorithm and hierarchical grouping.Three clusters of reference compounds with individual AOX-AINF-AAG-ASEN activity were identified through the organization of the molecular descriptors' matrix.The top cluster in the heatmap (marked in a red rectangle) revealed the most similar compounds, which included homologous nucleoside antibiotics, macrolides, anti-oxidants, aminoglycosides, and cardiac glycosides.The analysis conducted allowed us to use these molecules as reference compounds for structural analysis, which was performed using fingerprints.Table 4 summarizes the chemical structure and pharmacological activity of some compounds from cluster one. of MMPs with potential AOX-AINF-AAG-ASEN activity and phytochemicals with reported studied activities.Table 3 shows the top 10 of the most interconnected nodes.
3.1.3.Identification of Reference Compounds with Some Reported AOX-AINF-AAG-ASEN Activity through Hierarchical Analysis Figure 4 shows a heatmap of the main molecular similarities found among compounds with some AOX-AINF-AAG-ASEN properties.We analyzed the chemical space using the k-means algorithm and hierarchical grouping.Three clusters of reference compounds with individual AOX-AINF-AAG-ASEN activity were identified through the organization of the molecular descriptors' matrix.The top cluster in the heatmap (marked in a red rectangle) revealed the most similar compounds, which included homologous nucleoside antibiotics, macrolides, anti-oxidants, aminoglycosides, and cardiac glycosides.The analysis conducted allowed us to use these molecules as reference compounds for structural analysis, which was performed using fingerprints.Table 4 summarizes the chemical structure and pharmacological activity of some compounds from cluster one.

Identification of Drug-like Phytochemicals from MMPs
The literature search yielded 1025 MMPs, and 1382 phytochemicals were obtained from the PubChem server in SDF format.The SDF format is crucial for analyzing the compounds in DruLito software (version 1.0.0) and obtaining molecular descriptors.This search has yielded a result indicating that only 1373 out of 1382 molecules are available in the SDF format.We used Quantitative Estimation of Drug-likeness (QED) to identify phytochemicals with drug-like properties, a criterion for early drug discovery.Table 5 shows the number of phytochemicals analyzed and the number that passed the rule with a score greater than 0.5.We obtained only 718 drug-like compounds.

AOX-AINF-AAG-ASEN Activity
Given the potential toxic effect of some phytochemicals, we conducted a toxicity risk assessment of our drug-like phytochemicals.We removed those with tumorigenic, irritant, or mutagenic properties, resulting in 417 compounds.
The 417 drug-like phytochemicals identified through DataWarrior were subjected to a supplementary filter.We analyzed the fingerprints of 148 drug-like phytochemicals with PubChem-reported targets.

Identification of Drug-like Phytochemicals from MMPs
The literature search yielded 1025 MMPs, and 1382 phytochemicals were obtained from the PubChem server in SDF format.The SDF format is crucial for analyzing the compounds in DruLito software (version 1.0.0) and obtaining molecular descriptors.This search has yielded a result indicating that only 1373 out of 1382 molecules are available in the SDF format.We used Quantitative Estimation of Drug-likeness (QED) to identify phytochemicals with drug-like properties, a criterion for early drug discovery.Table 5 shows the number of phytochemicals analyzed and the number that passed the rule with a score greater than 0.5.We obtained only 718 drug-like compounds.Given the potential toxic effect of some phytochemicals, we conducted a toxicity risk assessment of our drug-like phytochemicals.We removed those with tumorigenic, irritant, or mutagenic properties, resulting in 417 compounds.
The 417 drug-like phytochemicals identified through DataWarrior were subjected to a supplementary filter.We analyzed the fingerprints of 148 drug-like phytochemicals with PubChem-reported targets.

Identification of Drug-like Phytochemicals via Fingerprint Analysis with Reference Compounds
We performed a comparison between the reference compounds we obtained from the top cluster of Figure 4 and 148 drug-like phytochemicals from MMPs.The reference compounds include paromomycin, magnesium ascorbate, calcium ascorbate, rapamycin, digoxin, geldanamycin, alvespimicyn, oleuropein, strophanthin k, natamycin, streptomycin, oubain, tunicamycin, ginsenoside Rb1, rapamycin, roxithromycin, timosaponin a-III, and azithromycin.By a fingerprint analysis using the Tanimoto index, the similarity or disparity among compounds was determined.According to the molecular distances and Ward's method, the molecules were classified into three groups by elbow method (Figure 5A). Figure 5B shows the similarity between phytochemicals and reference compounds in Ward's cluster map.K-means (Figure 5C) and silhouette method (Figure 5D) corroborated these results.Paromomycin, magnesium ascorbate, calcium ascorbate, rapamycin, geldanamycin, alvespimicyn, oleuropein, natamycin, streptomycin, tunicamycin, rifaximin, and roxithromycin all have structural similarities with 60 drug-like phytochemicals from MMPs.A total of 50 phytochemicals with drug-like qualities share similarities with digoxin and strophanthin.Ouabain and ginsenoside share a similar structure with 49 drug-like phytochemicals (see Supplementary Data).Drug-like phytochemicals from MMPs show structural similarities with several reference compounds.We depict an example of this analysis in Figure 5, using paromomycin as a reference compound.At the end of fingerprint analysis, we obtained 109 drug-like phytochemicals like reference compounds.

Identification of MMPs with One or More Potential AOX-AINF-AAG-ASEN Activities through Structural Network Analysis
We only considered MMPs that met these requirements: (1) Phytochemicals with QED value over 0.5 for considered drug-like phytochemicals.(2) Phytochemicals showing structural resemblance to reference compounds.(3) Phytochemicals with target genes reported in bioassays.(4) Phytochemicals lacking AOX-AINF-AAG-ASEN activity.From 109 drug-like phytochemicals with structural similarity to the reference compounds, we obtained only 24 that met the aforementioned criteria.

Effects of L. guatemalensis Benth on Carrageenan-Induced Paw Edema
To corroborate the results obtained in silico, and because of the accessibility of Lonchocarpus guatemalensis, an extract from the root of this plant was evaluated in a carrageenan-induced paw edema assay.The results can be seen in Figure 7. Carrageenan administration triggered a significant increase in the right paw volume of all pretreated groups.L. guatemalensis extract significantly inhibited the formation of the edema from 2 h up to 6 h at 31 and 100 mg/kg.The inhibition observed is comparable to that seen with the pre-treatment of indomethacin (10 mg/kg).The group treated with L. guatemalensis (10 mg/kg) did not show any statistically significant decrease in edema.The groups administered with vehicles, corn oil, or distilled water did not exhibit any impact on the rise of edema in the right paw.
Supplementary Materials provide information on the determination of the percentage of anti-inflammatory effects caused by the hexanic extract of L. guatemalensis.To find several new associations, we determined the XD-score ranking.The most connected genes mentioned above (ABCB1, ALOX12, APOBEC3G, ATM, BRCA1, CHRNA4, CHRNB2, FEN1, MCOLN1, OTUD3, PINK1, POLI, POLK, PRKN, SMAD3, SMPD1, SNCA, TDP1, UCHL1, USP10, USP17L5, USP28, USP30, USP7, and USP8) were used to perform this analysis.The most relevant signaling pathways or processes associated with this pharmacological target gene related to these MMPs are shown in Supplementary Data.

Effects of L. guatemalensis Benth on Carrageenan-Induced Paw Edema
To corroborate the results obtained in silico, and because of the accessibility of Lonchocarpus guatemalensis, an extract from the root of this plant was evaluated in a carrageenan-induced paw edema assay.The results can be seen in Figure 7. Carrageenan administration triggered a significant increase in the right paw volume of all pretreated groups.L. guatemalensis extract significantly inhibited the formation of the edema from 2 h up to 6 h at 31 and 100 mg/kg.The inhibition observed is comparable to that seen with the pretreatment of indomethacin (10 mg/kg).The group treated with L. guatemalensis (10 mg/Kg) did not show any statistically significant decrease in edema.The groups administered with vehicles, corn oil, or distilled water did not exhibit any impact on the rise of edema in the right paw.  on carrageenan-induced paw edema in mice.A statistically significant reduction in the percentage of edema can be observed from 1 to 6 h during treatment with L. guatemalensis at 30 and 100 mg/kg, which are similar to anti-inflammatory positive control indomethacin.Values are expressed as mean ± standard error of the mean (n = 6).Different symbols show significant (p < 0.05) differences between treatments according to ANOVA) post hoc Dunnett's test.# carrageenan vs. indomethacin.$ carrageenan vs. HE (10 mg/kg).∞ carrageenan vs. HE (31 mg/kg).* carrageenan vs. HE (100 mg/kg).

Discussion
The present study should rescue the knowledge about the truly native Mexican plants used medicinally.Motivating results show that only 35.51% of CONABIO's 1025 MMPs underwent any study of AOX-AINF-AAG-ASEN activity, and 64.48% lack scientific reports on it.To recognize MMPs with potential AOX-AINF-AAG-ASEN phytochemicals through network analysis and chemoinformatics, a six-phase workflow was devised.This method enhances the probability of desired biological activity when tested.
Phase I.In MMPs, the anti-oxidant and anti-inflammatory activities have been the most studied (414 plants).The interest in these properties is consistent with the literature, which shows that various plant compounds possess helpful properties such as anti-oxidative, anti-inflammatory, and anticancer, among others [55].MMPs such as Ximenia Americana, Manilkara zapota, and Curatella Americana have reported all the studied activities [56][57][58][59][60][61][62].
In preclinical medicinal research, these plants show promise in combating age-related, degenerative, and infectious diseases caused by inflammation, oxidative stress, and cellular senescence.Our networks of the MMPs show relevant nodes within the Asteraceae and Fabaceae families.Edible flowers acquired popularity for their polyphenol content, which has anti-oxidants and anti-inflammatory properties.They include Asteraceae [63] and Fabaceae [64] families.The polyphenols are present in several natural products like tea, chocolate, fruits, and vegetables.This provides anti-oxidant, cytotoxic, anti-inflammatory, antihypertensive, and anti-diabetic benefits [65].
Phase 3. Despite having 1025 MMPs, we found a great redundancy in the phytochemicals reported in these plants.A total of 1382 phytochemicals were found.Only 1373 were available on servers.Characterizing the phytochemicals and their pharmacological properties required a preliminary step because of their heterogeneous nature.Drug-like capacity using the QED criteria allowed us to establish if a compound is suitable for oral administration in an in vivo model.From the 1373 phytochemicals collected from the MMPs, only 52.3% (718 MMPs) met the requirements to be considered "drug-like compounds".Of these drug-like compounds, only 58.07%(417) were free of tumorigenic, irritant, or mutagenic properties.The importance of these findings lies because some compounds may be harmful and not well absorbed when taken orally.These issues encourage the optimization of the discarded molecules to find new potential drugs [72] or use its functional fragments for the design of evolutionary libraries for the design of drugs in the future [73].
Selecting drug-like phytochemicals with bioassay targets improved the chances of finding potential activity molecules.Of 417 drug-like phytochemicals, only 148 (35%) had reported targets.The selection based on the target reported excluded potential candidates with structural similarity.These excluded phytochemicals may be considered if we use predictive techniques, such as molecular docking.However, under our focus, we use the network with predicted target c bioassays to obtain more targets experimentally and thus guarantee a better selection of compounds with potential activity by network pharmacology or other methodologies.Using predictive data to choose candidates for structural networks may create biases instead of using experimental bioassays.These issues would generate unsuccessful candidates when evaluated in an experimental model.For this reason, it is of great importance to conduct additional bioassays to obtain more experimental targets and thus ensure a better selection of compounds with potential activity through the use of network pharmacology or other methodologies.
Phase 4. For structural similarity analysis, we build a dataset with reference compounds.To achieve this aim, we use Drug Bank Online, a comprehensive and cost-free online resource that comprises intricate data on medication.The included data involve drug targets, drug actions, and the interactions between FDA-approved medicines, besides experimental drugs awaiting FDA approval.Drug Bank's high-quality content and primary origin make it one of the world's most used reference drug resources [56].Reference molecules classified as steroids (6 molecules), macrolides (4 molecules), and macrocyclic compounds (4 molecules) are in Figure 4's top cluster.Most reference compounds have chemotherapeutic activity against bacteria, parasites, and cancer cells.Some of them are used to treat cardiovascular illnesses (3 molecules), and one compound is an immunosuppressor (Drugbank online).Most of these compounds are contained in the common functional lactonic group (63.15%).Lactones are outstanding exponents of secondary metabolites because of their remarkable biological activities and chemical architectures.δ-Lactones are the most abundant because of the high stability of their lactone rings δ-lactone moiety is common in several natural compounds with diverse biological activities, such as HIV protease inhibition, apoptosis induction (goniothalamin and rasfonin), antileukemic (dictyopyrone C), antitumor (pironetin), leishmanicidal, trypanocidal, antifungal (argentilactone), antibacterial, anti-inflammatory, and anticancer.This characteristic makes them highly valuable reference compounds for our research [74].The reference compounds that remain are used to treat cardiovascular illnesses.Various structures analyzed showed steroidal or steroid-related structures, including precursor steroids like triterpene [75].Antimicrobial, anticancer anti-oxidant, immunomodulatory, and anti-inflammatory effects [74,76] are related to their steroidal structure.This characteristic makes them highly valuable reference compounds for our research.
Hierarchical and K-means clustering techniques were used in fingerprint analysis due to their success in previous research [30,31,37].Our previous works were used for the search for senolytics [30,31].We applied it in this investigation to phytochemicals using compounds with different pharmacological functions.The verification of the clustering analysis was performed by various means, including the Ward method, Dunn index, and silhouette coefficient, all of which confirmed the existence of the same drug-like phytochemicals when compared to the corresponding reference drug.Of 148 phytochemicals obtained from MMPS, 109 showed structural similarity with reference compounds.
Phase 5.An examination of the network's nodes resulted in the identification of 24 compounds from a pool of 119 phytochemicals.It is noteworthy that most of these compounds are classified as alkaloids, among which are apparicine, vallesine, pseudo tropine, scopine, tropine, corytuberine, cryptocoryne, erysotrine, germanicol, laudanosine, liriodenine, magnocurarine, and scopoline.One criterion for selecting the phytochemicals is the similarity to reference compounds.Despite the predominant use of steroids, macrolides, and macrocyclic compounds as reference materials in this stage of the study, the structural analysis using fingerprints revealed that the phytochemicals with drug-like properties were mainly of an alkaloid nature.This trend may be because there is a greater representation of alkaloids, according to studies searching for phytochemicals in plants [77].Within the human or animal body, alkaloids are responsible for dynamic biological activities.Research findings provide evidence of the benefits of alkaloids, which could be applied in discovering and designing new analogs that could be of therapeutic use in various treatments [78].Alkaloids can provide feasible lead compounds for the discovery and development of anti-inflammatory drugs.Most of the alkaloids are simple in structure and easy to modify, providing a viable framework and pharmacophore options for the design of new anti-inflammatory drugs [79].Another criterion for the selection of phytochemicals was the target identification.With the advancement of synthetic, computational, and biological approaches, structural simplification will play an increasingly important role in drug discovery and contribute to improving the efficiency and success rate of drug development [80].Through the utilization of network analysis and chemoinformatics, Lonchocarpus guatemalensis, Vallesia glabra, Erythrina oaxacana, and Erythrina sousae have been identified as potential subjects for forthcoming pharmacological studies.The identification of candidate phytochemicals, including lonchocarpin, vallesine, and erysotrine, in these MMPs presents an intriguing opportunity for future chemotaxonomic research to examine their activity in several species of plants.
Our search for plants with potential drug-like phytochemicals was limited because of the exclusion of MMPs based on their reported active compounds or toxic effects.Phase 6.The evidence shows that many drugs act via the modulation of multiple targets rather than a single target [81].For this reason, the enrichment analysis was performed because of the limited number of targets of lonchocarpin, vallesine, and erysotrine found in Lonchocarpus guatemalensis, Vallesia glabra, Erythrina oaxacana, and Erythrina sousae.In our enrichment analysis, the targets identified are found in several cellular processes, infections, cancer, and age-related diseases, where oxidation, inflammation, senescence, and aging play an important role in their development.Four Mexican plants with medicinal potential for AOX-AINF-AAG-ASEN activities were identified through virtual screening and network analysis.
The application of these computational tools has played a pivotal role in the discovery and design of highly efficacious medications.The article highlights the crucial role that cheminformatics and network pharmacology play in the process of drug discovery.The fusion of computational techniques and extensive data analysis has brought about a paradigm shift in the selection of medicinal plants.This integration enables more efficient exploration of the chemical space and accurate prediction of both therapeutic and toxicological biological activities of phytochemicals.Our search methodology significantly highlights the essential characteristics necessary for the continued search and selection of medicinal plants based on phytochemicals.This approach transcends the boundaries of chemotaxonomic or ethnomedicinal studies, enabling the identification of potential candidates.Our model demonstrates its effectiveness in predicting molecular properties and activities, helping in the prioritization and optimization of phytochemicals.
Furthermore, clustering is one of the simplest and most popular unsupervised machine learning algorithms.K-means, supported by network pharmacology and fingerprint-based structural similarity methods, opens possibilities for computer-aided design and discovery of new molecules with the desired properties.
This convergence has great potential to speed up the drug discovery process, reduce costs, and increase success rates in the identification of new therapeutic agents.Continued research and development in this area of phytopharmacology will undoubtedly pave the way for more efficient and precise drug design strategies, ultimately benefiting patients and advancing the field of pharmaceutical sciences.
Phase 7. Of all the candidates, we selected L. guatemalensis because of its availability.The rest of the plants, such as Vallesia glabra, Erythina oaxacana, and Erythina souzae, are not in the flowering process.This study provides evidence that supports the anti-inflammatory effects of the hexane extract obtained from L. guatemalensis, specifically highlighting lonchocarpin as a significant constituent according to the existing literature.These results obtained on carrageenan-induced paw edema are preliminary, and further studies on the anti-inflammatory mechanisms of L. guatemalensis extracts, such as inflammation markers (determined in serum concentration, gene expression, etc.), are required.However, they are encouraged to be tested in other anti-oxidant or senescence models.This result should be noted that lonchocarpin demonstrates additional biological effects, such as its anti-oxidant and anticancer properties [82], which impact various targets involved in senescence processes [83].All these aspects reinforce our strategy to select plants with several therapeutic effects, proposing L guatemalensis as a potential phyto-drug to treat different pathologies.Inflammation is the most common feature of many chronic diseases and complications while playing a critical role in carcinogenesis [84].Several studies have shown that Nrf2 contributes to the early steps of the anti-inflammatory process by orchestrating the recruitment of inflammatory cells and regulating gene expression through the anti-oxidant response element (ARE).Therefore, lonchocarpin may be a potential agent for treating disorders that are associated with inflammation because it increases the expression of antioxidant enzymes, such as heme oxygenase-1 (HO-1), NAD(P)H: quinone oxidoreductase 1 (NQO1), and manganese superoxide dismutase (MnSOD), which are all under the control of Nrf2/anti-oxidant response element (ARE) signaling.Moreover, this phytochemical induces the nuclear translocation and DNA binding of Nrf2 to ARE as well as ARE-mediated transcriptional activities [85].Also, lonchocarpin inhibits Wnt/β-catenin signaling and suppresses colorectal cancer proliferation.The inflammatory response must be tightly regulated because uncontrolled inflammation may lead to tissue injury.Among the many signaling pathways activated, the canonical Wnt/β-catenin has been recently shown to play an important role in the expression of several inflammatory molecules [86].Also, Wnt/β-catenin mediates transcription and regulates oncogenic and immune responses [87] via the interaction with NF-κB, which is a key master of inflammation.Due to this reason, Wnt/β-catenin is a promising target for drug development [88].L. guatemalensis induces inhibition of inflammation at early times, reducing the formation of edema induced by carrageenan simply due to the biological properties of lonchocarpin and, in particular, the inhibition of Wnt/β-catenin [82] gives us the probability that this chalcone is responsible for the effect of this plant [89], just as its presence has been related to this activity in other plants [90].However, it is worth mentioning that other phytochemicals may also be participating in the anti-inflammatory activity.
The pharmaceutical industry has been confronting amplified drug development costs, intensified failure rates, augmented competition for established targets, and an escalated demand for first-in-class drug development predicated on new targets and pharmacological mechanisms [91].We hold the belief that exploring medicinal plants and their phytochemicals is crucial for drug development.For this research, the information collected was analyzed by introducing biological big data analysis-based systems and chemoinformatics.Owing to the complexity of the phytochemical mixtures existing in plants, a multi-component and multi-target approach is required, which differs from the approach used with conventional medicines [92].Qualitative and quantitative analysis of diverse phytochemicals and their association with systems biology is a crucial factor in the development of drugs from natural sources.
To know the poly-pharmacological effects of multiple components of medicinal plants, much time and effort are required through conventional experiments.It is, therefore, necessary to resort to integrative strategies to select and develop better treatments with lower costs.Efforts have been made to understand the disease by integrating conventional molecular biology experiments, bio-big data analysis, and artificial intelligence technologies, which are increasing rapidly [93].Future research should optimize phytochemical analysis and network pharmacology to predict multiple targets.Combinatorial strategies that target multiple mechanisms may offer better chances of achieving clinically meaningful treatment [94].Through implementing this integrative approach, we have a firm belief that our MMPs, chosen based on their phytochemical and pharmacological network components, are well-suitable candidates.The findings from this in vivo investigation of Lonchocarpus guatemalensis have confirmed the predicted anti-inflammatory activity observed in the in silico analysis, demonstrating promising results.However, additional in vitro studies with Lonchocarpus guatemalensis and experimental evaluation of other MMP candidates are necessary.We are currently conducting experimentation on these candidates, and we expect the results will be positive.

Figure 1 .
Figure 1.Workflow for identifying MMPs and drug-like phytochemicals related to AOX-AINF-AAG-ASEN activities.(1) Relevant MMPs by multi-effect reported.(2) Relevant MMPs without reported effect but with phytochemicals with reported activity.(3) Determination of drug-like phytochemicals from MMPs without AOX-AINF-AAG-ASEN background.(4) Structural comparison between drug-like phytochemicals and reference AOX-AINF-AAG-ASEN compounds.(5) MMPs that have potential drug-like phytochemicals can network with pharmacological targets.(6) Enrichment analysis of the pharmacological targets from the most relevant MMPs.

Figure 2 .
Figure 2. Graphical representation of the MMP dataset into a network.(A) The taxonomic characteristics of the MMPs are integrated into the structural network, which also includes information on

Figure 2 .
Figure 2. Graphical representation of the MMP dataset into a network.(A) The taxonomic characteristics of the MMPs are integrated into the structural network, which also includes information on the availability or lack of scientific support for AOX-AINF-AAG-ASEN properties.(B) Network with the most relevant node.The table summarizes the ten most connected nodes.

Figure 3 .
Figure 3. Structural network of MMPs that may have AOX-AINF-AAG-ASEN properties based on their chemical content.(A) MMP network lacks scientific evidence of AOX-AINF-AAG-ASEN activity but has proven phytochemicals.(B) The most relevant nodes of families, genera, and species

Figure 3 .
Figure 3. Structural network of MMPs that may have AOX-AINF-AAG-ASEN properties based on their chemical content.(A) MMP network lacks scientific evidence of AOX-AINF-AAG-ASEN activity but has proven phytochemicals.(B) The most relevant nodes of families, genera, and species of MMPs with potential AOX-AINF-AAG-ASEN activity and phytochemicals with reported studied activities.Table 3 shows the top 10 of the most interconnected nodes.

Table 4 .
The structure and pharmacological activity of the molecules from the top cluster, which possess some AOX-AINF-AAG-ASEN properties, are presented.

Figure 4 .
Figure 4. Hierarchical structural clustering of the AOX−AINF-AAG−ASEN compounds.The top cluster with the most similar compounds is shown and highlighted in a red rectangle.The cluster

Figure 5 .
Figure 5. Steps of the process of selecting some AOX−AINF−AAG−ASEN potential compounds involved comparing 149 drug-like phytochemicals with the reference paromomycin based on their structural similarity.(A) The dashed line represents the suitable number of clusters determined by the Elbow method.(B) Cluster analysis of Ward's method; the senolytic cluster is shaped in red and the senolytic molecule as 1. (C) Cluster plot using K-means showed the same molecules in the senolytic cluster marked in blue.(D) Silhouette cluster representation to corroborate the previously described cluster analysis methods.* Means matrix of distances used in clustering by Ward's method in subfigure B.

Figure 5 .
Figure 5. Steps of the process of selecting some AOX−AINF−AAG−ASEN potential compounds involved comparing 149 drug-like phytochemicals with the reference paromomycin based on their structural similarity.(A) The dashed line represents the suitable number of clusters determined by the Elbow method.(B) Cluster analysis of Ward's method; the senolytic cluster is shaped in red and the senolytic molecule as 1. (C) Cluster plot using K-means showed the same molecules in the senolytic cluster marked in blue.(D) Silhouette cluster representation to corroborate the previously described cluster analysis methods.* Means matrix of distances used in clustering by Ward's method in (B).

Biomolecules 2023, 13 , x 18 of 32 Figure 6 .Figure 6 .
Figure 6.MMP network lacking AOX-AINF-AAG-ASEN properties with their phytochemicals and biological targets.(A) The network was constructed using MMPs, their putative AOX-AINF-AAG-Figure 6. MMP network lacking AOX-AINF-AAG-ASEN properties with their phytochemicals and biological targets.(A) The network was constructed using MMPs, their putative AOX-AINF-AAG-ASEN compounds, and their biological targets reported in PubChem.The signaling pathways involved in these processes were many (see Methods).The network was designed and analyzed using Cytoscape software.(B) The structural network was built with the Cytohubba plugin for the top 10 most connected nodes.

Figure 7 .
Figure 7. Effects of Lonchocarpus guatemalensis (10, 31, and 100 mg/kg) and indomethacin (10 mg/kg) on carrageenan-induced paw edema in mice.A statistically significant reduction in the percentage of edema can be observed from 1 to 6 h during treatment with L. guatemalensis at 30 and 100 mg/kg,

Table 1 .
The quantity of MMPs with or without evidence of some AOX-AINF-AAG-ASEN activity.

Table 1 .
The quantity of MMPs with or without evidence of some AOX-AINF-AAG-ASEN activity.

Table 2 .
The most relevant MMPs with reported multi-effect and their phytochemicals.

Table 3 .
The most relevant potential MMPs without reported AOX-AINF-AAG-ASEN effect and their phytochemicals.

Table 4 .
The structure and pharmacological activity of the molecules from the top cluster, which possess some AOX-AINF-AAG-ASEN properties, are presented.

Table 4 .
The structure and pharmacological activity of the molecules from the top cluster, which possess some AOX-AINF-AAG-ASEN properties, are presented.

Table 5 .
Several drug-like molecules from the total of phytochemicals obtained from MMPs.
Number of Molecules from Number of Molecules with a Number of Molecules with a Score * Pharmacological activity was obtained by a literature review.The reference compounds had pharmacological effects on inflammation, oxidative stress, aging, and senescence.

Table 5 .
Several drug-like molecules from the total of phytochemicals obtained from MMPs.

Table 6 .
Relevant MPP candidates lack AOX-AINF-AAG-ASEN effects and have pharmacological targets.The pharmacological activity was obtained by a literature review (for more details, see Supplementary Materials). *