Glutamate-bound NMDARs arising from in vivo-like network activity extend spatio-temporal integration in a L5 cortical pyramidal cell model

Glutamate spillover drives robust all-or-none dendritic plateau potentials-an in silico investigation using models of striatal projection neurons

Plateau potentials are a critical feature of neuronal excitability, but their all-or-none behavior is not easily captured in modeling. In this study, we investigated models of plateau potentials in multi-compartment neuron models and found that including glutamate spillover provides robust all-or-none behavior. This result arises due to the prolonged duration of extrasynaptic glutamate. When glutamate spillover is not included, the all-or-none behavior is very sensitive to the steepness of the Mg²⁺ block. These results suggest a potentially significant role of glutamate spillover in plateau potential generation, providing a mechanism for robust all-or-none behavior across a wide range of slopes of the Mg²⁺ block curve. We also illustrate the importance of the all-or-none plateau potential behavior for nonlinear computation with regard to the nonlinear feature binding problem.

Spatial integration of dendrites in fast-spiking basket cells

Dendrites of fast-spiking basket cells (FS BCs) impact neural circuit functions in brain with both supralinear and sublinear integration strategies. Diverse spatial synaptic inputs and active properties of dendrites lead to distinct neuronal firing patterns. How the FS BCs with this bi-modal dendritic integration respond to different spatial dispersion of synaptic inputs remains unclear. In this study, we construct a multi-compartmental model of FS BC and analyze neuronal firings following simulated synaptic protocols from fully clustered to fully dispersed. Under these stimulation protocols, we find that supralinear dendrites dominate somatic firing of FS BC, while the preference for dispersing is due to sublinear dendrites. Moreover, we find that dendritic diameter and Ca2+-permeable AMPA conductance play an important role in it, while A-type K+ channel and NMDA conductance have little effect. The obtained results may give some implications for understanding dendritic computation.

Parvalbumin interneuron dendrites enhance gamma oscillations

Dendrites are essential determinants of the input-output relationship of single neurons, but their role in network computations is not well understood. Here, we use a combination of dendritic patch-clamp recordings and in silico modeling to determine how dendrites of parvalbumin (PV)-expressing basket cells contribute to network oscillations in the gamma frequency band. Simultaneous soma-dendrite recordings from PV basket cells in the dentate gyrus reveal that the slope, or gain, of the dendritic input-output relationship is exceptionally low, thereby reducing the cell's sensitivity to changes in its input. By simulating gamma oscillations in detailed network models, we demonstrate that the low gain is key to increase spike synchrony in PV basket cell assemblies when cells are driven by spatially and temporally heterogeneous synaptic inputs. These results highlight the role of inhibitory neuron dendrites in synchronized network oscillations.

Activity-dependent modulation of NMDA receptors by endogenous zinc shapes dendritic function in cortical neurons

NMDA receptors (NMDARs) have been proposed to control single-neuron computations in vivo. However, whether specific mechanisms regulate the function of such receptors and modulate input-output transformations performed by cortical neurons under in vivo-like conditions is understudied. Here, we report that in layer 2/3 pyramidal neurons (L2/3 PNs), repeated synaptic stimulation results in an activity-dependent decrease in NMDAR function by vesicular zinc. Such a mechanism shifts the threshold for dendritic non-linearities and strongly reduces LTP. Modulation of NMDARs is cell and pathway specific, being present selectively in L2/3-L2/3 connections but absent in inputs originating from L4 neurons. Numerical simulations highlight that activity-dependent modulation of NMDARs influences dendritic computations, endowing L2/3 PN dendrites with the ability to sustain non-linear integrations constant across different regimes of synaptic activity like those found in vivo. Our results unveil vesicular zinc as an important endogenous modulator of dendritic function in cortical PNs.

•

Vesicular zinc release downregulates function of synaptic NMDARs in cortical neurons

•

Zinc modulation of NMDARs is activity dependent, pathway and cell specific

•

Endogenous zinc controls dendritic non-linearities and synaptic plasticity in L2/3 PNs

•

Modulation of NMDARs normalizes dendritic function during ongoing synaptic activity

A general principle of dendritic constancy: A neuron's size- and shape-invariant excitability

