Alle in dieser Studie verwendeten Materialien, einschließlich Datenbanken, Software und experimentellen Reagenzien, sind in der Materialtabelle aufgeführt. Alle experimentellen Verfahren wurden vom Tierethikkomitee des Suzhou Dushu Lake Hospital (Genehmigungsnummer: 2410008) genehmigt.
Screening aktiver Komponenten und Zielerfassung von DCQD
Die aktiven Bestandteile der enthaltenen Kräuter Rhabarber, Mirabilit, Magnolienrinde und unreifer Bitterorange im Da-Cheng-Qi-Abkok (DCQD) wurden aus der traditionellen Datenbank der chinesischen Medizin (TCMSP) unter Verwendung folgender Kriterien geprüft: orale Bioverfügbarkeit (OB) ≥30 % und Arzneimittelähnlichkeit (DL) ≥0,18. Die entsprechenden Ziele der identifizierten aktiven Komponenten wurden dann gesammelt. Diese Schwellenwerte werden in Netzwerkpharmakologiestudien weit verbreitet verwendet, um Verbindungen mit günstigen pharmakokinetischen Eigenschaften und höherem Potenzial für arzneimittelähnliche Molekülezu filtern 34,35. Um die Zuverlässigkeit der Verbindungsauswahl zu verbessern und Verzerrungen durch das Screening in einer einzelnen Datenbank zu verringern, wurde eine Kreuzvalidierung mit der BATMAN-TCM-Datenbank durchgeführt. BATMAN-TCM ist eine integrative Datenbank, die bekannte und vorhergesagte Wechselwirkungen zwischen TCM-Inhaltsstoffen und Zielproteinen speichert und so die Erforschung pharmakologischer Mechanismen und der Arzneimittelentwicklung unterstützt. Die aktualisierte Version 2.0 bietet einen deutlich erweiterten TTI-Datensatz und erweiterte Funktionen und ist damit eine wertvolle Ressource zum Verständnis der molekularen Mechanismen der TCM und zur Entwicklung neuer Behandlungen für komplexe Krankheiten. Verbindungen ohne vorhergesagte Ziele oder mit BATMAN-TCM-Zielwerten unter 20 wurden ausgeschlossen. Zusätzlich wurden Verbindungen mit relativ niedrigen OB-Werten, aber vorhergesagten Zielen in BATMAN-TCM und gut dokumentierten pharmakologischen Aktivitäten (z. B. Anthraquinone und Flavonoide) durch Literaturauswertung beibehalten. Um potenzielle falsch-positive Komponenten weiter zu reduzieren, wurden endogene Hormone und neurotransmitterähnliche Substanzen (z. B. Progesteron und Serotonin) manuell ausgeschlossen. Der endgültige Verbindungssatz wurde für die spätere Zielvorhersage und Netzwerkkonstruktion verwendet. Der detaillierte Arbeitsablauf ist in ergänzender Abbildung 1 dargestellt. Die endgültige Liste der aktiven Komponenten ist in Ergänzungstabelle 1 bereitgestellt.
Die SMILES-Identifikatoren (eine lineare Notationsmethode, die ASCII-Strings verwendet, um die chemische Struktur von Molekülen zu beschreiben. Durch die Verwendung spezifischer Symbole und Regeln kodiert es Atome, chemische Bindungen und molekulare Topologie in eine kontinuierliche Zeichenfolge ohne Leerzeichen, was eine prägnante, eindeutige und computerlesbare Speicherung und Übertragung molekularer Informationen ermöglicht; jede aktive Komponente wurde aus der PubChem-Datenbank abgerufen und potenzielle Ziele wurden mit der Schweizer Zielvorhersage vorhergesagt, wobei der Organismusparameter auf 'Homo sapiens' gesetzt und eine Wahrscheinlichkeitsgrenze von >0 angewendet wurde. Die aus beiden Datenbanken gewonnenen Zielsätze wurden zusammengeführt und doppelte Einträge entfernt, um den endgültigen Kandidatenzielsatz für DCQD zu erzeugen.
Es ist wichtig zu beachten, dass Mirabilitum (Mangxiao, Na₂SO₄·10H₂O), ein Mineralstoff im DCQD, nicht in den TCMSP- oder ähnlichen Verbindungs-Ziel-Datenbanken enthalten ist. Folglich umfasst aufgrund dieser Datenbankbeschränkung das in dieser Studie aufgebaute Komponenten-Ziel-Netzwerk nicht die potenziellen Beiträge von Mirabilitum. Das Netzwerk basiert daher auf den phytochemischen Bestandteilen der anderen drei pflanzlichen Bestandteile: Rheum palmatum, Magnolia officinalis und Citrus aurantium.
Anreicherungsanalyse von DCQD-Zielen
DCQD-Ziele wurden zur Anreicherungsanalyse in die Metascape-Datenbank importiert, wobei die Art auf "Homo sapiens" beschränkt war. Genontologie (GO) und Kyoto-Encyclopedia of Gene and Genomes (KEGG) Anreicherungsanalysen wurden durchgeführt, wobei erstere biologische Prozesse (BP), zelluläre Komponenten (CC) und molekulare Funktionen (MF) umfassten. Die Top 6 der GO-Anreicherungsergebnisse und die Top 20 der KEGG-Anreicherungsergebnisse wurden ausgewählt. Die KEGG- und GO-Balkendiagramme wurden mit dem Online-Tool der Bioinformatik generiert. Die Daten wurden über das DOSE-Paket in R importiert, um das Krankheitsontologie-(DO)-Balkendiagramm zu erstellen.
Aufbau des Netzwerks "Herbal Medicine-Component-Target"
Die relevanten Komponenten von DCQD und ihre entsprechenden Ziele wurden in Cytoscape 3.10.3 importiert, um ein Netzwerkdiagramm zu erstellen, das die Beziehungen zwischen "Kräutermedizin, Bestandteil-Ziel" darstellt. Das CytoNCA-Plugin wurde verwendet, um Degree-Werte zu berechnen, und die fünf aktivsten Top-Komponenten wurden basierend auf dem Ranking gescreent.
Abruf von ARDS-Zielen aus öffentlichen Datenbanken
Mit dem Schlüsselwort "akutes Atemnotsyndrom" wurden ARDS-bezogene Ziele aus den OMIM- und GeneCards-Datenbanken abgerufen. Bei Genkarten wurden Ziele mit einem Relevanzwert >5 ausgewählt. Nach dem Zusammenführen und Entfernen von Duplikaten aus beiden Datenbanken wurde ein Kandidatenziel für ARDS ermittelt.
Bioinformatik-Analyse
Der sepsis-ARDS-bezogene Datensatz GSE32707 (GPL570-Plattform, n = 45), der RNA-seq-Daten von peripheren mononuklearen Blutzellen von 30 sepsisinduzierten ARDS-Patienten und 15 gesunden Kontrollpersonen umfasst, wurde aus der NCBI GEO-Datenbank heruntergeladen. Die Daten wurden mit dem GEOquery-Paket in R importiert, und die Normalisierung erfolgte mit dem Limma-Paket, um die Inter-Sample-Variabilität zu minimieren. Differentiell exprimierte Gene (DEGs) wurden anhand der Schwellenwerte |log2FC| identifiziert. > 1 und p < 0,05. Die Visualisierung erfolgte über Vulkandiagramme und Heatmaps. Für die Signalweganalyse wurden GO- und KEGG-Anreicherungsanalysen durchgeführt, wobei die Ergebnisse als Blasendiagramme und Barendiagramme präsentiert wurden.
Vorhersage potenzieller Ziele für DCQD in der ARDS-Behandlung
Ein Venn-Diagramm wurde mit dem bioinformatischen Online-Tool erstellt, um die DEGs, ARDS-bezogenen Ziele aus öffentlichen Datenbanken und DCQD-Arzneimittelziele zu vergleichen und so potenzielle therapeutische Ziele von DCQD für ARDS zu identifizieren.
Aufbau des Protein-Protein-Interaktionsnetzwerks (PPI) und Screening der Kernziele
Die sich überlappenden Ziele wurden in die STRING-Datenbank importiert, wobei die Spezies auf "Homo sapiens" gesetzt wurde und eine Interaktionswertschwelle von 0,4 erreicht wurde. Die PPI-Netzwerkdaten wurden im TSV-Format exportiert und mit Cytoscape 3.10.3 visualisiert. Das CytoHubba-Plugin wurde angewendet, um die Top 10 Ziele basierend auf MCC, MNC, Grad, EPC, Zwischenheit, Nähe, Radialität und Stresswerten zu screenen. Die Schnittstelle dieser Ziele ergab fünf Kernziele.
Aufbau des Netzwerks "Kräutermedizin-Komponent-Ziel-Krankheit"
Basierend auf den sich überlappenden Zielen wurden entsprechende DCQD-Komponenten identifiziert. Die relevanten Komponenten sowie ihre zugehörigen Ziele und Krankheitsverbindungen wurden in Cytoscape 3.10.3 importiert, um ein Netzwerkdiagramm zu erstellen, das die Beziehungen zwischen "Kräutermedizin, Komponente-Ziel-Krankheit" darstellt.
Anreicherungsanalyse der Kernziele und Aufbau des Netzwerks "Kräutermedizin-Komponente-Ziel-Krankheitsweg"
Um die potenziellen biologischen Funktionen und Schlüsselsignalwege von DCQD in der ARDS-Behandlung zu erforschen, wurden die Kernziele in Metascape (Spezies: "Homo sapiens") mit GO- und KEGG-Anreicherungsanalysen durchgeführt. Die Top 10 von GO und die Top 13 der KEGG-Anreicherungsergebnisse wurden ausgewählt. Die Visualisierung wurde mit Bioinformatik-Tools durchgeführt, um KEGG-Blasendiagramme, GO-Balkendiagramme und Sankey-Blasendiagramme zu erzeugen. Zusätzlich wurde auf Basis der Ergebnisse der funktionellen Anreicherung in Cytoscape 3.10.3 ein Netzwerkdiagramm erstellt, um die Beziehungen zwischen "Kräutermedizin-Komponente-Ziel-Krankheit-Weg" zu veranschaulichen.
Molekulare Andocken und Visualisierung
Die sekundären molekularen Strukturen der Verbindungen wurden aus der PubChem-Datenbank abgerufen, während die 3D-Strukturen der Kernziele aus der PDB-Datenbank (http://www.rcsb.org/) bezogen wurden. Wassermoleküle und Aminosäurereste wurden mit PyMOL entfernt. Zielproteine und -verbindungen wurden in AutoDock importiert, um das molekulare Andocken zur Berechnung der Bindungsenergien zu berechnen. Die Visualisierung erfolgte mit PyMOL. Mit dem Pheatmap-Paket in R wurde eine Bindungsenergie-Heatmap erzeugt, mit Kernzielen auf der x-Achse, Verbindungen auf der y-Achse und einem Farbgradient, der die Bindungsenergien darstellt. Zweidimensionale Interaktionsdiagramme der Andockergebnisse wurden mit Discovery Studio Visualizer 4.5 erstellt. Das molekulare Andocken wurde durchgeführt, um die Bindungsmachbarkeit zwischen Verbindungen und Zielgruppen zu bewerten, anstatt die tatsächliche inhibitorische Wirksamkeit vorherzusagen.
Konstruktion des ARDS-Mausmodells
Diese Studie wurde vom medizinischen Ethikausschuss des Suzhou Dushu Lake Hospital (2410008) genehmigt. Gesunde männliche C57BL/6J-Mäuse (n = 30; 20 ± 5 g) wurden gefunden. Die Mäuse wurden in eine Kontrollgruppe (n = 10), eine LPS-Gruppe (n = 10, LPS 10 mg/kg) und eine DCQD-Gruppe (n = 10, die DCQD 0,9 g/kg durch orale Gavage in Kombination mit LPS 10 mg/kg erhielt) unterteilt. Mäuse in der LPS-Gruppe erhielten eine intraperitoneale Injektion von LPS. Mäuse in der DCQD-Gruppe erhielten DCQD 0,9 g/kg unmittelbar nach der intraperitonealen Injektion von LPS per oraler Gavage. Unmittelbar vor und nach der Lungenextraktion wurden bei allen Mäusen Messungen durchgeführt. Anschließend wurden die Lungengewebe bei -80 °C gelagert.
Western-Blot-Analyse (WB)
Proteine wurden aus Lysaten des Mauslungengewebes isoliert, und ihre Konzentrationen wurden mit einem BCA-Proteinquantifizierungskit gemessen. Anschließend wurden die Proteine durch 10 Minuten Kochen bei 95 °C denaturiert. Die denaturierten Proteine wurden auf einem 10%igen SDS-PAGE-Gel getrennt, das mit einem Gel-Vorbereitungsset vorbereitet wurde, und dann auf eine PVDF-Membran übertragen. Um eine unspezifische Bindung zu verhindern, wurden die PVDF-Membranen bei Raumtemperatur für 5 Minuten mit einem proteinfreien, schnellen Blockpuffer blockiert. Nach dem Waschen mit PBST wurden die Membranen über Nacht bei 4 °C mit primären Antikörpern gegen MMP9, HSP90, MYC und ARG1 inkubiert. Nach einer weiteren Reinigung mit PBST wurden die Membranen 2 Stunden lang bei Zimmertemperatur mit HRP-konjugiertem Ziegen-Antikaninchen-IgG(H+L) inkubiert. Abschließend wurde das ultrasensitive Chemilumineszenz-Detektions-Kit-Substrat angewendet, und chemilumineszente Proteinbänder wurden mit einem BIO-RAD-Chemilumineszenz-Bildgebungssystem detektiert. Die Graustufenanalyse wurde mit ImageJ durchgeführt. Die in diesem Experiment verwendeten Antikörper sind in ergänzender Tabelle 2 aufgeführt.
Enzymgekoppelte Immunsorbent-Assay-Analyse (ELISA)
Zuerst wurde die bronchoalveoläre Lavageflüssigkeit (BALF) durch dreifache Lavage mit 2 ml Kochsalzlösung gesammelt. Die resultierenden Lavageflüssigkeiten wurden gesammelt und bei 10.000 g für 10 Minuten zentrifugiert, um Zellen zu entfernen. Der Supernatant wurde dann aspiriert und gefroren in flüssigem Stickstoff gelagert. Anschließend wurden die Konzentrationen der Zielzytokine gemäß den Anweisungen des Herstellers mit dem Mouse IL-6 ELISA Kit, Mouse IL-18 ELISA Kit und Mouse TNF-alpha ELISA Kit bestimmt. Anschließend wurde die Reaktionsplatte in einen Mikroplattenleser gesetzt, und der Wert der optischen Dichte (OD) jedes Bohrlochs wurde bei einer Wellenlänge von 450 nm gemessen, wobei der Blank-Kontrollbrunnen für Nullanpassung verwendet wurde. Abschließend wurde eine Standardkurve auf Basis der Konzentrationen und entsprechenden OD-Werte der Standards dargestellt. Die Konzentration des Zielfaktors in jeder Probe wurde gemäß der Standardkurvengleichung berechnet.
Statistische Analyse
Die Daten wurden mit GraphPad Prism 10.1.2 analysiert. Die Signifik der Unterschiede zwischen den Gruppen wurde statistisch durch eine Einweganalyse der Varianz oder den Schüler-t-Test bewertet. Ist der P-Wert kleiner als 0,05, gilt die Differenz als signifikant; ist der P-Wert kleiner als 0,01, gilt er als sehr bedeutend.