Hebbian plasticity in parallel synaptic pathways: A circuit mechanism for systems memory consolidation.

Systems memory consolidation involves the transfer of memories across brain regions and the transformation of memory content. For example, declarative memories that transiently depend on the hippocampal formation are transformed into long-term memory traces in neocortical networks, and procedural memories are transformed within cortico-striatal networks. These consolidation processes are thought to rely on replay and repetition of recently acquired memories, but the cellular and network mechanisms that mediate the changes of memories are poorly understood. Here, we suggest that systems memory consolidation could arise from Hebbian plasticity in networks with parallel synaptic pathways-two ubiquitous features of neural circuits in the brain. We explore this hypothesis in the context of hippocampus-dependent memories. Using computational models and mathematical analyses, we illustrate how memories are transferred across circuits and discuss why their representations could change. The analyses suggest that Hebbian plasticity mediates consolidation by transferring a linear approximation of a previously acquired memory into a parallel pathway. Our modelling results are further in quantitative agreement with lesion studies in rodents. Moreover, a hierarchical iteration of the mechanism yields power-law forgetting-as observed in psychophysical studies in humans. The predicted circuit mechanism thus bridges spatial scales from single cells to cortical areas and time scales from milliseconds to years.

Spikelets in pyramidal neurons: generating mechanisms, distinguishing properties, and functional implications.

Spikelets are small spike-like depolarizations that are found in somatic recordings of many neuron types. Spikelets have been assigned important functions, ranging from neuronal synchronization to the regulation of synaptic plasticity, which are specific to the particular mechanism of spikelet generation. As spikelets reflect spiking activity in neuronal compartments that are electrotonically distinct from the soma, four modes of spikelet generation can be envisaged: (1) dendritic spikes or (2) axonal action potentials occurring in a single cell as well as action potentials transmitted via (3) gap junctions or (4) ephaptic coupling in pairs of neurons. In one of the best studied neuron type, cortical pyramidal neurons, the origins and functions of spikelets are still unresolved; all four potential mechanisms have been proposed, but the experimental evidence remains ambiguous. Here we attempt to reconcile the scattered experimental findings in a coherent theoretical framework. We review in detail the various mechanisms that can give rise to spikelets. For each mechanism, we present the biophysical underpinnings as well as the resulting properties of spikelets and compare these predictions to experimental data from pyramidal neurons. We also discuss the functional implications of each mechanism. On the example of pyramidal neurons, we illustrate that several independent spikelet-generating mechanisms fulfilling vastly different functions might be operating in a single cell.

Soma-axon coupling configurations that enhance neuronal coincidence detection.

Coincidence detector neurons transmit timing information by responding preferentially to concurrent synaptic inputs. Principal cells of the medial superior olive (MSO) in the mammalian auditory brainstem are superb coincidence detectors. They encode sound source location with high temporal precision, distinguishing submillisecond timing differences among inputs. We investigate computationally how dynamic coupling between the input region (soma and dendrite) and the spike-generating output region (axon and axon initial segment) can enhance coincidence detection in MSO neurons. To do this, we formulate a two-compartment neuron model and characterize extensively coincidence detection sensitivity throughout a parameter space of coupling configurations. We focus on the interaction between coupling configuration and two currents that provide dynamic, voltage-gated, negative feedback in subthreshold voltage range: sodium current with rapid inactivation and low-threshold potassium current, IKLT. These currents reduce synaptic summation and can prevent spike generation unless inputs arrive with near simultaneity. We show that strong soma-to-axon coupling promotes the negative feedback effects of sodium inactivation and is, therefore, advantageous for coincidence detection. Furthermore, the feedforward combination of strong soma-to-axon coupling and weak axon-to-soma coupling enables spikes to be generated efficiently (few sodium channels needed) and with rapid recovery that enhances high-frequency coincidence detection. These observations detail the functional benefit of the strongly feedforward configuration that has been observed in physiological studies of MSO neurons. We find that IKLT further enhances coincidence detection sensitivity, but with effects that depend on coupling configuration. For instance, in models with weak soma-to-axon and weak axon-to-soma coupling, IKLT in the axon enhances coincidence detection more effectively than IKLT in the soma. By using a minimal model of soma-to-axon coupling, we connect structure, dynamics, and computation. Although we consider the particular case of MSO coincidence detectors, our method for creating and exploring a parameter space of two-compartment models can be applied to other neurons.

