Regulation of AMPA and NMDA receptor-mediated EPSPs in dendritic trees of thalamocortical cells.

Two main excitatory synapses are formed at the dendritic arbor of first-order nuclei thalamocortical (TC) neurons. Ascending sensory axons primarily establish contacts at large proximal dendrites, whereas descending corticothalamic fibers form synapses on thin distal dendrites. With the use of a multicomparment computational model based on fully reconstructed TC neurons from the ventroposterolateral nucleus of the cat, we compared local responses at the site of stimulation as well as somatic responses induced by both α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR)- and N-methyl-D-aspartate receptor (NMDAR)-mediated currents. We found that AMPAR-mediated responses, when synapses were located at proximal dendrites, induced a larger depolarization at the level of soma, whereas NMDAR-mediated responses were more efficient for synapses located at distal dendrites. The voltage transfer and transfer impedance were higher for NMDAR than for AMPAR activation at any location. For both types of synaptic current and for both input locations at the dendritic arbor, somatic responses were characterized by a low variability despite the large variability found in local responses in dendrites. The large neurons had overall smaller somatic responses than small neurons, but this relation was not found in local dendritic responses. We conclude that in TC cells, the dendritic location of small synaptic inputs does not play a major role in the amplitude of a somatic response, but the size of the neuron does. The variability of response amplitude between cells was much larger than the variability within cells. This suggests possible functional segregation of TC neurons of different size.

Efficacy of synaptic inhibition depends on multiple, dynamically interacting mechanisms implicated in chloride homeostasis.

Chloride homeostasis is a critical determinant of the strength and robustness of inhibition mediated by GABA(A) receptors (GABA(A)Rs). The impact of changes in steady state Cl(-) gradient is relatively straightforward to understand, but how dynamic interplay between Cl(-) influx, diffusion, extrusion and interaction with other ion species affects synaptic signaling remains uncertain. Here we used electrodiffusion modeling to investigate the nonlinear interactions between these processes. Results demonstrate that diffusion is crucial for redistributing intracellular Cl(-) load on a fast time scale, whereas Cl(-)extrusion controls steady state levels. Interaction between diffusion and extrusion can result in a somato-dendritic Cl(-) gradient even when KCC2 is distributed uniformly across the cell. Reducing KCC2 activity led to decreased efficacy of GABA(A)R-mediated inhibition, but increasing GABA(A)R input failed to fully compensate for this form of disinhibition because of activity-dependent accumulation of Cl(-). Furthermore, if spiking persisted despite the presence of GABA(A)R input, Cl(-) accumulation became accelerated because of the large Cl(-) driving force that occurs during spikes. The resulting positive feedback loop caused catastrophic failure of inhibition. Simulations also revealed other feedback loops, such as competition between Cl(-) and pH regulation. Several model predictions were tested and confirmed by [Cl(-)](i) imaging experiments. Our study has thus uncovered how Cl(-) regulation depends on a multiplicity of dynamically interacting mechanisms. Furthermore, the model revealed that enhancing KCC2 activity beyond normal levels did not negatively impact firing frequency or cause overt extracellular K(-) accumulation, demonstrating that enhancing KCC2 activity is a valid strategy for therapeutic intervention.

Analysis of morphological features of thalamocortical neurons from the ventroposterolateral nucleus of the cat.

Morphological features of the dendritic arborization can affect neuronal responses and thus the input-output function of a particular neuron. In this study, morphological data of eight fully reconstructed thalamocortical (TC) neurons from the ventroposterolateral (VPL) nucleus of adult cats have been analyzed. We examined several geometrical and topological parameters, which have been previously shown to have a high impact on the neuron firing pattern and propagation of signals in the dendritic tree. In addition to well-known morphological parameters such as number of dendritic trees (8.3 +/- 1.5) and number of branching points (80-120), we investigated the distribution of dendritic membrane area, branching points, geometrical ratio, asymmetry index, and mean path length for all subtrees of the TC neurons. We demonstrate that due to extensive branching in proximal and middle dendritic sections, the maximum value of the dendritic area distribution is reached at 120-160 mum from the soma. Our analysis reveals that TC neurons are highly branched cells and their dendritic branching pattern does not follow Rall's 3/2 power rule; average values at proximal vs. distal dendritic sections were different. We also found that the dendritic branching pattern of each subtree of the cell had a wide range in symmetry index, whereas the mean path length did not show a large variation through the dendritic arborizations.

