Summary
Il presente protocollo utilizza un pacchetto di simulazione biomolecolare e descrive l'approccio di dinamica molecolare (MD) per modellare la caspasi wild-type e le sue forme mutanti. Il metodo MD consente di valutare l'evoluzione dinamica della struttura delle caspasi e il potenziale effetto di mutazioni o modificazioni post-traduzionali.
Abstract
L'apoptosi è un tipo di morte cellulare programmata che elimina le cellule danneggiate e controlla lo sviluppo e l'omeostasi tissutale degli organismi multicellulari. Le caspasi, una famiglia di proteasi della cisteina, svolgono un ruolo chiave nell'inizio e nell'esecuzione dell'apoptosi. La maturazione delle caspasi e la loro attività sono messe a punto da modifiche post-traduzionali in modo altamente dinamico. Per valutare l'effetto dei cambiamenti post-traduzionali, i siti potenziali vengono regolarmente mutati con residui persistenti a qualsiasi modifica. Ad esempio, il residuo di serina viene sostituito con alanina o acido aspartico. Tuttavia, tali sostituzioni potrebbero alterare la conformazione del sito attivo della caspasi, portando a disturbi nell'attività catalitica e nelle funzioni cellulari. Inoltre, mutazioni di altri residui di amminoacidi situati in posizioni critiche potrebbero anche rompere la struttura e le funzioni delle caspasi e portare a perturbazioni dell'apoptosi. Per evitare le difficoltà di impiegare residui mutati, approcci di modellazione molecolare possono essere facilmente applicati per stimare il potenziale effetto delle sostituzioni aminoacidiche sulla struttura delle caspasi. Il presente protocollo consente la modellazione sia della caspasi wild-type che delle sue forme mutanti con il pacchetto di simulazione biomolecolare (Amber) e le strutture di supercomputer per testare l'effetto delle mutazioni sulla struttura e la funzione della proteina.
Introduction
L'apoptosi è uno dei processi cellulari più studiati che regolano la morfogenesi e l'omeostasi tissutale degli organismi pluricellulari. L'apoptosi può essere iniziata da una vasta gamma di stimoli esterni o interni, come l'attivazione dei recettori di morte, il disturbo nei segnali del ciclo cellulare, il danno al DNA, lo stress del reticolo endoplasmatico (ER) e varie infezioni batteriche e virali1. Le caspasi - attori chiave dell'apoptosi - sono convenzionalmente classificate in due gruppi: iniziatori (caspasi-2, caspasi-8, caspasi-9 e caspasi-10) ed effettori (caspasi-3, caspasi-6 e caspasi-7), a seconda della loro struttura di dominio e del posto nella cascata delle caspasi 2,3. Sui segnali di morte cellulare, le caspasi iniziatrici interagiscono con molecole adattatrici che facilitano la dimerizzazione indotta dalla prossimità e l'autoelaborazione per formare un enzima attivo. Le caspasi effettrici sono attivate attraverso la scissione da parte delle caspasi iniziatrici ed eseguono fasi di esecuzione a valle scindendo più substrati cellulari4.
La maturazione e la funzione delle caspasi iniziatrici ed effettrici sono regolate da un gran numero di diversi meccanismi intracellulari, tra i quali la modificazione post-traduzionale svolge un ruolo indispensabile nella modulazione della morte cellulare5. L'aggiunta di gruppi modificanti (fosforilazione, nitrosilazione, metilazione o acetilazione) o proteine (ubiquitinazione o SUMOilazione) modifica l'attività enzimatica delle caspasi o la conformazione e la stabilità proteica che regolano l'apoptosi. La mutagenesi sito-diretta è ampiamente applicata per studiare i potenziali siti di modificazione post-traduzionale e discernere il loro ruolo. Un sito di modifica putativo viene solitamente sostituito da un altro amminoacido, che non può essere ulteriormente modificato. Pertanto, la serina e la treonina potenzialmente fosforilate sono mutate in alanina e i siti di ubiquitinazione della lisina vengono sostituiti con arginina. Un'altra strategia include la sostituzione di un amminoacido che imita in particolare la modificazione post-traduzionale (ad esempio, glutammato e aspartato sono stati usati per imitare la serina o la treonina fosforilate)6. Tuttavia, alcune di queste sostituzioni situate nelle alte vicinanze di un sito attivo o in posizioni critiche potrebbero cambiare la struttura delle caspasi, disturbare l'attività catalitica e sopprimere la morte cellulare apoptotica7. Effetti simili potrebbero essere osservati nei casi di mutazioni missenso associate al tumore nei geni delle caspasi. Ad esempio, la mutazione associata al tumore della caspasi-6 - R259H - ha provocato cambiamenti conformazionali degli anelli nella tasca di legame del substrato, riducendo l'efficiente turnover catalitico dei substrati8. La sostituzione dell'aminoacido G325A nella caspasi-8 identificata nel carcinoma a cellule squamose della testa e del collo potrebbe ostacolare l'attività della caspasi-8, che ha portato alla modulazione della segnalazione del fattore nucleare-kB (NF-kB) e ha promosso la tumorigenesi9.
Per valutare il potenziale effetto delle sostituzioni aminoacidiche sulla struttura e la funzione delle caspasi, è possibile applicare la modellazione molecolare. L'approccio della dinamica molecolare (MD) è descritto in questo lavoro per modellare la caspasi wild-type e le sue forme mutanti usando il pacchetto di simulazione biomolecolare (Amber). Il metodo MD fornisce una visione dell'evoluzione dinamica della struttura proteica in seguito all'introduzione di mutazioni. Originariamente sviluppato dal gruppo di Peter Kollman, il pacchetto Amber è diventato uno degli strumenti software più popolari per le simulazioni biomolecolari10,11,12,13. Questo software è diviso in due parti: (1) AmberTools, una raccolta di programmi utilizzati abitualmente per la preparazione del sistema (assegnazione del tipo di atomo, aggiunta di idrogeni e molecole esplicite di acqua, ecc.) e l'analisi della traiettoria; e (2) Amber, che è incentrato sul programma di simulazione pmemd. AmberTools è un pacchetto gratuito (e un prerequisito per l'installazione di Amber stesso), mentre Amber è distribuito con una licenza separata e una struttura tariffaria. Simulazioni parallele su un supercomputer e/o utilizzando unità di elaborazione grafica (GPU) possono migliorare sostanzialmente le prestazioni per la ricerca scientifica della dinamica della struttura proteica14. Le ultime versioni del software disponibili sono AmberTools21 e Amber20, ma i protocolli descritti possono essere utilizzati anche con le versioni precedenti.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
1. Preparazione del sistema
NOTA: I modelli molecolari delle forme proteiche native e mutanti sono costruiti sulla base di un'appropriata struttura cristallina ottenuta dalla Protein Data Bank15,16.
- Per recuperare la struttura PDB selezionata, utilizzare l'elenco a discesa Scarica file e fare clic su Formato PDB. Rimuovere le osservazioni e i dati di connettività e inserire una scheda TER tra catene proteiche separate nel file PDB. Facoltativamente, impostare lo stato di protonazione dei residui di istidina sostituendo il nome del residuo HIS con HIE, HID o HIP (file supplementare 1).
NOTA: è ragionevole non rimuovere molecole d'acqua risolte cristallograficamente (se presenti nel file PDB). - Per preparare il modello iniziale, avviare il programma tleap (pacchetto AmberTools) e quindi immettere comandi che manipolano gli oggetti tleap .
- Avvia il programma: tleap. Caricare il campo di forza ff14SB per descrivere la proteina con la meccanica molecolare: > fonte leaprc.protein.ff14SB.
- Campi di forza di carico per molecole d'acqua e ioni atomici (Na+, Cl-): > sorgente leaprc.water.tip3p.
- Caricare il file PDB e costruire le coordinate per gli idrogeni, creando un oggetto chiamato mol: > mol = loadpdb protein.pdb.
- Verificare la presenza di incongruenze interne che potrebbero causare problemi: > controllare mol.
NOTA: i ponti disolfuro, se presenti, devono essere specificati manualmente. Sostituire i nomi delle cisteine legate covalentemente CYS con CYX e digitare il seguente comando in tleap per creare un legame tra gli atomi SG: legame mol.X.SG mol.Y.SG (X e Y sono numeri di residui di cisteina). - Creare una scatola di solventi attorno alla proteina (acqua TIP3P, distanza 12 Å): > solvatebox mol TIP3PBOX 12.0.
- Controllare la carica totale: > mol di carica e aggiungere controioni (Na+ o Cl-) per neutralizzare il sistema:
> addions mol Na+ 0. - Creare il file prmtop della topologia e il file inpcrd delle coordinate (input per l'esecuzione di una simulazione): > saveamberparm mol protein.prmtop protein.inpcrd. Esci dal programma tleap: > esci.
NOTA: i file di coordinate in formato ambra possono essere facilmente convertiti in formato PDB utilizzando lo strumento ambpdb (fare riferimento al Manuale dell'utente, vedere Tabella dei materiali).
2. Minimizzazione dell'energia
NOTA: La minimizzazione dell'energia è necessaria per rimuovere eventuali contatti difettosi e sovrapposizioni tra gli atomi nel sistema di partenza che portano all'instabilità durante l'esecuzione di MD.
- Eseguire il primo stadio della minimizzazione dell'energia (2.500 passi dell'algoritmo di discesa più ripida + 2.500 passi dell'algoritmo del gradiente coniugato) per ottimizzare le posizioni degli atomi di idrogeno aggiunti e delle molecole d'acqua mantenendo le coordinate proteiche fissate da restrizioni posizionali su atomi pesanti.
- Eseguire il programma pmemd come segue:
pmemd -O -i min1.in -p protein.prmtop -c protein.inpcrd -o protein_min1.out
-r protein_min1.rst -ref protein.inpcrd. - Seguire gli argomenti richiesti: -i dati di controllo FILE; -p FILE topologia molecolare, parametri del campo di forza, nomi degli atomi; -c FILE coordinate iniziali; -o FILE output del registro leggibile dall'utente; -r FILE coordinate finali; -ref FILE coordinate di riferimento per i vincoli di posizione.
- Specificare le opzioni nel file di input min1.in:
&cntrl
imin=1, maxcyc=5000, ncyc=2500,
cut=10.0, ntb=1,
ntc=1, ntf=1,
ntpr=10,
ntr=1,
restraintmask=':1-517 & !@H=',
restraint_wt=2,0
/
NOTA: Opzioni di input: imin = 1 eseguire la minimizzazione; maxcyc=5000 numero massimo di cicli di minimizzazione; Il metodo di minimizzazione NCYC=2500 passa dalla discesa più ripida al gradiente coniugato dopo i cicli NCYC ; cut=10.0 cutoff non legato (Å); ntb=1 vengono imposti limiti periodici, volume costante; ntc=1 non vengono applicati vincoli alle lunghezze di legame (per una migliore convergenza energetica); ntf=1 vengono calcolate tutte le interazioni; ntpr=10 Le informazioni sull'energia vengono stampate nel file out ogni passaggio di NTPR ; ntr=1 flag per trattenere atomi specificati usando un potenziale armonico; restraintmask=':1-517 & !@H=' specifica gli atomi trattenuti; restraint_wt=2,0 peso (kcal/mol∙Å2) per i vincoli di posizione.
- Eseguire il programma pmemd come segue:
- Eseguire la seconda fase della minimizzazione dell'energia (5.000 gradini di discesa più ripidi + 5.000 gradini di pendenza coniugati) senza vincoli per ottimizzare l'intero sistema.
- Utilizzare min2.in e protein_min1.rst come input: pmemd -O -i min2.in -p protein.prmtop
-c protein_min1.rst -o protein_min2.out -r protein_min2.rst. - Specificare le opzioni nel file di input min2.in:
&cntrl
imin=1, maxcyc=10000, ncyc=5000,
cut=10.0, ntb=1,
ntc=1, ntf=1,
ntpr=10
/
- Utilizzare min2.in e protein_min1.rst come input: pmemd -O -i min2.in -p protein.prmtop
3. Riscaldamento
NOTA: questa fase ha lo scopo di riscaldare il sistema da 0 K a 300 K. Le velocità iniziali sono assegnate agli atomi poiché il modello iniziale basato sul file PDB non contiene informazioni sulla velocità.
- Eseguire il processo di riscaldamento con restrizioni posizionali sugli atomi proteici (50 ps, volume costante).
- Utilizzare heat.in e protein_min2.rst come input: pmemd -O -i heat.in -p protein.prmtop
-c protein_min2.rst -o protein_heat.out
-r protein_heat.rst -x protein_heat.mdcrd -ref protein_min2.rst - Seguire gli argomenti richiesti: set di coordinate -x FILE salvati sulla traiettoria MD.
- Specificare le opzioni nel file di input heat.in:
&cntrl
imin=0, irest=0, ntx=1,
nstlim=25000, dt=0,002,
ntc=2, ntf=2,
cut=10.0, ntb=1,
ntpr=500, ntwx=500,
ntt=3, gamma_ln=2.0,
tempi=0.0, temp0=300.0,
ntr=1, restraintmask=':1-517',
restraint_wt=1,0,
nmropt=1
/
&wt TYPE='TEMP0',
istep1=0, istep2=25000,
valore1=0.1, valore2=300.0 /
&wt TYPE='END' /
NOTA: Opzioni di input: imin = 0 eseguire MD; irest=0 eseguire una nuova simulazione; ntx=1 nessuna velocità iniziale viene letta dal file rst ; nstlim=25000 numero di passi MD; dt=0,002 passo temporale (ps); ntc=2 legami che coinvolgono idrogeno sono vincolati usando l'algoritmo SHAKE; ntf=2 le forze per i legami vincolati non sono calcolate; Le coordinate ntwx=500 vengono scritte nel file MDCRD ogni passo di NTWX ; NTT=3 Regolazione della temperatura di Langevin; gamma_ln=2.0 frequenza di collisione per la dinamica di Langevin (ps−1); tempi=0.0 temperatura iniziale; temp0=300.0 temperatura di riferimento; nmropt=1 i parametri specificati in una &wt namelist vengono letti; TYPE='TEMP0' varia la temperatura target.
- Utilizzare heat.in e protein_min2.rst come input: pmemd -O -i heat.in -p protein.prmtop
4. Equilibrio
NOTA: Questa fase è necessaria per regolare la densità dell'acqua e ottenere lo stato di equilibrio della proteina.
- Eseguire l'equilibrio a 300 K senza alcuna restrizione (500 ps, pressione costante).
- Utilizzare equil.in e protein_heat.rst come input: pmemd -O -i equil.in -p protein.prmtop
-c protein_heat.rst -o protein_equil.out -r protein_equil.rst -x protein_equil.mdcrd. - Specificare le opzioni nel file di input equil.in:
&cntrl
imin=0, irest=1, ntx=5,
nstlim=250000,
dt=0,002,
ntc=2, ntf=2,
cut=10.0, ntb=2, ntp=1, taup=2.0,
ntpr=1000, ntwx=1000, ntwr=50000,
ntt=3, gamma_ln=2.0,
temp0=300.0
/
NOTA: Opzioni di input: irest=1 riavviare la simulazione; Le coordinate e le velocità NTX=5 vengono lette da un file RST generato in precedenza; ntb=2 vengono imposti limiti periodici, pressione costante; NTP=1 scala della pressione isotropa; taup=2.0 tempo di rilassamento della pressione (PS); NTWR=50000 Ogni passaggio NTWR , il file RST viene scritto, garantendo il ripristino da un arresto anomalo.
- Utilizzare equil.in e protein_heat.rst come input: pmemd -O -i equil.in -p protein.prmtop
5. Dinamica produttiva
- Dopo che l'equilibrio è stato raggiunto con successo, eseguire una simulazione MD di produzione (10 ns o più) a pressione costante e generare il file di traiettoria per la successiva analisi della struttura proteica.
- Utilizzare prod.in e protein_equil.rst come input: pmemd -O -i prod.in -p protein.prmtop
-c protein_equil.rst -o protein_prod.out -r protein_prod.rst -x protein_prod.mdcrd. - Specificare le opzioni nel file di input prod.in:
&cntrl
imin=0, irest=1, ntx=5,
nstlim=5000000,
dt=0,002,
ntc=2, ntf=2,
cut=10.0, ntb=2, ntp=1, taup=2.0,
ntpr=1000, ntwx=1000, ntwr=50000,
ntt=3, gamma_ln=2.0,
temp0=300.0, ig=-1
/
NOTA: la versione parallela di pmemd (pmemd. MPI) o la versione accelerata da GPU (pmemd.cuda) può essere utilizzata su cluster di computer e supercomputer. Una lunga simulazione MD può essere suddivisa in diversi segmenti ed eseguita in sequenza. Si consiglia vivamente di impostare ig=-1 (opzione seed casuale) per ogni riavvio della simulazione. I programmi ptraj o cppptraj possono essere utilizzati per l'analisi e l'elaborazione di traiettorie di coordinate. Per ulteriori informazioni, consultare il Manuale dell'utente del software.
- Utilizzare prod.in e protein_equil.rst come input: pmemd -O -i prod.in -p protein.prmtop
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
Il presente protocollo può essere facilmente applicato in studi di modificazione post-traduzionale di caspasi o mutazioni patogene. In questa sezione viene illustrato il flusso di lavoro di modellazione MD (Figura 1), che è stato utilizzato con successo nello studio della caspasi-27. Utilizzando la mutagenesi sito-diretta in vitro di potenziali siti di fosforilazione (da Ser/Thr ad Ala) e approcci biochimici, è stato dimostrato che la mutazione Ser384Ala impedisce l'elaborazione della caspasi-2 e blocca l'attività enzimatica e l'induzione della morte cellulare apoptotica (Figura supplementare 1). Tuttavia, né l'analisi della spettrometria di massa né la tecnologia Phos-Tag hanno confermato che Ser384 subisce fosforilazione7. Per capire come Ser384 regola l'attività della caspasi-2, è stata eseguita la modellazione MD della struttura proteica tridimensionale (Figura 1).
I modelli molecolari delle forme wild-type e mutanti di caspasi-2 sono stati costruiti utilizzando la struttura cristallina 1pyo (le strutture disponibili della caspasi-2 sono elencate nella Tabella supplementare 1). Il mutante Ser384Ala è stato creato rimuovendo l'atomo di O γ nel residuo Ser384 (in PDB, i residui Ser e Ala si distinguono per avere l'atomo Oγ). Le coordinate dell'idrogeno sono solitamente assenti nelle strutture cristalline. Pertanto, gli atomi di idrogeno sono stati aggiunti alla struttura proteica, e poi è stato risolto da uno strato spesso 12 Å di acqua TIP3P per consentire ulteriori simulazioni esplicite di solventi. Sono stati aggiunti ioni di sodio per neutralizzare il sistema. I modelli di partenza ottenuti di caspasi-2 sono stati sottoposti a minimizzazione dell'energia, equilibrio e successiva simulazione MD di 10 ns secondo il protocollo di modellazione MD, utilizzando i file di dati di controllo min1.in, min2.in, heat.in, equil.in e prod.in.
Nella struttura cristallina della caspasi-2, Ser384, insieme ad Arg219 e Arg378, sono coinvolti nella formazione della superficie della cavità nel sito attivo. Arg219 e Arg378 formano legami idrogeno con il gruppo carbossilico del substrato, mentre Ser384 non interviene nelle interazioni dirette con il substrato. Le simulazioni MD hanno confermato che la sostituzione di Ser384Ala non ha influenzato i residui catalitici Cys320 (nucleofilo) e His277 (base generale) ma ha indotto un importante cambiamento conformazionale in Arg378. Il suo gruppo guanidinio si rivolse verso il solvente di massa in modo che l'atomo diε N non potesse formare un legame idrogeno essenziale con il gruppo carbossilico del substrato (Figura 1). Nell'enzima nativo, si verifica un'interazione elettrostatica tra l'atomoO γ di Ser384 (carica negativa parziale) e il gruppo guanidinio Arg378 (carica positiva). Tuttavia, la sostituzione della serina apparentemente interrompe questa interazione. Pertanto, è stato dimostrato che la sostituzione Ser384Ala ha influenzato il riconoscimento del substrato da parte dei residui di arginina nel sito attivo della caspasi-2, compromettendo l'attività enzimatica e la capacità di innescare la morte cellulare apoptotica. Il meccanismo scoperto sembra evolutivamente conservato e comune per altri membri della famiglia delle caspasi. Pertanto, l'implementazione di simulazioni MD ha permesso di confermare i risultati biochimici e ottenere nuove informazioni sulla struttura molecolare del centro attivo delle caspasi.
Figura 1: Flusso di lavoro di modellazione MD utilizzato per studiare la struttura delle caspasi. La caspasi-2 wild-type e il suo mutante Ser384Ala sono stati studiati seguendo le istruzioni nella sezione protocollo. È stato rivelato che la sostituzione di Ser384Ala induce un importante cambiamento conformazionale nel residuo del sito attivo Arg378. Fare clic qui per visualizzare una versione ingrandita di questa figura.
Figura supplementare 1: Funzione della caspasi-2 nella morte cellulare. In risposta a stimoli estrinseci e intrinseci, la caspasi-2 wild-type può essere attivata da un meccanismo autoproteolitico. La caspasi-2 attiva fende Bid, che a sua volta promuove la permeabilizzazione della membrana esterna mitocondriale e la morte cellulare apoptotica. La mutazione Ser384Ala impedisce l'attivazione della caspasi-2 e l'induzione della morte cellulare. Clicca qui per scaricare questo file.
Tabella supplementare 1: Strutture di Caspasi-2 disponibili nella Protein Data Bank. Le strutture sono ordinate per data di rilascio. Clicca qui per scaricare questa tabella.
File supplementare 1: diversi passaggi nella modifica di un file PDB. Clicca qui per scaricare questo file.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
L'approccio MD descritto consente di modellare sia le forme wild-type che mutanti di caspasi utilizzando i pacchetti di simulazione biomolecolare. Qui vengono discusse diverse questioni importanti della metodologia. In primo luogo, una struttura cristallina rappresentativa della caspasi deve essere selezionata dalla Protein Data Bank. È importante sottolineare che entrambe le forme monomeriche e dimeriche di caspasi sono accettabili. Scegliere strutture ad alta risoluzione con un numero minimo di residui mancanti è una buona idea. Lo stato di protonazione di alcuni residui può essere impostato manualmente nel file PDB di input. Ad esempio, nella Figura 1, un atomo di idrogeno è stato attaccato all'atomo Nε2 (forma HIE) del residuo catalitico His277. In secondo luogo, una scatola di solventi deve essere creata attorno alla proteina per ulteriori simulazioni esplicite del solvente. La minimizzazione dell'energia è necessaria per ottimizzare le coordinate degli atomi di idrogeno aggiunti e delle molecole d'acqua. In terzo luogo, il riscaldamento e l'equilibratura devono essere eseguiti prima della simulazione della dinamica di produzione. Nella fase di riscaldamento, si dovrebbe assicurarsi che le cavità e le tasche leganti sulla superficie proteica si riempiano di molecole d'acqua (se necessario, è possibile aumentare il numero di passaggi MD). Nella fase di equilibrio, l'acquisizione della conformazione all'equilibrio della proteina deve essere confermata analizzando la deviazione quadratica media della radice degli atomi della spina dorsale dalle posizioni iniziali12. In quarto luogo, è necessario monitorare i cambiamenti conformazionali della struttura proteica durante la simulazione della dinamica di produzione. A tale scopo, si dovrebbero sovrapporre i frame di traiettoria alla struttura di partenza adattando gli atomi della spina dorsale e quindi analizzare i cambiamenti conformazionali utilizzando uno strumento di visualizzazione molecolare (VMD, PyMOL, ecc.) 17. Se necessario, è possibile aumentare il numero di passaggi MD e / o suddividere la simulazione in diversi segmenti eseguiti in sequenza. In quinto luogo, si consiglia di utilizzare la versione accelerata da GPU del programma di simulazione MD (pmemd.cuda), le cui prestazioni superano significativamente quelle ottenibili dalla tradizionale implementazione della CPU14,18.
Una possibile limitazione del metodo descritto è la disponibilità di alcune strutture caspasiche nella Protein Data Bank. Tuttavia, più di 900 voci PDB sono state trovate sotto la richiesta "caspasi", indicando che le strutture delle caspasi sono state studiate in dettaglio. Alcune difficoltà possono sorgere quando una struttura cristallina è ottenuta a bassa risoluzione (cioè ad un livello di dettaglio grossolano) o manca di alcuni frammenti importanti per la funzione caspasi.
In sintesi, gli approcci di modellazione MD si sono dimostrati efficaci nel prevedere i cambiamenti strutturali delle proteine dopo mutazioni dei loro geni, modifiche o interazioni intermolecolari. L'applicazione della simulazione MD nel campo della morte cellulare apre grandi opportunità per studiare i meccanismi molecolari di regolazione delle caspasi mediante modificazioni post-traduzionali. A questo proposito, la simulazione MD a lungo termine potrebbe valutare il comportamento dinamico delle regioni funzionali e il potenziale effetto delle sostituzioni aminoacidiche per chiarire le cause molecolari associate alla caspasi-2 mutata o modificata e ad altre caspasi. Ad esempio, mutazioni missenso dei geni iniziatori della caspasi (caspasi-2/caspasi-8/caspasi-9/caspasi-10) sono state rilevate nei tumori gastrointestinali e nei tumori del sistema nervoso centrale e periferico19. Insieme, la modellazione MD delle caspasi è un potente approccio in silico per comprendere i meccanismi che regolano la morte cellulare programmata e i disturbi associati e per lo sviluppo di nuove terapie efficaci.
Subscription Required. Please recommend JoVE to your librarian.
Disclosures
Gli autori non hanno conflitti di interesse da rivelare.
Acknowledgments
Questo lavoro è stato sostenuto da una sovvenzione della Russian Science Foundation (17-75-20102, lo sviluppo del protocollo). Gli esperimenti descritti nella sezione dei risultati rappresentativi (analisi della fosforilazione) sono stati sostenuti dalle società oncologiche di Stoccolma (181301) e svedese (190345).
Materials
| Name | Company | Catalog Number | Comments |
| Amber20 | University of California, San Francisco | Software for molecular dynamics simulation http://ambermd.org |
|
| AmberTools21 | University of California, San Francisco | Software for molecular modeling and analysis http://ambermd.org |
References
- Olsson, M., Zhivotovsky, B.
- Lavrik, I. N., Golks, A., Krammer, P. H. Caspase: Pharmacological manipulation of cell death. Journal of Clinical Investigation. 115 (10), 2665-2672 (2005).
- Degterev, A., Boyce, M., Yuan, J.
- Pop, C., Salvesen, G. S. Human caspases: Activation, specificity, and regulation. Journal of Biological Chemistry. 284 (33), 21777-21781 (2009).
- Zamaraev, A. V., Kopeina, G. S., Prokhorova, E. A., Zhivotovsky, B., Lavrik, I. N. Post-translational modification of caspases: The other side of apoptosis regulation. Trends in Cell Biology. 27 (5), 322-339 (2017).
- Pearlman, S. M., Serber, Z., Ferrell, J. E. A mechanism for the evolution of phosphorylation sites. Cell. 147 (4), 934-946 (2011).
- Zamaraev, A. V., et al. Requirement for Serine-384 in Caspase-2 processing and activity. Cell Death and Disease. 11 (10), 825 (2020).
- Dagbay, K. B., Hill, M. E., Barrett, E., Hardy, J. A. Tumor-associated mutations in caspase-6 negatively impact catalytic efficiency. Biochemistry. 56 (34), 4568-4577 (2017).
- Ando, M., et al. Cancer-associated missense mutations of caspase-8 activate nuclear factor-κB signaling. Cancer Science. 104 (8), 1002-1008 (2013).
- Case, D. A., et al. AMBER 2020. , University of California, San Francisco. (2020).
- Salomon-Ferrer, R., Case, D. A., Walker, R. C. An overview of the Amber biomolecular simulation package. WIREs Computational Molecular Science. 3 (2), 198-210 (2013).
- Roe, D. R., Cheatham, T. E. PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data. Journal of Chemical Theory and Computation. 9 (7), 3084-3095 (2013).
- Maier, J. A., et al. ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. Journal of Chemical Theory and Computation. 11 (8), 3696-3713 (2015).
- Salomon-Ferrer, R., Götz, A. W., Poole, D., Le Grand, S., Walker, R. C. Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh ewald. Journal of Chemical Theory and Computation. 9 (9), 3878-3888 (2013).
- Berman, H. M., et al.
- Burley, S. K., et al. RCSB Protein Data Bank: Powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences. Nucleic Acids Research. 49, 437-451 (2021).
- Martinez, X., et al. Molecular graphics: Bridging structural biologists and computer scientists. Structure. 27 (11), 1617-1623 (2019).
- Lee, T. S., et al. GPU-accelerated molecular dynamics and free energy methods in Amber18: Performance enhancements and new features. Journal of Chemical Information and Modeling. 58 (10), 2043-2050 (2018).
- Jäger, R., Zwacka, R. M. The enigmatic roles of caspases in tumor development. Cancers. 2 (4), 1952-1979 (2010).