Function and energy consumption constrain neuronal biophysics in a canonical computation: Coincidence detection.

Neural morphology and membrane properties vary greatly between cell types in the nervous system. The computations and local circuit connectivity that neurons support are thought to be the key factors constraining the cells' biophysical properties. Nevertheless, additional constraints can be expected to further shape neuronal design. Here, we focus on a particularly energy-intense system (as indicated by metabolic markers): principal neurons in the medial superior olive (MSO) nucleus of the auditory brainstem. Based on a modeling approach, we show that a trade-off between the level of performance of a functionally relevant computation and energy consumption predicts optimal ranges for cell morphology and membrane properties. The biophysical parameters appear most strongly constrained by functional needs, while energy use is minimized as long as function can be maintained. The key factors that determine model performance and energy consumption are 1) the saturation of the synaptic conductance input and 2) the temporal resolution of the postsynaptic signals as they reach the soma, which is largely determined by active membrane properties. MSO cells seem to operate close to pareto optimality, i.e., the trade-off boundary between performance and energy consumption that is formed by the set of optimal models. Good performance for drastically lower costs could in theory be achieved by small neurons without dendrites, as seen in the avian auditory system, pointing to additional constraints for mammalian MSO cells, including their circuit connectivity.

Dendrite-Specific Amplification of Weak Synaptic Input during Network Activity In Vivo.

Excitatory synaptic input reaches the soma of a cortical excitatory pyramidal neuron via anatomically segregated apical and basal dendrites. In vivo, dendritic inputs are integrated during depolarized network activity, but how network activity affects apical and basal inputs is not understood. Using subcellular two-photon stimulation of Channelrhodopsin2-expressing layer 2/3 pyramidal neurons in somatosensory cortex, nucleus-specific thalamic optogenetic stimulation, and paired recordings, we show that slow, depolarized network activity amplifies small-amplitude synaptic inputs targeted to basal dendrites but reduces the amplitude of all inputs from apical dendrites and the cell soma. Intracellular pharmacology and mathematical modeling suggests that the amplification of weak basal inputs is mediated by postsynaptic voltage-gated channels. Thus, network activity dynamically reconfigures the relative somatic contribution of apical and basal inputs and could act to enhance the detectability of weak synaptic inputs.

Extracellular waveforms reveal an axonal origin of spikelets in pyramidal neurons.

Spikelets are small spike-like membrane depolarizations measured at the soma whose origin in pyramidal neurons is still unresolved. We investigated the mechanism of spikelet generation using detailed models of pyramidal neurons. We simulated extracellular waveforms associated with action potentials and spikelets and compared these with experimental data obtained by Chorev and Brecht ( J Neurophysiol 108: 1584-1593, 2012) from hippocampal pyramidal neurons in vivo. We considered spikelets originating in the axon of a single cell as well as spikelets generated in two cells coupled with gap junctions. We found that in both cases the experimental data can be explained by an axonal origin of spikelets: in the single-cell case, action potentials are generated in the axon but fail to activate the soma. Such spikelets can be evoked by dendritic input. Alternatively, spikelets resulting from axoaxonal gap junction coupling with a large (greater than several hundred µm) distance between the somata of the coupled cells are also consistent with the data. Our results demonstrate that a cell firing a somatic spikelet generates a detectable extracellular potential that is different from the action potential-correlated extracellular waveform generated by the same cell and recorded at the same location. This, together with the absence of a refractory period between action potentials and spikelets, implies that spikelets and action potentials generated in one cell may easily get misclassified in extracellular recordings as two different cells, albeit they both constitute the output of a single pyramidal neuron. NEW & NOTEWORTHY We addressed the origin of spikelets, using compartmental models of pyramidal neurons. Comparing our simulation results with published extracellular spikelet recordings revealed an axonal origin of spikelets. Our results imply that action potential- and spikelet-associated extracellular waveforms may easily get misclassified as two different cells, albeit they both constitute the output of a single pyramidal cell.

h-Type Membrane Current Shapes the Local Field Potential from Populations of Pyramidal Neurons.

