Building a mathematical model of the brain.

Automatic leveraging of information in a hippocampal neuron database to generate mathematical models should help foster interactions between experimental and computational neuroscientists.

Hyperpolarization-Activated Cation Channels Shape the Spiking Frequency Preference of Human Cortical Layer 5 Pyramidal Neurons.

Discerning the contribution of specific ionic currents to complex neuronal dynamics is a difficult, but important, task. This challenge is exacerbated in the human setting, although the widely characterized uniqueness of the human brain compared with preclinical models necessitates the direct study of human neurons. Neuronal spiking frequency preference is of particular interest given its role in rhythm generation and signal transmission in cortical circuits. Here, we combine the frequency-dependent gain (FDG), a measure of spiking frequency preference, and novel in silico analyses to dissect the contributions of individual ionic currents to the suprathreshold features of human layer 5 (L5) neurons captured by the FDG. We confirm that a contemporary model of such a neuron, primarily constrained to capture subthreshold activity driven by the hyperpolarization-activated cyclic nucleotide gated (h-) current, replicates key features of the in vitro FDG both with and without h-current activity. With the model confirmed as a viable approximation of the biophysical features of interest, we applied new analysis techniques to quantify the activity of each modeled ionic current in the moments before spiking, revealing unique dynamics of the h-current. These findings motivated patch-clamp recordings in analogous rodent neurons to characterize their FDG, which confirmed that a biophysically detailed model of these neurons captures key interspecies differences in the FDG. These differences are correlated with distinct contributions of the h-current to neuronal activity. Together, this interdisciplinary and multispecies study provides new insights directly relating the dynamics of the h-current to suprathreshold spiking frequency preference in human L5 neurons.

Adaptive unscented Kalman filter for neuronal state and parameter estimation.

Data assimilation techniques for state and parameter estimation are frequently applied in the context of computational neuroscience. In this work, we show how an adaptive variant of the unscented Kalman filter (UKF) performs on the tracking of a conductance-based neuron model. Unlike standard recursive filter implementations, the robust adaptive unscented Kalman filter (RAUKF) jointly estimates the states and parameters of the neuronal model while adjusting noise covariance matrices online based on innovation and residual information. We benchmark the adaptive filter's performance against existing nonlinear Kalman filters and explore the sensitivity of the filter parameters to the system being modelled. To evaluate the robustness of the proposed solution, we simulate practical settings that challenge tracking performance, such as a model mismatch and measurement faults. Compared to standard variants of the Kalman filter the adaptive variant implemented here is more accurate and robust to faults.

Reduced oriens-lacunosum/moleculare cell model identifies biophysical current balances for in vivo theta frequency spiking resonance.

Conductance-based models have played an important role in the development of modern neuroscience. These mathematical models are powerful "tools" that enable theoretical explorations in experimentally untenable situations, and can lead to the development of novel hypotheses and predictions. With advances in cell imaging and computational power, multi-compartment models with morphological accuracy are becoming common practice. However, as more biological details are added, they make extensive explorations and analyses more challenging largely due to their huge computational expense. Here, we focus on oriens-lacunosum/moleculare (OLM) cell models. OLM cells can contribute to functionally relevant theta rhythms in the hippocampus by virtue of their ability to express spiking resonance at theta frequencies, but what characteristics underlie this is far from clear. We converted a previously developed detailed multi-compartment OLM cell model into a reduced single compartment model that retained biophysical fidelity with its underlying ion currents. We showed that the reduced OLM cell model can capture complex output that includes spiking resonance in in vivo-like scenarios as previously obtained with the multi-compartment model. Using the reduced model, we were able to greatly expand our in vivo-like scenarios. Applying spike-triggered average analyses, we were able to to determine that it is a combination of hyperpolarization-activated cation and muscarinic type potassium currents that specifically allow OLM cells to exhibit spiking resonance at theta frequencies. Further, we developed a robust Kalman Filtering (KF) method to estimate parameters of the reduced model in real-time. We showed that it may be possible to directly estimate conductance parameters from experiments since this KF method can reliably extract parameter values from model voltage recordings. Overall, our work showcases how the contribution of cellular biophysical current details could be determined and assessed for spiking resonance. As well, our work shows that it may be possible to directly extract these parameters from current clamp voltage recordings.

