Name: Explorando mutações da caspase e modificação pós-translacional por abordagens de modelagem molecular (Vídeo) | JoVE Journal
Uploaded: 2022-10-13T14:00:00.000+00:00
Duration: 5 min 56 s
Description: Explore simulações de dinâmica molecular usando Amber 20 para investigar mutações na caspase e modificações pós-traducionais. Este método auxilia no desenvolvimento de terapias direcionadas contra câncer.

ai_agent_IP

ai_agent

maintenance

instantTrialDurationMinutes

isBiopharma

enableNewMyJoveSectionDesign

educationFeatures

seoMetaApi

ai_copilot_institution_blacklist

ai_chat_institution_blacklist

freeAccessWindow

ssoRoleCapture

scheduledMaintenanceBanner

joveLabsFeatures

disable_autoplay_for_bots

extended_trial_duration

5 minutes

Este trabalho foi apoiado por uma bolsa da Fundação Russa de Ciência (17-75-20102, o desenvolvimento do protocolo). Os experimentos descritos na seção de resultados representativos (análise da fosforilação) foram apoiados pelas Sociedades de Câncer de Estocolmo (181301) e Suécia (190345).

1. Preparação do sistema
NOTA: Os modelos moleculares das formas proteicas nativas e mutantes são construídos com base em uma estrutura cristalina apropriada obtida do Protein Data Bank15,16.
<ol><li>Para recuperar a estrutura PDB selecionada, use a lista suspensa Baixar arquivos e clique em Formato PDB. Remova comentários e dados de conectividade e insira um cartão TER entre cadeias de proteínas separadas no arquivo PDB. Opcionalmente, defina o estado de protonação dos resíduos de histidina substituindo o nome do resíduo HIS por HIE, HID ou HIP (Arquivo Suplementar 1). NOTA: É razoável não remover moléculas de água cristalograficamente resolvidas (se presentes no arquivo PDB).</li>
	<li>Para preparar o modelo inicial, inicie o programa tleap (pacote AmberTools) e insira comandos que manipulam objetos tleap .<ol><li>Inicie o programa: tleap. Carregue o campo de força ff14SB para descrever a proteína com mecânica molecular: > fonte leaprc.protein.ff14SB.</li>
		<li>Campos de força de carga para moléculas de água e íons atômicos (Na+, Cl-): fonte de > leaprc.water.tip3p.</li>
		<li>Carregue o arquivo PDB e construa coordenadas para hidrogênios, criando um objeto chamado mol: > mol = loadpdb protein.pdb.</li>
		<li>Verifique se há inconsistências internas que possam causar problemas: > verifique mol. NOTA: Pontes de dissulfeto, se presentes, devem ser especificadas manualmente. Substitua os nomes das cisteínas ligadas covalentemente CYS por CYX e digite o seguinte comando em tleap para criar uma ligação entre átomos SG: ligação mol.X.SG mol.Y.SG (X e Y são números de resíduos de cisteína).</li>
		<li>Crie uma caixa de solvente em torno da proteína (água TIP3P, distância de 12 Å): > solvatebox mol TIP3PBOX 12.0.</li>
		<li>Verifique a carga total: > mol de carga e adicione contraíons (Na+ ou Cl-) para neutralizar o sistema: > adições mol Na+ 0.</li>
		<li>Crie o arquivo prmtop de topologia e o arquivo inpcrd de coordenadas (entradas para executar uma simulação): > saveamberparm mol protein.prmtop protein.inpcrd. Saia do programa tleap: > saia. NOTA: Os arquivos de coordenadas no formato âmbar podem ser facilmente convertidos para o formato PDB usando a ferramenta ambpdb (consulte o Manual do Usuário, consulte Tabela de Materiais).</li>
	</ol></li>
</ol>2. Minimização de energia
NOTA: A minimização de energia é necessária para remover quaisquer contatos defeituosos e sobreposições entre átomos no sistema de partida que levam à instabilidade ao executar o MD.
<ol><li>Realizar o primeiro estágio da minimização de energia (2.500 passos do algoritmo de descida mais íngreme + 2.500 passos do algoritmo de gradiente conjugado) para otimizar as posições de átomos de hidrogênio e moléculas de água adicionados, mantendo as coordenadas de proteína fixadas por restrições posicionais em átomos pesados.<ol><li>Execute o programa pmemd da seguinte maneira: pmemd -O -i min1.in -p protein.prmtop -c protein.inpcrd -o protein_min1.out -r protein_min1.rst -ref protein.inpcrd.</li>
		<li>Siga os argumentos necessários: -i FILE dados de controle; -p Topologia molecular FILE, parâmetros de campo de força, nomes de átomos; -c Coordenadas iniciais do FILE; -o FILE de saída de log legível pelo usuário; -r Coordenadas finais do FILE; -ref Coordenadas de referência FILE para restrições de posição.</li>
		<li>Especifique as opções no arquivo de entrada 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: Opções de entrada: imin=1 executar minimização; maxcyc=5000 número máximo de ciclos de minimização; O método de minimização ncyc=2500 é alternado da descida mais íngreme para o gradiente conjugado após ciclos ncyc ; corte=10,0 ponto de corte não colado (Å); ntb=1 são impostos limites periódicos, volume constante; ntc=1 nenhuma restrição é aplicada aos comprimentos de ligação (para melhor convergência de energia); ntf=1 todas as interações são calculadas; ntpr=10 informações de energia são impressas no arquivo de saída a cada passo ntpr ; ntr=1 sinalizador para conter átomos especificados usando um potencial harmônico; restraintmask=':1-517 & !@H=' especifica átomos contidos; restraint_wt=2,0 peso (kcal/mol∙Å2) para restrições posicionais.</li>
	</ol></li>
	<li>Realize o segundo estágio da minimização de energia (5.000 degraus de descida mais íngremes + 5.000 degraus de gradiente conjugado) sem restrições para otimizar todo o sistema.<ol><li>Use min2.in e protein_min1.rst como entradas: pmemd -O -i min2.in -p protein.prmtop -c protein_min1.rst -o protein_min2.out -r protein_min2.rst.</li>
		<li>Especifique as opções no arquivo de entrada min2.in: &cntrl imin=1, maxcyc=10000, ncyc=5000, cut=10.0, ntb=1, ntc=1, ntf=1, ntpr=10 /</li>
	</ol></li>
