O protocolo desenvolve um modelo tridimensional (3D) de um segmento dendrático com espinhos dendráticos para modelagem da plasticidade sináptica. A malha construída pode ser usada para modelagem computacional do tráfico de receptores AMPA na plasticidade sináptica de longo prazo usando o programa de software Blender com CellBlender e MCell.
A modelagem computacional da difusão e reação de espécies químicas em uma geometria tridimensional (3D) é um método fundamental para entender os mecanismos de plasticidade sináptica em espinhas dendríticas. Neste protocolo, a estrutura 3D detalhada dos dendritos e colunas dendríticas é modelada com malhas no blender de software com CellBlender. As regiões sinápticas e extrassinápticas são definidas na malha. Em seguida, o receptor sináptico e as moléculas de âncora sináptica são definidas com suas constantes de difusão. Finalmente, as reações químicas entre receptores sinápticos e âncoras sinápticas são incluídas e o modelo computacional é resolvido numericamente com o software MCell. Este método descreve o caminho espacial de cada molécula em uma estrutura geométrica 3D. Assim, é muito útil estudar o tráfico de receptores sinápticos dentro e fora das espinhas dendríticas durante a ocorrência de plasticidade sináptica. Uma limitação deste método é que o alto número de moléculas diminui a velocidade das simulações. A modelagem das espinhas dendríticas com este método permite o estudo da potencialização homossináptica e da depressão dentro de espinhas simples e plasticidade heterossináptica entre as espinhas dendríticas vizinhas.
A plasticidade sináptica tem sido associada ao aprendizado e à memória1. A plasticidade sináptica, como a potencialização a longo prazo (LTP) e a depressão de longo prazo (LTD), está associada, respectivamente, à inserção e remoção de receptores AMPA (AMPARs) dentro e fora da membrana sináptica2. As sinapses AMPAR estão localizadas no topo das pequenas estruturas de volume chamadas espinhos dendríticas3. Cada coluna contém uma região densa de proteína na membrana pós-sináptica chamada densidade pós-sináptica (PSD). Proteínas âncoras no PSD prendem AMPARs na região sináptica. Há poucas cópias de AMPARs dentro de uma única sinapse e o tráfico e reação de AMPARs com outras espécies em espinhas dendríticas é um processo estocástico2,4. Existem vários modelos compartimentais de tráfico de receptores sinápticos em lombadas dendríticas5,,6,,7,8. No entanto, há falta de modelos computacionais estocásticos do tráfico de AMPARs associados à plasticidade sináptica nas estruturas 3D dos dendritos e suas espinhas dendríticas.
A modelagem computacional é uma ferramenta útil para investigar os mecanismos subjacentes à dinâmica de sistemas complexos, como a retração-difusão de AMPARs em lombadas dendríticas durante a ocorrência de plasticidade sináptica9,,10,,11,,12. O modelo pode ser utilizado para visualizar cenários complexos, variando parâmetros sensíveis e fazendo previsões importantes em condições científicas envolvendo muitas variáveis difíceis ou impossíveis de controlar o experimental12,13. Definir o nível de detalhe de um modelo computacional é um passo fundamental na obtenção de informações precisas sobre o fenômeno modelado. Um modelo computacional ideal é um delicado equilíbrio entre complexidade e simplicidade para capturar as características essenciais dos fenômenos naturais sem ser computacionalmente proibitivo. Modelos computacionais muito detalhados podem ser caros para calcular. Por outro lado, sistemas pouco detalhados podem não ter os componentes fundamentais essenciais para capturar a dinâmica do fenômeno. Embora a modelagem 3D de lombadas dendríticas seja computacionalmente mais cara que 2D e 1D, existem condições, como em sistemas complexos com muitas variáveis não lineares reagindo e difundindo no tempo e espaço 3D, para o qual a modelagem em nível 3D é essencial para obter insights sobre o funcionamento do sistema. Além disso, a complexidade pode ser reduzida cuidadosamente para preservar as características essenciais de um modelo de dimensão inferior.
Em um sistema estocástico com poucas cópias de uma determinada espécie dentro de um pequeno volume, a dinâmica média do sistema se desvia da dinâmica média de uma grande população. Neste caso, é necessária a modelagem computacional estocástica de partículas de difusão de reação. Este trabalho introduz um método para modelagem estocástica de difusão de reação de algumas cópias de AMPARs em espinhas dendríticas 3D. O objetivo deste método é desenvolver um modelo computacional 3D de um segmento dendrático com espinhas dendríticas e suas sinapses para modelagem da plasticidade sináptica.
O método usa o software MCell para resolver o modelo numericamente, Blender para a construção de malhas 3D e CellBlender para criar e visualizar as simulações MCell, incluindo a reação espessora-difusão de moléculas em malhas 3D14,,15,,16. Blender é uma suíte para a criação de malhas e CellBlender é um complemento para o liquidificador de software base. MCell é um simulador de Monte Carlo para a difusão de reação de moléculas únicas17.
A lógica por trás do uso desse método consiste na modelagem da plasticidade sináptica para obter uma melhor compreensão desse fenômeno no ambiente microfisiológico das espinhas dendríticas14. Particularmente, este método permite a simulação de potencialização homossináptica, depressão homossináptica e plasticidade heterossináptica entre espinhas dendríticas14.
As características deste método incluem a modelagem da estrutura geométrica 3D do dendrite e suas sinapses, a difusão por caminhada aleatória e as reações químicas das moléculas envolvidas com a plasticidade sináptica. Este método oferece a vantagem de criar ambientes ricos para testar hipóteses e fazer previsões sobre o funcionamento de um sistema complexo não linear com um grande número de variáveis. Além disso, este método pode ser aplicado não apenas para estudar a plasticidade sináptica, mas também para estudar a difusão de reações estocásticas de moléculas em estruturas de malha 3D em geral.
Alternativamente, as malhas 3D de estruturas dendríticas podem ser construídas diretamente no Liquidificador a partir de reconstruções seriais de microscópio eletrônico18. Embora as malhas baseadas em reconstruções seriais forneçam estruturas 3D, o acesso aos dados experimentais nem sempre está disponível. Assim, a construção de malhas adaptadas a partir de estruturas geométricas básicas, conforme descrito no presente protocolo, proporciona flexibilidade para desenvolver segmentos dendríticos personalizados com colunas dendríticas.
Outro método computacional alternativo é a simulação em massa de reações bem mistas em um volume regular9,10,11,19,20,21,22. As simulações a granel são muito eficientes na resolução das reações de muitas espécies dentro de um único volume bem misturado23,mas a abordagem em massa é extremamente lenta para resolver a difusão de reação de moléculas dentro de muitos voxels bem misturados em uma malha 3D de alta resolução. Por outro lado, o presente método utilizando simulações MCell de difusão de reação de partículas individuais funciona eficientemente em malhas 3D de alta resolução15.
Antes de usar este método, deve-se perguntar se o fenômeno estudado requer uma abordagem de retração estocástica em uma malha 3D. Se o fenômeno tem poucas cópias (menos de 1.000) de pelo menos uma das espécies reagindo difundindo em uma estrutura geométrica complexa com pequenos compartimentos de volume, como espinhos dendráticos, então a modelagem estocástica de difusão de reação em malhas 3D é apropriada para a aplicação.
Existem várias etapas necessárias para a construção de um modelo computacional 3D de um segmento dendrático contendo espinhos dendráticos com plasticidade sináptica. As principais etapas são a instalação do software adequado para a construção do modelo, a construção de uma única coluna dendrítica a ser usada como modelo para criar múltiplas lombadas, e a criação de um segmento dendrático conectado com múltiplas espinhas dendríticas. O passo para modelagem da plasticidade sináptica consiste na inserção de âncoras na região do PSD e AMPARs no segmento dendrático e colunas dendríticas. Em seguida, são definidas reações cinéticas entre as âncoras localizadas no PSD e AMPARs para produzir espécies complexas de anchor-AMPAR que prendem os AMPARs na região sináptica. Respectivamente, o aumento e diminuição da afinidade entre as âncoras e as AMPARs sinápticas criam o processo de LTP e LTD.
Este artigo apresenta um método para a construção de malhas 3D para modelagem de processos de plasticidade sináptica de reação-difusão em um segmento dendrático com espinhas dendríticas. O modelo desenvolvido contém um segmento dendrítica com poucas espinhas dendríticas. A difusão lateral e a reação de AMPARs com âncoras sinápticas permitem a simulação da dinâmica basal. Os passos críticos do protocolo são o corte da esfera para a criação da parte superior da cabeça da coluna vertebral (<strong c…
The authors have nothing to disclose.
Este trabalho foi apoiado em parte pela Fundação de Ciência do Estado de São Paulo (FAPESP) #2015/50122-0 e IRTG-GRTK 1740/2, pela bolsa IBM/FAPESP #2016/18825-4, e pela concessão da FAPESP #2018/06504-4.
Blender | Blender Foundation | https://www.blender.org/ | |
CellBlender | University of Pittsburgh | https://mcell.org/ | |
Mcell | University of Pittsburgh | https://mcell.org/ |