Tackling recalcitrant Pseudomonas aeruginosa infections in critical illness via anti-virulence monotherapy.

Intestinal barrier derangement allows intestinal bacteria and their products to translocate to the systemic circulation. Pseudomonas aeruginosa (PA) superimposed infection in critically ill patients increases gut permeability and leads to gut-driven sepsis. PA infections are challenging due to multi-drug resistance (MDR), biofilms, and/or antibiotic tolerance. Inhibition of the quorum-sensing transcriptional regulator MvfR(PqsR) is a desirable anti-PA anti-virulence strategy as MvfR controls multiple acute and chronic virulence functions. Here we show that MvfR promotes intestinal permeability and report potent anti-MvfR compounds, the N-Aryl Malonamides (NAMs), resulting from extensive structure-activity-relationship studies and thorough assessment of the inhibition of MvfR-controlled virulence functions. This class of anti-virulence non-native ligand-based agents has a half-maximal inhibitory concentration in the nanomolar range and strong target engagement. Using a NAM lead in monotherapy protects murine intestinal barrier function, abolishes MvfR-regulated small molecules, ameliorates bacterial dissemination, and lowers inflammatory cytokines. This study demonstrates the importance of MvfR in PA-driven intestinal permeability. It underscores the utility of anti-MvfR agents in maintaining gut mucosal integrity, which should be part of any successful strategy to prevent/treat PA infections and associated gut-derived sepsis in critical illness settings. NAMs provide for the development of crucial preventive/therapeutic monotherapy options against untreatable MDR PA infections.

Corrigendum: Inhibitory Network Bistability Explains Increased Interneuronal Activity Prior to Seizure Onset.

[This corrects the article DOI: 10.3389/fncir.2019.00081.].

Deciphering how interneuron specific 3 cells control oriens lacunosum-moleculare cells to contribute to circuit function.

The wide diversity of inhibitory cells across the brain makes them suitable to contribute to network dynamics in specialized fashions. However, the contributions of a particular inhibitory cell type in a behaving animal are challenging to untangle as one needs to both record cellular activities and identify the cell type being recorded. Thus, using computational modeling and theory to predict and hypothesize cell-specific contributions is desirable. Here, we examine potential contributions of interneuron-specific 3 (I-S3) cells-an inhibitory interneuron found in CA1 hippocampus that only targets other inhibitory interneurons-during simulated θ rhythms. We use previously developed multicompartment models of oriens lacunosum-moleculare (OLM) cells, the main target of I-S3 cells, and explore how I-S3 cell inputs during in vitro and in vivo scenarios contribute to θ. We find that I-S3 cells suppress OLM cell spiking, rather than engender its spiking via postinhibitory rebound mechanisms, and contribute to θ frequency spike resonance during simulated in vivo scenarios. To elicit recruitment similar to in vitro experiments, inclusion of disinhibited pyramidal cell inputs is necessary, implying that I-S3 cell firing broadens the window for pyramidal cell disinhibition. Using in vivo virtual networks, we show that I-S3 cells contribute to a sharpening of OLM cell recruitment at θ frequencies. Furthermore, shifting the timing of I-S3 cell spiking due to external modulation shifts the timing of the OLM cell firing and thus disinhibitory windows. We propose a specialized contribution of I-S3 cells to create temporally precise coordination of modulation pathways.NEW & NOTEWORTHY How information is processed across different brain structures is an important question that relates to the different functions that the brain performs. Using modeling and theoretical analyses, we show that an inhibitory cell type that only inhibits other inhibitory cells can broaden the window for disinhibition of excitatory cells, manage input pathway switching, and modulate inhibitory cell spiking. This work contributes to the knowledge of how coordination between sensory and memory consolidation information can be attained.

