Soliton-like Waves in a Two-Dimensional Recurrent Spiking Neural Network with Weighted Spike-Timing-Dependent Plasticity

Technische Samenvatting: Soliton-achtige Golven in een Twee-Dimensionaal Recurrent Spiking Neuraal Netwerk met Gewogen Spike-Timing-Dependent Plasticiteit

Probleemstelling
Reizende golven van neurale activiteit zijn alomtegenwoordig in de cortex, maar de mechanismen waarmee zij ontstaan en zichzelf in recurrente circuits in stand houden, blijven onvolledig begrepen. Hoewel neurale veldmodellen reizende golven en solitaire "bumps" hebben voorspeld op basis van "Mexicaanse hoed"-connectiviteitsprofielen (kort bereikende excitatie gebalanceerd door lang bereikende inhibitie), abstraheren deze modellen doorgaans van individuele spikes en synaptische plasticiteit. Daarentegen is bekend dat Spike-Timing-Dependent Plasticiteit (STDP) receptieve velden vormgeeft en paden versterkt, maar het is onduidelijk of een recurrent excitatoir-inhibitoir (E/I) netwerk uitgerust met STDP autonoom stabiele, zelf-propagerende golfpakketjes (soliton-achtige structuren) kan genereren vanuit een uniforme staat, en of deze structuren specifieke botsingsdynamica vertonen zoals annihilatie.

Methodologie
De auteurs hebben een minimaal, biologisch plausibel computationeel model geconstrueerd van een twee-dimensionaal recurrent netwerk bestaande uit 200 × 100 paren excitatore en inhibitoire neuronen gerangschikt op een discreet rooster. Het model werkt in discrete tijdstappen met synchrone updates.

• Neuronmodel: Elk neuron integreert gewogen spikes van presynaptische buren. Een cruciaal kenmerk is divisieve normalisatie van de synaptische input, waarbij de totale excitatore drive wordt genormaliseerd door de som van de excitatore synaptische gewichten ($S_E$). Deze normalisatie wordt voorgesteld als een biologisch plausibele dendritische implementatie via lokale attenuatie bij binaire juncties. De drempelwaarde ($\theta$) ondergaat homeostatische adaptatie: deze daalt naar een minimum ($\theta_{min}$) tijdens stilte om "stall states" te voorkomen en stijgt naar een maximum ($\theta_{max}$) bij het vuren om spike-frequentie-adaptatie te implementeren.
• Connectiviteit: Het netwerk vertoont een geometrische asymmetrie: de excitatore-naar-excitatore ($E \to E$) verbindingen hebben een grotere radius ($r_{EE}=9$) dan de excitatore-naar-inhibitoire ($E \to I$) en inhibitoire-naar-excitatore ($I \to E$) verbindingen ($r_{EI}=r_{IE}=5$). Er bestaan geen $I \to I$ verbindingen. De initiële gewichten zijn uniform, waarbij de inhibitoire gewichten aanzienlijk hoger zijn ingesteld dan de excitatore gewichten om het kleinere aantal inhibitoire inputs te compenseren.
• Leerregel: Synaptische gewichten evolueren via Gewogen STDP (WSTDP), een multiplicatieve, temporeel asymmetrische leerregel. De regel zorgt ervoor dat gewichten begrensd blijven zonder harde afkapping (clipping), en convergeert naar een vast punt bepaald door de balans tussen de Long-Term Potentiation (LTP) en Long-Term Depression (LTD) gebeurtenissen.
• Stimulatie: Het netwerk wordt onderworpen aan herhaalde, gelokaliseerde stimulatie (het dwingen van een schijf van excitatore neuronen om te vuren) met een vaste periode. Experimenten omvatten enkelvoudige bron-stimulatie en dubbele bron-stimulatie met variërende fasen en frequenties om interacties tussen golven te bestuderen.

Belangrijkste Bijdragen en Resultaten