Reducing neuronal size results in less membrane and therefore lower input conductance. Smaller neurons are thus more excitable, as seen in their responses to somatic current injections. However, the impact of a neuron's size and shape on its voltage responses to dendritic synaptic activation is much less understood. Here we use analytical cable theory to predict voltage responses to distributed synaptic inputs in unbranched cables, showing that these are entirely independent of dendritic length. For a given synaptic density, neuronal responses depend only on the average dendritic diameter and intrinsic conductivity. This remains valid for a wide range of morphologies irrespective of their arborization complexity. Spiking models indicate that morphology-invariant numbers of spikes approximate the percentage of active synapses. In contrast to spike rate, spike times do depend on dendrite morphology. In summary, neuronal excitability in response to distributed synaptic inputs is largely unaffected by dendrite length or complexity.

Structure, Function, and Pharmacology of Glutamate Receptor Ion Channels

Many physiologic effects of l-glutamate, the major excitatory neurotransmitter in the mammalian central nervous system, are mediated via signaling by ionotropic glutamate receptors (iGluRs). These ligand-gated ion channels are critical to brain function and are centrally implicated in numerous psychiatric and neurologic disorders. There are different classes of iGluRs with a variety of receptor subtypes in each class that play distinct roles in neuronal functions. The diversity in iGluR subtypes, with their unique functional properties and physiologic roles, has motivated a large number of studies. Our understanding of receptor subtypes has advanced considerably since the first iGluR subunit gene was cloned in 1989, and the research focus has expanded to encompass facets of biology that have been recently discovered and to exploit experimental paradigms made possible by technological advances. Here, we review insights from more than 3 decades of iGluR studies with an emphasis on the progress that has occurred in the past decade. We cover structure, function, pharmacology, roles in neurophysiology, and therapeutic implications for all classes of receptors assembled from the subunits encoded by the 18 ionotropic glutamate receptor genes. SIGNIFICANCE STATEMENT: Glutamate receptors play important roles in virtually all aspects of brain function and are either involved in mediating some clinical features of neurological disease or represent a therapeutic target for treatment. Therefore, understanding the structure, function, and pharmacology of this class of receptors will advance our understanding of many aspects of brain function at molecular, cellular, and system levels and provide new opportunities to treat patients.

Active dendrites enable strong but sparse inputs to determine orientation selectivity

The dendrites of neocortical pyramidal neurons are excitable. However, it is unknown how synaptic inputs engage nonlinear dendritic mechanisms during sensory processing in vivo, and how they in turn influence action potential output. Here, we provide a quantitative account of the relationship between synaptic inputs, nonlinear dendritic events, and action potential output. We developed a detailed pyramidal neuron model constrained by in vivo dendritic recordings. We drive this model with realistic input patterns constrained by sensory responses measured in vivo and connectivity measured in vitro. We show mechanistically that under realistic conditions, dendritic Na+ and NMDA spikes are the major determinants of neuronal output in vivo. We demonstrate that these dendritic spikes can be triggered by a surprisingly small number of strong synaptic inputs, in some cases even by single synapses. We predict that dendritic excitability allows the 1% strongest synaptic inputs of a neuron to control the tuning of its output. Active dendrites therefore allow smaller subcircuits consisting of only a few strongly connected neurons to achieve selectivity for specific sensory features.

NMDA Receptors Enhance the Fidelity of Synaptic Integration

Excitatory synaptic transmission in many neurons is mediated by two coexpressed ionotropic glutamate receptor subtypes, AMPA and NMDA receptors, that differ in kinetics, ion selectivity, and voltage-sensitivity. AMPA receptors have fast kinetics and are voltage-insensitive, while NMDA receptors have slower kinetics and increased conductance at depolarized membrane potentials. Here, we report that the voltage dependency and kinetics of NMDA receptors act synergistically to stabilize synaptic integration of EPSPs across spatial and voltage domains. Simulations of synaptic integration in simplified and morphologically realistic dendritic trees revealed that the combined presence of AMPA and NMDA conductances reduce the variability of somatic responses to spatiotemporal patterns of excitatory synaptic input presented at different initial membrane potentials and/or in different dendritic domains. This moderating effect of the NMDA conductance on synaptic integration was robust across a wide range of AMPA-to-NMDA ratios, and results from synergistic interaction of NMDA kinetics (which reduces variability across membrane potential) and voltage dependence (which favors stabilization across dendritic location). When combined with AMPA conductance, the NMDA conductance compensates for voltage-dependent and impedance-dependent changes in synaptic driving force, and distance-dependent attenuation of synaptic potentials arriving at the axon, to increase the fidelity of synaptic integration and EPSP-spike coupling across both neuron state (i.e., initial membrane potential) and dendritic location of synaptic input. Thus, synaptic NMDA receptors convey advantages for synaptic integration that are independent of, but fully compatible with, their importance for coincidence detection and synaptic plasticity.

