A computational model of the ribbon synapse

A dynamical adaptation model of visual spatiotemporal processing in cones and horizontal cells

In this work, we introduce a phenomenological model for the cone-horizontal cell assembly, including spatial integration and formation of receptive field-like structures. The model extends our previous dynamical adaptation description with gain control accounting for processes in single cones, valid in severe nonlinear regimes. Here, a spatially extended feedback mechanism is introduced from horizontal cells to cones to account for experimental evidence, contributing thus to the development of a center-surround receptive field in cones and downstream bipolar cells. Feedback gain is defined on different spatial scales by weighting spatial filters: a short scale accounting for cone input to the feedback mechanism and a large scale driven by the syncytium characteristics of horizontal cells. A third spatial scale improves the description, mimicking neighboring cone-cone coupling. This overall spatial integration couples to temporal signal processing, thus obtaining a spatiotemporal model of outer retina responses capable of reproducing nonlinear features in both dimensions (space and time). The model was tested and validated using measurements on horizontal cells from different studies, with excellent performance. By its phenomenological nature, signal processing properties are inferred from model parameters. The model can be used in arrays of processing units with more complex incoming patterns of visual stimuli.

Model-based comparison of current flow in rod bipolar cells of healthy and early-stage degenerated retina

Retinal degenerative diseases, such as retinitis pigmentosa, are generally thought to initiate with the loss of photoreceptors, though recent work suggests that plasticity and remodeling occurs prior to photoreceptor cell loss. This degeneration subsequently leads to death of other retinal neurons, creating functional alterations and extensive remodeling of retinal networks. Retinal prosthetic devices stimulate the surviving retinal cells by applying external current using implanted electrodes. Although these devices restore partial vision, the quality of restored vision is limited. Further knowledge about the precise changes in degenerated retina as the disease progresses is essential to understand how current flows in retinas undergoing degenerative disease and to improve the performance of retinal prostheses. We developed computational models that describe current flow from rod photoreceptors to rod bipolar cells (RodBCs) in the healthy and early-stage degenerated retina. Morphologically accurate models of retinal cells with their synapses are constructed based on retinal connectome datasets, created using serial section transmission electron microscopy (TEM) images of 70 nm-thick slices of either healthy (RC1) or early-stage degenerated (RPC1) rabbit retina. The passive membrane and active ion currents of each cell are implemented using conductance-based models in the Neuron simulation environment. In response to photocurrent input at rod photoreceptors, the simulated membrane potential at RodBCs in early degenerate tissue is approximately 10-20 mV lower than that of RodBCs of that observed in wild type retina. Results presented here suggest that although RodBCs in RPC1 show early, altered morphology compared to RC1, the lower membrane potential is primarily a consequence of reduced rod photoreceptor input to RodBCs in the degenerated retina. Frequency response and step input analyses suggest that individual cell responses of RodBCs in either healthy or early-degenerated retina, prior to substantial photoreceptor cell loss, do not differ significantly.

On optimal coupling of the 'electronic photoreceptors' into the degenerate retina

Objective

To restore sight in atrophic age-related macular degeneration, the lost photoreceptors can be replaced with electronic implants, which replicate their two major functions: (1) converting light into an electric signal, and (2) transferring visual information to the secondary neurons in the retinal neural network—the bipolar cells (BC). We study the selectivity of BC activation by subretinal implants and dynamics of their response to pulsatile waveforms in order to optimize the electrical stimulation scheme such that retinal signal processing with 'electronic photoreceptors' remains as close to natural as possible.

Approach

A multicompartmental model of a BC was implemented to simulate responses of the voltage-gated calcium channels and subsequent synaptic vesicle release under continuous and pulsatile stimuli. We compared the predicted response under various frequencies, pulse durations, and alternating gratings to the corresponding experimental measurements. In addition, electric field was computed for various electrode configurations in a 3-d finite element model to assess the stimulation selectivity via spatial confinement of the field.

Main results

The modeled BC-mediated retinal responses were, in general, in good agreement with previously published experimental results. Kinetics of the calcium pumps and of the neurotransmitter release in ribbon synapses, which underpin the BC's temporal filtering and rectifying functions, allow mimicking the natural BC response with high frequency pulsatile stimulation, thereby preserving features of the retinal signal processing, such as flicker fusion, adaptation to static stimuli and non-linear summation of subunits in receptive field. Selectivity of the BC stimulation while avoiding direct activation of the downstream neurons (amacrine and ganglion cells—RGCs) is improved with local return electrodes.

Significance