Linking minimal and detailed models of CA1 microcircuits reveals how theta rhythms emerge and their frequencies controlled.

The wide variety of cell types and their biophysical complexities pose a challenge in our ability to understand oscillatory activities produced by cellular-based computational network models. This challenge stems from their high-dimensional and multiparametric natures. To overcome this, we implement a solution by linking minimal and detailed models of CA1 microcircuits that generate intrahippocampal (3-12 Hz) theta rhythms. We leverage insights from minimal models to guide explorations of more detailed models and obtain a cellular perspective of theta generation. Our findings distinguish the pyramidal cells as the theta rhythm initiators and reveal that their activity is regularized by the inhibitory cell populations, supporting a proposed hypothesis of an "inhibition-based tuning" mechanism. We find a strong correlation between input current to the pyramidal cells and the resulting local field potential theta frequency, indicating that intrinsic pyramidal cell properties underpin network frequency characteristics. This work provides a cellular-based foundation from which in vivo theta activities can be explored.

A Hypothesis for Theta Rhythm Frequency Control in CA1 Microcircuits.

Computational models of neural circuits with varying levels of biophysical detail have been generated in pursuit of an underlying mechanism explaining the ubiquitous hippocampal theta rhythm. However, within the theta rhythm are at least two types with distinct frequencies associated with different behavioral states, an aspect that must be considered in pursuit of these mechanistic explanations. Here, using our previously developed excitatory-inhibitory network models that generate theta rhythms, we investigate the robustness of theta generation to intrinsic neuronal variability by building a database of heterogeneous excitatory cells and implementing them in our microcircuit model. We specifically investigate the impact of three key "building block" features of the excitatory cell model that underlie our model design: these cells' rheobase, their capacity for post-inhibitory rebound, and their spike-frequency adaptation. We show that theta rhythms at various frequencies can arise dependent upon the combination of these building block features, and we find that the speed of these oscillations are dependent upon the excitatory cells' response to inhibitory drive, as encapsulated by their phase response curves. Taken together, these findings support a hypothesis for theta frequency control that includes two aspects: (i) an internal mechanism that stems from the building block features of excitatory cell dynamics; (ii) an external mechanism that we describe as "inhibition-based tuning" of excitatory cell firing. We propose that these mechanisms control theta rhythm frequencies and underlie their robustness.

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.

Alterations in Intrinsic and Synaptic Properties of Hippocampal CA1 VIP Interneurons During Aging.

Learning and memory deficits are hallmarks of the aging brain, with cortical neuronal circuits representing the main target in cognitive deterioration. While GABAergic inhibitory and disinhibitory circuits are critical in supporting cognitive processes, their roles in age-related cognitive decline remain largely unknown. Here, we examined the morphological and physiological properties of the hippocampal CA1 vasoactive intestinal peptide/calretinin-expressing (VIP+/CR+) type 3 interneuron-specific (I-S3) cells across mouse lifespan. Our data showed that while the number and morphological features of I-S3 cells remained unchanged, their firing and synaptic properties were significantly altered in old animals. In particular, the action potential duration and the level of steady-state depolarization were significantly increased in old animals in parallel with a significant decrease in the maximal firing frequency. Reducing the fast-delayed rectifier potassium or transient sodium conductances in I-S3 cell computational models could reproduce the age-related changes in I-S3 cell firing properties. However, experimental data revealed no difference in the activation properties of the Kv3.1 and A-type potassium currents, indicating that transient sodium together with other ion conductances may be responsible for the observed phenomena. Furthermore, I-S3 cells in aged mice received a stronger inhibitory drive due to concomitant increase in the amplitude and frequency of spontaneous inhibitory currents. These age-associated changes in the I-S3 cell properties occurred in parallel with an increased inhibition of their target interneurons and were associated with spatial memory deficits and increased anxiety. Taken together, these data indicate that VIP+/CR+ interneurons responsible for local circuit disinhibition survive during aging but exhibit significantly altered physiological properties, which may result in the increased inhibition of hippocampal interneurons and distorted mnemonic functions.