Modeling Reveals Human-Rodent Differences in H-Current Kinetics Influencing Resonance in Cortical Layer 5 Neurons

While our understanding of human neurons is often inferred from rodent data, inter-species differences between neurons can be captured by building cellular models specifically from human data. This includes understanding differences at the level of ion channels and their implications for human brain function. Thus, we here present a full spiking, biophysically detailed multi-compartment model of a human layer 5 (L5) cortical pyramidal cell. Model development was primarily based on morphological and electrophysiological data from the same human L5 neuron, avoiding confounds of experimental variability. Focus was placed on describing the behavior of the hyperpolarization-activated cation (h-) channel, given increasing interest in this channel due to its role in pacemaking and differentiating cell types. We ensured that the model exhibited post-inhibitory rebound spiking considering its relationship with the h-current, along with other general spiking characteristics. The model was validated against data not used in its development, which highlighted distinctly slower kinetics of the human h-current relative to the rodent setting. We linked the lack of subthreshold resonance observed in human L5 neurons to these human-specific h-current kinetics. This work shows that it is possible and necessary to build human-specific biophysical neuron models in order to understand human brain dynamics.

Illuminating dendritic function with computational models

Dendrites have always fascinated researchers: from the artistic drawings by Ramon y Cajal to the beautiful recordings of today, neuroscientists have been striving to unravel the mysteries of these structures. Theoretical work in the 1960s predicted important dendritic effects on neuronal processing, establishing computational modelling as a powerful technique for their investigation. Since then, modelling of dendrites has been instrumental in driving neuroscience research in a targeted manner, providing experimentally testable predictions that range from the subcellular level to the systems level, and their relevance extends to fields beyond neuroscience, such as machine learning and artificial intelligence. Validation of modelling predictions often requires - and drives - new technological advances, thus closing the loop with theory-driven experimentation that moves the field forward. This Review features the most important, to our understanding, contributions of modelling of dendritic computations, including those pending experimental verification, and highlights studies of successful interactions between the modelling and experimental neuroscience communities.

Electrical Compartmentalization in Neurons

The dendritic tree of neurons plays an important role in information processing in the brain. While it is thought that dendrites require independent subunits to perform most of their computations, it is still not understood how they compartmentalize into functional subunits. Here, we show how these subunits can be deduced from the properties of dendrites. We devised a formalism that links the dendritic arborization to an impedance-based tree graph and show how the topology of this graph reveals independent subunits. This analysis reveals that cooperativity between synapses decreases slowly with increasing electrical separation and thus that few independent subunits coexist. We nevertheless find that balanced inputs or shunting inhibition can modify this topology and increase the number and size of the subunits in a context-dependent manner. We also find that this dynamic recompartmentalization can enable branch-specific learning of stimulus features. Analysis of dendritic patch-clamp recording experiments confirmed our theoretical predictions.

Impact of functional synapse clusters on neuronal response selectivity

Clustering of functionally similar synapses in dendrites is thought to affect neuronal input-output transformation by triggering local nonlinearities. However, neither the in vivo impact of synaptic clusters on somatic membrane potential (sVm), nor the rules of cluster formation are elucidated. We develop a computational approach to measure the effect of functional synaptic clusters on sVm response of biophysical model CA1 and L2/3 pyramidal neurons to in vivo-like inputs. We demonstrate that small synaptic clusters appearing with random connectivity do not influence sVm. With structured connectivity,  ~10-20 synapses/cluster are optimal for clustering-based tuning via state-dependent mechanisms, but larger selectivity is achieved by 2-fold potentiation of the same synapses. We further show that without nonlinear amplification of the effect of random clusters, action potential-based, global plasticity rules cannot generate functional clustering. Our results suggest that clusters likely form via local synaptic interactions, and have to be moderately large to impact sVm responses.

Global and Multiplexed Dendritic Computations under In Vivo-like Conditions