If the retinal neural network is preserved to a large extent in age-related macular degeneration, selective stimulation of BCs with proper spatial and temporal modulation of the extracellular electric field may retain many features of the natural retinal signal processing and hence allow highly functional restoration of sight.

Sensory Processing at Ribbon Synapses in the Retina and the Cochlea

In recent years, sensory neuroscientists have made major efforts to dissect the structure and function of ribbon synapses which process sensory information in the eye and ear. This review aims to summarize our current understanding of two key aspects of ribbon synapses: 1) their mechanisms of exocytosis and endocytosis and 2) their molecular anatomy and physiology. Our comparison of ribbon synapses in the cochlea and the retina reveals convergent signaling mechanisms, as well as divergent strategies in different sensory systems.

Increasing Electrical Stimulation Efficacy in Degenerated Retina: Stimulus Waveform Design in a Multiscale Computational Model

A computational model of electrical stimulation of the retina is proposed for investigating current waveforms used in prosthetic devices for restoring partial vision lost to retinal degenerative diseases. The model framework combines a connectome-based neural network model characterized by accurate morphological and synaptic properties with an admittance method model of bulk tissue and prosthetic electronics. In this model, the retina was computationally "degenerated," considering cellular death and anatomical changes that occur early in disease, as well as altered neural behavior that develops throughout the neurodegeneration and is likely interfering with current attempts at restoring vision. A resulting analysis of stimulation range and threshold of ON ganglion cells within the retina that are either healthy or in beginning stages of degeneration is presented for currently used stimulation waveforms, and an asymmetric biphasic current stimulation for subduing spontaneous firing to allow increased control over ganglion cell firing patterns in degenerated retina is proposed. Results show that stimulation thresholds of retinal ganglion cells do not notably vary after beginning stages of retina degeneration. In addition, simulation of proposed asymmetric waveforms showed the ability to enhance the control of ganglion cell firing via electrical stimulation.

AMPA receptors at ribbon synapses in the mammalian retina: kinetic models and molecular identity

In chemical synapses, neurotransmitter molecules released from presynaptic vesicles activate populations of postsynaptic receptors that vary in functional properties depending on their subunit composition. Differential expression and localization of specific receptor subunits are thought to play fundamental roles in signal processing, but our understanding of how that expression is adapted to the signal processing in individual synapses and microcircuits is limited. At ribbon synapses, glutamate release is independent of action potentials and characterized by a high and rapidly changing rate of release. Adequately translating such presynaptic signals into postsynaptic electrical signals poses a considerable challenge for the receptor channels in these synapses. Here, we investigated the functional properties of AMPA receptors of AII amacrine cells in rat retina that receive input at spatially segregated ribbon synapses from OFF-cone and rod bipolar cells. Using patch-clamp recording from outside-out patches, we measured the concentration dependence of response amplitude and steady-state desensitization, the single-channel conductance and the maximum open probability. The GluA4 subunit seems critical for the functional properties of AMPA receptors in AII amacrines and immunocytochemical labeling suggested that GluA4 is located at synapses made by both OFF-cone bipolar cells and rod bipolar cells. Finally, we used a series of experimental observables to develop kinetic models for AII amacrine AMPA receptors and subsequently used the models to explore the behavior of the receptors and responses generated by glutamate concentration profiles mimicking those occurring in synapses. These models will facilitate future in silico modeling of synaptic signaling and processing in AII amacrine cells.

Hair cell-type dependent expression of basolateral ion channels shapes response dynamics in the frog utricle

The dynamics of vestibular afferent responses are thought to be strongly influenced by presynaptic properties. In this paper, by performing whole-cell perforated-patch experiments in the frog utricle, we characterized voltage-dependent currents and voltage responses to current steps and 0.3–100 Hz sinusoids. Current expression and voltage responses are strongly related to hair cell type. In particular, voltage responses of extrastriolar type eB (low pass, −3 dB corner at 52.5 ± 12.8 Hz) and striolar type F cells (resonant, tuned at 60 ± 46 Hz) agree with the dynamics (tonic and phasic, respectively) of the afferent fibers they contact. On the other hand, hair cell release (measured with single-sine membrane ΔCm measurements) was linearly related to Ca in both cell types, and therefore did not appear to contribute to dynamics differences. As a tool for quantifying the relative contribution of basolateral currents and other presynaptic factors to afferent dynamics, the recorded current, voltage and release data were used to build a NEURON model of the average extrastriolar type eB and striolar type F hair cell. The model contained all recorded conductances, a basic mechanosensitive hair bundle and a ribbon synapse sustained by stochastic voltage-dependent Ca channels, and could reproduce the recorded hair cell voltage responses. Simulated release obtained from eB-type and F-type models display significant differences in dynamics, supporting the idea that basolateral currents are able to contribute to afferent dynamics; however, release in type eB and F cell models does not reproduce tonic and phasic dynamics, mainly because of an excessive phase lag present in both cell types. This suggests the presence in vestibular hair cells of an additional, phase-advancing mechanism, in cascade with voltage modulation.