</ol>3. Aquecimento
NOTA: Esta etapa visa aquecer o sistema de 0 K a 300 K. As velocidades iniciais são atribuídas aos átomos, uma vez que o modelo inicial baseado no arquivo PDB não contém informações de velocidade.
<ol><li>Realizar o processo de aquecimento com restrições posicionais nos átomos de proteína (50 ps, volume constante).<ol><li>Use heat.in e protein_min2.rst como entradas: 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</li>
		<li>Siga os argumentos necessários: -x conjuntos de coordenadas FILE salvos ao longo da trajetória do MD.</li>
		<li>Especifique as opções no arquivo de entrada 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, máscara de retenção=':1-517', restraint_wt=1,0, nmropt=1 / &wt TYPE='TEMP0', istep1=0, istep2=25000, valor1=0,1, valor2=300,0 / &wt TYPE='FIM' / Observação : opções de entrada: imin = 0 executar MD; irest=0 execute uma nova simulação; ntx=1 nenhuma velocidade inicial é lida a partir do primeiro arquivo; nstlim=25000 número de passos de MD; dt=0,002 passo de tempo (ps); As ligações ntc=2 envolvendo hidrogênio são restringidas usando o algoritmo SHAKE; ntf=2 forças para as ligações restritas não são calculadas; ntwx=500 coordenadas são gravadas no arquivo mdcrd a cada passo ntwx ; ntt=3 Regulação da temperatura de Langevin; gamma_ln=2,0 frequência de colisão para a dinâmica de Langevin (ps−1); tempi=0,0 temperatura inicial; temperatura de referência temp0=300,0 ; nmropt=1 parâmetros especificados em uma lista de nomes &wt são lidos; TYPE='TEMP0' varia a temperatura alvo.</li>
	</ol></li>
</ol>4. Equilíbrio
NOTA: Esta etapa é necessária para ajustar a densidade da água e obter o estado de equilíbrio da proteína.
<ol><li>Realizar o equilíbrio a 300 K sem quaisquer restrições (500 ps, pressão constante).<ol><li>Use equil.in e protein_heat.rst como entradas: pmemd -O -i equil.in -p protein.prmtop -c protein_heat.rst -o protein_equil.out -r protein_equil.rst -x protein_equil.mdcrd.</li>
		<li>Especifique as opções no arquivo de entrada 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: Opções de entrada: irest=1 reinicie a simulação; As coordenadas e velocidades ntx=5 são lidas a partir de um arquivo rst gerado anteriormente; ntb=2 limites periódicos são impostos, pressão constante; ntp=1 escala de pressão isotrópica; taup=2,0 tempo de relaxamento da pressão (ps); ntwr=50000 a cada passo ntwr , o primeiro arquivo é gravado, garantindo a recuperação de uma falha.</li>
	</ol></li>
</ol>5. Dinâmica de produção
<ol><li>Após o equilíbrio ter sido alcançado com sucesso, realize uma simulação de MD de produção (10 ns ou mais) a pressão constante e gere o arquivo de trajetória para análise subsequente da estrutura da proteína.<ol><li>Use prod.in e protein_equil.rst como entradas: pmemd -O -i prod.in -p protein.prmtop -c protein_equil.rst -o protein_prod.out -r protein_prod.rst -x protein_prod.mdcrd.</li>
		<li>Especifique as opções no arquivo de entrada 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: A versão paralela do pmemd (pmemd. MPI) ou versão acelerada por GPU (pmemd.cuda) pode ser usada em clusters de computadores e supercomputadores. Uma longa simulação de MD pode ser dividida em vários segmentos e realizada sequencialmente. É altamente recomendável definir ig=-1 (opção de semente aleatória) para cada reinicialização da simulação. Os programas ptraj ou cppptraj podem ser usados para a análise e processamento de trajetórias de coordenadas. Para obter mais informações, consulte o Manual do Usuário do software.</li>
	</ol></li>
</ol>

