Die Wirkstoffe von RD
Insgesamt wurden 71 rohe chemische Inhaltsstoffe aus TCMSP mit „Rhizoma Drynariae“ als Stichwort abgerufen. Eine sekundäre Selektion erfolgte gemäß ADME-Kriterien (OB ≥ 30 %, DL ≥ 0,18), wobei schließlich 18 bioaktive Inhaltsstoffe mit günstiger oraler Bioverfügbarkeit und Arzneimitteltauglichkeit erhalten wurden. Tabelle 1 fasst detaillierte Grundinformationen zu diesen 18 Wirkstoffen zusammen. Diese 18 Verbindungen bilden die mehrkomponentige stoffliche Grundlage von RD gegen RCI.
Relevante Genziele von RD
In den Datenbanken TCMSP, ETCM und SymMap wurden nach Zielgenen gesucht, die zu den wirksamen Inhaltsstoffen von RD passen, und insgesamt 327 ursprüngliche Zielgeneinträge gesammelt. Nach Integration der Daten und Entfernung wiederholter Gensymbole wurden 213 nicht redundante potenzielle Zielgene von RD für die anschließende Schnittmengenanalyse gewonnen.
Relevante Genziele für RCI
In den Datenbanken Genecards und OMIM wurde mit dem Stichwort „rotator cuff injury“ gesucht, und es wurden 989 rohe, mit RCI assoziierte Gen-Einträge erhalten. Nach Datenintegration und Entfernung von Duplikaten wurden 870 eindeutige humane pathogene Gene identifiziert, die eng mit dem Fortschreiten einer RCI assoziiert sind.
Aufbau des PPI-Netzwerks
Die Schnittmengenanalyse zwischen 213 RD-Zielgenen und 870 RCI-pathogenen Genen ergab 71 überlappende Gene, die als potenzielle therapeutische Zielstrukturen von RD für die Behandlung von RCI definiert wurden. Ein Venn-Diagramm zur Darstellung der Schnittmengenbeziehung wurde mit der Software Venny 2.0 erstellt (Abbildung 2). Die 71 überlappenden Zielgene wurden in STRING importiert, um das PPI-Netzwerk von RD, das in RCI eingreift, aufzubauen (Abbildung 3). Die vollständigen PPI-Daten wurden in Cytoscape importiert, und cytoHubba identifizierte die 10 wichtigsten Hub-Gene mit den höchsten MCC-Scores: TNF, IL6, VEGFA, PTGS2, CASP3, AKT1, CCL2, PPARG, MAPK3, HMOX1. Ihre topologischen Informationen bezüglich des MCC-Werts sind in Tabelle 2 zusammengefasst und in Abbildung 4 visualisiert. Die 71 überschneidenden Gene repräsentieren die zentralen Zielstrukturen, die die therapeutische Wirkung von RD gegen RCI vermitteln.
GO- und KEGG-Anreicherungsanalyse
Um den molekularen Wirkmechanismus der RD-Kandidatenziele bei der Behandlung von RCI zu untersuchen, wurde mit dem clusterProfiler-Paket eine GO- und KEGG-Enrichment-Analyse an den überlappenden therapeutischen Zielstrukturen durchgeführt. Separate GO-Enrichment-Statistiken wurden für BP, CC und MF durchgeführt. Abbildungen 5A–B, 6A–B und 7A–B zeigen die 10 signifikant angereicherten Einträge für BP, CC und MF; in jedem Teilbild B ist die Genverteilung für jeden funktionellen Begriff dargestellt. Eine weitere KEGG-Enrichment-Analyse wurde durchgeführt, um die zentralen Signalwege aufzudecken, die an der RD-bedingten Intervention bei RCI beteiligt sind; die 30 am stärksten angereicherten Signalwege sind in Abbildung 8 visualisiert.
Molekulare Docking-Simulation
Die 10 wichtigsten Hub-Gene wurden als Proteinrezeptoren festgelegt, und alle 18 RD-Wirkstoffe wurden als Liganden für die Docking-Berechnung mittels AutoDock Vina und Open Babel verwendet. PyMOL visualisierte die Bindungskonformationen von Ligand und Rezeptor. Die Docking-Ergebnisse zeigten, dass die meisten Hub-Zielproteine stabile Bindungsstrukturen mit den RD-Wirkstoffen bildeten. Tabelle 3 und Abbildung 9 enthalten die vollständigen Docking-Parameter sowie dreidimensionale Bindungsdiagramme der 20 Komplexe mit der besten Bindungsaffinität.
In der vorliegenden Netzwerkpharmakologie-Analyse wurden 18 bioaktive Inhaltsstoffe von RD und 213 nicht redundante arzneimittelbezogene Zielstrukturen ermittelt. Nach dem Abgleich mit 870 RCI-assoziierten Genen wurden 71 gemeinsame therapeutische Kandidatenziele identifiziert. Die Konstruktion des PPI-Netzwerks und die cytoHubba-MCC-Auswertung ergaben zehn zentrale Gene: TNF, IL6, VEGFA, PTGS2, CASP3, AKT1, CCL2, PPARG, MAPK3 und HMOX1. Die GO-Enrichment-Analyse umfasste biologische Prozesse (BP), zelluläre Komponenten (CC) und molekulare Funktionen (MF), während die KEGG-Enrichment-Analyse mehrere angereicherte Signalwege identifizierte. Die molekulare Docking-Simulation bestätigte eine günstige Bindungsaffinität zwischen den wichtigsten aktiven Inhaltsstoffen von RD und diesen zentralen Zielproteinen.
VERFÜGBARKEIT VON DATEN:
Alle in dieser Studie generierten Roh- und verarbeiteten Daten, einschließlich Ziel-Listen, PPI-Netzwerkdateien, GO- und KEGG-Anreicherungsergebnisse sowie die molekulare Docking-Matrix, sind in der Science Data Bank hinterlegt. Der Zugangslink zu den Rohdaten lautet https://www.scidb.cn/anonymous/alVOVlJu.