In cortex, the local field potential (LFP) is thought to mainly stem from correlated synaptic input to populations of geometrically aligned neurons. Computer models of single cortical pyramidal neurons showed that subthreshold voltage-dependent membrane conductances can also shape the LFP signal, in particular the hyperpolarization-activated cation current (Ih; h-type). This ion channel is prominent in various types of pyramidal neurons, typically showing an increasing density gradient along the apical dendrites. Here, we investigate how Ih affects the LFP generated by a model of a population of cortical pyramidal neurons. We find that the LFP from populations of neurons that receive uncorrelated synaptic input can be well predicted by the LFP from single neurons. In this case, when input impinges on the distal dendrites, where most h-type channels are located, a strong resonance in the LFP was measured near the soma, whereas the opposite configuration does not reveal an Ih contribution to the LFP. Introducing correlations in the synaptic inputs to the pyramidal cells strongly amplifies the LFP, while maintaining the differential effects of Ih for distal dendritic versus perisomatic input. Previous theoretical work showed that input correlations do not amplify LFP power when neurons receive synaptic input uniformly across the cell. We find that this crucially depends on the membrane conductance distribution: the asymmetric distribution of Ih results in a strong amplification of the LFP when synaptic inputs to the cell population are correlated. In conclusion, we find that the h-type current is particularly suited to shape the LFP signal in cortical populations.SIGNIFICANCE STATEMENT The local field potential (LFP), the low-frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. While the cortical LFP is thought to mainly reflect synaptic inputs onto pyramidal neurons, little is known about the role of subthreshold active conductances in shaping the LFP. By means of biophysical modeling we obtain a comprehensive, qualitative understanding of how LFPs generated by populations of cortical pyramidal neurons depend on active subthreshold currents, and identify the key importance of the h-type channel. Our results show that LFPs can give information about the active properties of neurons and that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics.

Differential polarization of cortical pyramidal neuron dendrites through weak extracellular fields.

The rise of transcranial current stimulation (tCS) techniques have sparked an increasing interest in the effects of weak extracellular electric fields on neural activity. These fields modulate ongoing neural activity through polarization of the neuronal membrane. While the somatic polarization has been investigated experimentally, the frequency-dependent polarization of the dendritic trees in the presence of alternating (AC) fields has received little attention yet. Using a biophysically detailed model with experimentally constrained active conductances, we analyze the subthreshold response of cortical pyramidal cells to weak AC fields, as induced during tCS. We observe a strong frequency resonance around 10-20 Hz in the apical dendrites sensitivity to polarize in response to electric fields but not in the basal dendrites nor the soma. To disentangle the relative roles of the cell morphology and active and passive membrane properties in this resonance, we perform a thorough analysis using simplified models, e.g. a passive pyramidal neuron model, simple passive cables and reconstructed cell model with simplified ion channels. We attribute the origin of the resonance in the apical dendrites to (i) a locally increased sensitivity due to the morphology and to (ii) the high density of h-type channels. Our systematic study provides an improved understanding of the subthreshold response of cortical cells to weak electric fields and, importantly, allows for an improved design of tCS stimuli.

Spikelets in Pyramidal Neurons: Action Potentials Initiated in the Axon Initial Segment That Do Not Activate the Soma.

Spikelets are small spike-like depolarizations that can be measured in somatic intracellular recordings. Their origin in pyramidal neurons remains controversial. To explain spikelet generation, we propose a novel single-cell mechanism: somato-dendritic input generates action potentials at the axon initial segment that may fail to activate the soma and manifest as somatic spikelets. Using mathematical analysis and numerical simulations of compartmental neuron models, we identified four key factors controlling spikelet generation: (1) difference in firing threshold, (2) impedance mismatch, and (3) electrotonic separation between the soma and the axon initial segment, as well as (4) input amplitude. Because spikelets involve forward propagation of action potentials along the axon while they avoid full depolarization of the somato-dendritic compartments, we conjecture that this mode of operation saves energy and regulates dendritic plasticity while still allowing for a read-out of results of neuronal computations.

Active subthreshold dendritic conductances shape the local field potential.