Realistic modelling of receptor activation in hippocampal excitatory synapses: analysis of multivesicular release, release location, temperature and synaptic cross-talk.

Chemically mediated synaptic transmission results from fusion of synaptic vesicles with the presynaptic plasma membrane, subsequent release of the vesicular content into the cleft and binding to postsynaptic receptors. Previous modelling studies of excitatory neurotransmitter glutamate were based on simplified geometries failing to account for the biologically realistic synaptic environment, in particular, the presence of astrocytes, the geometry of extracellular space, and the neurotransmitter uptake mechanism. Using 3-dimensional reconstructions of hippocampal glutamatergic synapses including the surrounding astrocytic processes we have developed a biologically realistic model to analyse receptor activation in different conditions. We used the finite element method to simulate glutamate release, analyse glutamate diffusion following single and multiple vesicle release and binding at the postsynaptic site to AMPA and NMDA receptors. We demonstrate that: (1) the transmitter diffusion is highly temperature-sensitive; (2) release conditions and geometry more specifically affect AMPARs than NMDARs; (3) the sensitivities of AMPARs and NMDARs to simultaneous vesicular release are different; (4) in the case of multivesicle neurotransmitter release with variable delays, the binding of glutamate to AMPARs is additive up to 1 ms after the release, then becomes independent, but to NMDARs the binding is additive up to 33 ms; (5) the number of AMPARs varies more than the number of NMDRs in response to the input firing patterns; (6) the presence of astrocytes effectively blocks synaptic cross-talk; and (7) synaptic cross-talk, mediated by NMDARs but not AMPARs, is only possible after quasi-simultaneous multivesicular release at physiological temperature (35 degrees C) without intervening astrocytes, but not at 25 degrees C. Our simulations demonstrate the importance of temperature and ultrastructural synaptic environment in synaptic transmission and synaptic cross-talk.

Modeling thalamocortical cell: impact of ca channel distribution and cell geometry on firing pattern.

The influence of calcium channel distribution and geometry of the thalamocortical cell upon its tonic firing and the low threshold spike (LTS) generation was studied in a 3-compartment model, which represents soma, proximal and distal dendrites as well as in multi-compartment model using the morphology of a real reconstructed neuron. Using an uniform distribution of Ca(2+) channels, we determined the minimal number of low threshold voltage-activated calcium channels and their permeability required for the onset of LTS in response to a hyperpolarizing current pulse. In the 3-compartment model, we found that the channel distribution influences the firing pattern only in the range of 3% below the threshold value of total T-channel density. In the multi-compartmental model, the LTS could be generated by only 64% of unequally distributed T-channels compared to the minimal number of equally distributed T-channels. For a given channel density and injected current, the tonic firing frequency was found to be inversely proportional to the size of the cell. However, when the Ca(2+) channel density was elevated in soma or proximal dendrites, then the amplitude of LTS response and burst spike frequencies were determined by the ratio of total to threshold number of T-channels in the cell for a specific geometry.

Modeling extracellular field potentials and the frequency-filtering properties of extracellular space.

Extracellular local field potentials are usually modeled as arising from a set of current sources embedded in a homogeneous extracellular medium. Although this formalism can successfully model several properties of extracellular local field potentials, it does not account for their frequency-dependent attenuation with distance, a property essential to correctly model extracellular spikes. Here we derive expressions for the extracellular potential that include this frequency-dependent attenuation. We first show that, if the extracellular conductivity is nonhomogeneous, there is induction of nonhomogeneous charge densities that may result in a low-pass filter. We next derive a simplified model consisting of a punctual (or spherical) current source with spherically symmetric conductivity/permittivity gradients around the source. We analyze the effect of different radial profiles of conductivity and permittivity on the frequency-filtering behavior of this model. We show that this simple model generally displays low-pass filtering behavior, in which fast electrical events (such as Na(+)-mediated action potentials) attenuate very steeply with distance, whereas slower (K(+)-mediated) events propagate over larger distances in extracellular space, in qualitative agreement with experimental observations. This simple model can be used to obtain frequency-dependent extracellular field potentials without taking into account explicitly the complex folding of extracellular space.