Abbildung 1: Schematischer Arbeitsablauf zur Darstellung der in dieser Studie angewendeten bioinformatischen Pipeline. Beginnend links wurden Kandidaten bioaktiver Verbindungen und ihre entsprechenden RD-Ziele aus den Datenbanken TCMSP, ETCM und SymMap abgerufen. Auf der rechten Seite wurden Gene, die mit RCI assoziiert sind, aus GeneCards und OMIM extrahiert. Die gemeinsamen Gengesichtspunkte dieser beiden Datensätze wurden mithilfe eines Venn-Diagramms identifiziert und als potenzielle therapeutische Zielstrukturen definiert. Die gemeinsamen Zielgene wurden anschließend in zwei getrennten Zweigen weiterverarbeitet: Ein Zweig erstellte ein Protein-Protein-Interaktionsnetzwerk unter Verwendung von STRING und Cytoscape; der andere Zweig führte GO- und KEGG-Enrichment-Analysen mit dem R-Paket clusterProfiler durch. Die aus dem PPI-Netzwerk extrahierten Hub-Gene wurden anschließend durch molekulares Docking validiert, wobei die strukturelle Vorbereitung und Simulation mit PubChem, PDB, AutoDock Vina und PyMOL durchgeführt wurde. Abkürzungen: RD = Rhizoma Drynariae, RCI = rotator cuff injury, TCMSP = The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform, ETCM = The Encyclopedia of Traditional Chinese Medicine, OMIM = The Online Mendelian Inheritance in Man, GO = gene ontology, KEGG = Kyoto encyclopedia of genes and genomes. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 2: Mengendiagramm der überlappenden Zielstrukturen zwischen RD und RCI. Der blaue Kreis stellt die mit RCI assoziierten Zielstrukturen dar. Der gelbe Kreis stellt die potenziellen therapeutischen Zielstrukturen von RD dar. Der überlappende olivgrüne Bereich zeigt die gemeinsamen Zielstrukturen zwischen RD und RCI an. Die Zahlen und Prozentsätze in den einzelnen Bereichen geben die Anzahl und den Anteil eindeutiger oder gemeinsamer Zielstrukturen an. Abkürzungen: RD = Rhizoma Drynariae, RCI = Rotatorenmanschettenverletzung. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 3:PPI-Netzwerk, erstellt aus den 71 überlappenden Genzielen, die RD und RCI verbinden. Jeder Kreis im Netzwerk repräsentiert ein Zielgen, beschriftet mit seinem offiziellen Gensymbol. Linien, die zwei Kreise verbinden, veranschaulichen Protein-Protein-Assoziationen, die aus öffentlichen bioinformatischen Datenbanken gewonnen wurden. Die Knotengröße und die Farbsättigung spiegeln die topologische Bedeutung einzelner Gene innerhalb des Netzwerks wider. Gene mit größeren Kreisen und intensiveren Grüntönen weisen höhere topologische Werte auf und nehmen zentralere Positionen ein; kleinere, blassgelbe Kreise kennzeichnen Gene mit geringerer topologischer Bedeutung. Abkürzungen: RD = Rhizoma Drynariae, RCI = Rotator-Manschette-Verletzung, PPI = Protein-Protein-Interaktion. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 4: Die 10 wichtigsten Hub-Gene, identifiziert im PPI-Netzwerk mithilfe des cytoHubba-Plugins. Rechteckige Kästen stellen Hub-Gene dar, die jeweils mit ihrem offiziellen Gensymbol beschriftet sind. Schwarze Linien zwischen den Kästen zeigen Proteinwechselwirkungen zwischen diesen zentralen Zielstrukturen. Der Farbverlauf von hellrot über orange bis gelb veranschaulicht die Rangfolge der Hub-Gene entsprechend den von cytoHubba berechneten MCC-Scores; hellrote Kästen kennzeichnen Gene mit den höchsten MCC-Scores und somit der größten Netzwerkbedeutung, während gelbe Kästen Gene mit relativ niedrigeren Werten anzeigen. Abkürzungen: PPI = Protein-Protein-Interaktion, MCC = Maximal Clique Centrality.Bitte klicken Sie hier, um eine vergrößerte Ansicht dieser Abbildung einzusehen.