<ol><li>Olsson, M., Zhivotovsky, B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Caspases+and+cancer%5BTitle%5D)+AND+%22Cell+Death+and+Differentiation%22%5BJournal%5D)">Caspases and cancer.</a> Cell Death and Differentiation. 18 (9), 1441-1449 (2011).</li>
<li>Lavrik, I. N., Golks, A., Krammer, P. H. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Caspase%3A+Pharmacological+manipulation+of+cell+death%5BTitle%5D)+AND+%22Journal+of+Clinical+Investigation%22%5BJournal%5D)">Caspase: Pharmacological manipulation of cell death.</a> Journal of Clinical Investigation. 115 (10), 2665-2672 (2005).</li>
<li>Degterev, A., Boyce, M., Yuan, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+decade+of+caspases%5BTitle%5D)+AND+%22Oncogene%22%5BJournal%5D)">A decade of caspases.</a> Oncogene. 22 (53), 8543-8567 (2003).</li>
<li>Pop, C., Salvesen, G. S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Human+caspases%3A+Activation%2C+specificity%2C+and+regulation%5BTitle%5D)+AND+%22Journal+of+Biological+Chemistry%22%5BJournal%5D)">Human caspases: Activation, specificity, and regulation.</a> Journal of Biological Chemistry. 284 (33), 21777-21781 (2009).</li>
<li>Zamaraev, A. V., Kopeina, G. S., Prokhorova, E. A., Zhivotovsky, B., Lavrik, I. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Post-translational+modification+of+caspases%3A+The+other+side+of+apoptosis+regulation%5BTitle%5D)+AND+%22Trends+in+Cell+Biology%22%5BJournal%5D)">Post-translational modification of caspases: The other side of apoptosis regulation.</a> Trends in Cell Biology. 27 (5), 322-339 (2017).</li>
<li>Pearlman, S. M., Serber, Z., Ferrell, J. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+mechanism+for+the+evolution+of+phosphorylation+sites%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">A mechanism for the evolution of phosphorylation sites.</a> Cell. 147 (4), 934-946 (2011).</li>
<li>Zamaraev, A. V., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Requirement+for+Serine-384+in+Caspase-2+processing+and+activity%5BTitle%5D)+AND+%22Cell+Death+and+Disease%22%5BJournal%5D)">Requirement for Serine-384 in Caspase-2 processing and activity.</a> Cell Death and Disease. 11 (10), 825(2020).</li>
<li>Dagbay, K. B., Hill, M. E., Barrett, E., Hardy, J. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Tumor-associated+mutations+in+caspase-6+negatively+impact+catalytic+efficiency%5BTitle%5D)+AND+%22Biochemistry%22%5BJournal%5D)">Tumor-associated mutations in caspase-6 negatively impact catalytic efficiency.</a> Biochemistry. 56 (34), 4568-4577 (2017).</li>
<li>Ando, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Cancer-associated+missense+mutations+of+caspase-8+activate+nuclear+factor-%CE%BAB+signaling%5BTitle%5D)+AND+%22Cancer+Science%22%5BJournal%5D)">Cancer-associated missense mutations of caspase-8 activate nuclear factor-κB signaling.</a> Cancer Science. 104 (8), 1002-1008 (2013).</li>
<li>Case, D. A., et al. <a target="_blank" href="http://scholar.google.com/scholar?hl=en&safe=off&q=author%3aCase%2C+DA+%22AMBER+2020%22">AMBER 2020</a>. , University of California, San Francisco. (2020).</li>
<li>Salomon-Ferrer, R., Case, D. A., Walker, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((An+overview+of+the+Amber+biomolecular+simulation+package%5BTitle%5D)+AND+%22WIREs+Computational+Molecular+Science%22%5BJournal%5D)">An overview of the Amber biomolecular simulation package.</a> WIREs Computational Molecular Science. 3 (2), 198-210 (2013).</li>
<li>Roe, D. R., Cheatham, T. E. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((PTRAJ+and+CPPTRAJ%3A+Software+for+processing+and+analysis+of+molecular+dynamics+trajectory+data%5BTitle%5D)+AND+%22Journal+of+Chemical+Theory+and+Computation%22%5BJournal%5D)">PTRAJ and CPPTRAJ: Software for processing and analysis of molecular dynamics trajectory data.</a> Journal of Chemical Theory and Computation. 9 (7), 3084-3095 (2013).</li>
<li>Maier, J. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((ff14SB%3A+Improving+the+accuracy+of+protein+side+chain+and+backbone+parameters+from+ff99SB%5BTitle%5D)+AND+%22Journal+of+Chemical+Theory+and+Computation%22%5BJournal%5D)">ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB.</a> Journal of Chemical Theory and Computation. 11 (8), 3696-3713 (2015).</li>
<li>Salomon-Ferrer, R., Götz, A. W., Poole, D., Le Grand, S., Walker, R. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Routine+microsecond+molecular+dynamics+simulations+with+AMBER+on+GPUs.+2.+Explicit+solvent+particle+mesh+ewald%5BTitle%5D)+AND+%22Journal+of+Chemical+Theory+and+Computation%22%5BJournal%5D)">Routine microsecond molecular dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh ewald.</a> Journal of Chemical Theory and Computation. 9 (9), 3878-3888 (2013).</li>
<li>Berman, H. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+Protein+Data+Bank%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">The Protein Data Bank.</a> Nucleic Acids Research. 28 (1), 235-242 (2000).</li>
<li>Burley, S. K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((RCSB+Protein+Data+Bank%3A+Powerful+new+tools+for+exploring+3D+structures+of+biological+macromolecules+for+basic+and+applied+research+and+education+in+fundamental+biology%2C+biomedicine%2C+biotechnology%2C+bioengineering+and+energy+sciences%5BTitle%5D)+AND+%22Nucleic+Acids+Research%22%5BJournal%5D)">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.</a> Nucleic Acids Research. 49, 437-451 (2021).</li>
<li>Martinez, X., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Molecular+graphics%3A+Bridging+structural+biologists+and+computer+scientists%5BTitle%5D)+AND+%22Structure%22%5BJournal%5D)">Molecular graphics: Bridging structural biologists and computer scientists.</a> Structure. 27 (11), 1617-1623 (2019).</li>
<li>Lee, T. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((GPU-accelerated+molecular+dynamics+and+free+energy+methods+in+Amber18%3A+Performance+enhancements+and+new+features%5BTitle%5D)+AND+%22Journal+of+Chemical+Information+and+Modeling%22%5BJournal%5D)">GPU-accelerated molecular dynamics and free energy methods in Amber18: Performance enhancements and new features.</a> Journal of Chemical Information and Modeling. 58 (10), 2043-2050 (2018).</li>
<li>Jäger, R., Zwacka, R. M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((The+enigmatic+roles+of+caspases+in+tumor+development%5BTitle%5D)+AND+%22Cancers%22%5BJournal%5D)">The enigmatic roles of caspases in tumor development.</a> Cancers. 2 (4), 1952-1979 (2010).</li>
</ol>

