Voorspelling en screening van potentiële doelen van PLB
De informatie over de chemische structuur van PLB werd verkregen uit de PubChem-database (https://pubchem.ncbi.nlm.nih.gov/). Moleculaire targets van PLB werden voorspeld op basis van ligand-gebaseerde tweedimensionale en driedimensionale gelijkenis met behulp van de databanken 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) en HERB (http://herb.ac.cn/). Om de betrouwbaarheid en reproduceerbaarheid van de gegevens te waarborgen, werden database-specifieke filterdrempels als volgt toegepast: TCMSP behield targets met een orale biobeschikbaarheid (OB) ≥ 30% en drug-likeness (DL) ≥ 0,18; SwissTargetPrediction behield vermeldingen met een voorspellingswaarschijnlijkheid ≥ 0,5; PharmMapper selecteerde targets met een genormaliseerde fit-score ≥ 4,0; SEA includerde alleen targets met een E-waarde < 0,001 als significant verrijkt; TargetNet behield targets met een voorspellingswaarschijnlijkheid > 0,5; CTD includerde alleen targets met gecureerde bewijslast (afgeleid van chemische-geninteracties) en een interactiescore > 0,3; en HERB behield targets met een op literatuur gebaseerde score ≥ 0,4. Alle opgehaalde eiwitidentificatoren uit de bovengenoemde databanken werden gestandaardiseerd naar officiële humane HGNC-gensymbolen met behulp van de UniProt-database (https://www.uniprot.org/), waarbij de soort beperkt bleef tot Homo sapiens. Na het verwijderen van dubbele vermeldingen uit alle bronnen werden in totaal 500 unieke PLB-gerelateerde targets verkregen voor verdere analyse.
Ophalen van prostaatkanker-geassocieerde genen
Ziektegerelateerde genen voor prostaatkanker werden verkregen uit GeneCards (https://www.GeneCardss.org/), DrugBank (https://go.drugbank.com/), CTD en HERB. Er werden databasespecifieke inclusiecriteria toegepast om de betrouwbaarheid van de gegevens te waarborgen. GeneCardss behield genen met een relevantiescore ≥ 0,5, aangezien deze drempelwaarde genen omvat met matig tot sterk bewijs dat hen verbindt met de ziektevraag. DrugBank bevatte alleen vermeldingen met experimenteel bewijs (bijv. FDA-goedgekeurde of onderzoekende geneesmiddelen voor prostaatkanker) en sloot computationeel voorspelde of theoretische interacties uit. CTD behield alleen records met bewijsniveaus geclassificeerd als "marker" of "mechanisme" op basis van gecureerde chemische-gen-ziekte-interacties. HERB bevatte targets met een in de literatuur ondersteunde betrouwbaarheidsscore ≥ 0,4 om voldoende experimenteel bewijs of bewijs uit tekstmining te garanderen. Alle gensymbolen werden via de UniProt-database (https://www.uniprot.org/) gestandaardiseerd naar de HGNC- nomenclatuur, waarbij de soorten beperkt bleven tot Homo sapiens, en duplicaten werden verwijderd, wat resulteerde in 1.199 unieke targets gerelateerd aan prostaatkanker voor verdere analyse.
Identificatie van overlappende doelwitten en correctie van conflictannotaties
De 500 PLB-targets en 1.199 prostaatkanker-targets werden geïntersecteerd met behulp van gestandaardiseerde HGNC-symbolen in Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html), wat resulteerde in 151 overlappende kandidaat-targets. Om annotatieconflicten tussen databases op te lossen, werd een correctiepijplijn in drie stappen geïmplementeerd: (i) inconsistenties in gen-aliassen werden geünificeerd via de UniProt ID mapping tool; (ii) redundante paraloge genen werden uitgesloten met CD-HIT met een drempelwaarde voor sequentie-overeenkomst > 0.4; en (iii) tegenstrijdige functionele annotaties werden alleen behouden als ze werden ondersteund door ten minste twee onafhankelijke databases, terwijl unieke conflicterende beschrijvingen uit een enkele bron werden verworpen.
Netwerkconstructie en identificatie van hubs
Een eiwit-eiwitinteractienetwerk (PPI-netwerk) werd geconstrueerd door de overlappende targets in te voeren in de STRING-database (versie 12.0, https://cn.string-db.org/), waarbij de zoekopdracht werd beperkt tot Homo sapiens. Alleen interacties met een gecombineerde betrouwbaarheidsscore ≥ 0,700 werden behouden om een hoge betrouwbaarheid te waarborgen, en alle knopen zonder verbindingen werden uit het netwerk uitgesloten. Het geconstrueerde PPI-netwerk werd gevisualiseerd en topologisch geanalyseerd in Cytoscape (versie 3.10.2), waarna hub-genen werden bepaald met de cytoHubba-plugin, waarbij MCC (Maximal Clique Centrality) als primair rangschikkingsalgoritme en Degree als aanvullende crossvalidatiemetriek werden gebruikt.
Functionele verrijkingsanalyse
Functionele annotatie van de Gene Ontology (GO) en verrijkingsanalyses van Kyoto Encyclopedia of Genes and Genomes (KEGG)-paden voor de differentieel tot expressie gebrachte genen werden uitgevoerd met het Bioinformatics Online Analysis Platform (https://www.bioinformatics.com.cn, laatst geraadpleegd op 4 mei 2026). Het platform maakt gebruik van de exacte toets van Fisher om ruwe P-waarden te berekenen, en de Benjamini-Hochberg (BH)-methode werd gebruikt voor correctie op meervoudig testen om de false discovery rate (FDR) te beheersen. GO-termen (biologisch proces, cellulaire component en moleculaire functie) en KEGG-paden met een FDR < 0,05 werden beschouwd als significant verrijkt.
Moleculaire docking
Kristalstructuren van de geselecteerde doeleiwitten, die actieve conformaties vertegenwoordigen, werden verkregen uit de RCSB Protein Data Bank (https://www.rcsb.org/). De voorbereiding van de eiwitten werd uitgevoerd in Discovery Studio, inclusief het verwijderen van watermoleculen en heteroatomen, het verwijderen van redundante ketens, het herstellen van ontbrekende residuen, het toevoegen van waterstofatomen bij pH 7.4, het toewijzen van Gasteiger-ladingen en energieminimalisatie. De bindingspocket werd gedefinieerd op basis van de co-gekristalliseerde ligandcoördinaten, met gebruik van een grid dat alle kritische substraatbindende residuen omvatte. Moleculaire docking van PLB werd uitgevoerd met de SwissDock webserver (http://www.swissdock.ch/) via de Attracting Cavities-module. De vrije bindingsenergieën werden berekend om de affiniteit te beoordelen, en de conformatie met de laagste energie voor elk doeleiwit werd geselecteerd als de definitieve dockingpose. De resultaten van de moleculaire dockingsimulaties werden gevisualiseerd om deze te valideren.
Moleculaire dynamicasimulatie
Moleculaire dynamica (MD)-simulaties werden uitgevoerd met GROMACS 2022.2. Het Amber14SB-krachtveld werd gebruikt om het eiwit te beschrijven, waarbij het systeem werd gesolvateerd in TIP3P-water. Plumbagin-parameters, inclusief AM1-BCC partiële ladingen en GAFF2-atoomtypen, werden gegenereerd met Antechamber, gevolgd door topologieconversie met ACPYPE en toewijzing van Joung-Cheatham ionparameters. Elk eiwit-ligandcomplex werd ingebed in een afgeplatte dodecahedrische simulatiebox met een minimale afstand van eiwit tot box van 1,2 nm, gesolvateerd met TIP3P-watermoleculen en geneutraliseerd door toevoeging van 0,15 M Na⁺/Cl⁻. Na energieminimalisatie met het steepest descent-algoritme (Fmax < 1.000 kJ·mol⁻1·nm⁻1) onderging het systeem sequentiële NVT- en NPT-equilibratie van elk 200 ps bij 298 K. Productiesimulaties werden vervolgens uitgevoerd gedurende 200 ns onder NPT-condities met een integratietijdstap van 2 fs met behulp van het Verlet-cutoff-schema. Langeafstandselektrostatische interacties werden berekend met de particle mesh Ewald (PME)-methode, terwijl voor zowel elektrostatische als van der Waals-interacties een cutoff-afstand van 1,2 nm werd gehanteerd. Het LINCS-algoritme werd gebruikt om waterstofbevattende bindingen te beperken. De temperatuur werd op 298 K gehouden met de Nosé-Hoover-thermostaat en de druk op 1 bar met de Parrinello-Rahman-barostaat. Coördinaten werden elke 10 ps opgeslagen voor latere analyses. Trajectanalyse en visualisatie werden uitgevoerd met GROMACS-hulpprogramma's, VMD en PyMOL, terwijl MM-PBSA bindingvrije energieberekeningen, waar van toepassing, werden uitgevoerd met gmx_MMPBSA.
Equilibratiebeoordeling van moleculaire dynamiektrajecten en MM-PBSA vrije-energiebemonsteringsprotocol
Voor elk eiwit–plumbagine-complex werd een all-atom productie-moleculaire dynamics-simulatie van 200 ns uitgevoerd met GROMACS. Het tijdspunt van 100 ns werd ingesteld als de equilibratiegrens: het segment van 0–100 ns werd aangemerkt als de conformationele relaxatiefase, gedurende welke de eiwitbackbone en de ligandbindende pocket continue conformationele aanpassingen ondergingen; het segment van 100–200 ns werd geïdentificeerd als het thermodynamisch stabiele plateau, evidenced door de afwezigheid van unidirectionele drift in RMSD, gyratie-straal (Rg), residu-specifieke RMSF, begraven SASA van de ligand, intermoleculaire waterstofbruggen en ligand–receptor interactie-energieën, welke enkel geringe steady-state fluctuaties vertoonden. Alle kwantitatieve kinetische parameters en MM/PBSA-bindingsvrije energieën werden exclusief berekend uit de equilibriumfase van 100–200 ns, waaruit frames uniform elke 100 ps werden geëxtraheerd, wat resulteerde in 1.000 equilibrium-snapshots per systeem voor gmx_MMPBSA-input. De eerste 100 ns aan relaxatietrajecten werden verworpen om interferentie door conformationele drift in de vrije-energieberekeningen te elimineren.
Prognostische analyse van targetgenen bij prostaatkanker
Prognostische analyses werden uitgevoerd met prostaatkankerdatasets van The Cancer Genome Atlas (TCGA). RNA-sequencing-gegevens (STAR-counts) en de bijbehorende klinische informatie werden verkregen via het Genomic Data Commons-portaal (https://portal.gdc.cancer.gov). Genexpressietellingen werden omgezet naar transcripts per million (TPM) en genormaliseerd met behulp van log2(TPM + 1). Na uitsluiting van monsters met onvolledige klinische informatie werden 498 gevallen in de analyse opgenomen. Het mediane expressieniveau van elk gen werd gebruikt om patiënten te verdelen in groepen met een hoge en lage expressie. Kaplan-Meier-overlevingsanalyse met de log-ranktest en univariate Cox-proportional hazards-regressie werden uitgevoerd om de algehele overleving (OS) voor SRC en de progressievrije overleving (PFS) voor MAPK3 te evalueren, waarbij hazard ratio's (HR's) en 95% betrouwbaarheidsintervallen (CI's) werden gerapporteerd. De voorspellende prestaties werden verder beoordeeld door tijdafhankelijke receiver operating characteristic (ROC)-curves te genereren voor 1, 3 en 5 jaar. Statistische analyses werden uitgevoerd in R versie 4.0.3, waarbij P < 0,05 als statistisch significant werd beschouwd.