Ionic mechanisms underlying tonic and phasic firing behaviors in retinal ganglion cells: a model study

In the retina, the firing behaviors that ganglion cells exhibit when exposed to light stimuli are very important due to the significant roles they play in encoding the visual information. However, the detailed mechanisms, especially the intrinsic properties that generate and modulate these firing behaviors is not completely clear yet. In this study, 2 typical firing behaviors—i.e., tonic and phasic activities, which are widely observed in retinal ganglion cells (RGCs)—are investigated. A modified computational model was developed to explore the possible ionic mechanisms that underlie the generation of these 2 firing patterns. Computational results indicate that the generation of tonic and phasic activities may be attributed to the collective actions of 2 kinds of adaptation currents, i.e., an inactivating sodium current and a delayed-rectifier potassium current. The concentration of magnesium ions has crucial but differential effects in the modulation of tonic and phasic firings, when the model neuron is driven by N-methyl-D-aspartate (NMDA) -type synaptic input instead of constant current injections. The proposed model has robust features that account for the ionic mechanisms underlying the tonic and phasic firing behaviors, and it may also be used as a good candidate for modeling some other firing patterns in RGCs.

Light adaptation increases response latency of alpha ganglion cells via a threshold-like nonlinearity

Adaptation is an important process of sensory systems to adjust sensitivity to ensure the appropriate information encoding. Sensitivity and kinetics of retinal ganglion cell (RGC) responses have been studied extensively using a brief flash superimposed on different but steady backgrounds. However, it is still unclear if light adaptation exerts any effect on more complex response properties, such as response nonlinearity. In this study, we found that the latency of spike responses to a repeated flashing spot stimulation increased by 30 ms in the mouse ON α RGCs (An ON-type RGC is excited when a spot is turned on in the center of its receptive field). A single dimming event preceding the test flash on a steady adapting background could also produce similar effect in increasing latency of light responses. A simple computational model with a linear transformation of the light stimulus and a threshold-like nonlinearity could account for the experimental data. Moreover, the strength of the measured nonlinearity and the response latency were affected by the duration of light adaptation. The possible biological processes underlying this nonlinearity were explored. Voltage clamp recording revealed the presence of the increase in latency and threshold-like nonlinearity in the excitatory input of RGCs. However, no comparable nonlinearity was observed in the light responses of the ON cone bipolar cells. We further excluded GABAergic and glycinergic inhibition, N-methyl-D-aspartate receptor rectification and voltage-gated Na(+) channels as potential sources of this nonlinearity by pharmacological experiments. Our results indicate the bipolar cell terminals as the potential site of nonlinearity. Computational modeling constrained by experimental data supports that conclusion and suggests the voltage-sensitive Ca(++) channels and Ca(++)-dependent vesicle release in the bipolar cell terminals as mechanistic basis.

Amperometric resolution of a prespike stammer and evoked phases of fast release from retinal bipolar cells

The neurotransmitter glutamate is used by most neurons in the brain to activate a multitude of different types of glutamate receptors and transporters involved in fast and relatively slower signaling. Synaptic ribbons are large presynaptic structures found in neurons involved in vision, balance, and hearing, which use a large number of glutamate-filled synaptic vesicles to meet their signaling demands. To directly measure synaptic vesicle release events, the ribbon-type presynaptic terminals of goldfish retinal bipolar cells were coaxed to release a false transmitter that could be monitored with amperometry by placing the carbon fiber directly on the larger synaptic terminal. Spontaneous secretion events formed a unimodal charge distribution, but single spike properties were heterogeneous. Larger events rose exponentially without interruption (τ ∼ 30 μs), and smaller events exhibited a stammer in their rising phase that is interpreted as a brief pause in pore dilation, a characteristic commonly associated with large dense core granule fusion pores. These events were entirely Ca(2+)-dependent. Holding the cells at -60 mV halted spontaneous release; and when the voltage was stepped to >-40 mV, secretion ensued. When stepping the voltage to 0 mV, novel kinetic phases of vesicle recruitment were revealed. Approximately 14 vesicles were released per ribbon in two kinetic phases with time constants of 1.5 and 16 ms, which are proposed to represent different primed states within the population of docked vesicles.