Abbildung 5: BP-Anreicherungsanalyse der 71 gemeinsamen Zielgene von RD im Vergleich zu RCI. (A) GO-Anreicherungsanalyse potenzieller therapeutischer Zielstrukturen für BP. Es handelt sich um ein horizontales Balkendiagramm, das die 10 signifikant angereicherten GO-Biologische-Prozesse-Terme zeigt. Die horizontale X-Achse stellt die Genanzahl dar, die der Anzahl der Zielgene entspricht, die jedem GO-Term zugeordnet sind. Die vertikale Y-Achse listet einzelne angereicherte GO-Terme auf. Die kontinuierliche Farbleiste rechts kodiert -log10(p-Wert); wärmere Rottöne zeigen einen höheren -log10(p-Wert) und eine stärkere statistische Signifikanz an, während kühlere Violett- und Blautöne eine geringere Signifikanz anzeigen. (B) Beziehung zwischen potenziellen therapeutischen Zielstrukturen und BP. Panel B zeigt ein Chord-Diagramm, das die paarweisen Assoziationen zwischen angereicherten GO-Termen und den entsprechenden Zielgenen veranschaulicht. Die äußeren kreisförmigen Segmente sind in zwei Kategorien unterteilt: farbige Blöcke stehen für unterschiedliche GO-Biologische-Prozesse-Terme, und graue rechteckige Markierungen kennzeichnen einzelne Zielgen-Symbole. Gebogene Bänder innerhalb des Kreises verbinden GO-Terme mit ihren zugehörigen Zielgenen, wobei jedes Band ein Gen mit dem GO-Term verknüpft, an dem es beteiligt ist. Die Farbe jedes Bandes entspricht dem Farbton des entsprechenden GO-Terms aus dem äußeren Ring. Eine Legende unten ordnet jeder Farbe den jeweiligen GO-Term zu. Abkürzungen: RD = Rhizoma Drynariae, RCI = Rotatorenmanschettenverletzung, BP = biologischer Prozess, GO = Genontologie. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.