KEY POINTS: The local field potential (LFP), the low-frequency part of extracellular potentials recorded in neural tissue, is often used for probing neural circuit activity. Interpreting the LFP signal is difficult, however. While the cortical LFP is thought mainly to reflect synaptic inputs onto pyramidal neurons, little is known about the role of the various subthreshold active conductances in shaping the LFP. By means of biophysical modelling we obtain a comprehensive qualitative understanding of how the LFP generated by a single pyramidal neuron depends on the type and spatial distribution of active subthreshold currents. For pyramidal neurons, the h-type channels probably play a key role and can cause a distinct resonance in the LFP power spectrum. Our results show that the LFP signal can give information about the active properties of neurons and imply that preferred frequencies in the LFP can result from those cellular properties instead of, for example, network dynamics. ABSTRACT: The main contribution to the local field potential (LFP) is thought to stem from synaptic input to neurons and the ensuing subthreshold dendritic processing. The role of active dendritic conductances in shaping the LFP has received little attention, even though such ion channels are known to affect the subthreshold neuron dynamics. Here we used a modelling approach to investigate the effects of subthreshold dendritic conductances on the LFP. Using a biophysically detailed, experimentally constrained model of a cortical pyramidal neuron, we identified conditions under which subthreshold active conductances are a major factor in shaping the LFP. We found that, in particular, the hyperpolarization-activated inward current, Ih , can have a sizable effect and cause a resonance in the LFP power spectral density. To get a general, qualitative understanding of how any subthreshold active dendritic conductance and its cellular distribution can affect the LFP, we next performed a systematic study with a simplified model. We found that the effect on the LFP is most pronounced when (1) the synaptic drive to the cell is asymmetrically distributed (i.e. either basal or apical), (2) the active conductances are distributed non-uniformly with the highest channel densities near the synaptic input and (3) when the LFP is measured at the opposite pole of the cell relative to the synaptic input. In summary, we show that subthreshold active conductances can be strongly reflected in LFP signals, opening up the possibility that the LFP can be used to characterize the properties and cellular distributions of active conductances.

Subthreshold resonance properties contribute to the efficient coding of auditory spatial cues.

Neurons in the medial superior olive (MSO) and lateral superior olive (LSO) of the auditory brainstem code for sound-source location in the horizontal plane, extracting interaural time differences (ITDs) from the stimulus fine structure and interaural level differences (ILDs) from the stimulus envelope. Here, we demonstrate a postsynaptic gradient in temporal processing properties across the presumed tonotopic axis; neurons in the MSO and the low-frequency limb of the LSO exhibit fast intrinsic electrical resonances and low input impedances, consistent with their processing of ITDs in the temporal fine structure. Neurons in the high-frequency limb of the LSO show low-pass electrical properties, indicating they are better suited to extracting information from the slower, modulated envelopes of sounds. Using a modeling approach, we assess ITD and ILD sensitivity of the neural filters to natural sounds, demonstrating that the transformation in temporal processing along the tonotopic axis contributes to efficient extraction of auditory spatial cues.

Somatic versus dendritic resonance: differential filtering of inputs through non-uniform distributions of active conductances.

Synaptic inputs to neurons are processed in a frequency-dependent manner, with either low-pass or resonant response characteristics. These types of filtering play a key role in the frequency-specific information flow in neuronal networks. While the generation of resonance by specific ionic conductances is well investigated, less attention has been paid to the spatial distribution of the resonance-generating conductances across a neuron. In pyramidal neurons - one of the major excitatory cell-types in the mammalian brain - a steep gradient of resonance-generating h-conductances with a 60-fold increase towards distal dendrites has been demonstrated experimentally. Because the dendritic trees of these cells are large, spatial compartmentalization of resonant properties can be expected. Here, we use mathematical descriptions of spatially extended neurons to investigate the consequences of such a distal, dendritic localization of h-conductances for signal processing. While neurons with short dendrites do not exhibit a pronounced compartmentalization of resonance, i.e. the filter properties of dendrites and soma are similar, we find that neurons with longer dendrites ([Formula: see text] space constant) can show distinct filtering of dendritic and somatic inputs due to electrotonic segregation. Moreover, we show that for such neurons, experimental classification as resonant versus nonresonant can be misleading when based on somatic recordings, because for these morphologies a dendritic resonance could easily be undetectable when using somatic input. Nevertheless, noise-driven membrane-potential oscillations caused by dendritic resonance can propagate to the soma where they can be recorded, hence contrasting with the low-pass filtering at the soma. We conclude that non-uniform distributions of active conductances can underlie differential filtering of synaptic input in neurons with spatially extended dendrites, like pyramidal neurons, bearing relevance for the localization-dependent targeting of synaptic input pathways to these cells.

Homeostatic scaling of excitability in recurrent neural networks.

