# Network pharmacology refined with non-ubiquity and decoy-controlled molecular docking reveals insights into Moringa oleifera phytochemicals targeting insulin resistance

**Authors:** Armi Katrina Santos-Enriquez, Fabian M. Dayrit, Armando Jerome de Jesus, Nina Rosario L. Rojas

PMC · DOI: 10.3389/fbinf.2026.1756081 · Frontiers in Bioinformatics · 2026-03-10

## TL;DR

This study uses a refined network pharmacology approach to identify Moringa oleifera compounds that may help treat insulin resistance.

## Contribution

A novel network pharmacology method combining non-ubiquity filtering and decoy-controlled docking to identify Moringa oleifera phytochemicals targeting insulin resistance.

## Key findings

- Moringa oleifera phytochemicals like boldione, aurantiamide acetate, and flavonoids were predicted to target insulin resistance proteins.
- Network centrality measures identified key targets such as PPARα, PPARγ, and PTPN1.
- The approach filtered out ubiquitous compounds to highlight Moringa-specific effects on insulin resistance.

## Abstract

Moringa oleifera phytochemicals were predicted to target insulin resistance proteins using a modified network pharmacology and molecular docking approach. Two hundred ninety M. oleifera phytochemicals with their aglycones, acetylase and myrosinase degradation products were compiled from literature and phytochemical databases. Nine protein targets were identified from the intersection of gene lists with high relevance to insulin resistance from GeneCards and DisGeNET and the target genes predicted by reverse screening using Swiss Target Prediction: protein-tyrosine phosphatase 1B (PTPN1), 11-beta-hydroxysteroid dehydrogenase 1 (HSD11B1), peroxisome proliferator-activated receptor α (PPARα), peroxisome proliferator-activated receptor γ (PPARγ), PI3-kinase p85-alpha subunit (PIK3R1), insulin receptor (INSR), tumor necrosis factor α (TNF-α), endothelial nitric oxide synthase (eNOS) and hepatic lipase (LIPC). Binding affinities of phytochemicals with the targets were predicted using Autodock Vina. The predicted binding affinities were classified according to calculated thresholds using receiver operating characteristic (ROC) calculations of binding affinities of: (a) binders (annotated drugs and other molecules with known interaction with each target), and (b) decoys (molecules not expected to bind to a specific target). In addition, ubiquitous phytochemicals were filtered out to differentiate the effect on insulin resistance of M. oleifera from that of other plants. The resulting phytochemical-protein interaction network was visualized using Cytoscape. All mentioned targets, except hepatic lipase, were key targets based on various network centrality measures. Previous studies on murine models have shown that isothiocyanate-rich M. oleifera extracts ameliorate insulin resistance. Using our approach, the following phytochemicals, with predicted moderate bioavailability, high GI absorption, and probable binding with insulin resistance targets, are recommended for further in vivo or in vitro validation for insulin resistance activity: boldione (a steroid); aurantiamide acetate and aurantiamide (peptide derivatives); O-ethyl-[(3,4-dihydroxyphenyl)methyl] carbamothioate and O-methyl-N-[(4-hydroxyphenyl)methyl] carbamothioate (thiocarbamates); 4α,6α-dihydroxyeudesman-8β,12-olide (a sesquiterpenoid); sanleng acid and tianshic acid (fatty acid derivatives); 2′,5,5′,7-tetrahydroxyflavone; 2′,3,5,7-tetrahydroxyflavone; and 6-hydroxykaempferol (flavonoids). By combining network centrality measures of targets, using ROC-derived thresholds for docking energies, and considering ubiquity of phytochemicals, our refined network pharmacology approach may aid in discovering key bioactive phytochemicals as potential chemical markers for standardization and differentiation of an herbal drug.

Flowchart illustrates a bioinformatics pipeline identifying Moringa oleifera phytochemicals targeting insulin resistance, including literature mining, ADME screening, target prediction, binding affinity docking, classification using ROC, and final annotated network visualization.

## Linked entities

- **Genes:** PTPN1 (protein tyrosine phosphatase non-receptor type 1) [NCBI Gene 5770], HSD11B1 (hydroxysteroid 11-beta dehydrogenase 1) [NCBI Gene 3290], PPARA (peroxisome proliferator activated receptor alpha) [NCBI Gene 5465], PPARG (peroxisome proliferator activated receptor gamma) [NCBI Gene 5468], PIK3R1 (phosphoinositide-3-kinase regulatory subunit 1) [NCBI Gene 5295], INSR (insulin receptor) [NCBI Gene 3643], TNF (tumor necrosis factor) [NCBI Gene 7124], NOS3 (nitric oxide synthase 3) [NCBI Gene 4846], LIPC (lipase C, hepatic type) [NCBI Gene 3990]
- **Chemicals:** boldione (PubChem CID 13472), aurantiamide acetate (PubChem CID 9832120), aurantiamide (PubChem CID 185904), sanleng acid (PubChem CID 5321100), tianshic acid (PubChem CID 5321949), 2′,3,5,7-tetrahydroxyflavone (PubChem CID 5281610), 6-hydroxykaempferol (PubChem CID 5281638)
- **Species:** Moringa oleifera (taxon 3735)

## Full-text entities

- **Diseases:** insulin resistance (MESH:D007333)
- **Chemicals:** aglycones (MESH:C458179), aurantiamide (MESH:C574018), aurantiamide acetate (MESH:C011670), fatty acid (MESH:D005227), 2',3,5,7-tetrahydroxyflavone (-), thiocarbamates (MESH:D013859), steroid (MESH:D013256), isothiocyanate (MESH:C037152), flavonoids (MESH:D005419), boldione (MESH:C006573), sesquiterpenoid (MESH:D012717), peptide (MESH:D010455)
- **Species:** Mus musculus (house mouse, species) [taxon 10090]

## Full text

_Full body text omitted from this summary view._ Fetch the complete paper as Markdown: https://tomesphere.com/paper/PMC13008888/full.md

## Figures

14 figures with captions in the complete paper: https://tomesphere.com/paper/PMC13008888/full.md

## References

179 references — full list in the complete paper: https://tomesphere.com/paper/PMC13008888/full.md

---
Source: https://tomesphere.com/paper/PMC13008888