key: cord-0026748-ro8ejbcl
authors: Obaidullah, Ahmad J.; Alanazi, Mohammed M.; Alsaif, Nawaf A.; Alanazi, Ashwag S.; Albassam, Hussam; AZ, Alanazi; Alwassil, Osama I.; Alqahtani, Ali M.; Tareq, Abu Montakim
title: Network Pharmacology- and Molecular Docking-Based Identification of Potential Phytocompounds from Argyreia capitiformis in the Treatment of Inflammation
date: 2022-01-31
journal: Evid Based Complement Alternat Med
DOI: 10.1155/2022/8037488
sha: 68ef2a04dc5e582b3822b7559f4d21e3502a7064
doc_id: 26748
cord_uid: ro8ejbcl

The methanolic extract of Argyreia capitiformis stem was examined for anti-inflammatory activities following network pharmacology analysis and molecular docking study. Based on gas chromatography-mass spectrometry (GC-MS) analysis, 49 compounds were identified from the methanolic extract of A. capitiformis stem. A network pharmacology analysis was conducted against the identified compounds, and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene Ontology analysis of biological processes and molecular functions were performed. Six proteins (IL1R1, IRAK4, MYD88, TIRAP, TLR4, and TRAF6) were identified from the KEGG pathway analysis and subjected to molecular docking study. Additionally, six best ligand efficiency compounds and positive control (aspirin) from each protein were evaluated for their stability using the molecular dynamics simulation study. Our study suggested that IL1R1, IRAK4, MYD88, TIRAP, TLR4, and TRAF6 proteins may be targeted by compounds in the methanolic extract of A. capitiformis stem to provide anti-inflammatory effects.

Inflammation describes a biological process that occurs in tissues to protect the host against harmful stimuli, such as microorganisms and abnormal or damaged cells. Inflammation stimulates the immune system and regulates protective responses via immune cells, blood vessels, and molecular biological agents [1, 2] . Many chronic diseases, including cardiovascular and gastrointestinal illnesses, diabetes, rheumatism, and cancer, are associated with upregulated inflammation [3] . Chronic diseases represent a major human health concern according to the World Health Organization (WHO). e incidence of chronic inflammation-related disorders is expected to steadily increase in the United States (US) over the next 30 years. Approximately 125 million people in the US were diagnosed with chronic diseases in 2000, with 61 million (21%) having multiple conditions [4] [5] [6] . Typically, cellular and molecular mechanisms and interactions among various factors can efficiently limit the potential damage and prevent further infection during an acute inflammatory response, resulting in the eventual repair of cellular homeostasis following the resolution of acute inflammation. However, uncontrolled acute inflammation can develop into chronic inflammation, causing a range of chronic inflammatory diseases [7] . ree out of five people worldwide die from chronic inflammatory conditions, such as stroke, respiratory infections, cardiovascular diseases, cancers, diabetes, and obesity [4] [5] [6] . A growing interest in the use of medicinal plants has developed for the treatment and management of diseases in an effort to identify safer and more efficient anti-inflammatory agents for the prevention of inflammatory conditions rather than using synthetic anti-inflammatory drugs [8] .

Argyreia capitiformis (Poir.) Ooststr. is a member of the Convolvulaceae family of the Argyreia genus, which is not toxic and has medicinal and ornamental uses [9, 10] . Traditionally, a paste made from the leaves of A. capitiformis has been used as an effective treatment of bruising on the legs. A. capitiformis has also been used in traditional medicine as a purgative and to treat sexual debility and ear pain [10] [11] [12] . Several studies have been reported for the Argyreia species, including antioxidant, anti-inflammatory, immunomodulatory, and CNS activities with several bioactive compounds [13] [14] [15] . However, no such study has been evaluated for A. capitiformis except the recent study on the methanolic extract of A. capitiformis leaves that suppressed the nuclear factorkappa B (NF-κB) pathway and inhibited the lipopolysaccharide-induced production of nitric oxide and inducible nitric oxide synthase in RAW 264.7 cells, demonstrating anti-inflammatory activities [9] . However, the chemical compounds found in A. capitiformis that are responsible for these anti-inflammatory effects and the underlying mechanisms remain unknown. Further research is necessary to identify the potential lead compounds responsible for these biological anti-inflammatory outcomes.

Network pharmacology has become a widely accessible analysis method following the increased availability of biomedical data sets during the postgenomic period, supporting the growth of the fields of systems biology and polypharmacology [16] . Complex compound-gene and compound-protein interactions can be evaluated systematically to develop a prototype for efficient therapy. New therapeutic mechanisms may be discovered by network pharmacology analysis, which is oriented toward a "multigoals, multi-disease" paradigm rather than "one target, one drug" [17] [18] [19] . Network pharmacology represents an effective method for selecting and elucidating the synergistic effects among bioactive chemicals through the mechanistic exploration of effects on multiple disease pathways [19, 20] . Additionally, spectrometric and chromatographic technologies used in the initial evaluation of medicinal plants provide valuable information on bioactive activities that aid in the selection of biologically active species. Alkaloids, phenolic compounds, organic acids, esters, and amino acids are among the chemicals that GC-MS can detect quickly and accurately. us, in this investigation, GC-MS was used to detect and identify phytochemical constituents in A. capitiformis [21] [22] [23] [24] . e network pharmacology approach connects targeted genes with the effects of bioactive compounds; thus, the present study was designed to elucidate the anti-inflammatory effects of the methanolic extract of A. capitiformis stem using a network pharmacology approach. Bioactive compounds in the methanolic extract of A. capitiformis stem were identified for this study using gas chromatographymass spectrometry (GC-MS) analysis, followed by a molecular docking assay to investigate potential ligand-receptor interactions, including the assessment of binding affinity and stability.