Integration of Within-Cell Experimental Data With Multi-Compartmental Modeling Predicts H-Channel Densities and Distributions in Hippocampal OLM Cells.

Determining biophysical details of spatially extended neurons is a challenge that needs to be overcome if we are to understand the dynamics of brain function from cellular perspectives. Moreover, we now know that we should not average across recordings from many cells of a given cell type to obtain quantitative measures such as conductance since measures can vary multiple-fold for a given cell type. In this work we examine whether a tight combination of experimental and computational work can address this challenge. The oriens-lacunosum/moleculare (OLM) interneuron operates as a "gate" that controls incoming sensory and ongoing contextual information in the CA1 of the hippocampus, making it essential to understand how its biophysical properties contribute to memory function. OLM cells fire phase-locked to the prominent hippocampal theta rhythms, and we previously used computational models to show that OLM cells exhibit high or low theta spiking resonance frequencies that depend respectively on whether their dendrites have hyperpolarization-activated cation channels (h-channels) or not. However, whether OLM cells actually possess dendritic h-channels is unknown at present. We performed a set of whole-cell recordings of OLM cells from mouse hippocampus and constructed three multi-compartment models using morphological and electrophysiological parameters extracted from the same OLM cell, including per-cell pharmacologically isolated h-channel currents. We found that the models best matched experiments when h-channels were present in the dendrites of each of the three model cells created. This strongly suggests that h-channels must be present in OLM cell dendrites and are not localized to their somata. Importantly, this work shows that a tight integration of model and experiment can help tackle the challenge of characterizing biophysical details and distributions in spatially extended neurons. Full spiking models were built for two of the OLM cells, matching their current clamp cell-specific electrophysiological recordings. Overall, our work presents a technical advancement in modeling OLM cells. Our models are available to the community to use to gain insight into cellular dynamics underlying hippocampal function.

Neurostimulation stabilizes spiking neural networks by disrupting seizure-like oscillatory transitions.

An improved understanding of the mechanisms underlying neuromodulatory approaches to mitigate seizure onset is needed to identify clinical targets for the treatment of epilepsy. Using a Wilson-Cowan-motivated network of inhibitory and excitatory populations, we examined the role played by intrinsic and extrinsic stimuli on the network's predisposition to sudden transitions into oscillatory dynamics, similar to the transition to the seizure state. Our joint computational and mathematical analyses revealed that such stimuli, be they noisy or periodic in nature, exert a stabilizing influence on network responses, disrupting the development of such oscillations. Based on a combination of numerical simulations and mean-field analyses, our results suggest that high variance and/or high frequency stimulation waveforms can prevent multi-stability, a mathematical harbinger of sudden changes in network dynamics. By tuning the neurons' responses to input, stimuli stabilize network dynamics away from these transitions. Furthermore, our research shows that such stabilization of neural activity occurs through a selective recruitment of inhibitory cells, providing a theoretical undergird for the known key role these cells play in both the healthy and diseased brain. Taken together, these findings provide new vistas on neuromodulatory approaches to stabilize neural microcircuit activity.

Computationally going where experiments cannot: a dynamical assessment of dendritic ion channel currents during in vivo-like states.

Background: Despite technological advances, how specific cell types are involved in brain function remains shrouded in mystery. Further, little is known about the contribution of different ion channel currents to cell excitability across different neuronal subtypes and their dendritic compartments in vivo. The picture that we do have is largely based on somatic recordings performed in vitro. Uncovering dendritic ion channel current contributions in neuron subtypes that represent a minority of the neuronal population is not currently a feasible task using purely experimental means. Methods: We employ two morphologically-detailed multi-compartment models of a specific type of inhibitory interneuron, the oriens lacunosum moleculare (OLM) cell. The OLM cell is a well-studied cell type in CA1 hippocampus that is important in gating sensory and contextual information. We create in vivo-like states for these cellular models by including levels of synaptic bombardment that would occur in vivo. Using visualization tools and analyses we assess the ion channel current contribution profile across the different somatic and dendritic compartments of the models. Results: We identify changes in dendritic excitability, ion channel current contributions and co-activation patterns between in vitro and in vivo-like states. Primarily, we find that the relative timing between ion channel currents are mostly invariant between states, but exhibit changes in magnitudes and decreased propagation across dendritic compartments. We also find enhanced dendritic hyperpolarization-activated cyclic nucleotide-gated channel (h-channel) activation during in vivo-like states, which suggests that dendritically located h-channels are functionally important in altering signal propagation in the behaving animal. Conclusions: Overall, we have demonstrated, using computational modelling, the dynamical changes that can occur to ion channel mechanisms governing neuronal spiking in vitro and in vivo. In particular, we have shown that the magnitudes of some ion channel current contributions are differentially altered during in vivo-like states relative to in vitro.

