Synaptische plasticiteit is geassocieerd met leren en geheugen1. Synaptische plasticiteit, zoals langdurige potentiëring (LTP) en langdurige depressie (LTD), wordt geassocieerd met het inbrengen en verwijderen van AMPA-receptoren (AMPARs) in en uit het synaptische membraan2. De AMPAR synapsen bevinden zich bovenop de kleine volumestructuren genaamd dendritische stekels3. Elke wervelkolom bevat een eiwit dichte regio in de postsynaptische membraan genaamd de postsynaptische dichtheid (PSD). Veranker eiwitten bij de PSD trap AMPARs in het synaptische gebied. Er zijn weinig kopieën van AMPARs binnen een enkele synaps en de handel en reactie van AMPARs met andere soorten in dendritische stekels is een stochastisch proces2,4. Er zijn verschillende compartimenteringsmodellen van synaptische receptorhandel bij dendritische stekels5,6,7,8. Echter, er is een gebrek aan stochastische computationele modellen van de handel in AMPARs in verband met synaptische plasticiteit op de 3D-structuren van de dendrieten en hun dendritische stekels.
Computationele modellering is een nuttig hulpmiddel om de mechanismen te onderzoeken die ten grondslag liggen aan de dynamiek van complexe systemen, zoals de reactieverspreiding van AMPARs in dendritische stekels tijdens het optreden van synaptische plasticiteit9,10,11,12. Het model kan worden gebruikt om complexe scenario's te visualiseren, gevoelige parameters te variëren en belangrijke voorspellingen te doen in wetenschappelijke omstandigheden waarbij veel variabelen betrokken zijn die moeilijk of onmogelijk te controleren zijn12,13. Het definiëren van het detailniveau van een computationeel model is een fundamentele stap in het verkrijgen van nauwkeurige informatie over het gemodelleerde fenomeen. Een ideaal computationeel model is een delicaat evenwicht tussen complexiteit en eenvoud om de essentiële kenmerken van de natuurlijke fenomenen vast te leggen zonder computationeel onbetaalbaar te zijn. Computermodellen die te gedetailleerd zijn, kunnen duur zijn om te berekenen. Aan de andere kant kunnen systemen die slecht gedetailleerd zijn, de fundamentele componenten missen die essentieel zijn om de dynamiek van het fenomeen vast te leggen. Hoewel 3D-modellering van dendritische stekels rekenkundig duurder is dan 2D en 1D, zijn er voorwaarden, zoals in complexe systemen met veel niet-lineaire variabelen die reageren en verspreiden in tijd en 3D-ruimte, waarvoor modellering op 3D-niveau essentieel is om inzicht te krijgen in de werking van het systeem. Verder kan de complexiteit zorgvuldig worden verminderd om de essentiële kenmerken van een lager dimensionaal model te behouden.
In een stochastisch systeem met weinig exemplaren van een bepaalde soort binnen een klein volume, wijkt de gemiddelde dynamiek van het systeem af van de gemiddelde dynamiek van een grote populatie. In dit geval is de stochastische computationele modellering van reactieverspreidende deeltjes vereist. Dit werk introduceert een methode voor stochastische modellering reactie-diffusie van een paar exemplaren van AMPARs in 3D dendritische stekels. Het doel van deze methode is het ontwikkelen van een 3D computationeel model van een dendritisch segment met dendritische stekels en hun synapsen voor het modelleren van synaptische plasticiteit.
De methode maakt gebruik van de software MCell om het model numeriek op te lossen, Blender voor het construeren van 3D-mazen en CellBlender om de MCell-simulaties te maken en te visualiseren, inclusief de spatiotemporale reactieverspreiding van moleculen in 3D-mazen14,15,16. Blender is een suite voor het maken van mazen en CellBlender is een add-on voor de basissoftware Blender. MCell is een Monte Carlo simulator voor de reactie-verspreiding van enkele moleculen17.
De reden achter het gebruik van deze methode bestaat uit het modelleren van synaptische plasticiteit om een beter begrip van dit fenomeen te bereiken in de microfysiologische omgeving van de dendritische stekels14. Met name deze methode maakt de simulatie van homosynaptische potentiëring, homosynaptische depressie, en heterosynaptische plasticiteit tussen dendritische stekels14.
De kenmerken van deze methode zijn het modelleren van de 3D geometrische structuur van het dendriet en de synapsen, de verspreiding door willekeurige wandeling, en de chemische reacties van de moleculen die betrokken zijn bij synaptische plasticiteit. Deze methode biedt het voordeel van het creëren van rijke omgevingen om hypothesen te testen en voorspellingen te doen over de werking van een complex niet-lineair systeem met een groot aantal variabelen. Bovendien kan deze methode niet alleen worden toegepast voor het bestuderen van synaptische plasticiteit, maar ook voor het bestuderen van stochastische reactie-diffusie van moleculen in 3D-meshstructuren in het algemeen.
Als alternatief kunnen 3D-mazen van dendritische structuren direct in Blender worden geconstrueerd uit seriële reconstructies van elektronenmicroscoop18. Hoewel mazen op basis van seriële reconstructies 3D-structuren bieden, is toegang tot de experimentele gegevens niet altijd beschikbaar. Zo biedt de constructie van mazen aangepast aan basis geometrische structuren, zoals beschreven in het huidige protocol, flexibiliteit om aangepaste dendritische segmenten met dendritische stekels te ontwikkelen.
Een andere alternatieve rekenmethode is de bulksimulatie van goed gemengde reacties in een regelmatig volume9,10,11,19,20,21,22. De bulk simulaties zijn zeer efficiënt in het oplossen van de reacties van vele soorten binnen een enkele goed gemengde volume23, maar de bulk aanpak is uiterst traag om de reactie-diffusie van moleculen binnen veel goed gemengde voxels in een hoge resolutie 3D-mesh op te lossen. Anderzijds werkt de huidige methode met behulp van MCell-simulaties van reactieverspreiding van individuele deeltjes efficiënt in 3D-mazen met hoge resolutie15.
Alvorens deze methode te gebruiken, moet men zich afvragen of het onderzochte fenomeen een stochastische reactie-diffusiebenadering in een 3D-mesh vereist. Als het fenomeen weinig kopieën heeft (minder dan 1.000) van ten minste één van de reagerende soorten die zich verspreiden in een complexe geometrische structuur met kleine volumecompartimenten zoals dendritische stekels, dan is stochastische modellering van reactieverspreiding in 3D-mazen geschikt voor de toepassing.
Er zijn verschillende stappen nodig om een 3D-computationeel model te bouwen van een dendritisch segment met dendritische stekels met synaptische plasticiteit. De belangrijkste stappen zijn de installatie van de juiste software voor de bouw van het model, de bouw van een enkele dendritische wervelkolom te gebruiken als een sjabloon om meerdere stekels te creëren, en de oprichting van een dendritisch segment dat is verbonden met meerdere dendritische stekels. De stap voor het modelleren van synaptische plasticiteit bestaat uit het plaatsen van ankers in de PSD-regio en AMPARs in het dendritische segment en dendritische stekels. Vervolgens worden kinetische reacties tussen de ankers van de PSD en AMPARs gedefinieerd om complexe anker-AMPAR-soorten te produceren die de AMPARs in het synaptische gebied vangen. Respectievelijk de toename en afname van de affiniteit tussen de ankers en de synaptische AMPARs creëren het proces van LTP en LTD.