Neurons adjust their intrinsic excitability when experiencing a persistent change in synaptic drive. This process can prevent neural activity from moving into either a quiescent state or a saturated state in the face of ongoing plasticity, and is thought to promote stability of the network in which neurons reside. However, most neurons are embedded in recurrent networks, which require a delicate balance between excitation and inhibition to maintain network stability. This balance could be disrupted when neurons independently adjust their intrinsic excitability. Here, we study the functioning of activity-dependent homeostatic scaling of intrinsic excitability (HSE) in a recurrent neural network. Using both simulations of a recurrent network consisting of excitatory and inhibitory neurons that implement HSE, and a mean-field description of adapting excitatory and inhibitory populations, we show that the stability of such adapting networks critically depends on the relationship between the adaptation time scales of both neuron populations. In a stable adapting network, HSE can keep all neurons functioning within their dynamic range, while the network is undergoing several (patho)physiologically relevant types of plasticity, such as persistent changes in external drive, changes in connection strengths, or the loss of inhibitory cells from the network. However, HSE cannot prevent the unstable network dynamics that result when, due to such plasticity, recurrent excitation in the network becomes too strong compared to feedback inhibition. This suggests that keeping a neural network in a stable and functional state requires the coordination of distinct homeostatic mechanisms that operate not only by adjusting neural excitability, but also by controlling network connectivity.

Dynamic interaction of Ih and IK-LVA during trains of synaptic potentials in principal neurons of the medial superior olive.

In neurons of the medial superior olive (MSO), voltage-gated ion channels control the submillisecond time resolution of binaural coincidence detection, but little is known about their interplay during trains of synaptic activity that would be experienced during auditory stimuli. Here, using modeling and patch-clamp recordings from MSO principal neurons in gerbil brainstem slices, we examined interactions between two major currents controlling subthreshold synaptic integration: a low-voltage-activated potassium current (I(K-LVA)) and a hyperpolarization-activated cation current (I(h)). Both I(h) and I(K-LVA) contributed strongly to the resting membrane conductance and, during trains of simulated EPSPs, exhibited cumulative deactivation and inactivation, respectively. In current-clamp recordings, regular and irregular trains of simulated EPSCs increased input resistance up to 60%, effects that accumulated and decayed (after train) over hundreds of milliseconds. Surprisingly, the mean voltage and peaks of EPSPs increased by only a few millivolts during trains. Using a model of an MSO cell, we demonstrated that the nearly uniform response during modest depolarizing stimuli relied on changes in I(h) and I(K-LVA), such that their sum remained nearly constant over time. Experiments and modeling showed that, for simplified binaural stimuli (EPSC pairs in a noisy background), spike probability gradually increased in parallel with the increasing input resistance. Nevertheless, the interplay between I(h) and I(K-LVA) helps to maintain a nearly uniform shape of individual synaptic responses, and we show that the time resolution of synaptic coincidence detection can be maintained during trains if EPSC size gradually decreases (as in synaptic depression), counteracting slow increases in excitability.

Role of active dendritic conductances in subthreshold input integration.

Dendrites of many types of neurons contain voltage-dependent conductances that are active at subthreshold membrane potentials. To understand the computations neurons perform it is key to understand the role of active dendrites in the subthreshold processing of synaptic inputs. We examine systematically how active dendritic conductances affect the time course of postsynaptic potentials propagating along dendrites, and how they affect the interaction between such signals. Voltage-dependent currents can be classified into two types that have qualitatively different effects on subthreshold input responses: regenerative dendritic currents boost and broaden EPSPs, while restorative currents attenuate and narrow EPSPs. Importantly, the effects of active dendritic currents on EPSP shape increase as the EPSP travels along the dendrite. The effectiveness of active currents in modulating the EPSP shape is determined by their activation time constant: the faster it is, the stronger the effect on EPSP amplitude, while the largest effects on EPSP width occur when it is comparable to the membrane time constant. We finally demonstrate that the two current types can differentially improve precision and robustness of neural computations: restorative currents enhance coincidence detection of dendritic inputs, whereas direction selectivity to sequences of dendritic inputs is enhanced by regenerative dendritic currents.

Democracy-independence trade-off in oscillating dendrites and its implications for grid cells.