e stems of A. capitiformis were collected from the Sitakunda Eco-park, Chittagong, Bangladesh, in March 2020 and later identified by a taxonomist. e stems were subjected to air-drying and ground to a coarse powder. e powder (200 g) was soaked in methanol (1 L) for 7 days [25, 26] . Subsequently, the extract was filtered through Whatman filter paper and evaporated at 45°C. After the evaporation, 2.67 g of the black methanol extract yield (1.34%) was collected in an amber glass vial and refrigerated at 4°C until further use.

Analysis. An Agilent GC 7890A (Agilent Technologies Inc., Wilmington, DE, USA), combined with a triple-axis detector 5975 C single-quadrupole mass spectrometer, was used for GC-MS analysis. e chromatographic column was an Agilent HP 5MS column (30 m × 0.25 mm × 0.25 µm film thickness), using high-purity helium as the gas carrier at a flow rate of 1 mL per min. e injector temperature was 230°C, and the sample was injected using a splitless injector at 20 : 1. e temperature was set primarily to 40°C (held for 1 min), raised to 150°C at a rate of 5°C per min (held for 2 min), before being increased to 300°C at a rate of 5°C per min (held for 10 min). e temperature of the MS ion source was set to 150°C, and the temperature of the inlet line was set to 280°C. e scan range was set between 50 and 550 mass, with 70 eV electron energy and a 4-min solvent delay. Finally, by comparing the spectra against the NIST 2008 database (National Institute of Standard and Technology library), tentative compounds were identified. e total analysis time required for the sample was 65 min [27] .