Common Principles in Functional Organization of VIP/Calretinin Cell-Driven Disinhibitory Circuits Across Cortical Areas.

In the brain, there is a vast diversity of different structures, circuitries, cell types, and cellular genetic expression profiles. While this large diversity can often occlude a clear understanding of how the brain works, careful analyses of analogous studies performed across different brain areas can hint at commonalities in neuronal organization. This in turn can yield a fundamental understanding of necessary circuitry components that are crucial for how information is processed across the brain. In this review, we outline recent in vivo and in vitro studies that have been performed in different cortical areas to characterize the vasoactive intestinal polypeptide (VIP)- and/or calretinin (CR)-expressing cells that specialize in inhibiting GABAergic interneurons. In doing so, we make the case that, across cortical structures, interneuron-specific cells commonly specialize in the synaptic disinhibition of excitatory neurons, which can ungate the integration and plasticity of external inputs onto excitatory neurons. In line with this, activation of interneuron- specific cells enhances animal performance across a variety of behavioral tasks that involve learning, memory formation, and sensory discrimination, and may represent a key target for therapeutic interventions under different pathological conditions. As such, interneuron-specific cells across different cortical structures are an essential network component for information processing and normal brain function.

Synaptic Mechanisms Underlying the Network State-Dependent Recruitment of VIP-Expressing Interneurons in the CA1 Hippocampus.

Disinhibition is a widespread circuit mechanism for information selection and transfer. In the hippocampus, disinhibition of principal cells is provided by the interneuron-specific interneurons that express the vasoactive intestinal polypeptide (VIP-IS) and innervate selectively inhibitory interneurons. By combining optophysiological experiments with computational models, we determined the impact of synaptic inputs onto the network state-dependent recruitment of VIP-IS cells. We found that VIP-IS cells fire spikes in response to both the Schaffer collateral and the temporoammonic pathway activation. Moreover, by integrating their intrinsic and synaptic properties into computational models, we predicted recruitment of these cells between the rising phase and peak of theta oscillation and during ripples. Two-photon Ca2+-imaging in awake mice supported in part the theoretical predictions, revealing a significant speed modulation of VIP-IS cells and their preferential albeit delayed recruitment during theta-run epochs, with estimated firing at the rising phase and peak of the theta cycle. However, it also uncovered that VIP-IS cells are not activated during ripples. Thus, given the preferential theta-modulated firing of VIP-IS cells in awake hippocampus, we postulate that these cells may be important for information gating during spatial navigation and memory encoding.

Using computational models to predict in vivo synaptic inputs to interneuron specific 3 (IS3) cells of CA1 hippocampus that also allow their recruitment during rhythmic states.

Brain coding strategies are enabled by the balance of synaptic inputs that individual neurons receive as determined by the networks in which they reside. Inhibitory cell types contribute to brain function in distinct ways but recording from specific, inhibitory cell types during behaviour to determine their contributions is highly challenging. In particular, the in vivo activities of vasoactive intestinal peptide-expressing interneuron specific 3 (IS3) cells in the hippocampus that only target other inhibitory cells are unknown at present. We perform a massive, computational exploration of possible synaptic inputs to IS3 cells using multi-compartment models and optimized synaptic parameters. We find that asynchronous, in vivo-like states that are sensitive to additional theta-timed inputs (8 Hz) exist when excitatory and inhibitory synaptic conductances are approximately equally balanced and with low numbers of activated synapses receiving correlated inputs. Specifically, under these balanced conditions, the input resistance is larger with higher mean spike firing rates relative to other activated synaptic conditions investigated. Incoming theta-timed inputs result in strongly increased spectral power relative to baseline. Thus, using a generally applicable computational approach we predict the existence and features of background, balanced states in hippocampal circuits.