Os autores não têm conflitos de interesse a divulgar.

A abordagem MD descrita permite modelar as formas selvagens e mutantes de caspase usando os pacotes de simulação biomolecular. Várias questões importantes da metodologia são discutidas aqui. Primeiro, uma estrutura cristalina representativa da caspase precisa ser selecionada do Banco de Dados de Proteínas. É importante ressaltar que as formas monomérica e dimérica de caspase são aceitáveis. Escolher estruturas de alta resolução com um número mínimo de resíduos ausentes é uma boa ideia. O estado de protonação de alguns resíduos pode ser definido manualmente no arquivo PDB de entrada. Por exemplo, na Figura 1, um átomo de hidrogênio foi anexado ao átomo de Nε2 (forma HIE) do resíduo catalítico His277. Em segundo lugar, uma caixa de solvente deve ser criada em torno da proteína para simulações de solvente explícito adicionais. A minimização de energia é necessária para otimizar as coordenadas dos átomos de hidrogênio adicionados e moléculas de água. Em terceiro lugar, o aquecimento e o equilíbrio devem ser realizados antes da simulação da dinâmica de produção. Na fase de aquecimento, deve-se garantir que as cavidades e bolsas de ligação na superfície da proteína se encham de moléculas de água (se necessário, pode-se aumentar o número de etapas de MD). Na fase de equilíbrio, a aquisição da conformação de equilíbrio da proteína precisa ser confirmada pela análise do desvio raiz-média-quadrado de seus átomos da espinha dorsal em relação às posições iniciais12. Em quarto lugar, é necessário monitorar as mudanças conformacionais da estrutura da proteína durante a simulação da dinâmica de produção. Para este propósito, deve-se sobrepor quadros de trajetória na estrutura inicial, encaixando os átomos da espinha dorsal e, em seguida, analisar as mudanças conformacionais usando uma ferramenta de visualização molecular (VMD, PyMOL, etc.) 17. Se necessário, pode-se aumentar o número de etapas de DM e/ou dividir a simulação em vários segmentos realizados sequencialmente. Em quinto lugar, recomenda-se a utilização da versão acelerada por GPU do programa de simulação MD (pmemd.cuda), cujo desempenho excede significativamente o alcançável pela implementação tradicional da CPU14,18.
Uma possível limitação do método descrito é a disponibilidade de certas estruturas de caspase no Banco de Dados de Proteínas. No entanto, mais de 900 entradas do APO são encontradas no âmbito do pedido "caspase", indicando que as estruturas da caspase foram investigadas em pormenor. Algumas dificuldades podem surgir quando uma estrutura cristalina é obtida em baixa resolução (ou seja, em um nível grosseiro de detalhe) ou carece de alguns fragmentos importantes para a função caspase.
Em resumo, as abordagens de modelagem MD provaram ser eficazes na previsão de mudanças estruturais de proteínas após mutações de seus genes, modificações ou interações intermoleculares. A aplicação da simulação de DM no campo da morte celular abre grandes oportunidades para estudar os mecanismos moleculares de regulação da caspase por modificações pós-traducionais. A este respeito, a simulação de DM a longo prazo poderia avaliar o comportamento dinâmico de regiões funcionais e o efeito potencial das substituições de aminoácidos para elucidar as causas moleculares associadas à caspase-2 mutada ou modificada, bem como a outras caspases. Por exemplo, mutações missense de genes iniciadores da caspase (caspase-2/caspase-8/caspase-9/caspase-10) foram detectadas em cânceres gastrointestinais e em tumores dos sistemas nervoso central e periférico19. Juntos, a modelagem MD de caspases é uma poderosa abordagem in silico para entender os mecanismos que regulam a morte celular programada e distúrbios associados e para o desenvolvimento de novas terapias eficazes.