Abbildung 6: GO-CC-Terme, abgeleitet von den 71 Kandidatentargets von RD bei der Behandlung von RCI. (A) GO-Anreicherungsanalyse potenzieller therapeutischer Targetstrukturen für CC. (B) Beziehung zwischen potenziellen therapeutischen Targetstrukturen und CC. Panel A zeigt ein horizontales Balkendiagramm mit den 10 signifikant angereicherten GO-Zellkomponenteneinträgen. Die horizontale X-Achse quantifiziert die Anzahl der Zielgene, die jedem GO-Term zugeordnet sind, gekennzeichnet als „Counts“. Die Beschreibungen der angereicherten GO-Zellkomponenten sind entlang der vertikalen Y-Achse angeordnet. Die Farbleiste rechts ordnet die Farbvariation den -log10(pvalue)-Werten zu; helles Rot entspricht höheren -log10(pvalue)-Werten und stärkerer statistischer Signifikanz, während blaue Töne geringere Signifikanzniveaus darstellen. Panel B zeigt ein Chord-Diagramm, das Verbindungen zwischen angereicherten GO-Zellkomponententermen und ihren assoziierten Zielgenen veranschaulicht. Segmente am äußeren Ring sind in zwei Gruppen unterteilt: farbige Blöcke stehen für unterschiedliche GO-Zellkomponententerme, und graue Markierungen bezeichnen individuelle Gen-Symbole der Zielgene. Geschwungene Bänder innerhalb des Kreises verbinden jedes Gen mit seinem zugehörigen GO-Term, wobei jedes Band die Farbe seines entsprechenden GO-Terms aus dem äußeren Ring übernimmt. Eine Legende unten erklärt die Farbcodierung aller GO-Zellkomponentenkategorien. Abkürzungen: RD = Rhizoma Drynariae, RCI = Rotatorenmanschettenverletzung, CC = Zellkomponente, GO = Genontologie Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 7: Analysen der Anreicherung molekularer Funktionen (MF) anhand der 71 überlappenden Zielgene von RD bei RCI-Behandlung. (A) GO-Anreicherungsanalyse potenzieller therapeutischer Zielstrukturen für molekulare Funktionen (MF). (B) Zusammenhang zwischen potenziellen therapeutischen Zielstrukturen und molekularen Funktionen (MF). Panel A zeigt ein horizontales Balkendiagramm, das die 10 signifikant angereicherten GO-Terme für molekulare Funktionen auflistet. Die Werte auf der horizontalen X-Achse repräsentieren die Genanzahl, die der Gesamtanzahl der Zielgene entspricht, die jeweils einem Eintrag für molekulare Funktionen zugeordnet sind. Auf der vertikalen Y-Achse sind die einzelnen Beschreibungen der GO-Terme für molekulare Funktionen aufgeführt. Die farbige Gradientenleiste rechts veranschaulicht die Farbunterschiede als -log10(p-Wert)-Maßstab; hellrot steht für höhere -log10(p-Wert)-Werte und stärkere statistische Signifikanz, während blaue Töne eine geringere Signifikanz anzeigen. Panel B zeigt ein Sehnendiagramm (Chord-Diagramm), das die Beziehungen zwischen angereicherten GO-Kategorien für molekulare Funktionen und den dazugehörigen Zielgenen darstellt. Die Segmente im äußeren Ring enthalten zwei Arten von Beschriftungen: farbige Blöcke kennzeichnen unterschiedliche GO-Terme für molekulare Funktionen, und kleine graue rechteckige Markierungen zeigen die einzelnen Symbole der Zielgene. Geschwungene Bänder, die vom äußeren Ring nach innen verlaufen, verbinden jedes Gen mit dem zugehörigen Funktionsbegriff; dabei übernimmt jedes Band die Farbe des entsprechenden GO-Terms. Eine Legende unten erklärt die Farbcodierung aller im Diagramm dargestellten GO-Kategorien für molekulare Funktionen. Abkürzungen: RD = Rhizoma Drynariae, RCI = Rotatorenmanschettenverletzung, MF = molekulare Funktion, GO = Genontologie Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.

Abbildung 8: KEGG-Pfad-Anreicherungsanalyse der 71 gemeinsamen Zielgene von RD im Vergleich zu RCI. Dieses Punktdiagramm veranschaulicht die am stärksten angereicherten KEGG-Pfade. Die vertikale Y-Achse listet die Namen der einzelnen angereicherten Pfade auf. Die horizontale X-Achse stellt das Genverhältnis (GeneRatio) dar, definiert als das Verhältnis der Zielgene, die einem Pfad zugeordnet wurden, zur Gesamtanzahl der Gene in diesem Pfad. Jeder Kreispunkt steht für einen angereicherten Pfad. Die Punktgröße entspricht der Genanzahl, wobei größere Kreise auf mehr Zielgene, die dem Pfad zugeordnet sind, hinweisen. Die Punktfarbe spiegelt den p-Wert wider: Rot deutet auf niedrigere p-Werte und stärkere statistische Signifikanz hin, während tiefes Blau höhere p-Werte und schwächere statistische Signifikanz anzeigt. Zwei seitliche Legenden definieren jeweils die Skala für die Genanzahl und die Farbkodierung des p-Werts. Abkürzungen: RD = Rhizoma Drynariae, RCI = Rotatorenmanschettenverletzung, KEGG = Kyoto Encyclopedia of Genes and Genomes. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.