Inhibitory Network Bistability Explains Increased Interneuronal Activity Prior to Seizure Onset.

Recent experimental literature has revealed that GABAergic interneurons exhibit increased activity prior to seizure onset, alongside additional evidence that such activity is synchronous and may arise abruptly. These findings have led some to hypothesize that this interneuronal activity may serve a causal role in driving the sudden change in brain activity that heralds seizure onset. However, the mechanisms predisposing an inhibitory network toward increased activity, specifically prior to ictogenesis, without a permanent change to inputs to the system remain unknown. We address this question by comparing simulated inhibitory networks containing control interneurons and networks containing hyperexcitable interneurons modeled to mimic treatment with 4-Aminopyridine (4-AP), an agent commonly used to model seizures in vivo and in vitro. Our in silico study demonstrates that model inhibitory networks with 4-AP interneurons are more prone than their control counterparts to exist in a bistable state in which asynchronously firing networks can abruptly transition into synchrony driven by a brief perturbation. This transition into synchrony brings about a corresponding increase in overall firing rate. We further show that perturbations driving this transition could arise in vivo from background excitatory synaptic activity in the cortex. Thus, we propose that bistability explains the increase in interneuron activity observed experimentally prior to seizure via a transition from incoherent to coherent dynamics. Moreover, bistability explains why inhibitory networks containing hyperexcitable interneurons are more vulnerable to this change in dynamics, and how such networks can undergo a transition without a permanent change in the drive. We note that while our comparisons are between networks of control and ictogenic neurons, the conclusions drawn specifically relate to the unusual dynamics that arise prior to seizure, and not seizure onset itself. However, providing a mechanistic explanation for this phenomenon specifically in a pro-ictogenic setting generates experimentally testable hypotheses regarding the role of inhibitory neurons in pre-ictal neural dynamics, and motivates further computational research into mechanisms underlying a newly hypothesized multi-step pathway to seizure initiated by inhibition.

Deciphering the Contribution of Oriens-Lacunosum/Moleculare (OLM) Cells to Intrinsic θ Rhythms Using Biophysical Local Field Potential (LFP) Models.

Oscillations in local field potentials (LFPs) are prevalent and contribute to brain function. An understanding of the cellular correlates and pathways affecting LFPs is needed, but many overlapping pathways in vivo make this difficult to achieve. A prevalent LFP rhythm in the hippocampus associated with memory processing and spatial navigation is the θ (3-12 Hz) oscillation. θ rhythms emerge intrinsically in an in vitro whole hippocampus preparation and this reduced preparation makes it possible to assess the contribution of different cell types to LFP generation. We focus on oriens-lacunosum/moleculare (OLM) cells as a major class of interneurons in the hippocampus. OLM cells can influence pyramidal (PYR) cells through two distinct pathways: by direct inhibition of PYR cell distal dendrites, and by indirect disinhibition of PYR cell proximal dendrites. We use previous inhibitory network models and build biophysical LFP models using volume conductor theory. We examine the effect of OLM cells to ongoing intrinsic LFP θ rhythms by directly comparing our model LFP features with experiment. We find that OLM cell inputs regulate the robustness of LFP responses without affecting their average power and that this robust response depends on coactivation of distal inhibition and basal excitation. We use our models to estimate the spatial extent of the region generating LFP θ rhythms, leading us to predict that about 22,000 PYR cells participate in intrinsic θ generation. Besides obtaining an understanding of OLM cell contributions to intrinsic LFP θ rhythms, our work can help decipher cellular correlates of in vivo LFPs.