Tutti gli esperimenti sono stati approvati dal Consiglio di Revisione Istituzionale dell'Università di Medicina Cinese di Nanchino (Codice di Protocollo: 202501A018) e sono stati condotti in conformità con la Guida per la Cura e l'Uso degli Animali da Laboratorio degli Istituti Nazionali di Salute (NIH). Tutte le procedure sperimentali rispettavano le linee guida ARRIVE.
Farmaci sperimentali
La soluzione campione YQHXD conteneva 20 g di radice di membranaceus di Astragalus , 20 g di radice di Angelica sinensis , 15 g di rizoma Ligusticum chuanxiong , 6 g di radice di Stephania tetrandra , 15 g di frutto di Chaenomeles speciosa e 10 g di semi maturi di Sinapis alba , con un peso totale grezzo di 86 g per prescrizione. Tutti i materiali erboristici sono stati autenticati secondo la Farmacopea Cinese (edizione 2020) dal Dipartimento di Farmacia dell'Ospedale Ruikang affiliato all'Università di Medicina Cinese del Guangxi. La soluzione campione YQHXD è stata preparata utilizzando un metodo di decoction acqueosa. È stata aggiunta acqua distillata sterile a 8 volte il volume totale di grezza del farmaco, e le erbe sono state immerse per 30 minuti fino a essere completamente sommerse. La miscela è stata poi bollita e decoctata a fuoco basso per 30 minuti prima della filtrazione. I residui sono stati estratti nuovamente con lo stesso volume di acqua distillata sterile seguendo la stessa procedura. I filtrati delle due estrazioni sono stati combinati e concentrati a 55,5 mL a pressione ridotta tramite evaporazione rotativa per ottenere la soluzione di stocco YQHXD. In base al peso totale del farmaco grezzo, la concentrazione finale del farmaco grezzo era di 1,55 g/mL, con una resa di estrazione di circa il 64,5%. Il dosaggio somministrato ai ratti sperimentali è stato convertito dalla dose clinica adulta utilizzando il metodo di normalizzazione della superficie corporea. Il coefficiente di conversione tra umani e ratti era di circa 6,3; Pertanto, la dose per ratti (mg/kg) è stata calcolata come la dose umana (mg/kg) × 6,3. Il gruppo ad alta dose ha ricevuto 15,48 g/kg/giorno come dose di riferimento, mentre i gruppi a dose media e bassa hanno ricevuto rispettivamente 7,74 g/kg/giorno e 3,87 g/kg/giorno. Per un volume di gavaga di 10 mL/kg, le corrispondenti concentrazioni operative erano rispettivamente 1,548, 0,774 e 0,387 g/mL per i gruppi ad alta, media e bassa dose. La soluzione standard è stata diluita di conseguenza per raggiungere il dosaggio target per ciascun gruppo sperimentale prima della somministrazione. Dopo la preparazione, i decocti venivano raffreddati, allicitati in condizioni sterili e conservati a −40 °C fino all'uso.
Analisi farmacologica delle reti
Previsione delle proteine target dei principi attivi
Il 24 novembre 2025, composti attivi sono stati recuperati dal database TCMSP utilizzando i termini di ricerca "Huangqi, Angelica, Chuanxiong, Stephania, Mugua e Jiezi." I composti sono stati selezionati in base alla biodisponibilità orale (OB ≥ 30%) e alla somiglianza del farmaco (DL ≥ 0,18)9. La motivazione di queste soglie è la seguente: il database TCMSP integra dati di composti provenienti da fonti autorevole come PubChem, ChEMBL e ChemSpider e fornisce parametri relativi all'ADME, inclusa la biodisponibilità orale. La biodisponibilità orale si riferisce alla proporzione di una dose somministrata per via orale che raggiunge la circolazione sistemica senza variazioni, e un alto livello di OB è essenziale per identificare molecole bioattive. La somiglianza farmacologica è un indicatore utilizzato per valutare la farmacogabilità dei composti e per prevedere la loro solubilità e stabilità. In questo studio, è stata applicata una soglia di OB più rigorosa (≥30%), superiore al valore raccomandato dal database (≥20%), per identificare composti candidati con assorbimento orale favorevole e profili farmacocinetici. Nel frattempo, DL è stato utilizzato per stimare il potenziale farmacologico dei composti, e è stato adottato un cutoff di DL ≥ 0,18 basato sul valore medio di DL osservato nei composti della medicina tradizionale cinese. I composti attivi selezionati sono stati recuperati in formato SMILES dal database PubChem e importati nel database SwissTargetPrediction per la previsione dei bersagli, impostando la specie su Homo sapiens. Le proteine target sono state standardizzate utilizzando simboli genici ufficiali, e sono stati mantenuti solo bersagli con valori di probabilità superiori a 0 per migliorare l'affidabilità della previsione. I composti senza obiettivi previsti sono stati esclusi dall'analisi di rete successiva. Dopo aver unito tutti i target, le voci duplicate venivano rimosse tramite software di fogli di calcolo per garantire che ogni target fosse rappresentato solo una volta nel dataset finale.
Raccolta di obiettivi chiave
I geni correlati alla malattia sono stati recuperati dai database GeneCards10, OMIM11 e TTD12 utilizzando la parola chiave "degenerazione del disco lombare". I bersagli genici raccolti sono stati fusi e le voci duplicate sono state rimosse per ottenere un insieme finale di bersagli associati alla malattia.
Potenziali obiettivi di YQHXD per il trattamento della LDD
La piattaforma Venny 2.1.0 veniva utilizzata per determinare i potenziali obiettivi di YQHXD in LDD. È stata quindi generata una rete per mostrare le associazioni tra i componenti attivi e i loro obiettivi previsti. I componenti attivi principali sono stati selezionati in base ai loro valori di grado nella rete.
Costruzione di reti di componenti farmaceutici e bersagli di intersezione di malattie
I principi attivi, le proteine bersaglio e i geni correlati alla LDD dalla formula TCM YQHXD sono stati importati nel softwareCytoscape 13 per costruire un diagramma di rete farmaco-componente-bersaglio-malattia.
Costruzione della rete PPI e identificazione dei principali target principali
Bersagli intersecanti sono stati inviati alla piattaformaSTRING 14 per l'analisi, con l'organismo impostato su Homo sapiens. Per costruire la rete di interazione proteina-proteina (PPI), la soglia del punteggio di interazione è stata impostata a >0,4, i nodi disconnessi sono stati nascosti e tutte le altre impostazioni sono state lasciate come predefinite. I dati risultanti venivano importati nel software, dove lo strumento NetworkAnalyzer veniva utilizzato per valutare le proprietà topologiche e calcolare i gradi dei nodi. I primi cinque nodi, classificati per grado, sono stati identificati come obiettivi chiave dell'hub.
Analisi dell'arricchimento dei target principali
I target intersecanti sono stati caricati nel databaseDAVID 15 per l'annotazione funzionale Gene Ontology (GO) e l'analisi di arricchimento dei percorsi della Kyoto Encyclopedia of Genes and Genomes (KEGG). Il background dell'organismo è stato impostato come Homo sapiens, con il cutoff di arricchimento definito secondo i parametri DAVID predefiniti, inclusi una soglia EASE < 0,1 e una soglia di conteggio ≥ 2. L'universo genico si basava sul background predefinito del genoma completo fornito da DAVID. I test multipli sono stati corretti utilizzando il metodo Benjamini-Hochberg per controllare il tasso di false scoperte (FDR), e i risultati DAVID sono stati adeguati su questa base. I percorsi arricchiti sono stati ordinati per valore di conteggio in ordine decrescente. I grafici a barre di arricchimento GO e KEGG sono stati prodotti utilizzando una piattaforma di bioinformatica online.
Aggancio molecolare
Per prevedere le interazioni tra i componenti attivi chiave di YQHXD e i bersagli principali di LDD, in questo studio è stata effettuata un'analisi di docking molecolare. I primi cinque composti attivi chiave, classificati secondo i valori topologici della rete e la rilevanza del bersaglio, sono stati selezionati come liganti, mentre le corrispondenti proteine target principali sono state utilizzate come recettori. I file SDF dei ligandi sono stati recuperati da PubChem, e le strutture tridimensionali delle proteine target sono state scaricate dal Protein Data Bank (PDB). L'analisi dell'attracco è stata effettuata sulla piattaforma CB-Dock2 utilizzando il suo protocollo di attracco cieco completamente automatizzato. Dopo aver caricato le strutture proteiche, la piattaforma rimuoveva automaticamente gli eteroatomi, inclusi piccole molecole, ioni e molecole d'acqua, e successivamente riparava gli atomi mancanti. È stata utilizzata la modalità Auto Blind Docking. In questa modalità, CB-Dock2 identifica potenziali sacche di legame proteico-superficie con un algoritmo di rilevamento delle cavità basato sulla curvatura e poi esegue calcoli di attracco utilizzando il motore Vina16 integrato AutoDock. Durante il processo di aggancio, le proteine sono state trattate come strutture rigide, mentre i liganti sono stati considerati flessibili e tutti i parametri di aggancio sono stati impostati ai valori predefiniti della piattaforma. I risultati dell'agganciamento sono stati sottoposti a visualizzazione tridimensionale e analisi di interazione utilizzando gli strumenti integrati di CB-Dock217, e sono state calcolate le energie di legame per stimare le affinità di legame previste tra i composti attivi e le proteine bersaglio.
Simulazione della dinamica molecolare
Sono state effettuate simulazioni di dinamica molecolare per il complesso proteina-ligando target utilizzando GROMACS 2022.2. La proteina è stata modellata con il campo di forza Amber14SB, e il sistema solvente è stato rappresentato utilizzando il modello idrico TIP3P. Le cariche atomiche dei ligandi venivano calcolate utilizzando il metodo AM1-BCC implementato in Antechamber, seguiti dall'assegnazione dei tipi di atomi GAFF2. I file topologici risultanti furono convertiti in un formato compatibile con GROMACS utilizzando ACPYPE. I parametri ionici coerenti con il modello idrico TIP3P sono stati assegnati in base al set di parametri Joung-Cheatham. Il complesso di partenza è stato inserito in una scatola di simulazione dodecaedra troncata con almeno 1,2 nm tra la superficie della proteina e il bordo della scatola. Sono state aggiunte molecole d'acqua TIP3P per solventare il sistema, e sono stati introdotti ioni Na+ e Cl− per neutralizzare la carica del sistema e impostare la forza ionica a 0,15 M. La minimizzazione dell'energia è stata eseguita con l'algoritmo di discesa più ripida fino a quando la forza massima è stata inferiore a 1.000 kJ·mol−1·nm−1. Il sistema veniva quindi bilanciato per 100 ps sotto l'ensemble NVT a 298 K e 1 bar. Dopo l'equilibrio, veniva eseguita una simulazione di produzione di 100 ns sotto l'insieme NPT con un passo temporale di 2 fs. Le interazioni non legate venivano gestite con il sistema di corte Verlet. Le interazioni elettrostatiche venivano calcolate utilizzando il metodo Ewald a rete particellare (PME), e veniva utilizzato un cutoff di 1,2 nm sia per le interazioni di van der Waals che per quelle di Coulomb. I legami contenenti idrogeno erano vincolati dall'algoritmo LINCS. La temperatura veniva controllata con il termostato Nosé-Hoover, e la pressione con il barostato Parrinello-Rahman, mantenendo il sistema a 298 K e 1 bar. Le traiettorie venivano salvate ogni 10 ps. La dinamica strutturale e le interazioni intermolecolari sono state analizzate utilizzando gli strumenti integrati di GROMACS, VMD e PyMOL. L'energia libera di legame è stata stimata con gmx_MMPBSA dove applicabile.
Studi sugli animali
Allevamento degli animali
Settanta maschi di ratti Sprague-Dawley (SD) di età compresa tra 3 e 8 settimane sono stati assegnati casualmente a sette gruppi sperimentali (n = 10 per gruppo) utilizzando numeri casuali generati al computer, assicurando che non ci siano differenze statisticamente significative nel peso corporeo di base tra i gruppi. Tutti gli animali sono stati alloggiati in condizioni controllate (22 ± 1 °C, 55 ± 10% di umidità relativa, 12 ore di ciclo chiaro/buio) con libero accesso a cibo e acqua. Le posizioni della gabbia sono state randomizzate ad ogni cambio di gabbia per minimizzare il bias ambientale. Dopo una settimana di acclimatazione, i ratti dei gruppi modello hanno subito una puntura caudale IVD per indurre la degenerazione. È stata somministrata un'anestesia di infiltrazione locale con bupivacaina preoperatoriamente per ridurre il dolore. Gli animali sono stati monitorati attentamente durante l'intervento chirurgico e durante il recupero post-operatorio. Le osservazioni postoperatorie includevano tempi di recupero, ripristino dell'attività spontanea, assunzione di cibo e acqua e cambiamenti di peso corporeo. Una settimana dopo l'istituzione del modello, due ratti di ciascun gruppo sono stati selezionati casualmente per l'imaging radiografico digitale (DR) post-operatorio e la colorazione di ematossilina ed eosina (HE). L'analisi radiografica ha mostrato una riduzione dell'altezza della IVD caudale, mentre la colorazione con HE ha rivelato una marcata riduzione delle cellule notocordiali del nucleo pulposo (NP) e una disorganizzazione della struttura lamelare dell'anulusfibrosus 18, confermando il successo dell'istituzione del modello. Otto settimane dopo l'intervento, i ratti sono stati soppressi per inalazione di un'overdose di isoflurano, e sono stati raccolti campioni di tessuto caudal e sangue. Dopo la validazione preliminare, sono stati utilizzati tre campioni indipendenti per gruppo che soddisfacevano i criteri di qualità predefiniti per le analisi finali di istologia, ELISA, PCR e western blot. Le sezioni istologiche sono state valutate in condizioni di cieco singolo, con vetrine etichettate solo tramite numeri identificabili degli animali. Per ELISA, PCR e western blot, è stata applicata la codifica cieca durante l'aliquotazione e l'etichettatura del campione, e l'ordine di test di tutti i campioni è stato randomizzato prima degli esperimenti. I dati grezzi, inclusi i valori di assorbanza, valori Ct e intensità delle bande, venivano decodificati in identità di gruppo da uno statistico non coinvolto nel raggruppamento sperimentale o nell'elaborazione del campione e successivamente analizzati dal personale statistico. Tutte le analisi statistiche sono state effettuate da statistici non coinvolti nel raggruppamento sperimentale o nell'elaborazione del campione. Eventuali eventi avversi, inclusi infezioni della ferita, deiscenza, sanguinamento o letargia grave, venivano immediatamente registrati e valutati. Le condizioni lievi sono state attentamente monitorate con cure di supporto rafforzate. Sono stati applicati parametri umani se i ratti hanno mostrato perdita di peso corporea superiore al 20% rispetto al valore di base, dolore intenso non alleviato dagli analgesici, perdita di mobilità, incapacità di accedere a cibo o acqua, infezione, necrosi, altre gravi complicazioni post-operatorie, difficoltà respiratoria o qualsiasi condizione potenzialmente letale. In tali casi, gli animali venivano soppressi umanamente per inalazione di un'overdose di isoflurano.
Modellazione sperimentale
I ratti sono stati anestesiati con isoflurano inalato e posizionati a terra prono. Le IVD di coda Co6/7, Co7/8 e Co8/9 erano localizzate sotto guida fluoroscopica a raggi X. Un ago da 21 G veniva avanzato verticalmente nella IVD fino a raggiungere il lato opposto della pelle, dopodiché la punta dell'ago veniva ritirata verso il centro della coda. Un secondo ago serviva come riferimento posizionale. Per indurre una degenerazione persistente della IVD, l'ago veniva ruotato di 360°, mantenuto in posizione per 30 secondi e poi rimosso. Nel gruppo Sham, la profondità della puntura attraverso la pelle e il muscolo era di circa 2 mm, senza danneggiare la flebo. Negli altri gruppi, l'ago veniva inserito 5 mm verso il centro della IVD, evitando i vasi sanguigni. Dopo l'inserimento, l'ago veniva ruotato di 360° e lasciato in posizione per 30 secondi. Il sito della puntura è stato disinfettato, segnato e protetto con garze sterili. Una settimana dopo la modellazione, le immagini a raggi X della colonna caudale hanno mostrato una diminuzione dell'altezza della IVD. La colorazione con HE ha rivelato una significativa riduzione delle cellule notocordali NP e una disorganizzazione della struttura della lamella annulare, confermando l'istituzione del modello di successo. Dopo una modellazione di successo, ai gruppi che ricevevano trattamenti con MTC e al gruppo di controllo positivo (PC) sono state somministrate le dosi appropriate, calcolate in base al rapporto sulla superficie corporea tra umani e ratti.
Raggruppamento e intervento
In questo esperimento sono stati utilizzati in totale 70 maschi di ratti SD di età 3–8 settimane, suddivisi in sette gruppi:
(1) Gruppo di controllo (Ctrl): I ratti di questo gruppo sono stati assegnati casualmente nello stesso momento e nelle stesse condizioni dei gruppi sperimentali. Venivano ospitati in gabbie individuali in una struttura per animali a livello SPF, con accesso gratuito al cibo standard e all'acqua potabile.
(2) Gruppo finto: I ratti sono stati anestetizzati inalando isoflurano. Dopo l'anestesia, la pelle della coda veniva disinfestata con povidone-iodio e i ratti venivano posizionati a terra. Le IVD Co6/7 e Co7/8 erano localizzate sotto guida fluoroscopica a raggi X. L'ago è stato avanzato di circa 2 mm attraverso la pelle e il muscolo senza danneggiare la flebo IVD. L'acqua distillata (10 mL/kg) veniva somministrata una volta al giorno tramite gavage per 8 settimane.
(3) Gruppo modello: I ratti sono stati anestetizzati con isoflurano inalato. Dopo l'anestesia, la pelle della coda veniva disinfestata con povidone-iodio e i ratti venivano posizionati a terra. Sotto guida fluoroscopica a raggi X, sono stati individuati i IVD Co6/7, Co7/8 e Co8/9. Un ago da 21 G veniva avanzato verticalmente nella IVD fino a raggiungere il lato opposto della pelle, e la punta dell'ago veniva poi ritratta verso il centro della coda.
Un secondo ago veniva usato come riferimento posizionale. L'ago veniva ruotato di 360° e mantenuto fermo per 30 secondi per indurre la degenerazione continua della IVD prima della rimozione. L'acqua distillata (10 mL/kg) veniva somministrata quotidianamente tramite gavage per 8 settimane, e il recupero dei ratti veniva monitorato attentamente dopo l'intervento.
(4) gruppo PC: Il gruppo di controllo positivo è stato sottoposto alla stessa procedura di modellazione del gruppo modello e successivamente ha ricevuto una gara orale di celecoxib (20 mg/kg/giorno), preparata come sospensione in carbossimetilcellulosa al 0,5% di sodio (CMC-Na), per 8 settimane consecutive, con un volume gavage di 10 mL/kg di peso corporeo.
(5) YQHXD gruppo a bassa dose [YQHXD(L)]: La modellazione è stata eseguita come nel gruppo modello. I ratti hanno ricevuto YQHXD a 3,87 g/kg/giorno tramite gavage per 8 settimane, con un volume gavage di 10 mL/kg e una concentrazione operativa di 0,387 g/mL.
(6) YQHXD gruppo a dosi medie [YQHXD(M)]: La modellazione è stata eseguita come nel gruppo modello. I ratti hanno ricevuto YQHXD a 7,74 g/kg/giorno per gavage per 8 settimane, con un volume gavage di 10 mL/kg e una concentrazione operativa di 0,774 g/mL.
(7) YQHXD gruppo ad alta dose [YQHXD(H)]: La modellazione è stata eseguita come nel gruppo modello. I ratti hanno ricevuto YQHXD a 15,48 g/kg/giorno per gavage per 8 settimane, con un volume gavage di 10 mL/kg e una concentrazione operativa di 1,548 g/mL.
Cambiamenti patologici
Il tessuto caudale IVD del ratto è stato isolato chirurgicamente e sono stati osservati cambiamenti nella IVD e nei tessuti muscolari circostanti. Il tessuto raccolto veniva poi lavorato in sezioni patologiche, con i campioni preparati in sezioni longitudinali per l'osservazione morfologica. La IVD vertebrale caudale, fissata in formaldeide per 96 ore, è stata rapidamente decalcificata usando una soluzione decalcificante, disidratata tramite una serie alcolica gradata e eliminata con xilene. Dopo l'immersione nella paraffina, sono state preparate sezioni di paraffina da 5 μm. Dopo la deparaffinizzazione con xilene e la disidratazione attraverso una serie di gradienti alcolici, le sezioni sono state colorate con ematossilina ed eosina. Furono poi disidratati, liberati e montati. I campioni sono stati osservati al microscopio ottico per esaminare i cambiamenti morfologici nella struttura caudale IVD tra i vari gruppi. Successivamente fu eseguita la colorazione tricroma di Masson. Le sezioni sono state colorate sequencialmente con soluzione di dicromato di potassio ed ematossilina per 3 minuti ciascuna, seguite da differenziazione con etanolo acido, e poi colorate con Picrosirius Red per 10 minuti. Dopo il trattamento con acido fosfotongstico per 10 minuti, le sezioni sono state colorate con blu di toluidina per 5 minuti. Dopo la disidratazione, le sezioni sono state pulite con xilene, montate e osservate al microscopio per esaminare i cambiamenti nella IVD caudale.
Rilevamento ELISA delle proteine bersaglio sierie
Otto settimane dopo l'intervento, i ratti sono stati anestetizzati con isoflurano inalato e è stato prelevato sangue aortico addominale. I campioni di sangue periferico sono stati conservati a temperatura ambiente per 1 ora e centrifugati a 1.400 × g per 10 minuti. Il siero veniva separato e conservato a −80 °C fino all'uso. I livelli sierici di AKT1, STAT3 e Bcl-2 sono stati quantificati tramite un saggio immunoenzimatico legato (ELISA). I test venivano eseguiti secondo i protocolli dei produttori. L'assorbanza è stata rilevata a 450 nm e le concentrazioni sono state determinate da curve standard. I range di rilevamento per AKT1, STAT3 e BCL-2 erano rispettivamente 0,625–40 ng/mL, 31,25–1000 pg/mL e 0,16–10 ng/mL, con sensibilità di 0,1 ng/mL, 10 pg/mL e 0,1 ng/mL. Tutti i campioni di siero sono stati analizzati senza diluizione e tutti i risultati rientravano nei corrispondenti intervalli standard di curva. Tutti i campioni e gli standard sono stati analizzati in duplicato e la normalizzazione è stata effettuata sulla base di volumi sierici uguali.
Rilevamento PCR dei livelli di espressione genica
I tessuti del disco intervertebrale sono stati sezionati dalle vertebre caudali del ratto e sono stati raccolti circa 20 mg da ciascun campione. I tessuti sono stati polverizzati in azoto liquido e lisati in 500 μL di tampone di lisi non enzimatica per 15 minuti. Veniva quindi estratto l'RNA totale e la sua concentrazione e purezza venivano misurate con uno spettrofotometro. Campioni con valori A260/A 280 di 1,8–2,0 sono stati accettati per l'analisi e le concentrazioni di RNA sono state regolate a 500 ng/μL. La trascrizione inversa è stata effettuata secondo il protocollo del produttore. La PCR quantitativa (qPCR) è stata eseguita su un sistema PCR in tempo reale per misurare l'espressione di mRNA di Akt1, Bcl-2 e Stat3 . Le sequenze di primer sono fornite nella Tabella 1. Il programma PCR consisteva in una denaturazione iniziale a 95 °C per 30 secondi, seguita da 40 cicli a 95 °C per 5 secondi e 60 °C per 30 secondi. Dopo l'amplificazione è stata condotta un'analisi della curva di fusione per verificare la specificità del prodotto. Il Gapdh è stato utilizzato come gene di riferimento e tutti i campioni sono stati analizzati in triplice. L'efficienza dell'amplificazione degli innesco era del 90%–110%. L'espressione genica relativa è stata calcolata utilizzando il metodo 2−ΔΔCt .
Western blotting per i livelli di espressione proteica
La proteina totale è stata isolata dai tessuti caudali del disco intervertebrale di ciascun gruppo. I tessuti sono stati lisati in un tampone RIPA integrato con inibitori della proteasi. La concentrazione proteica è stata misurata utilizzando un kit di analisi proteica BCA e regolata a 2 μg/μL. Quantità proteiche uguali (15 μg per campione) sono state separate da SDS-PAGE al 10% e trasferite alle membrane di difluoruro di polivinilidene (PVDF) a corrente costante di 300 mA per 90 minuti utilizzando un sistema di trasferimento umido. Le membrane sono state poi bloccate con latte scremo al 5% per 1 ora a temperatura ambiente e incubate durante la notte a 4 °C con anticorpi primari contro AKT1 (1:5.000), STAT3 (1:6.000), Bcl-2 (1:1.500) e β-actin (1:4.000). Dopo il lavaggio con TBST, le membrane sono state incubate con anticorpi secondari coniugati con perossidasi di rafano (HRP) a una diluizione di 1:5.000 per 1 ora a temperatura ambiente. Le bande proteiche sono state rilevate utilizzando un reagente a chemiluminescenza potenziata (ECL) e le immagini sono state acquisite con un sistema di imaging a gel in condizioni di esposizione identiche. Le intensità delle bande sono state misurate utilizzando ImageJ e normalizzate in β-actin per determinare l'espressione proteica relativa. Ogni esperimento è stato condotto indipendentemente in triplice, con ogni campione ottenuto da un animale diverso.
Statistiche
I risultati sono espressi come media ± deviazione standard (SD). Le differenze tra i gruppi sono state valutate tramite analisi unidirezionale della varianza (ANOVA), seguite dal test post hoc della differenza minima significativa (LSD). La significatività statistica è stata definita come P < 0,05.