Abbildung 9: Die am höchsten bewerteten Docking-Komplexe. (A) (PTGS2, Eriodictyol, -9,7 kcal/mol). (B) (PTGS2, Luteolin, -9,6 kcal/mol). (C) (PTGS2, Catechin, -9,5 kcal/mol). (D) (PTGS2, Naringenin, -9,5 kcal/mol). (E) (PTGS2, Aureusidin, -9,4 kcal/mol). (F) (PTGS2, Kaempferol, -9,2 kcal/mol). (G) (PTGS2, Chroman, -9,1 kcal/mol). (H) (TNF, Catechin, -9,1 kcal/mol). (I) (PTGS2, Digallat, -9 kcal/mol). (J) (TNF, Digallat, -9 kcal/mol). (K) (PPARG, Chroman, -8,9 kcal/mol). (L) (PPARG, Naringenin, -8,9 kcal/mol). (M) (TNF, Aureusidin, -8,9 kcal/mol). (N) (TNF, Kaempferol, -8,8 kcal/mol). (O) (MAPK3, Chroman, -8,7 kcal/mol). (P) (MAPK3, Naringenin, -8,7 kcal/mol). (Q) (PPARG, Aureusidin, -8,7 kcal/mol). (R) (PPARG, Kaempferol, -8,7 kcal/mol). (S) (PTGS2, 22_Stigmasten_3_one, -8,7 kcal/mol). (T) (MAPK3, Luteolin, -8,6 kcal/mol). Die Panels A bis T zeigen jeweils die berechnete Bindungsgeometrie einzelner kleiner Moleküle innerhalb der Bindungstasche der entsprechenden Proteintargets. Cyanfarbene Stabmodelle stellen die Aminosäurereste dar, die die Proteinbindungstasche auskleiden. Kleine Molekül-Liganden sind als magenta- oder grüne Stäbchen dargestellt, um unterschiedliche chemische Strukturen zu unterscheiden. Gelbe gestrichelte Linien markieren intermolekulare Wasserstoffbrückenbindungen zwischen Liganden und Taschenresten, wobei numerische Anmerkungen die Abstände der Wasserstoffbrücken in Ångström (Å) angeben. Die durchscheinende weiße Oberfläche skizziert die gesamte dreidimensionale Topologie jedes Proteinrezeptors. Alle Berechnungen der Bindungspositionen wurden mittels starrem Rezeptor-Moleküldocking mit der AutoDock-Vina-Software durchgeführt. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
| Mol-ID | Molekülname | orale Bioverfügbarkeit (OB) (%) | Arzneimittelähnlichkeit (DL) |
| MOL001040 | (2R)-5,7-Dihydroxy-2-(4-hydroxyphenyl)chroman-4-on | 42,36 | 0,21 |
| MOL001978 | Aureusidin | 53,42 | 0,24 |
| MOL002914 | Eriodyctiol (Flavanon) | 41,35 | 0,24 |
| MOL000449 | Stigmasterol | 43,83 | 0,76 |
| MOL000358 | Beta-Sitosterol | 36,91 | 0,75 |
| MOL000422 | Kaempferol | 41,88 | 0,24 |
| MOL004328 | Naringenin | 59,29 | 0,21 |
| MOL000492 | (+)-Catechin | 54,83 | 0,24 |
| MOL005190 | Eriodictyol | 71,79 | 0,24 |
| MOL000569 | Digallat | 61,85 | 0,26 |
| MOL000006 | Luteolin | 36,16 | 0,25 |
| MOL009061 | 22-Stigmasten-3-on | 39,25 | 0,76 |
| MOL009063 | Cyclolaudenolacetat | 41,66 | 0,79 |
| MOL009075 | Cycloartenon | 40,57 | 0,79 |
| MOL009076 | Cyclolaudenol | 39,05 | 0,79 |
| MOL009078 | Davalliosid A_qt | 62,65 | 0,51 |
| MOL009087 | Mariosid_qt | 70,79 | 0,19 |
| MOL009091 | Xanthogalenol | 41,08 | 0,32 |
Tabelle 1: Pharmakokinetische Profile von achtzehn Kandidaten bioaktiver Inhaltsstoffe aus RD. Verbindungen wurden aus der TCMSP-Datenbank abgerufen und anhand von zwei ADME-bezogenen Schwellenwerten gefiltert: OB ≥ 30 % und DL ≥ 0,18. Die Spalten enthalten die TCMSP-Zugangsnummer (Mol ID), den Namen der Verbindung, den berechneten OB-Wert (%) sowie den DL-Score jedes untersuchten Inhaltsstoffs. Abkürzungen: RD = Rhizoma Drynariae, RCI = Rotatorenmanschettenverletzung, ADME = Absorption, Verteilung, Metabolismus, Exkretion, TCMSP = The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform, OB = orale Bioverfügbarkeit, DL = Arzneimittelähnlichkeit.
| Rang | Gen | Name | Punktzahl |
| 1 | TNF | Tumor-Nekrose-Faktor | 328291218366500 |
| 2 | IL6 | Interleukin 6 | 328291215913592 |
| 3 | VEGFA | Wachstumsfaktor A für vaskuläre endotheliale Zellen | 328291215896284 |
| 4 | PTGS2 | Prostaglandin-Endoperoxid-Synthase 2 | 328288656326402 |
| 5 | CASP3 | Caspase 3 | 328288180323600 |
| 6 | AKT1 | AKT-Serin/Threonin-Kinase 1 | 328287745700516 |
| 7 | CCL2 | C-C-Motiv-Chemokin-Ligand 2 | 328161521686320 |
| 8 | PPARG | Peroxisom-Proliferator-aktivierter Rezeptor Gamma | 325302724852296 |
| 9 | MAPK3 | mitogenaktivierte Proteinkinase 3 | 308779166166218 |
| 10 | HMOX1 | Häm-Oxygenase 1 | 299171945318400 |
Tabelle 2: Die 10 wichtigsten Hub-Gene, identifiziert aus dem PPI-Netzwerk. Die Gene wurden nach MCC-Scores geordnet, die mit dem cytoHubba-Plugin in Cytoscape berechnet wurden. MCC ist ein topologisches Maß von cytoHubba, bei dem höhere numerische Werte eine größere topologische Bedeutung innerhalb des Netzwerks widerspiegeln. Die Spalten zeigen die Rangfolge, das offizielle Gen-Symbol, den vollständigen Gen-Namen und den entsprechenden rohen MCC-Score für jedes Hub-Gen an. Abkürzungen: PPI = Protein-Protein-Interaktion, MCC = Maximal Clique Centrality.
| Rezeptoren | Liganden | Freie Energie (kcal/mol) | Entsprechende Seriennummern in Abbildung 9 |
| PTGS2 | eriodictyol | -9.7 | A |
| PTGS2 | luteolin | -9.6 | B |
| PTGS2 | catechin | -9.5 | C |
| PTGS2 | naringenin | -9.5 | D |
| PTGS2 | Aureusidin | -9.4 | E |
| PTGS2 | kaempferol | -9.2 | F |
| PTGS2 | chroman | -9.1 | G |
| TNF | catechin | -9.1 | H |
| PTGS2 | digallate | -9 | I |
| TNF | digallate | -9 | J |
| PPARG | chroman | -8.9 | K |
| PPARG | naringenin | -8.9 | L |
| TNF | Aureusidin | -8.9 | M |
| TNF | kaempferol | -8.8 | N |
| MAPK3 | chroman | -8.7 | O |
| MAPK3 | naringenin | -8.7 | P |
| PPARG | Aureusidin | -8.7 | Q |
| PPARG | kaempferol | -8.7 | R |
| PTGS2 | 22_Stigmasten_3_one | -8.7 | S |
| MAPK3 | luteolin | -8.6 | T |
Tabelle 3: Ergebnisse der molekularen Docking-Bindungsaffinität für die 20 ausgewählten Rezeptor-Ligand-Komplexe. Jede Zeile enthält das Rezeptorprotein, den kleinen molekularen Liganden, die berechnete Bindungsenergie (kcal/mol) sowie die entsprechende Unterdiagramm-Bezeichnung in Abbildung 9. Je negativer die Werte der freien Energie sind, desto stärker ist die Bindungsstabilität zwischen der Proteintasche und der Ligandverbindung. Abkürzungen: PTGS2 = Prostaglandin-Endoperoxid-Synthase 2, TNF = Tumornekrosefaktor, PPARG = Proliferator-aktivierter Rezeptor gamma, MAPK3 = mitogenaktivierter Proteinkinase 3.