Large-Scale Mechanistic Models of Brain Circuits with Biophysically and Morphologically Detailed Neurons.

Understanding the brain requires studying its multiscale interactions from molecules to networks. The increasing availability of large-scale datasets detailing brain circuit composition, connectivity, and activity is transforming neuroscience. However, integrating and interpreting this data remains challenging. Concurrently, advances in supercomputing and sophisticated modeling tools now enable the development of highly detailed, large-scale biophysical circuit models. These mechanistic multiscale models offer a method to systematically integrate experimental data, facilitating investigations into brain structure, function, and disease. This review, based on a Society for Neuroscience 2024 MiniSymposium, aims to disseminate recent advances in large-scale mechanistic modeling to the broader community. It highlights (1) examples of current models for various brain regions developed through experimental data integration; (2) their predictive capabilities regarding cellular and circuit mechanisms underlying experimental recordings (e.g., membrane voltage, spikes, local-field potential, electroencephalography/magnetoencephalography) and brain function; and (3) their use in simulating biomarkers for brain diseases like epilepsy, depression, schizophrenia, and Parkinson's, aiding in understanding their biophysical underpinnings and developing novel treatments. The review showcases state-of-the-art models covering hippocampus, somatosensory, visual, motor, auditory cortical, and thalamic circuits across species. These models predict neural activity at multiple scales and provide insights into the biophysical mechanisms underlying sensation, motor behavior, brain signals, neural coding, disease, pharmacological interventions, and neural stimulation. Collaboration with experimental neuroscientists and clinicians is essential for the development and validation of these models, particularly as datasets grow. Hence, this review aims to foster interest in detailed brain circuit models, leading to cross-disciplinary collaborations that accelerate brain research.

Theta-gamma phase amplitude coupling in a hippocampal CA1 microcircuit.

Phase amplitude coupling (PAC) between slow and fast oscillations is found throughout the brain and plays important functional roles. Its neural origin remains unclear. Experimental findings are often puzzling and sometimes contradictory. Most computational models rely on pairs of pacemaker neurons or neural populations tuned at different frequencies to produce PAC. Here, using a data-driven model of a hippocampal microcircuit, we demonstrate that PAC can naturally emerge from a single feedback mechanism involving an inhibitory and excitatory neuron population, which interplay to generate theta frequency periodic bursts of higher frequency gamma. The model suggests the conditions under which a CA1 microcircuit can operate to elicit theta-gamma PAC, and highlights the modulatory role of OLM and PVBC cells, recurrent connectivity, and short term synaptic plasticity. Surprisingly, the results suggest the experimentally testable prediction that the generation of the slow population oscillation requires the fast one and cannot occur without it.

An Adaptive Generalized Leaky Integrate-and-Fire Model for Hippocampal CA1 Pyramidal Neurons and Interneurons.

Full-scale morphologically and biophysically realistic model networks, aiming at modeling multiple brain areas, provide an invaluable tool to make significant scientific advances from in-silico experiments on cognitive functions to digital twin implementations. Due to the current technical limitations of supercomputer systems in terms of computational power and memory requirements, these networks must be implemented using (at least) simplified neurons. A class of models which achieve a reasonable compromise between accuracy and computational efficiency is given by generalized leaky integrate-and fire models complemented by suitable initial and update conditions. However, we found that these models cannot reproduce the complex and highly variable firing dynamics exhibited by neurons in several brain regions, such as the hippocampus. In this work, we propose an adaptive generalized leaky integrate-and-fire model for hippocampal CA1 neurons and interneurons, in which the nonlinear nature of the firing dynamics is successfully reproduced by linear ordinary differential equations equipped with nonlinear and more realistic initial and update conditions after each spike event, which strictly depends on the external stimulation current. A mathematical analysis of the equilibria stability as well as the monotonicity properties of the analytical solution for the membrane potential allowed (i) to determine general constraints on model parameters, reducing the computational cost of an optimization procedure based on spike times in response to a set of constant currents injections; (ii) to identify additional constraints to quantitatively reproduce and predict the experimental traces from 85 neurons and interneurons in response to any stimulation protocol using constant and piecewise constant current injections. Finally, this approach allows to easily implement a procedure to create infinite copies of neurons with mathematically controlled firing properties, statistically indistinguishable from experiments, to better reproduce the full range and variability of the firing scenarios observed in a real network.

NeuronAlg: An Innovative Neuronal Computational Model for Immunofluorescence Image Segmentation.