A apoptose é um dos processos celulares mais amplamente estudados que regulam a morfogênese e a homeostase tecidual de organismos multicelulares. A apoptose pode ser iniciada por uma ampla gama de estímulos externos ou internos, como a ativação de receptores de morte, distúrbios nos sinais do ciclo celular, danos ao DNA, estresse do retículo endoplasmático (RE) e várias infecções bacterianas e virais1. As caspases - principais atores apoptóticos - são convencionalmente classificadas em dois grupos: iniciadores (caspase-2, caspase-8, caspase-9 e caspase-10) e efetores (caspase-3, caspase-6 e caspase-7), dependendo de sua estrutura de domínio e do lugar na cascata da caspase 2,3. Após os sinais de morte celular, as caspases iniciadoras interagem com moléculas adaptadoras que facilitam a dimerização induzida pela proximidade e o autoprocessamento para formar uma enzima ativa. As caspases efetoras são ativadas através de clivagem por caspases iniciadoras e realizam etapas de execução a jusante por clivagem múltiplos substratos celulares4.
A maturação e a função das caspases iniciadora e efetora são reguladas por um grande número de diferentes mecanismos intracelulares, dentre os quais a modificação pós-translacional desempenha um papel indispensável na modulação da morte celular5. A adição de grupos modificadores (fosforilação, nitrosilação, metilação ou acetilação) ou proteínas (ubiquitinação ou SUMOilação) altera a atividade enzimática das caspases ou a conformação e estabilidade proteica que regulam a apoptose. A mutagênese dirigida ao site é amplamente aplicada para investigar os possíveis locais de modificação pós-translacional e discernir seu papel. Um suposto local de modificação é geralmente substituído por outro aminoácido, que não pode ser modificado ainda mais. Assim, serina e treonina potencialmente fosforiladas são mutadas para alanina, e os locais de ubiquitinação de lisina são substituídos por arginina. Outra estratégia inclui a substituição de um aminoácido que imita particularmente a modificação pós-translacional (por exemplo, glutamato e aspartato têm sido usados para imitar serina fosforilada ou treonina)6. No entanto, algumas dessas substituições localizadas na alta vizinhança de um sítio ativo ou em posições críticas podem alterar a estrutura da caspase, perturbar a atividade catalítica e suprimir a morte celular apoptótica7. Efeitos semelhantes podem ser observados em casos de mutações missense associadas a tumores em genes caspase. Por exemplo, a mutação tumoral associada à caspase-6 - R259H - resultou em alterações conformacionais de alças na bolsa de ligação ao substrato, reduzindo o eficiente turnover catalítico dos substratos8. A substituição do aminoácido G325A na caspase-8 identificada no carcinoma espinocelular de cabeça e pescoço poderia dificultar a atividade da caspase-8, o que levou à modulação da sinalização do fator nuclear kB (NF-kB) e promoveu a tumorigênese9.
Para avaliar o efeito potencial das substituições de aminoácidos na estrutura e função da caspase, a modelagem molecular pode ser aplicada. A abordagem da dinâmica molecular (MD) é descrita neste trabalho para modelar a caspase do tipo selvagem e suas formas mutantes usando o pacote de simulação biomolecular (Amber). O método MD fornece uma visão da evolução dinâmica da estrutura da proteína após a introdução de mutações. Originalmente desenvolvido pelo grupo de Peter Kollman, o pacote Amber tornou-se uma das ferramentas de software mais populares para simulações biomoleculares10,11,12,13. Este software é dividido em duas partes: (1) AmberTools, uma coleção de programas rotineiramente utilizados para a preparação do sistema (atribuição de tipo de átomo, adição de hidrogênios e moléculas de água explícita, etc.) e análise de trajetória; e (2) Âmbar, que é centrado em torno do programa de simulação pmemd. O AmberTools é um pacote gratuito (e um pré-requisito para instalar o próprio Amber), enquanto o Amber é distribuído com uma estrutura separada de licenças e taxas. Simulações paralelas em um supercomputador e/ou utilizando unidades de processamento gráfico (GPUs) podem melhorar substancialmente o desempenho para a pesquisa científica da dinâmica da estrutura de proteínas14. As versões de software mais recentes disponíveis são AmberTools21 e Amber20, mas os protocolos descritos também podem ser usados com as versões anteriores.

<table><tbody><tr><td>Âmbar 20</td><td>Universidade da Califórnia, São Francisco</td><td></td><td>Software para simulação de dinâmica molecular http://ambermd.org</td></tr><tr><td>Ferramentas Âmbar 21</td><td>Universidade da Califórnia, São Francisco</td><td></td><td>Software para modelagem e análise molecular http://ambermd.org</td></tr></tbody></table>

explorando mutações da caspase e modificação pós-translacional por abordagens de modelagem molecular

Exploring Caspase Mutations and Post-Translational Modification by Molecular Modeling Approaches