Dendrites integrate inputs nonlinearly, but it is unclear how these nonlinearities contribute to the overall input-output transformation of single neurons. We developed statistically principled methods using a hierarchical cascade of linear-nonlinear subunits (hLN) to model the dynamically evolving somatic response of neurons receiving complex, in vivo-like spatiotemporal synaptic input patterns. We used the hLN to predict the somatic membrane potential of an in vivo-validated detailed biophysical model of a L2/3 pyramidal cell. Linear input integration with a single global dendritic nonlinearity achieved above 90% prediction accuracy. A novel hLN motif, input multiplexing into parallel processing channels, could improve predictions as much as conventionally used additional layers of local nonlinearities. We obtained similar results in two other cell types. This approach provides a data-driven characterization of a key component of cortical circuit computations: the input-output transformation of neurons during in vivo-like conditions.

Inhibition enhances spatially-specific calcium encoding of synaptic input patterns in a biologically constrained model

Synaptic plasticity, which underlies learning and memory, depends on calcium elevation in neurons, but the precise relationship between calcium and spatiotemporal patterns of synaptic inputs is unclear. Here, we develop a biologically realistic computational model of striatal spiny projection neurons with sophisticated calcium dynamics, based on data from rodents of both sexes, to investigate how spatiotemporally clustered and distributed excitatory and inhibitory inputs affect spine calcium. We demonstrate that coordinated excitatory synaptic inputs evoke enhanced calcium elevation specific to stimulated spines, with lower but physiologically relevant calcium elevation in nearby non-stimulated spines. Results further show a novel and important function of inhibition-to enhance the difference in calcium between stimulated and non-stimulated spines. These findings suggest that spine calcium dynamics encode synaptic input patterns and may serve as a signal for both stimulus-specific potentiation and heterosynaptic depression, maintaining balanced activity in a dendritic branch while inducing pattern-specific plasticity.

Human Cortical Pyramidal Neurons: From Spines to Spikes via Models

We present detailed models of pyramidal cells from human neocortex, including models on their excitatory synapses, dendritic spines, dendritic NMDA- and somatic/axonal Na⁺ spikes that provided new insights into signal processing and computational capabilities of these principal cells. Six human layer 2 and layer 3 pyramidal cells (HL2/L3 PCs) were modeled, integrating detailed anatomical and physiological data from both fresh and postmortem tissues from human temporal cortex. The models predicted particularly large AMPA- and NMDA-conductances per synaptic contact (0.88 and 1.31 nS, respectively) and a steep dependence of the NMDA-conductance on voltage. These estimates were based on intracellular recordings from synaptically-connected HL2/L3 pairs, combined with extra-cellular current injections and use of synaptic blockers, and the assumption of five contacts per synaptic connection. A large dataset of high-resolution reconstructed HL2/L3 dendritic spines provided estimates for the EPSPs at the spine head (12.7 ± 4.6 mV), spine base (9.7 ± 5.0 mV), and soma (0.3 ± 0.1 mV), and for the spine neck resistance (50–80 MΩ). Matching the shape and firing pattern of experimental somatic Na⁺-spikes provided estimates for the density of the somatic/axonal excitable membrane ion channels, predicting that 134 ± 28 simultaneously activated HL2/L3-HL2/L3 synapses are required for generating (with 50% probability) a somatic Na⁺ spike. Dendritic NMDA spikes were triggered in the model when 20 ± 10 excitatory spinous synapses were simultaneously activated on individual dendritic branches. The particularly large number of basal dendrites in HL2/L3 PCs and the distinctive cable elongation of their terminals imply that ~25 NMDA-spikes could be generated independently and simultaneously in these cells, as compared to ~14 in L2/3 PCs from the rat somatosensory cortex. These multi-sites non-linear signals, together with the large (~30,000) excitatory synapses/cell, equip human L2/L3 PCs with enhanced computational capabilities. Our study provides the most comprehensive model of any human neuron to-date demonstrating the biophysical and computational distinctiveness of human cortical neurons.

Cell-type-specific inhibition of the dendritic plateau potential in striatal spiny projection neurons

Striatal spiny projection neurons (SPNs) receive convergent excitatory synaptic inputs from the cortex and thalamus. Activation of spatially clustered and temporally synchronized excitatory inputs at the distal dendrites could trigger plateau potentials in SPNs. Such supralinear synaptic integration is crucial for dendritic computation. However, how plateau potentials interact with subsequent excitatory and inhibitory synaptic inputs remains unknown. By combining computational simulation, two-photon imaging, optogenetics, and dual-color uncaging of glutamate and GABA, we demonstrate that plateau potentials can broaden the spatiotemporal window for integrating excitatory inputs and promote spiking. The temporal window of spiking can be delicately controlled by GABAergic inhibition in a cell-type-specific manner. This subtle inhibitory control of plateau potential depends on the location and kinetics of the GABAergic inputs and is achieved by the balance between relief and reestablishment of NMDA receptor Mg²⁺ block. These findings represent a mechanism for controlling spatiotemporal synaptic integration in SPNs.