Background: Image analysis applications in digital pathology include various methods for segmenting regions of interest. Their identification is one of the most complex steps and therefore of great interest for the study of robust methods that do not necessarily rely on a machine learning (ML) approach. Method: A fully automatic and optimized segmentation process for different datasets is a prerequisite for classifying and diagnosing indirect immunofluorescence (IIF) raw data. This study describes a deterministic computational neuroscience approach for identifying cells and nuclei. It is very different from the conventional neural network approaches but has an equivalent quantitative and qualitative performance, and it is also robust against adversative noise. The method is robust, based on formally correct functions, and does not suffer from having to be tuned on specific data sets. Results: This work demonstrates the robustness of the method against variability of parameters, such as image size, mode, and signal-to-noise ratio. We validated the method on three datasets (Neuroblastoma, NucleusSegData, and ISBI 2009 Dataset) using images annotated by independent medical doctors. Conclusions: The definition of deterministic and formally correct methods, from a functional and structural point of view, guarantees the achievement of optimized and functionally correct results. The excellent performance of our deterministic method (NeuronalAlg) in segmenting cells and nuclei from fluorescence images was measured with quantitative indicators and compared with those achieved by three published ML approaches.

Membrane electrical properties of mouse hippocampal CA1 pyramidal neurons during strong inputs.

In this work, we highlight an electrophysiological feature often observed in recordings from mouse CA1 pyramidal cells that has so far been ignored by experimentalists and modelers. It consists of a large and dynamic increase in the depolarization baseline (i.e., the minimum value of the membrane potential between successive action potentials during a sustained input) in response to strong somatic current injections. Such an increase can directly affect neurotransmitter release properties and, more generally, the efficacy of synaptic transmission. However, it cannot be explained by any currently available conductance-based computational model. Here we present a model addressing this issue, demonstrating that experimental recordings can be reproduced by assuming that an input current modifies, in a time-dependent manner, the electrical and permeability properties of the neuron membrane by shifting the ionic reversal potentials and channel kinetics. For this reason, we propose that any detailed model of ion channel kinetics for neurons exhibiting this characteristic should be adapted to correctly represent the response and the synaptic integration process during strong and sustained inputs.

Reconstruction of the Hippocampus.

The hippocampus is a widely studied brain region thought to play an important role in higher cognitive functions such as learning, memory, and navigation. The amount of data on this region increases every day and delineates a complex and fragmented picture, but an integrated understanding of hippocampal function remains elusive. Computational methods can help to move the research forward, and reconstructing a full-scale model of the hippocampus is a challenging yet feasible task that the research community should undertake.In this chapter, we present strategies for reconstructing a large-scale model of the hippocampus. Based on a previously published approach to reconstruct and simulate brain tissue, which is also explained in Chap. 10 , we discuss the characteristics of the hippocampus in the light of its special anatomical and physiological features, data availability, and existing large-scale hippocampus models. A large-scale model of the hippocampus is a compound model of several building blocks: ion channels, morphologies, single cell models, connections, synapses. We discuss each of those building blocks separately and discuss how to merge them back and simulate the resulting network model.

Effects of Ih and TASK-like shunting current on dendritic impedance in layer 5 pyramidal-tract neurons.

Pyramidal neurons in neocortex have complex input-output relationships that depend on their morphologies, ion channel distributions, and the nature of their inputs, but which cannot be replicated by simple integrate-and-fire models. The impedance properties of their dendritic arbors, such as resonance and phase shift, shape neuronal responses to synaptic inputs and provide intraneuronal functional maps reflecting their intrinsic dynamics and excitability. Experimental studies of dendritic impedance have shown that neocortical pyramidal tract neurons exhibit distance-dependent changes in resonance and impedance phase with respect to the soma. We, therefore, investigated how well several biophysically detailed multicompartment models of neocortical layer 5 pyramidal tract neurons reproduce the location-dependent impedance profiles observed experimentally. Each model tested here exhibited location-dependent impedance profiles, but most captured either the observed impedance amplitude or phase, not both. The only model that captured features from both incorporates hyperpolarization-activated cyclic nucleotide-gated (HCN) channels and a shunting current, such as that produced by Twik-related acid-sensitive K+ (TASK) channels. TASK-like channel density in this model was proportional to local HCN channel density. We found that although this shunting current alone is insufficient to produce resonance or realistic phase response, it modulates all features of dendritic impedance, including resonance frequencies, resonance strength, synchronous frequencies, and total inductive phase. We also explored how the interaction of HCN channel current (Ih) and a TASK-like shunting current shape synaptic potentials and produce degeneracy in dendritic impedance profiles, wherein different combinations of Ih and shunting current can produce the same impedance profile.NEW & NOTEWORTHY We simulated chirp current stimulation in the apical dendrites of 5 biophysically detailed multicompartment models of neocortical pyramidal tract neurons and found that a combination of HCN channels and TASK-like channels produced the best fit to experimental measurements of dendritic impedance. We then explored how HCN and TASK-like channels can shape the dendritic impedance as well as the voltage response to synaptic currents.