A apoptose é um tipo de morte celular programada que elimina células danificadas e controla o desenvolvimento e a homeostase tecidual de organismos multicelulares. Caspases, uma família de proteases de cisteína, desempenham um papel fundamental na iniciação e execução da apoptose. A maturação das caspases e sua atividade é ajustada por modificações pós-traducionais de forma altamente dinâmica. Para avaliar o efeito das alterações pós-translacionais, os locais potenciais são rotineiramente mutados com resíduos persistentes a quaisquer modificações. Por exemplo, o resíduo de serina é substituído por alanina ou ácido aspártico. No entanto, tais substituições poderiam alterar a conformação do sítio ativo da caspase, levando a distúrbios na atividade catalítica e nas funções celulares. Além disso, mutações de outros resíduos de aminoácidos localizados em posições críticas também podem quebrar a estrutura e as funções das caspases e levar à perturbação da apoptose. Para evitar as dificuldades de empregar resíduos mutantes, abordagens de modelagem molecular podem ser prontamente aplicadas para estimar o efeito potencial das substituições de aminoácidos na estrutura da caspase. O presente protocolo permite a modelagem tanto da caspase do tipo selvagem quanto de suas formas mutantes com o pacote de simulação biomolecular (Amber) e instalações de supercomputador para testar o efeito de mutações na estrutura e função da proteína.

O presente protocolo pode ser prontamente aplicado em estudos de modificação pós-translacional de caspases ou mutações patogênicas. Nesta seção, é ilustrado o fluxo de trabalho de modelagem de DM (Figura 1), que tem sido utilizado com sucesso no estudo da caspase-27. Usando mutagênese in vitro direcionada ao local de potenciais sítios de fosforilação (Ser/Thr a Ala) e abordagens bioquímicas, demonstrou-se que a mutação Ser384Ala impediu o processamento da caspase-2 e bloqueou a atividade enzimática e a indução da morte celular apoptótica (Figura 1 Suplementar). No entanto, nem a análise de espectrometria de massa nem a tecnologia Phos-Tag confirmaram que o Ser384 sofre fosforilação7. Para compreender como o Ser384 regula a atividade da caspase-2, foi realizada a modelagem MD da estrutura proteica tridimensional (Figura 1).
Os modelos moleculares das formas selvagens e mutantes da caspase-2 foram construídos utilizando a estrutura cristalina 1pyo (as estruturas disponíveis da caspase-2 estão listadas na Tabela Suplementar 1). O mutante Ser384Ala foi criado removendo o átomo de O γ no resíduo Ser384 (em PDB, os resíduos Ser e Ala distinguem-se por ter o átomo Oγ). As coordenadas de hidrogênio geralmente estão ausentes nas estruturas cristalinas. Portanto, átomos de hidrogênio foram adicionados à estrutura da proteína e, em seguida, ela foi solvida por uma camada de 12 Å de espessura de água TIP3P para permitir simulações de solventes explícitos. Íons de sódio foram adicionados para neutralizar o sistema. Os modelos iniciais obtidos de caspase-2 foram submetidos à minimização de energia, equilíbrio e subsequente simulação de MD de 10 ns de acordo com o protocolo de modelagem MD, utilizando os arquivos de dados de min1.in, min2.in, heat.in, equil.in e controle de prod.in.
Na estrutura cristalina da caspase-2, Ser384, juntamente com Arg219 e Arg378, estão envolvidos na formação da superfície da cavidade no sítio ativo. Arg219 e Arg378 formam ligações de hidrogênio com o grupo carboxila do substrato, enquanto Ser384 não intervém em interações diretas com o substrato. Simulações de DM confirmaram que a substituição de Ser384Ala não afetou os resíduos catalíticos Cys320 (nucleófilo) e His277 (base geral), mas induziu uma grande alteração conformacional em Arg378. Seu grupo guanidínio voltou-se para o solvente a granel de modo que o átomode N ε não poderia formar uma ligação essencial de hidrogênio com o grupo carboxila do substrato (Figura 1). Na enzima nativa, ocorre uma interação eletrostática entre o átomo Oγ de Ser384 (carga parcial negativa) e o grupo guanidínio Arg378 (carga positiva). No entanto, a substituição de serina aparentemente interrompe essa interação. Assim, demonstrou-se que a substituição do Ser384Ala afetou o reconhecimento do substrato pelos resíduos de arginina no sítio ativo da caspase-2, prejudicando a atividade enzimática e a capacidade de desencadear a morte celular apoptótica. O mecanismo descoberto parece evolutivamente conservado e comum para outros membros da família caspase. Assim, a implementação de simulações de DM permitiu confirmar os resultados bioquímicos e obter novos insights sobre a estrutura molecular do centro ativo das caspases.
<img alt="Figure 1" class="xfigimg" src="/files/ftp_upload/64206/64206fig1v2.jpg"> Figura 1: Fluxo de trabalho de modelagem MD usado para estudar a estrutura da caspase. A caspase-2 do tipo selvagem e seu mutante Ser384Ala foram investigados seguindo as instruções na seção de protocolo. Foi revelado que a substituição do Ser384Ala induz uma importante alteração conformacional no resíduo do sítio ativo Arg378. <a href="https://www.jove.com/files/ftp_upload/64206/64206fig01large.jpg" target="_blank">Por favor, clique aqui para ver uma versão maior desta figura.</a>
Figura 1 Suplementar: Função da caspase-2 na morte celular. Em resposta a estímulos extrínsecos e intrínsecos, a caspase-2 do tipo selvagem pode ser ativada por um mecanismo autoproteolítico. A caspase-2 ativa cliva Bid, que, por sua vez, promove a permeabilização da membrana externa mitocondrial e a morte celular apoptótica. A mutação Ser384Ala impede a ativação da caspase-2 e a indução da morte celular. <a href="https://www.jove.com/files/ftp_upload/64206/Supplementary%20Figure%201%20fin.jpg" target="_blank">Clique aqui para baixar este arquivo.</a>
Tabela Suplementar 1: Estruturas de Caspase-2 disponíveis no Banco de Dados de Proteínas. As estruturas são classificadas por data de lançamento. <a href="https://www.jove.com/files/ftp_upload/64206/Supplementary%20Table%201.xlsx" target="_blank">Por favor, clique aqui para baixar esta Tabela.</a>
Arquivo Suplementar 1: Diferentes etapas na edição de um arquivo PDB. <a href="https://www.jove.com/files/ftp_upload/64206/Supplementary%20File%201.docx" target="_blank">Clique aqui para baixar este arquivo.</a>