Modelling intrinsic electrophysiological properties of ON and OFF retinal ganglion cells

ON and OFF retinal ganglion cells (RGCs) display differences in their intrinsic electrophysiology: OFF cells maintain spontaneous activity in the absence of any input, exhibit subthreshold membrane potential oscillations, rebound excitation and burst firing; ON cells require excitatory input to drive their activity and display none of the aforementioned phenomena. The goal of this study was to identify and characterize ionic currents that explain these intrinsic electrophysiological differences between ON and OFF RGCs. A mathematical model of the electrophysiological properties of ON and OFF RGCs was constructed and validated using published patch-clamp data from isolated intact mouse retina. The model incorporates three ionic currents hypothesized to play a role in generating behaviors that are different between ON and OFF RGCs. These currents are persistent Na( + ), I (NaP), hyperpolarization-activated, I (h), and low voltage activated Ca(2 + ), I (T), currents. Using computer simulations of Hodgkin-Huxley type neuron with a single compartment model we found two distinct sets of I (NaP), I (h), I (T) conductances that correspond to ON and OFF RGCs populations. Simulations indicated that special properties of I (T) explain the differences in intrinsic electrophysiology between ON and OFF RGCs examined here. The modelling shows that the maximum conductance of I (T) is higher in OFF than in ON cells, in agreement with recent experimental data.

The mossy fiber bouton: the "common" or the "unique" synapse?

Synapses are the key elements for signal processing and plasticity in the brain. They are composed of nearly the same structural subelements, an apposition zone including a pre- and postsynaptic density, a cleft and a pool of vesicles. It is, however, their actual composition that determines their different behavior in synaptic transmission and plasticity. Here, we describe and discuss the structural factors underlying the unique functional properties of the hippocampal mossy fiber (MF) synapse. Two membrane specializations, active zones (AZs; transmitter release sites), and puncta adherentia (PA), putative adhesion complexes were found. On average, individual boutons had ∼20 AZs with a mean surface area of 0.1 μm² and a short distance of 0.45 μm between individual AZs. Mossy fiber boutons (MFBs) and their target structures were isolated from each other by astrocytes, but fine glial processes never reached the AZs. Therefore, two structural factors are likely to promote synaptic cross-talk: the short distance and the absence of fine glial processes between individual AZs. Thus, synaptic crosstalk may contribute to the high efficacy of hippocampal MF synapses. On average, an adult bouton contained ∼16,000 synaptic vesicles; ∼600 vesicles were located within 60 nm from the AZ, ∼4000 between 60 nm and 200 nm, and the remaining beyond 200 nm, suggesting large readily releasable, recycling, and reserve pools. Thus, the size of the three pools together with the number and distribution of AZs underlie the unique extent of synaptic efficacy and plasticity of the hippocampal MF synapse.

A biophysical model for modulation frequency encoding in the cochlear nucleus

Encoding of amplitude modulated (AM) acoustical signals is one of the most compelling tasks for the mammalian auditory system: environmental sounds, after being filtered and transduced by the cochlea, become narrowband AM signals. Despite much experimental work dedicated to the comprehension of auditory system extraction and encoding of AM information, the neural mechanisms underlying this remarkable feature are far from being understood (Joris et al., 2004). One of the most accepted theories for this processing is the existence of a periodotopic organization (based on temporal information) across the more studied tonotopic axis (Frisina et al., 1990b). In this work, we will review some recent advances in the study of the mechanisms involved in neural processing of AM sounds, and propose an integrated model that runs from the external ear, through the cochlea and the auditory nerve, up to a sub-circuit of the cochlear nucleus (the first processing unit in the central auditory system). We will show that varying the amount of inhibition in our model we can obtain a range of best modulation frequencies (BMF) in some principal cells of the cochlear nucleus. This could be a basis for a synchronicity based, low-level periodotopic organization.

A computational study on the role of gap junctions and rod Ih conductance in the enhancement of the dynamic range of the retina