Dendritic democracy and independence have been characterized for near-instantaneous processing of synaptic inputs. However, a wide class of neuronal computations requires input integration on long timescales. As a paradigmatic example, entorhinal grid fields have been thought to be generated by the democratic summation of independent dendritic oscillations performing direction-selective path integration. We analyzed how multiple dendritic oscillators embedded in the same neuron integrate inputs separately and determine somatic membrane voltage jointly. We found that the interaction of dendritic oscillations leads to phase locking, which sets an upper limit on the timescale for independent input integration. Factors that increase this timescale also decrease the influence that the dendritic oscillations exert on somatic voltage. In entorhinal stellate cells, interdendritic coupling dominates and causes these cells to act as single oscillators. Our results suggest a fundamental trade-off between local and global processing in dendritic trees integrating ongoing signals.

The role of ongoing dendritic oscillations in single-neuron dynamics.

The dendritic tree contributes significantly to the elementary computations a neuron performs while converting its synaptic inputs into action potential output. Traditionally, these computations have been characterized as both temporally and spatially localized. Under this localist account, neurons compute near-instantaneous mappings from their current input to their current output, brought about by somatic summation of dendritic contributions that are generated in functionally segregated compartments. However, recent evidence about the presence of oscillations in dendrites suggests a qualitatively different mode of operation: the instantaneous phase of such oscillations can depend on a long history of inputs, and under appropriate conditions, even dendritic oscillators that are remote may interact through synchronization. Here, we develop a mathematical framework to analyze the interactions of local dendritic oscillations and the way these interactions influence single cell computations. Combining weakly coupled oscillator methods with cable theoretic arguments, we derive phase-locking states for multiple oscillating dendritic compartments. We characterize how the phase-locking properties depend on key parameters of the oscillating dendrite: the electrotonic properties of the (active) dendritic segment, and the intrinsic properties of the dendritic oscillators. As a direct consequence, we show how input to the dendrites can modulate phase-locking behavior and hence global dendritic coherence. In turn, dendritic coherence is able to gate the integration and propagation of synaptic signals to the soma, ultimately leading to an effective control of somatic spike generation. Our results suggest that dendritic oscillations enable the dendritic tree to operate on more global temporal and spatial scales than previously thought; notably that local dendritic activity may be a mechanism for generating on-going whole-cell voltage oscillations.

Different levels of Ih determine distinct temporal integration in bursting and regular-spiking neurons in rat subiculum.

Pyramidal neurons in the subiculum typically display either bursting or regular-spiking behaviour. Although this classification into two neuronal classes is well described, it is unknown how these two classes of neurons contribute to the integration of input to the subiculum. Here, we report that bursting neurons possess a hyperpolarization-activated cation current (I(h)) that is two-fold larger (conductance, 5.3 +/- 0.5 nS) than in regular-spiking neurons (2.2 +/- 0.6 nS), whereas I(h) exhibits similar voltage-dependent and kinetic properties in both classes of neurons. Bursting and regular-spiking neurons display similar morphology. The difference in I(h) between the two classes of neurons is not responsible for the distinct firing patterns, as neither pharmacological blockade of I(h) nor enhancement of I(h) using a dynamic clamp affects the qualitative firing patterns. Instead, the difference in I(h) between bursting and regular-spiking neurons determines the temporal integration of evoked synaptic input from the CA1 area. In response to stimulation at 50 Hz, bursting neurons, with a large I(h), show approximately 50% less temporal summation than regular-spiking neurons. The amount of temporal summation in both neuronal classes is equal after pharmacological blockade of I(h). A computer simulation model of a subicular neuron with the properties of either a bursting or a regular-spiking neuron confirmed the pivotal role of I(h) in temporal integration of synaptic input. These data suggest that in the subicular network, bursting neurons are better suited to discriminate the content of high-frequency input, such as that occurring during gamma oscillations, than regular-spiking neurons.

The effect of dendritic topology on firing patterns in model neurons.

Neuronal firing patterns are influenced by both membrane properties and dendritic morphology. Distinguishing two sources of morphological variability-metrics and topology-we investigate the extent to which model neurons that have the same metrical and membrane properties can still produce different firing patterns as a result of differences in dendritic topology. Within a set of dendritic trees that have the same number of terminal segments and the same total dendritic length, we show that firing frequency strongly correlates with topology as expressed by the mean dendritic path length. The effect of dendritic topology on firing frequency is bigger for trees with equal segment diameters than for trees whose segment diameters obey Rall's 3/2 power law. If active dendritic channels are present, dendritic topology influences not only firing frequency but also type of firing (regular, bursting).