O presente protocolo utiliza um pacote de simulação biomolecular e descreve a abordagem de dinâmica molecular (MD) para modelar a caspase do tipo selvagem e suas formas mutantes. O método MD permite avaliar a evolução dinâmica da estrutura da caspase e o efeito potencial de mutações ou modificações pós-traducionais.

Explorando mutações da caspase e modificação pós-translacional por abordagens de modelagem molecular

Apoptosis is a type of programmed cell death that eliminates damaged cells and controls the development and tissue homeostasis of multicellular organisms. ...

desempenham um papel fundamental na iniciação e execução usando modelagem molecular. Estudamos o efeito de modificações e mutações pós-tradução na estrutura e função da caspase. Descrevemos a abordagem da dinâmica molecular que dá uma visão da evolução da proteína, após a introdução de modificações estruturais no nível atômico.Tal abordagem é um significado particular para a compreensão dos mecanismos que regulam a morte celular do programa e distúrbios associados e o desenvolvimento de novas terapias eficazes. Demonstramos capacidades de modelagem molecular com uma das ferramentas mais populares, Amber 20. Mas o protocolo apresentado também pode ser usado com as versões formais deste software.O vídeo a seguir descreve instruções passo a passo sobre a preparação da estrutura da Caspase para simulação e a realização de uma investigação ascilica deste tipo selvagem de proteínas e formas modificadas. Para recuperar o banco de dados de proteína selecionado ou a estrutura P D B, use a lista suspensa de arquivos de download e clique no formato P D B. Remova observações e dados de conectividade e insira um cartão T E R entre cadeias de proteínas separadas no arquivo P D B.Para preparar o modelo inicial, inicie o programa T-LEAP do pacote Amber Tools. Em seguida, carregue o campo de força F F 14 S B para descrever a proteína com mecânica molecular e os parâmetros para moléculas de água e íons atômicos, como sódio e cloreto, inserindo os comandos indicados na linha de comando. Em seguida, carregue o arquivo P D B e construa coordenadas para hidrogênios criando um objeto chamado mol.Verifique se há inconsistências internas que possam causar problemas usando o comando check mol. Crie uma caixa de solvente ao redor da proteína. Em seguida, verifique a carga total inserindo mol de carga e adicione íons contadores para neutralizar o sistema.Crie o arquivo de topologia P R M top e a coordenada I N P C R D digitando o comando indicado. Uma vez feito, feche o programa T-lEAP. Realizar o primeiro estágio da minimização de energia para otimizar as posições dos átomos de hidrogênio e moléculas de água adicionados, mantendo as coordenadas de proteína fixadas por restrições posicionais em átomos pesados.Execute o programa P M E M D inserindo o comando indicado. Siga os argumentos necessários enquanto controlo os dados. P topologia molecular, parâmetros de campo de força e nomes de átomos.Coordenadas iniciais C. O usuário legível log saída R coordenadas finais. REF, coordenadas de referência para restrições de posição.Em seguida, realize o segundo estágio da minimização de energia sem restrições para otimizar todo o sistema usando as entradas do comando indicado. Esta etapa visa aquecer o sistema de zero a 300 kelvins. Realizar o processo de aquecimento com restrições posicionais nos átomos de proteína com 50 picossegundos em volume constante, usando o comando dado como entradas.Siga o argumento necessário como X conjuntos de coordenadas salvos ao longo da trajetória da dinâmica molecular. O próximo estágio é necessário para ajustar a densidade da água e obter o estado de equilíbrio da proteína. Efectuar o equilíbrio a 300 kelvins durante 500 picossegundos a pressão constante, sem quaisquer restrições, utilizando o comando indicado após o equilíbrio ter sido alcançado com êxito.Realizar uma simulação de dinâmica molecular de produção por 10 nanossegundos ou mais a pressão constante e gerar o arquivo de trajetória para posterior análise da estrutura proteica usando o comando indicado. A versão paralela do programa, P M E M D M P I ou G P U versão acelerada P M E M D cuda pode ser usada em clusters de computadores e supercomputadores. Uma longa simulação de dinâmica molecular pode ser dividida em vários segmentos e realizada sequencialmente.Nesta análise, a Caspase dois do tipo selvagem e seu mutante Serina-384 alanina foram investigados seguindo o fluxo de trabalho de modelagem de dinâmica molecular. A substituição de serina-384 alanina induziu importante alteração conformacional no resíduo do sítio ativo Arginina-378. Além disso, foi demonstrado que a substituição de serina-384 alanina afetou o reconhecimento do substrato pelos resíduos de arginina no sítio ativo, prejudicando a atividade de base fundida.A mutagênese direcionada ao local e os testes bioquímicos demonstraram que a mutação Serina-384 alanina bloqueou a atividade enzimática no processamento da Caspase dois e suprimiu a apoptose das células cancerígenas. A abordagem descrita pode ser usada para avaliar o efeito de mutações de aminoácidos e modificações pós-traducionais na caspase dois, bem como em outras caspase envolvidas em várias doenças cancerígenas.