Predicting brain organization with a computational model: 50-year perspective on lateral inhibition and oscillatory gating by dendrodendritic synapses.

The first compartmental computer models of brain neurons using the Rall method predicted novel and unexpected dendrodendritic interactions between mitral and granule cells in the olfactory bulb. We review the models from a 50-year perspective on the work that has challenged, supported, and extended the original proposal that these interactions mediate both lateral inhibition and oscillatory activity, essential steps in the neural basis of olfactory processing and perception. We highlight strategies behind the neurophysiological experiments and the Rall methods that enhance the ability of detailed compartmental modeling to give counterintuitive predictions that lead to deeper insights into neural organization at the synaptic and circuit level. The application of these methods to mechanisms of neurogenesis and plasticity are exciting challenges for the future.

Data-driven integration of hippocampal CA1 synaptic physiology in silico.

The anatomy and physiology of monosynaptic connections in rodent hippocampal CA1 have been extensively studied in recent decades. Yet, the resulting knowledge remains disparate and difficult to reconcile. Here, we present a data-driven approach to integrate the current state-of-the-art knowledge on the synaptic anatomy and physiology of rodent hippocampal CA1, including axo-dendritic innervation patterns, number of synapses per connection, quantal conductances, neurotransmitter release probability, and short-term plasticity into a single coherent resource. First, we undertook an extensive literature review of paired recordings of hippocampal neurons and compiled experimental data on their synaptic anatomy and physiology. The data collected in this manner is sparse and inhomogeneous due to the diversity of experimental techniques used by different groups, which necessitates the need for an integrative framework to unify these data. To this end, we extended a previously developed workflow for the neocortex to constrain a unifying in silico reconstruction of the synaptic physiology of CA1 connections. Our work identifies gaps in the existing knowledge and provides a complementary resource toward a more complete quantification of synaptic anatomy and physiology in the rodent hippocampal CA1 region.

A kinetic model for Brain-Derived Neurotrophic Factor mediated spike timing-dependent LTP.

Across the mammalian nervous system, neurotrophins control synaptic plasticity, neuromodulation, and neuronal growth. The neurotrophin Brain-Derived Neurotrophic Factor (BDNF) is known to promote structural and functional synaptic plasticity in the hippocampus, the cerebral cortex, and many other brain areas. In recent years, a wealth of data has been accumulated revealing the paramount importance of BDNF for neuronal function. BDNF signaling gives rise to multiple complex signaling pathways that mediate neuronal survival and differentiation during development, and formation of new memories. These different roles of BDNF for neuronal function have essential consequences if BDNF signaling in the brain is reduced. Thus, BDNF knock-out mice or mice that are deficient in BDNF receptor signaling via TrkB and p75 receptors show deficits in neuronal development, synaptic plasticity, and memory formation. Accordingly, BDNF signaling dysfunctions are associated with many neurological and neurodegenerative conditions including Alzheimer's and Huntington's disease. However, despite the widespread implications of BDNF-dependent signaling in synaptic plasticity in healthy and pathological conditions, the interplay of the involved different biochemical pathways at the synaptic level remained mostly unknown. In this paper, we investigated the role of BDNF/TrkB signaling in spike-timing dependent plasticity (STDP) in rodent hippocampus CA1 pyramidal cells, by implementing the first subcellular model of BDNF regulated, spike timing-dependent long-term potentiation (t-LTP). The model is based on previously published experimental findings on STDP and accounts for the observed magnitude, time course, stimulation pattern and BDNF-dependence of t-LTP. It allows interpreting the main experimental findings concerning specific biomolecular processes, and it can be expanded to take into account more detailed biochemical reactions. The results point out a few predictions on how to enhance LTP induction in such a way to rescue or improve cognitive functions under pathological conditions.