1. Emergentie van Dissipatieve Solitonen: Startend vanaf een ruimtelijk uniforme connectiviteit, induceert herhaalde gelokaliseerde stimulatie de spontane vorming van stabiele, zelf-propagerende golfpakketjes. Deze pakketjes behouden een constante ruimtelijke profiel (een heldere excitatore ring gevolgd door een inhibitoire staart) en propageren met een constante snelheid. Cruciaal is dat deze golven, eenmaal gevestigd, ook blijven propageren nadat de externe stimulatie is verwijderd, wat aantoont dat het zelf-onderhoudende structuren zijn.
2. Rol van Geometrische Asymmetrie en WSTDP: De emergentie van deze golven vereist twee condities:
   • Geometrische Asymmetrie: De excitatore verbindingsradius moet groter zijn dan de inhibitoire radius om een functioneel Mexicaanse-hoed-profiel te creëren.
   • WSTDP: De leerregel "graveert" de richting van propagatie in het synaptische gewichtslandschap. Synapsen in de richting van de propagatie ondergaan netto LTP, terwijl synapsen in de tegenovergestelde richting netto LTD ondergaan. Dit creëert een directionele bias die achterwaartse propagatie voorkomt en de golf stabiliseert.
3. Botsingsdynamica en Annihilatie: Wanneer twee golfpakketjes frontaal botsen, annihileren ze in plaats van door elkaar heen te bewegen. Dit gedrag is kenmerkend voor dissipatieve solitonen. De annihilatie wordt veroorzaakt door de inhibitoire staart van de ene golf die de excitatore front van de andere onderdrukt.
4. Semi-persistente Ruimtelijke Grenzen: In dubbele bron-experimenten vormt het botsingspunt een "semi-persistente" grens, gecodeerd in de synaptische gewichten.
   • Fase-encodering: Als bronnen in fase zijn, vormt de grens zich op het geometrische middelpunt. Als er een faseverschuiving is, verschuift de grens naar de achterlopende bron.
   • Frequentiecompetitie: Als bronnen verschillende frequenties hebben, overschrijft de snellere bron progressief de synaptische gewichten van de tragere bron, waardoor de grens verschuift. Een voldoende groot frequentievoordeel leidt tot een volledige territoriale overname.
   • Robuustheid: Deze grenzen blijven bestaan zelfs nadat één bron is uitgeschakeld, hoewel hun stabiliteit wordt gemoduleerd door de adaptieve drempelparameters.
5. Noodzaak van Sterke Initiële Inhibitie: Het model toont aan dat sterke initiële inhibitoire gewichten een geometrische noodzaak zijn. Het verminderen van deze gewichten leidt tot globale synchrone oscillaties (een epileptische-achtige regime) in plaats van gelokaliseerde solitonen, wat de cruciale rol van inhibitie benadrukt bij het inkaderen van de excitatore golffront.

Betekenis en Claims
Het artikel claimt een minimaal computationeel kader te bieden voor het bestuderen van de emergentie van corticale reizende golven, de afbakening van activiteitszones en ruimtelijk geheugen, uitsluitend afgeleid van lokale plasticiteitsregels.

• Mechanisme van Emergentie: De studie suggereert dat complexe spatiotemporele dynamica, zoals stabiele reizende golven en botsingsannihilatie, kunnen voortkomen uit eenvoudige lokale regels (WSTDP en homeostatische drempeladaptatie), mits er een geometrische asymmetrie in de connectiviteit aanwezig is.
• Ruimtelijk Geheugen: De auteurs stellen voor dat het synaptische gewichtslandschap in dergelijke netwerken een vorm van ruimtelijk geheugen implementeert. Het netwerk "onthoudt" de relatieve timing (fase) en snelheid (frequentie) van voorgaande inputs via de positie en stabiliteit van de annihilatiegrenzen.
• Directionele Stabiliteit: Het WSTDP-mechanisme maakt niet alleen propagatie mogelijk, maar stabiliseert deze ook actief door onstabiele reverse modi te onderdrukken, wat ongecontroleerde oscillaties voorkomt.
• Biologische Plausibiliteit: Het model biedt een potentiële verklaring voor hoe corticale circuits onderscheidende activiteitsdomeinen kunnen segregeren en behouden, en hoe reizende golven kunnen worden gemoduleerd door concurrerende inputs, alles binnen een kader van lokale, spike-gebaseerde plasticiteit.

De auteurs merken op dat hoewel het model minimaal is (discrete tijd, binaire spikes), de robuustheid van de soliton-emergentie suggereert dat deze fenomenen fundamentele eigenschappen zijn van recurrente E/I-netwerken met lokale plasticiteit, in plaats van artefacten van specifieke parametervariaties. Toekomstig werk wordt gesuggereerd om asynchrone propagatie en gestructureerde synaptische afferenten te verkennen.