Recent works suggest that one of the roles of gap junctions in sensory systems is to enhance their dynamic range by avoiding early saturation in the first processing stages. In this work, we use a minimal conductance-based model of the ON rod pathways in the vertebrate retina to study the effects of electrical synaptic coupling via gap junctions among rods and among AII amacrine cells on the dynamic range of the retina. The model is also used to study the effects of the maximum conductance of rod hyperpolarization activated current I_h on the dynamic range of the retina, allowing a study of the interrelations between this intrinsic membrane parameter with those two retina connectivity characteristics. Our results show that for realistic values of I_h conductance the dynamic range is enhanced by rod-rod coupling, and that AII-AII coupling is less relevant to dynamic range amplification in comparison with receptor coupling. Furthermore, a plot of the retina output response versus input intensity for the optimal parameter configuration is well fitted by a power law with exponent . The results are consistent with predictions of more theoretical works and suggest that the earliest expression of gap junctions along the rod pathways, together with appropriate values of rod I_h conductance, has the highest impact on vertebrate retina dynamic range enhancement.

Mechanisms contributing to tonic release at the cone photoreceptor ribbon synapse

Time-resolved capacitance measurements in combination with fluorescence measurements of internal calcium suggested three kinetic components of release in acutely isolated cone photoreceptors of the tiger salamander. A 45-fF releasable pool, corresponding to about 1,000 vesicles, was identified. This pool could be depleted with a time constant of a few hundred milliseconds and its recovery from depletion was quite rapid (tau approximately 1 s). The fusion of vesicles in this pool was blocked by low-millimolar EGTA. Endocytosis was sufficiently slow that it is likely that refilling of the releasable pool occurred from preformed vesicles. A second, slower component of release (tau(depletion) approximately 3 s) was identified that was approximately twice the size of the releasable pool. This pool may serve as a first reserve pool that replenishes the releasable pool. Computer simulations indicate that the properties of the releasable and first reserve pools are sufficient to maintain synaptic signaling for several seconds in the face of near-maximal stimulations and in the absence of other sources of vesicles. Along with lower rates of depletion, additional mechanisms, such as replenishment from distal reserve pools and the fast recycling of vesicles, may further contribute to the maintenance of graded, tonic release from cone photoreceptors.

The morphology of excitatory central synapses: from structure to function

Synapses are the key elements for signal transduction and plasticity in the brain. For a better understanding of the functional signal cascades underlying synaptic transmission, a quantitative morphological analysis of the pre- and postsynaptic structures that represent morphological correlates for synaptic transmission is important. In particular, realistic values of the number, distribution, and geometry of synaptic contacts and the organization of the pool of synaptic vesicles provide important constraints for realistic models and numerical simulations of those parameters of synaptic transmission that, at present, are still not accessible to experiment. Although all synapses are composed of almost the same structural elements, the composition of these elements within a given synapse and the microcircuit in which they are embedded are the deciding factors determining its function. One possible way to analyze these structures is by computer-assisted three-dimensional reconstructions of synapses and their subsequent quantitative analysis based on ultrathin serial sections. The present review summarizes and discusses the morphology of five central excitatory synapses that are quantitatively well described: (1) a giant synapse, the so-called Calyx of Held, in the medial nucleus of the trapezoid body in the auditory brain stem, (2) the mossy fiber terminal establishing synapses with multiple cerebellar granule cell dendrites, (3) the mossy fiber bouton in the hippocampus predominantly terminating on proximal dendrites of CA3 pyramidal neurons, (4) the climbing fiber-Purkinje cell synapse in the cerebellum, and (5) cortical input synapses on the basal dendrites of layer 5 pyramidal cells. The detailed morphological description of these synaptic structures may help to define the morphological correlates of the functional parameters of synaptic transmission, such as the readily releasable pool of synaptic vesicles, of release, and of the variability of quantal size and might therefore explain the existing differences in the function between individual synapses embedded in different microcircuits.

Form and Function of on-off Amacrine Cells in the Amphibian Retina

on-off amacrine cells were studied with whole cell recording techniques and intracellular staining methods using intact retina-eyecup preparations of the tiger salamander ( Ambystoma tigrinum) and the mudpuppy ( Necturus maculosus). Morphological characterization of these cells included three-dimensional reconstruction methods based on serial optical sections obtained with a confocal microscope. Some cells had their detailed morphology digitized with a computer-assisted tracing system and converted to compartmental models for computer simulations. The dendrites of on-off amacrine cells have spines and numerous varicosities. Physiological recordings confirmed that on-off amacrine cells generate both large- and small-amplitude impulses attributed, respectively, to somatic and dendritic generation sites. Using a multichannel model for impulse generation, computer simulations were carried out to evaluate how impulses are likely to propagate throughout these structures. We conclude that the on-off amacrine cell is organized with multifocal dendritic impulse generating sites and that both dendritic and somatic impulse activity contribute to the functional repertoire of these interneurons: locally generated dendritic impulses can provide regional activation, while somatic impulse activity results in rapid activation of the entire dendritic tree.