The physiological variability of channel density in hippocampal CA1 pyramidal cells and interneurons explored using a unified data-driven modeling workflow.

Every neuron is part of a network, exerting its function by transforming multiple spatiotemporal synaptic input patterns into a single spiking output. This function is specified by the particular shape and passive electrical properties of the neuronal membrane, and the composition and spatial distribution of ion channels across its processes. For a variety of physiological or pathological reasons, the intrinsic input/output function may change during a neuron's lifetime. This process results in high variability in the peak specific conductance of ion channels in individual neurons. The mechanisms responsible for this variability are not well understood, although there are clear indications from experiments and modeling that degeneracy and correlation among multiple channels may be involved. Here, we studied this issue in biophysical models of hippocampal CA1 pyramidal neurons and interneurons. Using a unified data-driven simulation workflow and starting from a set of experimental recordings and morphological reconstructions obtained from rats, we built and analyzed several ensembles of morphologically and biophysically accurate single cell models with intrinsic electrophysiological properties consistent with experimental findings. The results suggest that the set of conductances expressed in any given hippocampal neuron may be considered as belonging to two groups: one subset is responsible for the major characteristics of the firing behavior in each population and the other is responsible for a robust degeneracy. Analysis of the model neurons suggests several experimentally testable predictions related to the combination and relative proportion of the different conductances that should be expressed on the membrane of different types of neurons for them to fulfill their role in the hippocampus circuitry.

Synaptic clusters function as odor operators in the olfactory bulb.

How the olfactory bulb organizes and processes odor inputs through fundamental operations of its microcircuits is largely unknown. To gain new insight we focus on odor-activated synaptic clusters related to individual glomeruli, which we call glomerular units. Using a 3D model of mitral and granule cell interactions supported by experimental findings, combined with a matrix-based representation of glomerular operations, we identify the mechanisms for forming one or more glomerular units in response to a given odor, how and to what extent the glomerular units interfere or interact with each other during learning, their computational role within the olfactory bulb microcircuit, and how their actions can be formalized into a theoretical framework in which the olfactory bulb can be considered to contain "odor operators" unique to each individual. The results provide new and specific theoretical and experimentally testable predictions.

Optimizing computer models of corticospinal neurons to replicate in vitro dynamics.

Corticospinal neurons (SPI), thick-tufted pyramidal neurons in motor cortex layer 5B that project caudally via the medullary pyramids, display distinct class-specific electrophysiological properties in vitro: strong sag with hyperpolarization, lack of adaptation, and a nearly linear frequency-current (F-I) relationship. We used our electrophysiological data to produce a pair of large archives of SPI neuron computer models in two model classes: 1) detailed models with full reconstruction; and 2) simplified models with six compartments. We used a PRAXIS and an evolutionary multiobjective optimization (EMO) in sequence to determine ion channel conductances. EMO selected good models from each of the two model classes to form the two model archives. Archived models showed tradeoffs across fitness functions. For example, parameters that produced excellent F-I fit produced a less-optimal fit for interspike voltage trajectory. Because of these tradeoffs, there was no single best model but rather models that would be best for particular usages for either single neuron or network explorations. Further exploration of exemplar models with strong F-I fit demonstrated that both the detailed and simple models produced excellent matches to the experimental data. Although dendritic ion identities and densities cannot yet be fully determined experimentally, we explored the consequences of a demonstrated proximal to distal density gradient of Ih, demonstrating that this would lead to a gradient of resonance properties with increased resonant frequencies more distally. We suggest that this dynamical feature could serve to make the cell particularly responsive to major frequency bands that differ by cortical layer. NEW & NOTEWORTHY: We developed models of motor cortex corticospinal neurons that replicate in vitro dynamics, including hyperpolarization-induced sag and realistic firing patterns. Models demonstrated resonance in response to synaptic stimulation, with resonance frequency increasing in apical dendrites with increasing distance from soma, matching the increasing oscillation frequencies spanning deep to superficial cortical layers. This gradient may enable specific corticospinal neuron dendrites to entrain to relevant oscillations in different cortical layers, contributing to appropriate motor output commands.

The possible consequences for cognitive functions of external electric fields at power line frequency on hippocampal CA1 pyramidal neurons.