Basal tree complexity shapes functional pathways in the prefrontal cortex

While the morphology of basal dendritic trees in cortical pyramidal neurons varies, the functional implications of this diversity are just starting to emerge. In layer 5 pyramidal neurons of the prefrontal cortex, for example, increased basal tree complexity determines the recruitment of these neurons into functional circuits. Here, we use a modeling approach to investigate whether and how the morphology of the basal tree mediates the functional output of neurons. We implemented 57 basal tree morphologies of layer 5 prefrontal pyramidal neurons of the rat and identified morphological types that were characterized by different response features, forming distinct functional types. These types were robust to a wide range of manipulations (distribution of active ionic mechanisms, NMDA conductance, somatic and apical tree morphology, or the number of activated synapses) and supported different temporal coding schemes at both the single neuron and the microcircuit level. We predict that the basal tree morphological diversity among neurons of the same class mediates their segregation into distinct functional pathways. Extension of our approach/findings to other cortical areas and/or layers or under pathological conditions may provide a generalized role of the basal trees for neuronal function.NEW & NOTEWORTHY Our results suggest that the segregation of neurons to different functional types based on their basal tree morphology is in large part independent of the distribution of active ionic mechanisms, NMDA conductance, somatic and apical tree morphology, and the number of activated synapses; different functional types support distinct temporal coding schemes. This can be exploited to create networks with diverse coding characteristics, thus contributing to the functional heterogeneity within the same layer and area.

Dendritic Spikes in Sensory Perception

What is the function of dendritic spikes? One might argue that they provide conditions for neuronal plasticity or that they are essential for neural computation. However, despite a long history of dendritic research, the physiological relevance of dendritic spikes in brain function remains unknown. This could stem from the fact that most studies on dendrites have been performed in vitro. Fortunately, the emergence of novel techniques such as improved two-photon microscopy, genetically encoded calcium indicators (GECIs), and optogenetic tools has provided the means for vital breakthroughs in in vivo dendritic research. These technologies enable the investigation of the functions of dendritic spikes in behaving animals, and thus, help uncover the causal relationship between dendritic spikes, and sensory information processing and synaptic plasticity. Understanding the roles of dendritic spikes in brain function would provide mechanistic insight into the relationship between the brain and the mind. In this review article, we summarize the results of studies on dendritic spikes from a historical perspective and discuss the recent advances in our understanding of the role of dendritic spikes in sensory perception.

Nanoscale diffusion in the synaptic cleft and beyond measured with time-resolved fluorescence anisotropy imaging

Neural activity relies on molecular diffusion within nanoscopic spaces outside and inside nerve cells, such as synaptic clefts or dendritic spines. Measuring diffusion on this small scale in situ has not hitherto been possible, yet this knowledge is critical for understanding the dynamics of molecular events and electric currents that shape physiological signals throughout the brain. Here we advance time-resolved fluorescence anisotropy imaging combined with two-photon excitation microscopy to map nanoscale diffusivity in ex vivo brain slices. We find that in the brain interstitial gaps small molecules move on average ~30% slower than in a free medium whereas inside neuronal dendrites this retardation is ~70%. In the synaptic cleft free nanodiffusion is decelerated by ~46%. These quantities provide previously unattainable basic constrains for the receptor actions of released neurotransmitters, the electrical conductance of the brain interstitial space and the limiting rate of molecular interactions or conformational changes in the synaptic microenvironment.

Dendritic integration: 60 years of progress

Understanding how individual neurons integrate the thousands of synaptic inputs they receive is critical to understanding how the brain works. Modeling studies in silico and experimental work in vitro, dating back more than half a century, have revealed that neurons can perform a variety of different passive and active forms of synaptic integration on their inputs. But how are synaptic inputs integrated in the intact brain? With the development of new techniques, this question has recently received substantial attention, with new findings suggesting that many of the forms of synaptic integration observed in vitro also occur in vivo, including in awake animals. Here we review six decades of progress, which collectively highlights the complex ways that single neurons integrate their inputs, emphasizing the critical role of dendrites in information processing in the brain.