e network pharmacology analysis was performed using the STITCH platform (http:// stitch.embl.de/) to identify putative associations between the identified compounds and target genes. Multiple compound targets were identified using the Homo sapiens genome [27] [28] [29] . Multiple functional nodes and edges were identified in the network. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed on the components identified in the network to obtain a biological interpretation of 2

Evidence-Based Complementary and Alternative Medicine the vast list of potential targets and to identify potential antiinflammatory pathways that are targeted. e STRING (search tool for retrieval of interacting genes) database (https://string-db.org), which includes predicted protein-protein interactions (PPIs), was used to predict functional protein interactions [30] .

2.4.1. Protein Selection and Preparation. Six targeted proteins were selected from the KEGG analysis: interleukin 1 (IL-1) receptor type 1 (IL1R1; PDB: 1ITB) [31] , IL-1 receptorassociated kinase 4 (IRAK4; PDB: 6EGA) [32] , myeloid differentiation factor 88 (MYD88; PDB: 4EO7) [33] , TIR domain-containing adaptor protein (TIRAP; PDB: 4FZ5) [34] , Toll-like receptor 4 (TLR4; PDB: 3FXI) [35] , and tumor necrosis factor (TNF) receptor-associated factor 6 (TRAF6; PDB: 3HCT) [36] . Protein structures were retrieved in .pdb format from the Protein Data Bank (http://www.rcsb.org/ pdb). e proteins were prepared using Schrödinger (Maestro v11.1), utilizing the force field OPLS3 [27] .

We selected 47 compounds identified from A. capitiformis, according to the qualitative GC-MS analysis, which we submitted to the molecular docking study. e selected compounds were retrieved in .sdf format from the PubChem database. In addition, aspirin (CID: 2244) was utilized in this study as a positive antiinflammatory control. e three-dimensional (3D) structures of the selected compounds were constructed in Schrödinger using LigPrep (Maestro v11.1), utilizing the force field OPLS3.

To create receptor grids and execute a molecular docking analysis, Glide (Schrödinger, Maestro v11.1) was used. e grids were generated in Glide with the default settings and the OPLS3 force field. A cubic box with a boundary box (14Å × 14Å × 14Å) was specified for the receptors. All docking studies were conducted using Glide's standard precision (SP) and flexible docking modes, and the lowest docking score for each ligand was recorded.

e free energies of binding (ΔG; kcal/mol) for each ligand and the target receptors were calculated using the Schrödinger software package Prime/MM-GBSA module (OPLS3) [37, 38] . e ligand efficiency (LE) was assessed for each ligand by obtaining the ratio of ΔG to the number of heavy atoms (NHA): LE � −(ΔG)/NHA [39] .

e molecular dynamics simulation study was conducted in YASARA dynamics by the aid of AMBER14 force field [40, 41] . e docked complexes were optimized and cleaned, and hydrogen bond network system was oriented. e cubic simulation cell was created where the TIP3P solvation model was used with periodic boundary conditions [42] . e simulation system was neutralized at 310 K temperature, pH 7.4, and 0.9% NaCl. e initial energy minimization was conducted by steepest grained algorithms by simulating annealing methods. e long-range electrostatic interactions were calculated by the particle mesh Ewald (PME) methods with a cutoff radius of 8.0Å [43, 44] . e simulation time step was set as 2.0 fs. e simulation trajectories were saved after 100 ps and finally run for 20 ns by following the constant pressure and Berendsen thermostat [45] . e simulation trajectories were used to calculate the root-meansquare deviation, solvent accessible surface area, radius of gyration, and hydrogen bond [46] [47] [48] [49] [50] [51] [52] [53] [54] .

In this work, methanol was utilized as the solvent for extraction, resulting in a 1.34% yield. e GC-MS analysis of the A. capitiformis stem methanolic extract revealed 49 compounds with different retention times and peak areas (Table 1 and Figure S1 ). e methanolic extract contained the following identified compounds: stigmast-4en-3-one (20.78%, RT: 58.161); hexadeca-2,6,10,14-tetraen-1-ol, 3,7,11,16-tetramethyl-, (E,E,E)-(18.36%, RT: 51.477); ursa-9(11),12-dien-3-one (10.35%, RT: 58.55); ursa-9(11),12-dien-3-ol (6.89%, RT: 55.077); 2- 

A KEGG pathway analysis of potential target genes (IL1R1, IRAK4, MYD88, TIRAP, TLR4, and TRAF6) revealed signaling pathways related with anti-inflammatory effects (see Table 2 ). Figure S2 shows the PPI network with 12 proteins (IL1R1, IRAK2, MYD88, IRAK4, IRAK2, MYD88, TIRAP,  TRAF6, TLR3, TLR4, TLR5 , and TLR6), all of which have anti-inflammatory effects. Tables S1 and S2 represent the biological processes and molecular functions of the genes interacting with the compounds, respectively. Proteinprotein interaction (PPI) status of 10 proteins with coexpression is demonstrated in Table S3 .

A total of 47 compounds docking results are presented in Tables 3-8 , which show the findings. Aspirin has been used as a positive control for this study.

is study's findings reveal that binding energies of most of the ligands to receptors are negative, as later validated by MM-GBSA analysis. e proteins and ligands' molecular interaction is presented in Supplementary Materials (Tables S4-S9 ). To better understand the docking score, we have studied the ligand efficiency, which demonstrated excellent support for molecular docking scores. e molecular dynamics simulation of the targeted receptors (IL1R1, IRAK4, MYD88, TIRAP, TLR4, and TRAF6) against the best stable compounds is presented in Figures 1-6 . 

Plants have been used throughout the history of traditional medicine to induce a variety of biological effects, and extensive pharmaceutical resources have recently been devoted to the identification and investigation of new remedies, including those derived from plants. A critical issue encountered by researchers who perform phytoscience is that a single plant can harbor a wide range of bioactive chemicals [55, 56] . e pharmaceutical industry relies on phytochemicals to develop new drugs and therapeutic agents. Finding natural bioactive components is the first step in developing novel drugs. Screening plant extracts for therapeutically active chemicals is a novel strategy. It is important to know that plants have a lot of different types of phytochemicals, which have a lot of different biological properties.

ese include antioxidant, anti-inflammatory, antidiarrhea, antiulcer, and anticancer activities [57] [58] [59] .

Determining which compounds are responsible for the biological activities associated with plant materials can help understand toxicities, determine suitable doses, and identify ideal methods for compound extraction. e successful acquisition of components from plant materials depends primarily on the solvent used during the extraction process [60, 61] . e stem methanolic extract was analyzed by GC-MS analysis, and the results indicated 49 different chemicals with varying retention times and peak areas.

e network pharmacology analysis was performed to evaluate potential interactions between the identified chemical compounds and proteins, followed by multiple comparisons to determine the number of genes responsible for the anti-inflammatory effects, which was represented by the STITCH platform ( Figure 7 ). Table 2 shows the results of KEGG pathway analysis performed on potential target Evidence-Based Complementary and Alternative Medicine genes, which identified signaling pathways associated with anti-inflammatory actions. Analyses of biological processes and molecular functions, which identified proteins and pathways with significant values, were also performed (Table S1) . A comparison of the compound-gene interaction network and the protein-protein interaction network was performed to reveal biological and molecular functions (Table S2) . From a molecular and functional perspective, these findings can assist in understanding the computational rules of compounds that have the potential to treat diseases. For example, protein kinase C (PKC) binding is known to treat inflammatory diseases [62] , and the present results showed that binding with TIRAP, UGT1A10, and UGT1A7 was significant (p � 0.0447) and similar to that for PKC. Also, retinoic acid binding has a great impact on the immune system and the inflammatory response. Our findings identified the genes UGT1A1, UGT1A3, UGT1A7, UGT1A8, and UGT1A9 as being the most significant interactions for the compounds in our extract (p � 2.69E − 06). Although a remarkable amount of functional and structural data has been compiled for each identified protein, our knowledge regarding protein-protein relationships remain scattered. e purpose of the STRING database is the collection, scoring, integration, and complementation with computational predictions for all public sources of PPIs [63] [64] [65] [66] . In the present study, 10 proteins (IL1R1, IRAK4,  IRAK2, MYD88, TIRAP, TRAF6, TLR3, TLR4, TLR5 , and TLR6) were used to analyze a PPI network ( Figure S2 ), which downstream. TAK1 activates the IKKβ-to-IκB-α-to-NF-κB pathway, inducing the transcription of proinflammatory genes. TAK1 also influences gene expression by activating mitogen-activated protein kinase (MAPK) cascades. TRAF6, βRIP1, and TAK1 are activated by the binding of TLR with TRIF, resulting in the activation of MAPK, interferon regulatory factor 3 (IRF3), NF-κB, and interferon-transcription activation. e TRIF pathway also promotes the release of proinflammatory cytokines, although to a lesser degree than the MYD88 pathway. TLRs are a class of specific receptors that are key players in mounting an effective innate immune response to infection [67] [68] [69] . Subsequently, gene co-occurrence can be used to identify gene families whose patterns exhibit similarity across genomes. ree types of analyses have been used to examine genomes (neighborhood, fusion, gene co-occurrence) based on the systemic comparison of all-against-all genomes to evaluate the impacts of historical genome restructurings, genetic gains and losses, and gene fusion [70, 71] . In the present study, 100% sequence conservation was observed for the ten proteins (IL1R1, IRAK4, IRAK2, MYD88, TIRAP,  TRAF6, TLR3, TLR4, TLR5 , and TLR6) selected for the PPI network analysis, as shown in Figure S3 . Additionally, gene coexpression was also studied in the present study, as summarized in Table S3 and Figure S4 . e coexpression pathway is predicated by performing gene-by-gene correlation testing across many gene expression databases. [63, 72] . e present findings suggested co-expression among the ten selected proteins (IL1R1,  IRAK4, IRAK2, MYD88, TIRAP, TRAF6, TLR3, TLR4,  TLR5 , and TLR6) in Homo sapiens. In addition, coexpression (transferred) of three more genes was observed in Gallus gallus, Mus musculus, and Rattus norvegicus. Structure-based drug discovery is gradually becoming a key technique for facilitating the rapid and cost-effective discovery and optimization of lead compounds. e use of a rational, structure-based drug design strategy is more efficient than conventional drug development techniques because this approach seeks to understand the molecular basis of diseases and incorporates information regarding the biological target's 3D structure during the drug design process [37] . A molecular docking study was incorporated into the present study to predict the complex structure formed by ligand-protein binding and analyze the ligand's conformational space within the protein-binding site. A score function is then used for each docking analysis to assess the free energy of the interaction between the protein and ligand [37, 73, 74] . Additionally, the LE is calculated, which can enrich docking functions and allow for the coordination between docking outcomes and experimental results. Critical information regarding a molecule's properties, such as the NHA, can then be combined into a single table [39] . IL-1 controls a range of innate immune pathways, making it a key regulator of inflammation [75] . Two IL-1 cell surface receptors and a decoy receptor have been identified, including IL1R1 and IL1R2. First, IL-1 binds with IL1R1, inducing the formation of a heterodimer between IL1R1 and either IL-1RAcP or IL1R3, followed by IL-1 receptor-associated kinase (IRAK) and MyD88. e inflammatory response induced by IL1R1 occurs when IL1R1 binds with either the IL-1α or IL-1β ligands, whereas the T-lymphocytes, fibroblastic cells, epithelial cells, and endothelial cells have been indicated [76] [77] [78] . In the present study, selected six target proteins based on the network pharmacology analysis for the molecular docking study: IL1R1 (PDB: 1ITB), IRAK4 (PDB: 6EGA), and MYD88 (PDB: 4EO7) was used for the docking study. For IL1R1 (PDB: 1ITB), the majority of the 47 tested compounds demonstrated good docking scores, except for 12 compounds (Table 3) . Among the 33 compounds with good docking scores, catechol (3.02), phenol (2.94), 4(1H)-pyrimidinone, 6-hydroxy-(2.79), hydroquinone (2.70), 3-furaldehyde (2.51), and sulcatone (2.51) showed the best ligand efficiencies, compared with the LE value of 1.63 for the positive control aspirin (Table S4 and Figure S5 ). e compounds with the best ligand efficiencies were found to interact with LYS-93 by H-bond (distances >4Å), and catechol interacted via two H-bonds. Similar findings were also observed for the positive control aspirin. According to an earlier study, the noncontiguous binding epitope containing LYS-93 was identified for IL-1β [79] .

IRAK4 (PDB: 6EGA) also exhibited similar findings in the docking experiment, with only 12 compounds that did not display any interactions. Catechol (4.45), phenol (4.35), p-vinylguaiacol (4.15), methylcyclohexane (4.08), 3-furaldehyde (4.06), and 4-methyl-1,5-heptadiene (3.94) showed Table 4 . e LE value of aspirin was 2.39. e molecular interactions between the compounds with the best LE values and IRAK4 (PDB: 6EGA) are presented in Table S5 and Figure S6 . ASP-329 formed H-bonds with catechol (4.02, 3.48Å), p-vinylguaiacol (4.05, 4.06Å), phenol (3.49Å), and aspirin (4.24Å), whereas no H-bond interactions were observed for 4-methyl-1,5-heptadiene or methylcyclohexane. Additionally, 3-furaldehyde interacted with MET-265 (4.0Å) via one H-bond. Remarkably, ASP-329, MET-265, and LYS-213 were all reported to interact with IRAK4 in a previous study [32] . Additionally, MYD88 (PDB: 4EO7) was also studied as a potential target of the anti-inflammatory effects of the methanolic extract of A. capitiformis stems. In the molecular docking experiment, 29 compounds were found to interact with MYD88 (Table 5) Table S6 and Figure S7 . ASP-156 was reported to form H-bonds with coniferol, 2-(2-hydroxy-2-phenylethyl)-3,5,6-trimethylpyrazine, and (Z,E)-farnesol, whereas hydroquinone interacted with ASP-156 through the formation of hydrophobic bonds (3.51Å). No H-bond interactions were reported for 4-methyl-1,5-heptadiene. e other compound, trans-13-docosenamide, interacted with ARG-160 (5.21Å) and GLU-159 (4.36Å) via H-bonds.

TIRAP, also known as MYD88 adapter-like (Mal), is an important link between MYD88 and the receptor complex formed by TLR2 and TLR4 activation following bacterial infection [32] . MYD88 activates IRAK1 and IRAK4 and eventually activates TRAF6, causing the prototypic inflammatory transcription factor NF-κB to translocate to the nucleus [80, 81] . TIRAP is the second adaptor that has been identified to mediate NF-κB activation through the activation of TLR4 and TLR2 signaling pathways [82] [83] [84] [85] . e molecular docking of TIRAP (PDB: 4FZ5) was examined against the 45 compounds identified in the methanolic extract of A. capitiformis stems, of which 35 demonstrated interactions and 10 did not (Table 6 ). e best interacting compounds were phenol (3.06), 4-ethylresorcinol (2.96), catechol (2.84), methylcyclohexane (2.72), hydroquinone (2.71), and 3-furaldehyde (2.50), compared with the positive control aspirin, which exhibited an LE of 1.82 ( Figure S8 ). e molecular interactions between the compounds with the best LE values and TIRAP (PDB : 4FZ5) are presented in Table S7 . LEU-107 formed H-bond with 3-furaldehyde, catechol, phenol, and aspirin, and a hydrophobic interaction was also observed for catechol, methylcyclohexane, and phenol. Another protein residue, LYS-210, interacted with 3-furaldehyde, 4-ethylresorcinol, hydroquinone, and phenol via hydrophobic interactions. In a previous study, LYS-210 was reported to form a hydrophobic interaction with NF-κB [86] . In addition, TLR4 (PDB: 3FXI) was compared against 45 compounds found in the methanolic extract of A. capitiformis stems, although 17 compounds did not interact with TLR4 (Table 7) . Phenol (4.59), 4-methyl-1,5heptadiene (4.45), catechol (4.39), sulcatone (4.17), hydroquinone (4.01), and coniferol (4.003) displayed the highest LE values, whereas aspirin showed an LE value of 2.75. e molecular interactions between the compounds with the best LE values and TLR4 (PDB: 3FXI) are presented in Table S8 and Figure S9 . SER-441 was found to form H-bond interactions with phenol and sulcatone, whereas hydrophobic interactions were exhibited for phenol, hydroquinone, catechol, and coniferol. VAL-82 interacted via hydrophobic interactions with coniferol, 4methyl-1,5-heptadiene, catechol, phenol, and hydroquinone. In a previous study, SER-441 was reported to interact strongly with apigenin-7-O-glucoside [87] . Similar interactions were also observed for ILE-80 and VAL-82 [88] .

Finally, TRAF6 (PDB: 3HCT) was assessed against the 45 chemicals found in the methanolic extract of A. capitiformis stems, of which 34 were identified as interacting (Table 8) Table S9 and Figure S10 . ARG-6 formed H-bond interactions with coniferol, 2-(2-hydroxy-2-phenylethyl)-3,5,6trimethylpyrazine, and aspirin, whereas hydrophobic interactions were identified for 2-(2-hydroxy-2-phenylethyl)-3,5,6-trimethylpyrazine and sulcatone. GLN-54 also formed H-bond interactions with multiple compounds, including 2-(2-hydroxy-2-phenylethyl)-3,5,6-trimethylpyrazine and 3-furaldehyde. e root-mean-square deviations from the C-alpha atoms from the docked complexes IL1R1 (PDB: 1ITB) proteins are illustrated in Figure 1 . e complexes' RMSD values were calculated to find out the deviations among the simulation complexes and structural stability. Figure 1 (a) demonstrates that the complexes had similar RMSD profiles and did not fluctuate much in the simulation trajectories. e RMSD profile of the complexes reached the steady state after 5 ns and maintained the structural stability till the last periods of the simulations, which defines the stability of the complexes. e SASA of the complexes was analyzed to find out the change in the surface area. e higher SASA represents the expansion of the surface area of the protein, whereas the lower SASA indicates the truncated nature of the complexes. Figure 1(b) shows that the complexes were in a stable state in SASA. e radius of gyration defines the labiality and mobility of the complexes, where Figure 1(c) indicates the lower deviations. e hydrogen bond patterning follows a similar stable profile (Figure 1(d) ). e molecular dynamics simulation of IRAK4 (PDB: 6EGA) is presented in Figure 2 . e RMSD of the complexes had a stable trend in the RMSD profile for all the complexes except catechol. e higher RMSD of these complexes defines the higher flexibility of these compounds in the simulating environments (Figure 2(a) ). e SASA profile of the complexes was stable, did not fluctuate much, and had a steady trend in SASA (Figure 2(b) ).

is SASA profile correlates with the stable and rigid profile of the complexes (Figure 2(c) ). e radius of gyration and hydrogen bond Evidence-Based Complementary and Alternative Medicine pattern systems were similar and did not change too much in the simulations (Figure 2(d) ). e RMSD of the MYD88 (4EO7) protein complexes had a lower level of fluctuations, and lower deviations were observed across the compounds. All compounds had a lower RMSD than 2.5Å in whole simulation periods (Figure 3(a) ). e SASA of the MYD88 (4EO7) complexes were stable, but the complex of trans-13docosenamide had a lower SASA than the other complexes.

is SASA profile of trans-13-docosenamide defines the MYD88 (4EO7) experienced the condensed conformation upon binding with the corresponding ligands (Figure 3(b) ). e 2-(2-Hydroxy-2-phenylethyl)-3,5,6-trimethylpyrazine complexes had a higher Rg than other complexes, which defines the complexes' flexible nature than other complexes (Figure 3(c) ). ese complexes also had a higher SASA in Figure 3 (b), which depicts the changes in the confirmations than other complexes in simulated environments. e hydrogenbonding pattern of the complexes for MYD88 (4EO7) protein was found stable and did not change too much in simulations (Figure 3(d) ).

e RMSD profile of TIRAP (PDB: 4FZ5) complexes is illustrated in Figure 4 . All complexes from TIRAP (PDB: 4FZ5) protein had an initial rise of RMSD due to a higher degree of flexibility but stabilized subsequently after 5 ns times. e complex 4-ethylresorcinol had a higher RMSD than the other complexes, which might be responsible for the more remarkable conformational changes and the flexibility of the compounds (Figure 4(a) ). e SASA of the complexes had a stable and similar trend for all the compounds. But the complex phenol had a lower SASA profile at the last phase of SASA, which defines the complexes' truncated nature in simulations (Figure 4(b) ). e radius of gyration profile of the complexes had a lower trend, which correlates with the less flexibility of the complexes (Figure 4(c) ). e hydrogen bond patterning of the complexes had a stable profile in Figure 4(d) .

e molecular dynamics simulation study of the TLR4 (PDB: 3FXI) and complexes was done to analyze the structural deviation in the docked structure. e root-meansquare deviations of all complexes are illustrated in Figure 5 (a). e RMSD value of the complexes initially followed the upper trend from the beginning. is might be happening due to the higher flexibility level. But all the complexes from TLR4 (PDB: 3FXI) had a stable profile after 10 ns times and followed a similar trend until the last periods, demonstrating structural stability. e SASA of the TLR4 (PDB: 3FXI) complexes had lowered the degree of the deviations from the beginning and followed lower fluctuations, which define no changes in the surface area of the complexes ( Figure 5(b) ). e Rg and hydrogen bond patterns of the simulation systems were stable and did not change too much, which correlates with the structural stability (Figures 5(c) and 5(d)). e docked complexes from the TRAF6 (PDB: 3HCT) protein and their simulation descriptors are illustrated in Figure 6 . e RMSD profile of the TRAF6 (PDB: 3HCT) complexes defines that the phenol and aspirin had a higher level of RMSD than other complexes, which correlates with the comparative higher degree of the deviations of the complexes. All complexes had reached a stable state after 5 ns of the simulation times. Moreover, the complexes exhibit RMSD lower than 2.5Å, which defines the complexes with a higher rigidity degree (Figure 6(a) ). e SASA profile of the complexes had a lower deviation, as illustrated in Figure 6 (b). e phenol had a lower SASA value than all the compounds, indicating the truncated nature of the protein complexes compared with the others. Moreover, the radius of gyration from Figure 6 (c) demonstrates that the complexes had a lower degree of deviations, and no significant higher fluctuations were observed.

is Rg profile indicates the complexes had lower mobility and flexibility during the simulation times. e hydrogen bond pattern of the complexes was stable during the simulations (Figure 6(d) ). 

e network pharmacology analysis revealed key pathways involved in the anti-inflammatory activities induced by the chemical compounds found in the methanolic extract of A. capitiformis stem. Six inflammation pathways were obtained from the KEGG pathway analysis (IL1R1, IRAK4, MYD88, TIRAP, TLR4, and TRAF6), and molecular docking studies of these pathways revealed that the identified chemical compounds had strong binding affinities with these pathway components.

e current study discovered that 3-furaldehyde, phenol, catechol, and hydroquinone were effective anti-inflammatory compounds found within the methanolic extract of A. capitiformis and played an important part in the inflammation pathway by targeting these six proteins. Additional in vitro and in vivo experiments will be helpful to validate and optimize the findings of this study.

GC-MS: Gas chromatography-mass spectrometry KEGG: Kyoto Encyclopedia of Genes and Genomes WHO:

World Health Organization NF-κB: Nuclear factor-kappa B STRING: Search Tool for Retrieval of Interacting Genes PPIs:

Protein-protein interactions IL: Interleukins PKC: Protein kinase C MAPK: Mitogen-activated protein kinase IRF:

Interferon regulatory factor IRAK:

Receptor-associated kinase TLR:

Toll-like receptor 3D: ree-dimensional PDB:

Protein Data Bank RMSD: Root-mean-square deviation of the complexes SASA:

Solvent accessible surface area Rg:

Radius of gyration.

Data Availability e data are available within the manuscript and also accessible from the corresponding author upon request.

e authors declare no conflicts of interest.

Ahmad J. Obaidullah, Mohammed M. Alanazi, Nawaf A. Alsaif, Ashwag S. Alanazi, Hussam Albassam, Osama I. Alwassil, and Alanazi AZ conducted the study and performed the corresponding data analysis, prepared figures, and drafted the manuscript. Abu Montakim Tareq, Mohammed M. Alanazi, and Ahmad J. Obaidullah supervised the study, performed the data analysis, prepared the figure, and contributed to writing with Shafi Mahmud and Ali M. Alqahtani. Ahmad J. Obaidullah, Mohammed M. Alanazi, and Abu Montakim Tareq were responsible for writing, review, and editing of the manuscript. All the authors approved the final version of the manuscript. Table S1 : biological processes of the genes targeted by the compounds; Table S2 : molecular functions of the genes interacting with the compounds; Table S3 : protein-protein interaction (PPI) status of 10 proteins with coexpression; Table S4 : molecular interactions with the best ligand efficiencies for compounds binding with IL1R1 (PDB: 1ITB); Table S5 : molecular interactions with the best ligand efficiencies for compounds binding with IRAK4 (PDB: 6EGA); Table S6 : molecular interactions with the best ligand efficiencies for compounds binding with MYD88 (4EO7); Table  S7 : molecular interactions with the best ligand efficiencies for compounds binding with TIRAP (PDB: 4FZ5); Table S8 : molecular interactions with the best ligand efficiencies for compounds binding with TLR4 (PDB: 3FXI); and Table S9 : molecular interactions with the best ligand efficiencies for compounds binding with TRAF6 (PDB: 3HCT). Figure S1 : GC-MS chromatograph of the methanolic extract of Argyreia capitiformis stem; Figure S2 : protein-protein interaction (PPI) network of 10 proteins; Figure S3 : gene cooccurrence of 10 proteins for Homo sapiens; Figure S4 : gene coexpression of 10 proteins in PPI network; and Figure S5 : two-dimensional representations of the best ligand efficiencies for compounds binding with IL1R1 (PDB: 1ITB). 3-Furaldehyde (A), 4(1H)-pyrimidinone, 6-hydroxy-(B), catechol (C), hydroquinone (D), phenol (E), sulcatone (F), and aspirin (G) are shown. Figure S6 : two-dimensional representations of the best ligand efficiencies for compounds binding to IRAK4 (PDB: 6EGA). 3-Furaldehyde (A), 4methyl-1,5-heptadiene (B), catechol (C), methylcyclohexane (D), p-vinylguaiacol (E), phenol (F), and aspirin (G) are shown. Figure S7 : two-dimensional representations of the best ligand efficiencies for compounds binding to MYD88 (4EO7). (Z,E)-farnesol (A), 2-(2-hydroxy-2-phenylethyl)-3,5,6-trimethylpyrazine (B), 4-methyl-1,5-heptadiene (C), coniferol (D), hydroquinone (E), trans-13-docosenamide (F), and aspirin (G) are shown. Figure S8 : two-dimensional representations of the best ligand efficiencies for compounds binding to TIRAP (PDB: 4FZ5). 3-Furaldehyde (A), 4ethylresorcinol (B), catechol (C), hydroquinone (D), methylcyclohexane (E), phenol (F), and aspirin (G) are shown. Figure S9 : two-dimensional representations of the best ligand efficiencies for compounds against TLR4 (PDB: 3FXI). 4-Methyl-1,5-heptadiene (A), catechol (B), coniferol (C), hydroquinone (D), phenol (E), sulcatone (F), and aspirin (G) are shown. Figure S10 : two-dimensional representations of the best ligand efficiencies for compounds against TRAF6 (PDB: 3HCT). 2-(2-Hydroxy-2-phenylethyl)-3,5,6-trimethylpyrazine (A), 3-furaldehyde (B), coniferol (C), coumaran (D), phenol (E), sulcatone (F), and aspirin (G) are shown. (Supplementary Materials) 

Alternative activation of macrophages: an immunologic functional perspective

Antagonistic crosstalk between NF-κB and SIRT1 in the regulation of inflammation and metabolic disorders

Inflammatory mechanisms: the molecular basis of inflammation and disease

Mitochondrial dysfunction and multiple sclerosis

Effects of short-and long-term exposures to particulate matter on inflammatory marker levels in the general population

e role of the radiologist in determining disease severity in inflammatory bowel diseases

Triptolide attenuates inflammatory response in membranous glomerulo-nephritis rat via downregulation of NF-κB signaling pathway

Medicinal plants with anti-inflammatory activities from selected countries and regions of Africa

Evaluation of anti-inflammatory effect by regulating NF-κB pathway of Argyreia capitata (Vahl) Choisy extract in LPS-induced RAW 264.7 macrophages

Ethnomedicinal plants used in sexual disorder in Balangir and Deogarh districts, Odisha, India

Medicinal uses, pharmacology, and phytochemistry of Convolvulaceae plants with central nervous system efficacies: a systematic review

Phytochemical and pharmacological studies on Convolvulus fatmensis Ktze

Argyreia speciosa (Linn. f.) sweet: a comprehensive review

A literature review on Argyreia nervosa (burm. F.) bojer

Investigations into the immunomodulatory activity of Argyreia speciosa

Network pharmacological identification of active compounds and potential actions of Erxian decoction in alleviating menopause-related symptoms

Traditional Chinese medicine network pharmacology: theory, methodology and application

Network pharmacology-based approaches capture essence of Chinese herbal medicines

Network pharmacology-based identification of potential targets of the flower of Trollius chinensis Bunge acting on anti-inflammatory effects

Multiscale modeling of inflammation-induced tumorigenesis reveals competing oncogenic and on coprotective roles for inflammation

Comparative analysis of EPA and DHA in fish oil nutritional capsules by GC-MS

Recent trends in the application of chromatographic techniques in the analysis of luteolin and its derivatives

Antioxidant, biomolecule oxidation protective activities of Nardostachys jatamansi DC and its phytochemical analysis by RP-HPLC and GC-MS

GC-MS analysis of phytoconstituents from Amomum nilgiricum and molecular docking interactions of bioactive serverogenin acetate with target proteins

Chemical profiles, pharmacological properties, and in silico studies provide new insights on Cycas pectinata

Effects of organic extracts of six Bangladeshi plants on in vitro thrombolysis and cytotoxicity

Deeper insights on Cnesmone javanica blume leaves extract: chemical profiles, biological attributes, network pharmacology and molecular docking

Network pharmacology-based assessment to elucidate the molecular mechanism of anti-diabetic action of Tinospora cordifolia

Pharmacological potential and chemical characterization of bridelia ferruginea benth-a native tropical african medicinal plant

STRING v10: protein-protein interaction networks, integrated over the tree of life

Crystal structure of the type-I interleukin-1 receptor complexed with interleukin-1β

Conformational flexibility and inhibitor binding to unphosphorylated interleukin-1 receptor-associated kinase 4 (IRAK4)

Molecular mechanisms for the subversion of MyD88 signaling by TcpC from virulent uropathogenic Escherichia coli

X-ray crystallographic structure of TIR-domain from the human TIR-domain containing adaptor protein/MyD88-adaptor-like protein (TIRAP/MAL)

e structural basis of lipopolysaccharide recognition by the TLR4-MD-2 complex

E2 interaction and dimerization in the crystal structure of TRAF6

Structure-based virtual screening for drug discovery: principles, applications and recent advances

Assessing the performance of the MM/PBSA and MM/GBSA methods. 6. Capability to predict protein-protein binding free energies and re-rank binding poses generated by protein-protein docking

Drug efficiency indices for improvement of molecular docking scoring functions

Development and testing of a general amber force field

YASARA: a tool to obtain structural guidance in biocatalytic investigations

Structure and dynamics of TIP3P, TIP4P, and TIP5P water near smooth and atomistic walls of different hydroaffinity

An implementation of the smooth particle mesh ewald method on GPU hardware

A smooth particle mesh Ewald method

New ways to boost molecular dynamics simulations

Virtual screening and molecular dynamics simulation study of plantderived compounds to identify potential inhibitors of main protease from SARS-CoV-2

Extensive immunoinformatics study for the prediction of novel peptide-based epitope vaccine with docking confirmation against envelope protein of Chikungunya virus: a computational biology approach

Comparative molecular investigation of the potential inhibitors against SARS-CoV-2 main protease: a molecular docking study

Fermentation optimization of cellulase production from sugarcane bagasse by Bacillus pseudomycoides and molecular modeling study of cellulase

Drug repurposing approach against novel coronavirus disease (COVID-19) through virtual screening targeting SARS-CoV-2 main protease

In silico design of an epitope-based vaccine against choline binding protein A of Streptococcus pneumoniae

A molecular modelling approach for identifying antiviral selenium-containing heterocyclic compounds that inhibit the main protease of SARS-CoV-2: an in silico investigation

Molecular docking and dynamics study of natural compound for potential inhibition of main protease of SARS-CoV-2

Chemical profiles and pharmacological properties with in silico studies on elatostema papillosum wedd

Historical review of medicinal plants' usage

Modern Phytomedicine: Turning Medical Plants into Drugs

GC-MS analysis and anti-microbial activity of Psidium guajava (leaves) grown in Malva region of India

GC-MS analysis of phytocompounds of leaf and stem of Marsilea quadrifolia (L.)

Screening and GC-MS profiling of ethanolic extract of Tylophora pauciflora

From traditional Chinese medicine to rational cancer therapy

Phytomedicine in otorhinolaryngology and pulmonology: clinical trials with herbal remedies

Protein kinase C and tolllike receptor signaling

STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets

Functional coverage of the human genome by existing structures, structural genomics targets, and homology models

Towards a knowledge-based human protein atlas

UniProt: e universal protein knowledgebase

Structure and function of Toll receptors and their ligands

Signaling to NF-κB by toll-like receptors

e role of pattern-recognition receptors in innate immunity: update on Toll-like receptors

Predicting protein function by genomic context: quantitative evaluation and qualitative inferences

Computational prediction of protein-protein interactions

NCBI GEO: archive for functional genomics data sets-update

Virtual screening strategies in drug discovery: a critical review

Structure-based virtual screening for drug discovery: a problem-centric review

Immunological and inflammatory functions of the interleukin-1 family

Mouse neutrophils express the decoy type 2 interleukin-1 receptor (IL-1R2) constitutively and in acute inflammatory conditions

IL-1 receptor type 2 suppresses collagen-induced arthritis by inhibiting IL-1 signal on macrophages

Mass spectrometric analysis of the endogenous type I interleukin-1 (IL-1) receptor signaling complex formed after IL-1 binding identifies IL-1RAcP, MyD88, and IRAK-4 as the stable components

Design and synthesis of small molecule interleukin-1 receptor antagonists based on a benzene template

e interleukin-1 receptor/Toll-like receptor superfamily: 10 years of progress

Pathogen recognition and innate immunity

Mal (MyD88-adapter-like) is required for toll-like receptor-4 signal transduction

TIRAP: an adapter molecule in the Toll signaling pathway

e adaptor molecule TIRAP provides signalling specificity for toll-like receptors

Essential role for TIRAP in activation of the signalling cascade shared by TLR2 and TLR4

Kaempferol is an anti-inflammatory compound with activity towards NF-κB pathway proteins

Anthocyanins from Hibiscus syriacus L. attenuate LPS-induced inflammation by inhibiting the TLR4-mediated NF-κB signaling pathway

TLR4-Targeting therapeutics: structural basis and computer-aided drug discovery approaches