exploring-caspase-mutations-post-translational-modification-molecular

Research

JoVE Journal

Biology

Biologia

O infarto do miocárdio (IM) não é apenas seguido por comprometimento da função cardíaca, mas também por apoptose na amígdala, uma região do cérebro envolvida nas consequências comportamentais do IM. Este protocolo descreve como induzir IM, coletar tecido da amígdala e medir a atividade da caspase-3, um marcador de apoptose.

Caspase-3 Atividade no Rat Amygdala Medido pelo espectrofluorometria Após Infarto do Miocárdio

Para detectar células saudáveis em animais inteiros que contêm baixos níveis de actividade de caspase, o biossensor altamente sensível designado CaspaseTracker foi gerado para Drosophila. biossensor actividade dependente de caspase é detectado em células saudáveis ​​de vida longa ao longo dos órgãos internos dos animais adultos criados em condições optimizadas na ausência de estímulos de morte.

In Vivo Biosensor Tracks apoptótica não-Caspase Atividade em Drosophila

Este protocolo descreve caspase Bimolecular complementação de fluorescência (BiFC); um método baseado em imagens que pode ser usado para visualizar induzida proximidade das caspases de iniciador, que é o primeiro passo para sua ativação.

Iluminando os caminhos para a ativação de Caspase usando fluorescência Bimolecular complementação

Este protocolo descreve o fluxo de trabalho para obter macrófagos derivados de monócitos (MDM) a partir de amostras de sangue humano, um método simples para introduzir eficientemente os repórteres de complementação de fluorescência bimolecular inflamatória (BiFC) em MDM humano sem comprometer a viabilidade e o comportamento celular, e uma abordagem baseada em imagens para medir a ativação inflamatória de caspase em células vivas.

Visualização de Caspases Inflamatórias Induzidas proximidade em macrófagos derivados de monócitos humanos

A infecção pelo vírus da gripe A (IAV) ativa as caspases que clivam o hospedeiro e as proteínas virais, que, por sua vez, têm funções pró e antivirais. Empregando-se inibidores, interferência de RNA, mutagênese sítio-dirigida e técnicas de western blotting e RT-qPCR, foram identificadas caspases em células de mamíferos infectadas que clivam a cortactina do hospedeiro e as histonas desacetilases.

Identificando Caspases e seus motivos que clivam proteínas durante a infecção pelo vírus da gripe A

Apresenta-se aqui um protocolo imuno-histoquímico detalhado para identificar, validar e direcionar caspases funcionalmente relevantes em tecidos complexos.

Quantificação da Caspase-9 Imunocorada em Tecido Retiniano

Este protocolo descreve um método abrangente para avaliar a ativação da caspase (caspase-1, caspase-3, caspase-7, caspase-8, caspase-9 e caspase-11) em resposta a modelos in vitro e in vivo (em camundongos) de infecção, insultos estéreis e câncer para determinar o início de vias de morte celular, como piroptose, apoptose, necroptose e PANoptose.

Avaliação da ativação da caspase para avaliar a morte de células imunes inatas

Aqui apresentamos um protocolo para expressar e purificar Drosophila caspases recombinantes Dronc e Drice, e seu uso em ensaios de clivagem in vitro .

In vitro Ensaios de Clivagem usando Caspases de Drosophila Recombinante Purificada para Triagem de Substrato

O presente protocolo descreve dois métodos para medir a atividade da caspase através de um substrato fluorogênico usando citometria de fluxo ou um espectrofluorômetro.

Medindo a atividade da caspase usando um ensaio fluorométrico ou citometria de fluxo

Aqui, descrevemos um protocolo para a otimização e parametrização de resíduos de aminoácidos modificados com espécies carbonilas reativas, adaptáveis a sistemas proteicos. As etapas do protocolo incluem projeto e otimização de estrutura, atribuições de carga, construção de parâmetros e preparação de sistemas de proteínas.

Sintetizando Aminoácidos Modificados com Carbonilas Reativas in Silico para Avaliar Efeitos Estruturais Usando Simulações de Dinâmica Molecular 

Explore simulações de dinâmica molecular usando Amber 20 para investigar mutações na caspase e modificações pós-traducionais. Este método auxilia no desenvolvimento de terapias direcionadas contra câncer.

Explorando mutações da caspase e modificação pós-translacional por abordagens de modelagem molecular

Explorar mais vídeos