Effects of Neural Morphology and Input Distribution on Synaptic Processing by Global and Focal NMDA-Spikes

Cortical neurons can respond to glutamatergic stimulation with regenerative N-Methyl-D-aspartic acid (NMDA)-spikes. NMDA-spikes were initially thought to depend on clustered synaptic activation. Recent work had shown however a new variety of a global NMDA-spike, which can be generated by randomly distributed inputs. Very little is known about the factors that influence the generation of these global NMDA-spikes, as well the potentially distinct rules of synaptic integration and the computational significance conferred by the two types of NMDA-spikes. Here I show that the input resistance (R_IN) plays a major role in influencing spike initiation; while the classical, focal NMDA-spike depended upon the local (dendritic) R_IN, the threshold of global NMDA-spike generation was set by the somatic R_IN. As cellular morphology can exert a large influence on R_IN, morphologically distinct neuron types can have dissimilar rules for NMDA-spikes generation. For example, cortical neurons in superficial layers were found to be generally prone to global NMDA-spike generation. In contrast, electric properties of cortical layer 5b cells clearly favor focal NMDA-spikes. These differences can translate into diverse synaptic integration rules for the different classes of cortical cells; simulated superficial layers neurons were found to exhibit strong synaptic interactions between different dendritic branches, giving rise to a single integrative compartment mediated by the global NMDA-spike. In these cells, efficiency of postsynaptic activation was relatively little dependent on synaptic distribution. By contrast, layer 5b neurons were capable of true multi-unit computation involving independent integrative compartments formed by clustered synaptic input which could trigger focal NMDA-spikes. In a sharp contrast to superficial layers neurons, randomly distributed synaptic inputs were not very effective in driving firing the layer 5b cells, indicating a possibility for different computation performed by these important cortical neurons.

Contribution of sublinear and supralinear dendritic integration to neuronal computations

Nonlinear dendritic integration is thought to increase the computational ability of neurons. Most studies focus on how supralinear summation of excitatory synaptic responses arising from clustered inputs within single dendrites result in the enhancement of neuronal firing, enabling simple computations such as feature detection. Recent reports have shown that sublinear summation is also a prominent dendritic operation, extending the range of subthreshold input-output (sI/O) transformations conferred by dendrites. Like supralinear operations, sublinear dendritic operations also increase the repertoire of neuronal computations, but feature extraction requires different synaptic connectivity strategies for each of these operations. In this article we will review the experimental and theoretical findings describing the biophysical determinants of the three primary classes of dendritic operations: linear, sublinear, and supralinear. We then review a Boolean algebra-based analysis of simplified neuron models, which provides insight into how dendritic operations influence neuronal computations. We highlight how neuronal computations are critically dependent on the interplay of dendritic properties (morphology and voltage-gated channel expression), spiking threshold and distribution of synaptic inputs carrying particular sensory features. Finally, we describe how global (scattered) and local (clustered) integration strategies permit the implementation of similar classes of computations, one example being the object feature binding problem.

LEMS: a language for expressing complex biological models in concise and hierarchical form and its use in underpinning NeuroML 2

Computational models are increasingly important for studying complex neurophysiological systems. As scientific tools, it is essential that such models can be reproduced and critically evaluated by a range of scientists. However, published models are currently implemented using a diverse set of modeling approaches, simulation tools, and computer languages making them inaccessible and difficult to reproduce. Models also typically contain concepts that are tightly linked to domain-specific simulators, or depend on knowledge that is described exclusively in text-based documentation. To address these issues we have developed a compact, hierarchical, XML-based language called LEMS (Low Entropy Model Specification), that can define the structure and dynamics of a wide range of biological models in a fully machine readable format. We describe how LEMS underpins the latest version of NeuroML and show that this framework can define models of ion channels, synapses, neurons and networks. Unit handling, often a source of error when reusing models, is built into the core of the language by specifying physical quantities in models in terms of the base dimensions. We show how LEMS, together with the open source Java and Python based libraries we have developed, facilitates the generation of scripts for multiple neuronal simulators and provides a route for simulator free code generation. We establish that LEMS can be used to define models from systems biology and map them to neuroscience-domain specific simulators, enabling models to be shared between these traditionally separate disciplines. LEMS and NeuroML 2 provide a new, comprehensive framework for defining computational models of neuronal and other biological systems in a machine readable format, making them more reproducible and increasing the transparency and accessibility of their underlying structure and properties.