Previsione e screening di potenziali bersagli della PLB
Le informazioni sulla struttura chimica del PLB sono state ottenute dal database PubChem (https://pubchem.ncbi.nlm.nih.gov/). I bersagli molecolari del PLB sono stati predetti in base alla similarità bidimensionale e tridimensionale di tipo ligando-dipendente, utilizzando i database 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) e HERB (http://herb.ac.cn/). Per garantire affidabilità e riproducibilità dei dati, sono state applicate soglie di filtraggio specifiche per ciascun database: TCMSP ha mantenuto i bersagli con una biodisponibilità orale (OB) ≥ 30% e una somiglianza ai farmaci (DL) ≥ 0,18; SwissTargetPrediction ha conservato le voci con una probabilità di predizione ≥ 0,5; PharmMapper ha selezionato i bersagli con un punteggio di adattamento normalizzato ≥ 4,0; SEA ha incluso soltanto i bersagli con un valore E < 0,001 come significativamente arricchiti; TargetNet ha mantenuto i bersagli con una probabilità di predizione > 0,5; CTD ha incluso solo i bersagli con evidenze curate (dedotte da interazioni chimico-geniche) e un punteggio di interazione > 0,3; HERB ha conservato i bersagli con un punteggio supportato dalla letteratura ≥ 0,4. Tutti gli identificatori proteici recuperati dai suddetti database sono stati standardizzati in simboli genici ufficiali umani HGNC utilizzando il database UniProt (https://www.uniprot.org/), limitando la specie a Homo sapiens. Dopo l'eliminazione delle voci duplicate provenienti da tutte le fonti, sono stati ottenuti complessivamente 500 bersagli univoci associati al PLB per le analisi successive.
Recupero di geni associati al cancro della prostata
I geni associati alla malattia per il cancro alla prostata sono stati ottenuti da GeneCards (https://www.GeneCardss.org/), DrugBank (https://go.drugbank.com/), CTD e HERB. Sono stati applicati criteri di inclusione specifici per ciascun database per garantire l'affidabilità dei dati. GeneCardss ha mantenuto i geni con un punteggio di rilevanza ≥ 0,5, poiché questa soglia include geni con prove da moderate a forti che li collegano alla malattia in esame. DrugBank ha incluso solo voci con evidenze sperimentali (ad esempio farmaci approvati dalla FDA o in fase di sperimentazione per il cancro alla prostata), escludendo interazioni previste computazionalmente o teoriche. CTD ha mantenuto solo i record con livelli di evidenza classificati come «marcatore» o «meccanismo», basati su interazioni chimico-gene-malattia curate. HERB ha incluso i bersagli con un punteggio di confidenza supportato dalla letteratura ≥ 0,4 per assicurare evidenze sperimentali o derivate da text-mining sufficienti. Tutti i simboli genici sono stati uniformati alla nomenclatura HGNC tramite il database UniProt (https://www.uniprot.org/), limitatamente alla specie Homo sapiens, e sono stati rimossi i duplicati, ottenendo 1.199 bersagli unici associati al cancro alla prostata per le analisi successive.
Identificazione di bersagli sovrapposti e correzione dell'annotazione dei conflitti
I 500 bersagli PLB e i 1.199 bersagli del cancro alla prostata sono stati incrociati utilizzando simboli HGNC standard in Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html), ottenendo 151 bersagli candidati sovrapposti. Per risolvere i conflitti di annotazione tra database, è stata implementata una procedura di correzione in tre fasi: (i) le incongruenze relative agli alias genici sono state unificate mediante lo strumento di mappatura degli ID UniProt; (ii) i geni ridondanti paraloghi sono stati esclusi utilizzando CD-HIT con una soglia di similarità di sequenza > 0,4; e (iii) le annotazioni funzionali contrastanti sono state mantenute solo se supportate da almeno due database indipendenti, mentre le descrizioni contrastanti univoche provenienti da un'unica fonte sono state scartate.
Costruzione della rete e identificazione degli hub
Una rete di interazioni proteina-proteina (PPI) è stata costruita inviando le proteine sovrapposte al database STRING (versione 12.0, https://cn.string-db.org/), limitando la ricerca a Homo sapiens. Sono state mantenute solo le interazioni con un punteggio di confidenza combinato ≥ 0,700 per garantire un'elevata affidabilità, ed esclusi dalla rete tutti i nodi privi di connessioni. La rete PPI costruita è stata visualizzata e analizzata dal punto di vista topologico in Cytoscape (versione 3.10.2); successivamente, i geni hub sono stati identificati utilizzando il plugin cytoHubba, adottando MCC (Centralità della Clique Massimale) come algoritmo principale di classificazione e Degree come metrica ausiliaria di validazione incrociata.
Analisi di arricchimento funzionale
Le analisi di annotazione funzionale Gene Ontology (GO) e di arricchimento dei percorsi metabolici Kyoto Encyclopedia of Genes and Genomes (KEGG) per i geni differenzialmente espressi sono state effettuate utilizzando la piattaforma online di analisi bioinformatica (https://www.bioinformatics.com.cn, ultima consultazione il 4 maggio 2026). La piattaforma impiega il test esatto di Fisher per calcolare i valori P grezzi e il metodo di Benjamini-Hochberg (BH) è stato utilizzato per la correzione dei test multipli al fine di controllare il tasso di falsa scoperta (FDR). I termini GO (processo biologico, componente cellulare e funzione molecolare) e i percorsi KEGG con un FDR < 0,05 sono stati considerati significativamente arricchiti.
Docking molecolare
Le strutture cristalline delle proteine bersaglio selezionate, rappresentanti conformazioni attive, sono state ottenute dal RCSB Protein Data Bank (https://www.rcsb.org/). La preparazione delle proteine è stata eseguita in Discovery Studio, inclusa la rimozione di molecole d'acqua e eteroatomi, l'eliminazione di catene ridondanti, la riparazione di residui mancanti, l'aggiunta di atomi di idrogeno a pH 7,4, l'assegnazione di cariche di Gasteiger e la minimizzazione dell'energia. La tasca di legame è stata definita in base alle coordinate del ligando co-cristallizzato, utilizzando una griglia che comprendeva tutti i residui critici per il legame con il substrato. Il docking molecolare di PLB è stato effettuato utilizzando il server web SwissDock (http://www.swissdock.ch/) tramite il modulo Attracting Cavities. Le energie libere di legame sono state calcolate per valutare l'affinità, e la conformazione a energia più bassa per ciascun bersaglio è stata selezionata come posa finale di docking. I risultati delle simulazioni di docking molecolare sono stati visualizzati per verificarne la validità.
Simulazione di dinamica molecolare
Le simulazioni di dinamica molecolare (MD) sono state eseguite utilizzando GROMACS 2022.2. Il campo di forza Amber14SB è stato utilizzato per descrivere la proteina, con il sistema solvatato in acqua TIP3P. I parametri della plumbagina, inclusi le cariche parziali AM1-BCC e i tipi atomici GAFF2, sono stati generati utilizzando Antechamber, seguiti dalla conversione della topologia con ACPYPE e dall'assegnazione dei parametri ionici Joung-Cheatham. Ciascun complesso proteina-ligando è stato inserito in una scatola di simulazione tronco-dodecaedrica con una distanza minima tra proteina e bordo della scatola di 1,2 nm, solvatato con molecole d'acqua TIP3P e neutralizzato aggiungendo Na⁺/Cl⁻ 0,15 M. Dopo la minimizzazione dell'energia mediante l'algoritmo di massima pendenza (Fmax < 1.000 kJ·mol⁻1·nm⁻1), il sistema ha subito un'equilibrazione sequenziale NVT e NPT per 200 ps ciascuna a 298 K. Le simulazioni produttive sono state quindi eseguite per 200 ns in condizioni NPT con un passo di integrazione di 2 fs, utilizzando lo schema di cutoff Verlet. Le interazioni elettrostatiche a lungo raggio sono state calcolate con il metodo particle mesh Ewald (PME), mentre sia le interazioni elettrostatiche che quelle di van der Waals hanno impiegato una distanza di cutoff di 1,2 nm. L'algoritmo LINCS è stato utilizzato per vincolare i legami contenenti idrogeno. La temperatura è stata mantenuta a 298 K mediante il termostato di Nosé-Hoover e la pressione a 1 bar mediante il barostato di Parrinello-Rahman. Le coordinate sono state salvate ogni 10 ps per analisi successive. L'analisi e la visualizzazione delle traiettorie sono state effettuate utilizzando strumenti di GROMACS, VMD e PyMOL, mentre i calcoli dell'energia libera di legame MM-PBSA sono stati eseguiti con gmx_MMPBSA quando appropriato.
Valutazione dell'equilibrazione delle traiettorie di dinamica molecolare e protocollo di campionamento dell'energia libera MM-PBSA
Per ogni complesso proteina-plumbagina, è stata eseguita una simulazione molecolare dinamica all-atom di produzione della durata di 200 ns utilizzando GROMACS. Il punto temporale di 100 ns è stato impostato come limite di equilibrazione: il segmento da 0 a 100 ns è stato definito come fase di rilassamento conformazionale, durante la quale la struttura secondaria della proteina e la tasca di legame del ligando hanno subito continui aggiustamenti conformazionali; il segmento da 100 a 200 ns è stato identificato come plateau termodinamicamente stabile, come dimostrato dall'assenza di deriva unidirezionale negli RMSD, nel raggio di girazione (Rg), negli RMSF a livello di residuo, nella SASA sepolta del ligando, nei legami idrogeno intermolecolari e nelle energie di interazione ligando-recettore, che hanno mostrato soltanto piccole fluttuazioni in regime stazionario. Tutti i parametri cinetici quantitativi e le energie libere di legame MM/PBSA sono stati calcolati esclusivamente a partire dalla fase di equilibrio compresa tra 100 e 200 ns, da cui sono stati estratti uniformemente frame ogni 100 ps, ottenendo 1.000 istantanee di equilibrio per sistema come input per gmx_MMPBSA. Le prime 100 ns delle traiettorie di rilassamento sono state scartate per eliminare interferenze dovute alla deriva conformazionale nei calcoli dell'energia libera.
Analisi prognostica dei geni bersaglio nel cancro della prostata
Le analisi prognostiche sono state condotte utilizzando set di dati sul cancro alla prostata provenienti da The Cancer Genome Atlas (TCGA). I dati di sequenziamento dell'RNA (STAR-counts) e le corrispondenti informazioni cliniche sono stati ottenuti dal portale Genomic Data Commons (https://portal.gdc.cancer.gov). I conteggi dell'espressione genica sono stati convertiti in trascritti per milione (TPM) e normalizzati utilizzando log2(TPM + 1). Dopo l'esclusione dei campioni con informazioni cliniche incomplete, l'analisi ha incluso 498 casi. Il livello mediano di espressione di ciascun gene è stato utilizzato per suddividere i pazienti in gruppi ad alta e bassa espressione. Sono state eseguite analisi di sopravvivenza di Kaplan-Meier con il test del log-rank e regressioni univariate di Cox per valutare la sopravvivenza complessiva (OS) per SRC e la sopravvivenza libera da progressione (PFS) per MAPK3, riportando i rapporti di rischio (HRs) e gli intervalli di confidenza al 95% (CIs). La performance predittiva è stata ulteriormente valutata generando curve caratteristiche operanti del ricevitore (ROC) dipendenti dal tempo a 1, 3 e 5 anni. Le analisi statistiche sono state effettuate in R versione 4.0.3, considerando statisticamente significativi valori di P < 0,05.