The possible effects on cognitive processes of external electric fields, such as those generated by power line pillars and household appliances are of increasing public concern. They are difficult to study experimentally, and the relatively scarce and contradictory evidence make it difficult to clearly assess these effects. In this study, we investigate how, why and to what extent external perturbations of the intrinsic neuronal activity, such as those that can be caused by generation, transmission and use of electrical energy can affect neuronal activity during cognitive processes. For this purpose, we used a morphologically and biophysically realistic three-dimensional model of CA1 pyramidal neurons. The simulation findings suggest that an electric field oscillating at power lines frequency, and environmentally measured strength, can significantly alter both the average firing rate and temporal spike distribution properties of a hippocampal CA1 pyramidal neuron. This effect strongly depends on the specific and instantaneous relative spatial location of the neuron with respect to the field, and on the synaptic input properties. The model makes experimentally testable predictions on the possible functional consequences for normal hippocampal functions such as object recognition and spatial navigation. The results suggest that, although EF effects on cognitive processes may be difficult to occur in everyday life, their functional consequences deserve some consideration, especially when they constitute a systematic presence in living environments.

Twenty years of ModelDB and beyond: building essential modeling tools for the future of neuroscience.

Neuron modeling may be said to have originated with the Hodgkin and Huxley action potential model in 1952 and Rall's models of integrative activity of dendrites in 1964. Over the ensuing decades, these approaches have led to a massive development of increasingly accurate and complex data-based models of neurons and neuronal circuits. ModelDB was founded in 1996 to support this new field and enhance the scientific credibility and utility of computational neuroscience models by providing a convenient venue for sharing them. It has grown to include over 1100 published models covering more than 130 research topics. It is actively curated and developed to help researchers discover and understand models of interest. ModelDB also provides mechanisms to assist running models both locally and remotely, and has a graphical tool that enables users to explore the anatomical and biophysical properties that are represented in a model. Each of its capabilities is undergoing continued refinement and improvement in response to user experience. Large research groups (Allen Brain Institute, EU Human Brain Project, etc.) are emerging that collect data across multiple scales and integrate that data into many complex models, presenting new challenges of scale. We end by predicting a future for neuroscience increasingly fueled by new technology and high performance computation, and increasingly in need of comprehensive user-friendly databases such as ModelDB to provide the means to integrate the data for deeper insights into brain function in health and disease.

Early-onset epileptic encephalopathy caused by gain-of-function mutations in the voltage sensor of Kv7.2 and Kv7.3 potassium channel subunits.

Mutations in Kv7.2 (KCNQ2) and Kv7.3 (KCNQ3) genes, encoding for voltage-gated K(+) channel subunits underlying the neuronal M-current, have been associated with a wide spectrum of early-onset epileptic disorders ranging from benign familial neonatal seizures to severe epileptic encephalopathies. The aim of the present work has been to investigate the molecular mechanisms of channel dysfunction caused by voltage-sensing domain mutations in Kv7.2 (R144Q, R201C, and R201H) or Kv7.3 (R230C) recently found in patients with epileptic encephalopathies and/or intellectual disability. Electrophysiological studies in mammalian cells transfected with human Kv7.2 and/or Kv7.3 cDNAs revealed that each of these four mutations stabilized the activated state of the channel, thereby producing gain-of-function effects, which are opposite to the loss-of-function effects produced by previously found mutations. Multistate structural modeling revealed that the R201 residue in Kv7.2, corresponding to R230 in Kv7.3, stabilized the resting and nearby voltage-sensing domain states by forming an intricate network of electrostatic interactions with neighboring negatively charged residues, a result also confirmed by disulfide trapping experiments. Using a realistic model of a feedforward inhibitory microcircuit in the hippocampal CA1 region, an increased excitability of pyramidal neurons was found upon incorporation of the experimentally defined parameters for mutant M-current, suggesting that changes in network interactions rather than in intrinsic cell properties may be responsible for the neuronal hyperexcitability by these gain-of-function mutations. Together, the present results suggest that gain-of-function mutations in Kv7.2/3 currents may cause human epilepsy with a severe clinical course, thus revealing a previously unexplored level of complexity in disease pathogenetic mechanisms.

Genotype-phenotype correlations in neonatal epilepsies caused by mutations in the voltage sensor of K(v)7.2 potassium channel subunits.

