Vorhersage und Screening potenzieller Zielstrukturen von PLB
Die chemische Strukturinformation von PLB wurde aus der PubChem-Datenbank (https://pubchem.ncbi.nlm.nih.gov/) bezogen. Molekulare Zielstrukturen von PLB wurden basierend auf ligandenbasierten zweidimensionalen und dreidimensionalen Ähnlichkeiten mittels der Datenbanken SwissTarget (https://swisstargetprediction.ch/index.php), SEA (https://sea.bkslab.org/), TargetNet (http://targetnet.scbdd.com/calcnet/index/), PharmMapper (https://www.lilab-ecust.cn/pharmmapper/index.html), Comparative Toxicogenomics Database (CTD, https://ctdbase.org/), Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, https://www.tcmsp-e.com/index.php) und HERB (http://herb.ac.cn/) vorhergesagt. Um Datenzuverlässigkeit und Reproduzierbarkeit sicherzustellen, wurden datenbankspezifische Filtergrenzwerte wie folgt angewendet: TCMSP behielt Zielstrukturen mit einer oralen Bioverfügbarkeit (OB) ≥ 30 % und Drug-Likeness (DL) ≥ 0,18 bei; SwissTargetPrediction behielt Einträge mit einer Vorhersagewahrscheinlichkeit ≥ 0,5 bei; PharmMapper wählte Zielstrukturen mit einem normalisierten Fit-Score ≥ 4,0 aus; SEA schloss nur Zielstrukturen mit einem E-Wert < 0,001 als signifikant angereichert ein; TargetNet behielt Zielstrukturen mit einer Vorhersagewahrscheinlichkeit > 0,5 bei; CTD schloss nur Zielstrukturen mit kuratierten Belegen (abgeleitet aus Chemikalie–Gen-Interaktionen) und einem Interaktionsscore > 0,3 ein; und HERB behielt Zielstrukturen mit einem literaturgestützten Score ≥ 0,4 bei. Alle abgerufenen Proteinidentifikatoren aus den oben genannten Datenbanken wurden mithilfe der UniProt-Datenbank (https://www.uniprot.org/) auf offizielle humane HGNC-Gen-Symbole standardisiert, wobei die Spezies auf Homo sapiens beschränkt wurde. Nach Entfernung doppelter Einträge aus allen Quellen wurden insgesamt 500 eindeutige PLB-bezogene Zielstrukturen für die nachfolgende Analyse erhalten.
Abruf von mit Prostatakrebs assoziierten Genen
Krankheitsbezogene Gene für Prostatakarzinom wurden aus GeneCards (https://www.GeneCardss.org/), DrugBank (https://go.drugbank.com/), CTD und HERB bezogen. Es wurden datenbankspezifische Einschlusskriterien angewandt, um die Datenzuverlässigkeit sicherzustellen. GeneCards behielt Gene mit einem Relevanzscore ≥ 0,5 bei, da diese Schwelle Gene mit mittlerer bis starker Evidenz, die sie mit der Krankheitsanfrage verbinden, erfasst. DrugBank schloss nur Einträge mit experimenteller Evidenz ein (z. B. FDA-zugelassene oder in der Prüfung befindliche Wirkstoffe gegen Prostatakarzinom) und nahm keine computergestützt vorhergesagten oder theoretischen Wechselwirkungen auf. CTD behielt ausschließlich Datensätze bei, deren Evidenzniveau als „Marker“ oder „Mechanismus“ klassifiziert war, basierend auf kurierten chemisch-gen-krankheitsbezogenen Interaktionen. HERB schloss Zielstrukturen mit einem literaturbasierten Konfidenzscore ≥ 0,4 ein, um ausreichende experimentelle oder textbasierte Evidenz sicherzustellen. Alle Gensymbole wurden über die UniProt-Datenbank (https://www.uniprot.org/) der HGNC-Nomenklatur angepasst, wobei die Art auf Homo sapiens beschränkt wurde; Duplikate wurden entfernt, wodurch 1.199 eindeutige, mit Prostatakarzinom assoziierte Zielstrukturen für die nachfolgende Analyse verblieben.
Identifizierung überlappender Zielstrukturen und Korrektur von Konfliktanmerkungen
Die 500 PLB-Ziele und 1.199 Prostatakrebs-Ziele wurden mittels standardisierter HGNC-Symbole in Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html) übereinandergelegt, wodurch sich 151 überlappende Kandidatenziele ergaben. Zur Auflösung von Querdatenbank-Anmerkungs-Konflikten wurde eine dreistufige Korrektur-Pipeline implementiert: (i) Gen-Alias-Inkonsistenzen wurden über das UniProt ID-Mapping-Tool vereinheitlicht; (ii) paraloge redundante Gene wurden mithilfe von CD-HIT mit einer Sequenzähnlichkeits-Schwelle > 0,4 ausgeschlossen; und (iii) widersprüchliche funktionelle Annotationen wurden nur beibehalten, wenn sie durch mindestens zwei unabhängige Datenbanken gestützt wurden, während eindeutige, konfliktbehaftete Beschreibungen aus einer einzigen Quelle verworfen wurden.
Netzwerkkonstruktion und Identifizierung von zentralen Knotenpunkten
Ein Protein-Protein-Interaktionsnetzwerk (PPI) wurde erstellt, indem die überlappenden Zielstrukturen in die STRING-Datenbank (Version 12.0, https://cn.string-db.org/) eingegeben wurden, wobei die Suche auf Homo sapiens beschränkt blieb. Es wurden nur Interaktionen mit einem kombinierten Konfidenzscore ≥ 0,700 berücksichtigt, um eine hohe Zuverlässigkeit sicherzustellen; alle Knoten ohne Verbindungen wurden aus dem Netzwerk ausgeschlossen. Das erstellte PPI-Netzwerk wurde in Cytoscape (Version 3.10.2) visualisiert und hinsichtlich seiner topologischen Eigenschaften analysiert. Anschließend wurden die zentralen Gene (Hub-Gene) mithilfe des cytoHubba-Plugins ermittelt, wobei MCC (Maximale Cliquen-Zentralität) als primärer Rangalgorithmus und Degree als zusätzliche Validierungsmetrik verwendet wurden.
Funktionsbereicherungsanalyse
Die funktionelle Annotation nach Gene Ontology (GO) und die Analyse der Kyoto Encyclopedia of Genes and Genomes (KEGG)-Pfadanreicherung für die differentiell exprimierten Gene wurden mit der Bioinformatics Online Analysis Platform (https://www.bioinformatics.com.cn, letzter Zugriff am 4. Mai 2026) durchgeführt. Die Plattform verwendet den exakten Fisher-Test zur Berechnung der unkorrigierten P-Werte, und die Korrektur für Mehrfachvergleiche erfolgte mittels der Benjamini-Hochberg-Methode (BH), um die Falsch-Entdeckungs-Rate (FDR) zu kontrollieren. GO-Terme (biologischer Prozess, zelluläre Komponente und molekulare Funktion) sowie KEGG-Pfade mit einem FDR < 0,05 galten als signifikant angereichert.
Molekulares Docking
Kristallstrukturen der ausgewählten Zielproteine, die aktive Konformationen repräsentieren, wurden aus der RCSB Protein Data Bank (https://www.rcsb.org/) bezogen. Die Proteinaufbereitung erfolgte in Discovery Studio, einschließlich der Entfernung von Wassermolekülen und Heteroatomen, des Löschens redundanter Ketten, der Reparatur fehlender Reste, der Addition von Wasserstoffatomen bei pH 7,4, der Zuweisung von Gasteiger-Ladungen sowie der Energie-Minimierung. Die Bindungstasche wurde basierend auf den Koordinaten des ko-kristallisierten Liganden definiert, unter Verwendung eines Rasters, das alle wichtigen Substrat-bindenden Reste umfasste. Das molekulare Docking von PLB wurde mittels des SwissDock-Web-Servers (http://www.swissdock.ch/) über das Modul „Attracting Cavities“ durchgeführt. Die freien Bindungsenergien wurden zur Bewertung der Affinität berechnet, und die Konformation mit der niedrigsten Energie für jedes Target wurde als endgültige Docking-Position ausgewählt. Die Ergebnisse der molekularen Docking-Simulationen wurden zur Validierung visualisiert.
Simulation der Moleküldynamik
Molekulardynamik-(MD-)Simulationen wurden mit GROMACS 2022.2 durchgeführt. Für die Beschreibung des Proteins kam das Amber14SB-Kraftfeld zum Einsatz, wobei das System in TIP3P-Wasser solvatisiert wurde. Plumbagin-Parameter, einschließlich AM1-BCC-Teilladungen und GAFF2-Atomtypen, wurden mittels Antechamber generiert, gefolgt von einer Topologiekonvertierung mit ACPYPE und der Zuweisung der Joung-Cheatham-Ionenparameter. Jeder Protein-Ligand-Komplex wurde in eine abgeschnittene dodekaedrische Simulationsbox eingebettet, wobei ein minimaler Abstand zwischen Protein und Boxwand von 1,2 nm eingehalten wurde, anschließend mit TIP3P-Wassermolekülen solvatisiert und durch Zugabe von 0,15 M Na⁺/Cl⁻ neutralisiert. Nach einer Energieminimierung mittels des steilsten Abstiegs-Algorithmus (Fmax < 1.000 kJ·mol⁻1·nm⁻1) erfolgte eine sequenzielle NVT- und NPT-Equilibration jeweils über 200 ps bei 298 K. Anschließend wurden Produktionsläufe über 200 ns unter NPT-Bedingungen mit einem Integrationszeitschritt von 2 fs unter Verwendung des Verlet-Cutoff-Schemas durchgeführt. Langreichweitige elektrostatische Wechselwirkungen wurden mit der Particle-Mesh-Ewald-(PME-)Methode berechnet, während sowohl elektrostatische als auch van-der-Waals-Wechselwirkungen einen Cutoff-Abstand von 1,2 nm verwendeten. Der LINCS-Algorithmus wurde zur Einschränkung von bindenden Atomen mit Wasserstoff eingesetzt. Die Temperatur wurde mittels Nosé-Hoover-Thermostat bei 298 K gehalten, der Druck mittels Parrinello-Rahman-Barostat bei 1 bar. Koordinaten wurden alle 10 ps zur späteren Analyse gespeichert. Die Trajektorienanalyse und Visualisierung erfolgte mit GROMACS-Programmen, VMD und PyMOL, während MM-PBSA-Bindungsenergieberechnungen gegebenenfalls mit gmx_MMPBSA durchgeführt wurden.
Auswertung der Gleichgewichtseinstellung von Molekulardynamik-Trajektorien und MM-PBSA-Freie-Energie-Stichprobenprotokoll
Für jeden Protein–Plumbagin-Komplex wurde eine 200 ns lange all-atom Produktionssimulation mittels molekularer Dynamik mit GROMACS durchgeführt. Der Zeitpunkt von 100 ns wurde als Gleichgewichtsgrenze festgelegt: Der Abschnitt von 0–100 ns wurde als konformationelle Relaxationsphase definiert, in der das Proteinrückgrat und die Ligandenbindungsstelle kontinuierliche konformationelle Anpassungen durchliefen; der Abschnitt von 100–200 ns wurde als thermodynamisch stabiles Plateau identifiziert, wie durch das Fehlen eines unidirektionalen Drifts in RMSD, Trägheitsradius (Rg), restenspezifischem RMSF, verborgener SASA des Liganden, intermolekularen Wasserstoffbrückenbindungen sowie Ligand–Rezeptor-Wechselwirkungsenergien belegt wurde, die lediglich geringfügige Schwankungen im stationären Zustand zeigten. Alle quantitativen kinetischen Parameter und MM/PBSA-Bindungsenergien wurden ausschließlich aus der Gleichgewichtsphase von 100–200 ns berechnet, aus der alle 100 ps gleichmäßig Frames extrahiert wurden, was 1.000 Gleichgewichtsaufnahmen pro System für die Eingabe in gmx_MMPBSA ergab. Die ersten 100 ns der Relaxationstrajektorien wurden verworfen, um Störungen durch konformationellen Drift bei den Freienergieberechnungen auszuschließen.
Prognostische Analyse der Zielgene beim Prostatakarzinom
Prognoseanalysen wurden unter Verwendung von Prostatakrebs-Datensätzen aus dem Cancer Genome Atlas (TCGA) durchgeführt. RNA-Sequenzierungsdaten (STAR-Counts) und die entsprechenden klinischen Informationen wurden vom Genomic Data Commons-Portal (https://portal.gdc.cancer.gov) bezogen. Genexpressionszahlen wurden in Transkripte pro Million (TPM) umgerechnet und mittels log2(TPM + 1) normalisiert. Nach Ausschluss von Proben mit unvollständigen klinischen Angaben wurden 498 Fälle in die Analyse eingeschlossen. Der mediane Expressionswert jedes Gens wurde verwendet, um die Patienten in Gruppen mit hoher und niedriger Expression einzuteilen. Die Kaplan-Meier-Überlebensanalyse mit dem Log-Rang-Test und die univariate Cox-Regression nach dem proportionalen Hazard-Modell wurden durchgeführt, um das Gesamtüberleben (OS) für SRC und das progressionsfreie Überleben (PFS) für MAPK3 zu bewerten, wobei Hazard-Ratios (HRs) und 95-%-Konfidenzintervalle (CIs) angegeben wurden. Die prädiktive Leistung wurde weiterhin durch Erstellung zeitabhängiger ROC-Kurven (Receiver Operating Characteristic) nach 1, 3 und 5 Jahren überprüft. Die statistischen Analysen wurden in R-Version 4.0.3 durchgeführt, wobei ein P < 0,05 als statistisch signifikant angesehen wurde.