Multiscale modeling of presynaptic dynamics from molecular to mesoscale

The architecture of invaginating rod synapses slows glutamate diffusion and shapes synaptic responses

Synapses of retinal rod photoreceptors involve deep invaginations occupied by second-order rod bipolar cell (RBP) and horizontal cell (HC) dendrites. Synaptic vesicles are released into this invagination at multiple sites beneath an elongated presynaptic ribbon. To study the impact of this architecture on glutamate diffusion and receptor activity, we reconstructed four rod terminals and their postsynaptic dendrites from serial electron micrographs of the mouse retina. We incorporated these structures into anatomically realistic Monte Carlo simulations of neurotransmitter diffusion and receptor activation. By comparing passive diffusion of glutamate in realistic structures with geometrically simplified models, we found that glutamate exits anatomically realistic synapses 10-fold more slowly than previously predicted. Constraining simulations with physiological data, we modeled activity of EAAT5 glutamate transporters in rods, AMPA receptors on HC dendrites, and metabotropic glutamate receptors (mGluR6) on RBP dendrites. Simulations suggested that ∼3,000 EAAT5 populate rod membranes. While uptake by surrounding glial Müller cells retrieves most glutamate released by rods, binding and uptake by EAAT5 influence RBP kinetics. Glutamate persistence allows mGluR6 on RBP dendrites to integrate the stream of vesicles released by rods in darkness. Glutamate's tortuous diffusional path confers quantal variability, as release from nearby ribbon sites exerts larger effects on RBP and HC receptors than release from more distant sites. Temporal integration supports slower sustained release rates, but additional quantal variability can impede postsynaptic detection of changes in release produced by rod light responses. These results show an example of the profound impact that synaptic architecture can have on postsynaptic responses.

Integrated multi-omic characterizations of the synapse reveal RNA processing factors and ubiquitin ligases associated with neurodevelopmental disorders

The molecular composition of the excitatory synapse is incompletely defined due to its dynamic nature across developmental stages and neuronal populations. To address this gap, we apply proteomic mass spectrometry to characterize the synapse in multiple biological models, including the fetal human brain and human induced pluripotent stem cell (hiPSC)-derived neurons. To prioritize the identified proteins, we develop an orthogonal multi-omic screen of genomic, transcriptomic, interactomic, and structural data. This data-driven framework identifies proteins with key molecular features intrinsic to the synapse, including characteristic patterns of biophysical interactions and cross-tissue expression. The multi-omic analysis captures synaptic proteins across developmental stages and experimental systems, including 493 synaptic candidates supported by proteomics. We further investigate three such proteins that are associated with neurodevelopmental disorders-Cullin 3 (CUL3), DEAD-box helicase 3 X-linked (DDX3X), and Y-box binding protein-1 (YBX1)-by mapping their networks of physically interacting synapse proteins or transcripts. Our study demonstrates the potential of an integrated multi-omic approach to more comprehensively resolve the synaptic architecture.

Geometric tortuosity at invaginating rod synapses slows glutamate diffusion and shapes synaptic responses: insights from anatomically realistic Monte Carlo simulations

At the first synapse in the vertebrate retina, rod photoreceptor terminals form deep invaginations occupied by multiple second-order rod bipolar and horizontal cell (RBP and HC) dendrites. Synaptic vesicles are released into this invagination at multiple sites beneath an elongated presynaptic ribbon. We investigated the impact of this complex architecture on the diffusion of synaptic glutamate and activity of postsynaptic receptors. We obtained serial electron micrographs of mouse retina and reconstructed four rod terminals along with their postsynaptic RBP and HC dendrites. We incorporated these structures into an anatomically realistic Monte Carlo simulation of neurotransmitter diffusion and receptor activation. We compared passive diffusion of glutamate in these realistic structures to existing, geometrically simplified models of the synapse and found that glutamate exits anatomically realistic synapses ten times more slowly than previously predicted. By comparing simulations with electrophysiological recordings, we modeled synaptic activation of EAAT5 glutamate transporters in rods, AMPA receptors on HC dendrites, and metabotropic glutamate receptors (mGluR6) on RRBP dendrites. Our simulations suggested that ~3,000 EAAT5 transporters populate the rod presynaptic membrane and that, while uptake by surrounding glial Müller cells retrieves much of the glutamate released by rods, binding and uptake by EAAT5 influences RBP response kinetics. The long lifetime of glutamate within the cleft allows mGluR6 on RBP dendrites to temporally integrate the steady stream of vesicles released at this synapse in darkness. Glutamate's tortuous diffusional path through realistic synaptic geometry confers quantal variability, as release from nearby ribbon sites exerts larger effects on RBP and HC receptors than release from more distant sites. While greater integration may allow slower sustained release rates, added quantal variability complicates the challenging task of detecting brief decreases in release produced by rod light responses at scotopic threshold.

Aβ-mediated synaptic glutamate dynamics and calcium dynamics in astrocytes associated with Alzheimer's disease

The accumulation of amyloid β peptide A β is assumed to be one of the main causes of Alzheimer's disease AD . There is increasing evidence that astrocytes are the primary targets of Aβ. Aβ can cause abnormal synaptic glutamate, aberrant extrasynaptic glutamate, and astrocytic calcium dysregulation through astrocyte glutamate transporters facing the synaptic cleft (GLT-syn), astrocyte glutamate transporters facing the extrasynaptic space (GLT-ess), metabotropic glutamate receptors in astrocytes (mGluR), N-methyl-D-aspartate receptors in astrocytes (NMDAR), and glutamatergic gliotransmitter release (Glio-Rel). However, it is difficult to experimentally identify the extent to which each pathway affects synaptic glutamate, extrasynaptic glutamate, and astrocytic calcium signaling. Motivated by these findings, this work established a concise mathematical model of astrocyteCa 2 + dynamics, including the above Aβ-mediated glutamate-related multiple pathways. The model results presented the extent to which five mechanisms acted upon by Aβ affect synaptic glutamate, extrasynaptic glutamate, and astrocytic intracellularCa 2 + signals. We found that GLT-syn is the main pathway through which Aβ affects synaptic glutamate. GLT-ess and Glio-Rel are the main pathways through which A β affects extrasynaptic glutamate. GLT-syn, mGluR, and NMDAR are the main pathways through which Aβ affects astrocytic intracellularCa 2 + signals. Additionally, we discovered a strong, monotonically increasing relationship between the mean glutamate concentration and the meanCa 2 + oscillation amplitude (or frequency). Our results may have therapeutic implications for slowing cell death induced by the combination of glutamate imbalance andCa 2 + dysregulation in AD.

CellVis2: a conference on visualizing the molecular cell

In January 2024, a targeted conference, 'CellVis2', was held at Scripps Research in La Jolla, USA, the second in a series designed to explore the promise, practices, roadblocks, and prospects of creating, visualizing, sharing, and communicating physical representations of entire biological cells at scales down to the atom.

Making time and space for calcium control of neuron activity

Calcium directly controls or indirectly regulates numerous functions that are critical for neuronal network activity. Intracellular calcium concentration is tightly regulated by numerous molecular mechanisms because spatial domains and temporal dynamics (not just peak amplitude) are critical for calcium control of synaptic plasticity and ion channel activation, which in turn determine neuron spiking activity. The computational models investigating calcium control are valuable because experiments achieving high spatial and temporal resolution simultaneously are technically unfeasible. Simulations of calcium nanodomains reveal that specific calcium sources can couple to specific calcium targets, providing a mechanism to determine the direction of synaptic plasticity. Cooperativity of calcium domains opposes specificity, suggesting that the dendritic branch might be the preferred computational unit of the neuron.