Mutations in the K(V)7.2 gene encoding for voltage-dependent K(+) channel subunits cause neonatal epilepsies with wide phenotypic heterogeneity. Two mutations affecting the same positively charged residue in the S4 domain of K(V)7.2 have been found in children affected with benign familial neonatal seizures (R213W mutation) or with neonatal epileptic encephalopathy with severe pharmacoresistant seizures and neurocognitive delay, suppression-burst pattern at EEG, and distinct neuroradiological features (R213Q mutation). To examine the molecular basis for this strikingly different phenotype, we studied the functional characteristics of mutant channels by using electrophysiological techniques, computational modeling, and homology modeling. Functional studies revealed that, in homomeric or heteromeric configuration with K(V)7.2 and/or K(V)7.3 subunits, both mutations markedly destabilized the open state, causing a dramatic decrease in channel voltage sensitivity. These functional changes were (i) more pronounced for channels incorporating R213Q- than R213W-carrying K(V)7.2 subunits; (ii) proportional to the number of mutant subunits incorporated; and (iii) fully restored by the neuronal K(v)7 activator retigabine. Homology modeling confirmed a critical role for the R213 residue in stabilizing the activated voltage sensor configuration. Modeling experiments in CA1 hippocampal pyramidal cells revealed that both mutations increased cell firing frequency, with the R213Q mutation prompting more dramatic functional changes compared with the R213W mutation. These results suggest that the clinical disease severity may be related to the extent of the mutation-induced functional K(+) channel impairment, and set the preclinical basis for the potential use of K(v)7 openers as a targeted anticonvulsant therapy to improve developmental outcome in neonates with K(V)7.2 encephalopathy.

Effect of the initial synaptic state on the probability to induce long-term potentiation and depression.

Long-term potentiation (LTP) and long-term depression (LTD) are the two major forms of long-lasting synaptic plasticity in the mammalian neurons, and are directly related to higher brain functions such as learning and memory. Experimentally, they are characterized by a change in the strength of a synaptic connection induced by repetitive and properly patterned stimulation protocols. Although many important details of the molecular events leading to LTP and LTD are known, experimenters often report problems in using standard induction protocols to obtain consistent results, especially for LTD in vivo. We hypothesize that a possible source of confusion in interpreting the results, from any given experiment on synaptic plasticity, can be the intrinsic limitation of the experimental techniques, which cannot take into account the actual state and peak conductance of the synapses before the conditioning protocol. In this article, we investigate the possibility that the same experimental protocol may result in different consequences (e.g., LTD instead of LTP), according to the initial conditions of the stimulated synapses, and can generate confusing results. Using biophysical models of synaptic plasticity and hippocampal CA1 pyramidal neurons, we study how, why, and to what extent the phenomena observed at the soma after induction of LTP/LTD reflects the actual (local) synaptic state. The model and the results suggest a physiologically plausible explanation for why LTD induction is experimentally difficult to obtain. They also suggest experimentally testable predictions on the stimulation protocols that may be more effective.

Sparse coding and lateral inhibition arising from balanced and unbalanced dendrodendritic excitation and inhibition.

The precise mechanism by which synaptic excitation and inhibition interact with each other in odor coding through the unique dendrodendritic synaptic microcircuits present in olfactory bulb is unknown. Here a scaled-up model of the mitral-granule cell network in the rodent olfactory bulb is used to analyze dendrodendritic processing of experimentally determined odor patterns. We found that the interaction between excitation and inhibition is responsible for two fundamental computational mechanisms: (1) a balanced excitation/inhibition in strongly activated mitral cells, leading to a sparse representation of odorant input, and (2) an unbalanced excitation/inhibition (inhibition dominated) in surrounding weakly activated mitral cells, leading to lateral inhibition. These results suggest how both mechanisms can carry information about the input patterns, with optimal level of synaptic excitation and inhibition producing the highest level of sparseness and decorrelation in the network response. The results suggest how the learning process, through the emergent development of these mechanisms, can enhance odor representation of olfactory bulb.

Functional impact of dendritic branch-point morphology.

Cortical pyramidal cells store multiple features of complex synaptic input in individual dendritic branches and independently regulate the coupling between dendritic and somatic spikes. Branch points in apical trees exhibit wide ranges of sizes and shapes, and the large diameter ratio between trunk and oblique dendrites exacerbates impedance mismatch. The morphological diversity of dendritic bifurcations could thus locally tune neuronal excitability and signal integration. However, these aspects have never been investigated. Here, we first quantified the morphological variability of branch points from two-photon images of rat CA1 pyramidal neurons. We then investigated the geometrical features affecting spike initiation, propagation, and timing with a computational model validated by glutamate uncaging experiments. The results suggest that even subtle membrane readjustments at branch points could drastically alter the ability of synaptic input to generate, propagate, and time action potentials.