Induction and modulation of persistent activity in a layer V PFC microcircuit model

Agranular frontal cortical microcircuit underlying cognitive control in macaques

The error-related negativity and an N2-component recorded over medial frontal cortex index core functions of cognitive control. While they are known to originate from agranular frontal areas, the underlying microcircuit mechanisms remain elusive. Most insights about microcircuit function have been derived from variations of the so-called canonical microcircuit model. These microcircuit architectures are based extensively on studies from granular sensory cortical areas in monkeys, cats, and rodents. However, evidence has shown striking cytoarchitectonic differences across species and differences in the functional relationships across cortical layers in agranular compared to granular sensory areas. In this minireview, we outline a tentative microcircuit model underlying cognitive control in the agranular frontal cortex of primates. The model incorporates the main GABAergic interneuron subclasses with specific laminar arrangements and target regions on pyramidal cells. We emphasize the role of layer 5 pyramidal cells in error and conflict detection. We offer several specific questions necessary for creating a specific intrinsic microcircuit model of the agranular frontal cortex.

The effect of alterations of schizophrenia-associated genes on gamma band oscillations

Abnormalities in the synchronized oscillatory activity of neurons in general and, specifically in the gamma band, might play a crucial role in the pathophysiology of schizophrenia. While these changes in oscillatory activity have traditionally been linked to alterations at the synaptic level, we demonstrate here, using computational modeling, that common genetic variants of ion channels can contribute strongly to this effect. Our model of primary auditory cortex highlights multiple schizophrenia-associated genetic variants that reduce gamma power in an auditory steady-state response task. Furthermore, we show that combinations of several of these schizophrenia-associated variants can produce similar effects as the more traditionally considered synaptic changes. Overall, our study provides a mechanistic link between schizophrenia-associated common genetic variants, as identified by genome-wide association studies, and one of the most robust neurophysiological endophenotypes of schizophrenia.

Inhibitory control of frontal metastability sets the temporal signature of cognition

Cortical dynamics are organized over multiple anatomical and temporal scales. The mechanistic origin of the temporal organization and its contribution to cognition remain unknown. Here we demonstrate the cause of this organization by studying a specific temporal signature (time constant and latency) of neural activity. In monkey frontal areas, recorded during flexible decisions, temporal signatures display specific area-dependent ranges, as well as anatomical and cell-type distributions. Moreover, temporal signatures are functionally adapted to behaviorally relevant timescales. Fine-grained biophysical network models, constrained to account for experimentally observed temporal signatures, reveal that after-hyperpolarization potassium and inhibitory GABA-B conductances critically determine areas' specificity. They mechanistically account for temporal signatures by organizing activity into metastable states, with inhibition controlling state stability and transitions. As predicted by models, state durations non-linearly scale with temporal signatures in monkey, matching behavioral timescales. Thus, local inhibitory-controlled metastability constitutes the dynamical core specifying the temporal organization of cognitive functions in frontal areas.

Hippocampal-Cortical Memory Trace Transfer and Reactivation Through Cell-Specific Stimulus and Spontaneous Background Noise

The hippocampus plays important roles in memory formation and retrieval through sharp-wave-ripples. Recent studies have shown that certain neuron populations in the prefrontal cortex (PFC) exhibit coordinated reactivations during awake ripple events. These experimental findings suggest that the awake ripple is an important biomarker, through which the hippocampus interacts with the neocortex to assist memory formation and retrieval. However, the computational mechanisms of this ripple based hippocampal-cortical coordination are still not clear due to the lack of unified models that include both the hippocampal and cortical networks. In this work, using a coupled biophysical model of both CA1 and PFC, we investigate possible mechanisms of hippocampal-cortical memory trace transfer and the conditions that assist reactivation of the transferred memory traces in the PFC. To validate our model, we first show that the local field potentials generated in the hippocampus and PFC exhibit ripple range activities that are consistent with the recent experimental studies. Then we demonstrate that during ripples, sequence replays can successfully transfer the information stored in the hippocampus to the PFC recurrent networks. We investigate possible mechanisms of memory retrieval in PFC networks. Our results suggest that the stored memory traces in the PFC network can be retrieved through two different mechanisms, namely the cell-specific input representing external stimuli and non-specific spontaneous background noise representing spontaneous memory recall events. Importantly, in both cases, the memory reactivation quality is robust to network connection loss. Finally, we find out that the quality of sequence reactivations is enhanced by both increased number of SWRs and an optimal background noise intensity, which tunes the excitability of neurons to a proper level. Our study presents a mechanistic explanation for the memory trace transfer from the hippocampus to neocortex through ripple coupling in awake states and reports two different mechanisms by which the stored memory traces can be successfully retrieved.

Alterations in Schizophrenia-Associated Genes Can Lead to Increased Power in Delta Oscillations

Genome-wide association studies have implicated many ion channels in schizophrenia pathophysiology. Although the functions of these channels are relatively well characterized by single-cell studies, the contributions of common variation in these channels to neurophysiological biomarkers and symptoms of schizophrenia remain elusive. Here, using computational modeling, we show that a common biomarker of schizophrenia, namely, an increase in delta-oscillation power, may be a direct consequence of altered expression or kinetics of voltage-gated ion channels or calcium transporters. Our model of a circuit of layer V pyramidal cells highlights multiple types of schizophrenia-related variants that contribute to altered dynamics in the delta-frequency band. Moreover, our model predicts that the same membrane mechanisms that increase the layer V pyramidal cell network gain and response to delta-frequency oscillations may also cause a deficit in a single-cell correlate of the prepulse inhibition, which is a behavioral biomarker highly associated with schizophrenia.

Biophysical Psychiatry-How Computational Neuroscience Can Help to Understand the Complex Mechanisms of Mental Disorders

The brain is the most complex of human organs, and the pathophysiology underlying abnormal brain function in psychiatric disorders is largely unknown. Despite the rapid development of diagnostic tools and treatments in most areas of medicine, our understanding of mental disorders and their treatment has made limited progress during the last decades. While recent advances in genetics and neuroscience have a large potential, the complexity and multidimensionality of the brain processes hinder the discovery of disease mechanisms that would link genetic findings to clinical symptoms and behavior. This applies also to schizophrenia, for which genome-wide association studies have identified a large number of genetic risk loci, spanning hundreds of genes with diverse functionalities. Importantly, the multitude of the associated variants and their prevalence in the healthy population limit the potential of a reductionist functional genetics approach as a stand-alone solution to discover the disease pathology. In this review, we outline the key concepts of a "biophysical psychiatry," an approach that employs large-scale mechanistic, biophysics-founded computational modelling to increase transdisciplinary understanding of the pathophysiology and strive toward robust predictions. We discuss recent scientific advances that allow a synthesis of previously disparate fields of psychiatry, neurophysiology, functional genomics, and computational modelling to tackle open questions regarding the pathophysiology of heritable mental disorders. We argue that the complexity of the increasing amount of genetic data exceeds the capabilities of classical experimental assays and requires computational approaches. Biophysical psychiatry, based on modelling diseased brain networks using existing and future knowledge of basic genetic, biochemical, and functional properties on a single neuron to a microcircuit level, may allow a leap forward in deriving interpretable biomarkers and move the field toward novel treatment options.

Subcellular Localization and Activity of TRPM4 in Medial Prefrontal Cortex Layer 2/3

TRPM4 is a Ca²⁺-activated non-selective cationic channel that conducts monovalent cations. TRPM4 has been proposed to contribute to burst firing and sustained activity in several brain regions, however, the cellular and subcellular pattern of TRPM4 expression in medial prefrontal cortex (mPFC) during postnatal development has not been elucidated. Here, we use multiplex immunofluorescence labeling of brain sections to characterize the postnatal developmental expression of TRPM4 in the mouse mPFC. We also performed electrophysiological recordings to correlate the expression of TRPM4 immunoreactivity with the presence of TRPM4-like currents. We found that TRPM4 is expressed from the first postnatal day, with expression increasing up to postnatal day 35. Additionally, in perforated patch clamp experiments, we found that TRPM4-like currents were active at resting membrane potentials at all postnatal ages studied. Moreover, TRPM4 is expressed in both pyramidal neurons and interneurons. TRPM4 expression is localized in the soma and proximal dendrites, but not in the axon initial segment of pyramidal neurons. This subcellular localization is consistent with a reduction in the basal current only when we locally perfused 9-Phenanthrol in the soma, but not upon perfusion in the medial or distal dendrites. Our results show a specific localization of TRPM4 expression in neurons in the mPFC and that a 9-Phenanthrol sensitive current is active at resting membrane potential, suggesting specific functional roles in mPFC neurons during postnatal development and in adulthood.

T2N as a new tool for robust electrophysiological modeling demonstrated for mature and adult-born dentate granule cells

Compartmental models are the theoretical tool of choice for understanding single neuron computations. However, many models are incomplete, built ad hoc and require tuning for each novel condition rendering them of limited usability. Here, we present T2N, a powerful interface to control NEURON with Matlab and TREES toolbox, which supports generating models stable over a broad range of reconstructed and synthetic morphologies. We illustrate this for a novel, highly detailed active model of dentate granule cells (GCs) replicating a wide palette of experiments from various labs. By implementing known differences in ion channel composition and morphology, our model reproduces data from mouse or rat, mature or adult-born GCs as well as pharmacological interventions and epileptic conditions. This work sets a new benchmark for detailed compartmental modeling. T2N is suitable for creating robust models useful for large-scale networks that could lead to novel predictions. We discuss possible T2N application in degeneracy studies.

On Rhythms in Neuronal Networks with Recurrent Excitation

We investigate rhythms in networks of neurons with recurrent excitation, that is, with excitatory cells exciting each other. Recurrent excitation can sustain activity even when the cells in the network are driven below threshold, too weak to fire on their own. This sort of "reverberating" activity is often thought to be the basis of working memory. Recurrent excitation can also lead to "runaway" transitions, sudden transitions to high-frequency firing; this may be related to epileptic seizures. Not all fundamental questions about these phenomena have been answered with clarity in the literature. We focus on three questions here: (1) How much recurrent excitation is needed to sustain reverberating activity? How does the answer depend on parameters? (2) Is there a positive minimum frequency of reverberating activity, a positive "onset frequency"? How does it depend on parameters? (3) When do runaway transitions occur? For reduced models, we give mathematical answers to these questions. We also examine computationally to which extent our findings are reflected in the behavior of biophysically more realistic model networks. Our main results can be summarized as follows. (1) Reverberating activity can be fueled by extremely weak slow recurrent excitation, but only by sufficiently strong fast recurrent excitation. (2) The onset of reverberating activity, as recurrent excitation is strengthened or external drive is raised, occurs at a positive frequency. It is faster when the external drive is weaker (and the recurrent excitation stronger). It is slower when the recurrent excitation has a longer decay time constant. (3) Runaway transitions occur only with fast, not with slow, recurrent excitation. We also demonstrate that the relation between reverberating activity fueled by recurrent excitation and runaway transitions can be visualized in an instructive way by a (generalized) cusp catastrophe surface.

Pleiotropic effects of schizophrenia-associated genetic variants in neuron firing and cardiac pacemaking revealed by computational modeling

Schizophrenia patients have an increased risk of cardiac dysfunction. A possible factor underlying this comorbidity are the common variants in the large set of genes that have recently been discovered in genome-wide association studies (GWASs) as risk genes of schizophrenia. Many of these genes control the cell electrogenesis and calcium homeostasis. We applied biophysically detailed models of layer V pyramidal cells and sinoatrial node cells to study the contribution of schizophrenia-associated genes on cellular excitability. By including data from functional genomics literature to simulate the effects of common variants of these genes, we showed that variants of voltage-gated Na⁺ channel or hyperpolarization-activated cation channel-encoding genes cause qualitatively similar effects on layer V pyramidal cell and sinoatrial node cell excitability. By contrast, variants of Ca²⁺ channel or transporter-encoding genes mostly have opposite effects on cellular excitability in the two cell types. We also show that the variants may crucially affect the propagation of the cardiac action potential in the sinus node. These results may help explain some of the cardiac comorbidity in schizophrenia, and may facilitate generation of effective antipsychotic medications without cardiac side-effects such as arrhythmia.

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.

Hetereogeneity in Neuronal Intrinsic Properties: A Possible Mechanism for Hub-Like Properties of the Rat Anterior Cingulate Cortex during Network Activity

The anterior cingulate cortex (ACC) is vital for a range of brain functions requiring cognitive control and has highly divergent inputs and outputs, thus manifesting as a hub in connectomic analyses. Studies show diverse functional interactions within the ACC are associated with network oscillations in the β (20–30 Hz) and γ (30-80 Hz) frequency range. Oscillations permit dynamic routing of information within cortex, a function that depends on bandpass filter–like behavior to selectively respond to specific inputs. However, a putative hub region such as ACC needs to be able to combine inputs from multiple sources rather than select a single input at the expense of others. To address this potential functional dichotomy, we modeled local ACC network dynamics in the rat in vitro. Modal peak oscillation frequencies in the β- and γ-frequency band corresponded to GABA_Aergic synaptic kinetics as seen in other regions; however, the intrinsic properties of ACC principal neurons were highly diverse. Computational modeling predicted that this neuronal response diversity broadened the bandwidth for filtering rhythmic inputs and supported combination—rather than selection—of different frequencies within the canonical γ and β electroencephalograph bands. These findings suggest that oscillating neuronal populations can support either response selection (routing) or combination, depending on the interplay between the kinetics of synaptic inhibition and the degree of heterogeneity of principal cell intrinsic conductances.

Calcium regulation of HCN channels supports persistent activity in a multiscale model of neocortex

Neuronal persistent activity has been primarily assessed in terms of electrical mechanisms, without attention to the complex array of molecular events that also control cell excitability. We developed a multiscale neocortical model proceeding from the molecular to the network level to assess the contributions of calcium (Ca(2+)) regulation of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels in providing additional and complementary support of continuing activation in the network. The network contained 776 compartmental neurons arranged in the cortical layers, connected using synapses containing AMPA/NMDA/GABAA/GABAB receptors. Metabotropic glutamate receptors (mGluR) produced inositol triphosphate (IP3) which caused the release of Ca(2+) from endoplasmic reticulum (ER) stores, with reuptake by sarco/ER Ca(2+)-ATP-ase pumps (SERCA), and influence on HCN channels. Stimulus-induced depolarization led to Ca(2+) influx via NMDA and voltage-gated Ca(2+) channels (VGCCs). After a delay, mGluR activation led to ER Ca(2+) release via IP3 receptors. These factors increased HCN channel conductance and produced firing lasting for ∼1min. The model displayed inter-scale synergies among synaptic weights, excitation/inhibition balance, firing rates, membrane depolarization, Ca(2+) levels, regulation of HCN channels, and induction of persistent activity. The interaction between inhibition and Ca(2+) at the HCN channel nexus determined a limited range of inhibition strengths for which intracellular Ca(2+) could prepare population-specific persistent activity. Interactions between metabotropic and ionotropic inputs to the neuron demonstrated how multiple pathways could contribute in a complementary manner to persistent activity. Such redundancy and complementarity via multiple pathways is a critical feature of biological systems. Mediation of activation at different time scales, and through different pathways, would be expected to protect against disruption, in this case providing stability for persistent activity.

A robust in vivo-like persistent firing supported by a hybrid of intracellular and synaptic mechanisms

Persistent firing is believed to support short-term information retention in the brain. Established hypotheses make use of the recurrent synaptic connectivity to support persistent firing. However, this mechanism is known to suffer from a lack of robustness. On the other hand, persistent firing can be supported by an intrinsic cellular mechanism in multiple brain areas. However, the consequences of having both the intrinsic and the synaptic mechanisms (a hybrid model) on persistent firing remain largely unknown. The goal of this study is to investigate whether a hybrid neural network model with these two mechanisms has advantages over a conventional recurrent network based model. Our computer simulations were based on in vitro recordings obtained from hippocampal CA3 pyramidal cells under cholinergic receptor activation. Calcium activated non-specific cationic (CAN) current supported persistent firing in the Hodgkin-Huxley style cellular models. Our results suggest that the hybrid model supports persistent firing within a physiological frequency range over a wide range of different parameters, eliminating parameter sensitivity issues generally recognized in network based persistent firing. In addition, persistent firing in the hybrid model is substantially more robust against distracting inputs, can coexist with theta frequency oscillations, and supports pattern completion.

Synchronous rhythmic interaction enhances children's perceived similarity and closeness towards each other

Inter-personal synchronization is important for performing many cooperative tasks. Notably, synchrony has also been shown to have considerable positive social influences, possibly mediated by synchrony-induced changes in social attitude such as an increased sense of similarity and affiliation between interacting individuals. This effect has been demonstrated in adults, but it is unknown whether synchrony might have a similar impact on the social attitudes of children. We thus set to directly examine the influence of synchronous rhythmic interaction on perceived similarity and closeness in pairs of 8-9 year old children. We found that children who had participated in a synchronous interaction regarded their interacting partner as more similar and closer to themselves than children who had not interacted at all or who had taken part in an asynchronous interaction. These findings reveal that synchronous interaction can positively alter social attitudes between interacting children, suggesting a potential mechanism by which synchrony may enhance positive social interaction through attitudinal shift.

SPIKY: a graphical user interface for monitoring spike train synchrony

Techniques for recording large-scale neuronal spiking activity are developing very fast. This leads to an increasing demand for algorithms capable of analyzing large amounts of experimental spike train data. One of the most crucial and demanding tasks is the identification of similarity patterns with a very high temporal resolution and across different spatial scales. To address this task, in recent years three time-resolved measures of spike train synchrony have been proposed, the ISI-distance, the SPIKE-distance, and event synchronization. The Matlab source codes for calculating and visualizing these measures have been made publicly available. However, due to the many different possible representations of the results the use of these codes is rather complicated and their application requires some basic knowledge of Matlab. Thus it became desirable to provide a more user-friendly and interactive interface. Here we address this need and present SPIKY, a graphical user interface that facilitates the application of time-resolved measures of spike train synchrony to both simulated and real data. SPIKY includes implementations of the ISI-distance, the SPIKE-distance, and the SPIKE-synchronization (an improved and simplified extension of event synchronization) that have been optimized with respect to computation speed and memory demand. It also comprises a spike train generator and an event detector that makes it capable of analyzing continuous data. Finally, the SPIKY package includes additional complementary programs aimed at the analysis of large numbers of datasets and the estimation of significance levels.

Synaptic clustering within dendrites: an emerging theory of memory formation

It is generally accepted that complex memories are stored in distributed representations throughout the brain, however the mechanisms underlying these representations are not understood. Here, we review recent findings regarding the subcellular mechanisms implicated in memory formation, which provide evidence for a dendrite-centered theory of memory. Plasticity-related phenomena which affect synaptic properties, such as synaptic tagging and capture, synaptic clustering, branch strength potentiation and spinogenesis provide the foundation for a model of memory storage that relies heavily on processes operating at the dendrite level. The emerging picture suggests that clusters of functionally related synapses may serve as key computational and memory storage units in the brain. We discuss both experimental evidence and theoretical models that support this hypothesis and explore its advantages for neuronal function.

Dendritic Excitability and Gain Control in Recurrent Cortical Microcircuits

Layer 5 thick tufted pyramidal cells (TTCs) in the neocortex are particularly electrically complex, owing to their highly excitable dendrites. The interplay between dendritic nonlinearities and recurrent cortical microcircuit activity in shaping network response is largely unknown. We simulated detailed conductance-based models of TTCs forming recurrent microcircuits that were interconnected as found experimentally; the network was embedded in a realistic background synaptic activity. TTCs microcircuits significantly amplified brief thalamocortical inputs; this cortical gain was mediated by back-propagation activated N-methyl-d-aspartate depolarizations and dendritic back-propagation-activated Ca²⁺ spike firing, ignited by the coincidence of thalamic-activated somatic spike and local dendritic synaptic inputs, originating from the cortical microcircuit. Surprisingly, dendritic nonlinearities in TTCs microcircuits linearly multiplied thalamic inputs—amplifying them while maintaining input selectivity. Our findings indicate that dendritic nonlinearities are pivotal in controlling the gain and the computational functions of TTCs microcircuits, which serve as a dominant output source for the neocortex.

Dendritic nonlinearities reduce network size requirements and mediate ON and OFF states of persistent activity in a PFC microcircuit model

Technological advances have unraveled the existence of small clusters of co-active neurons in the neocortex. The functional implications of these microcircuits are in large part unexplored. Using a heavily constrained biophysical model of a L5 PFC microcircuit, we recently showed that these structures act as tunable modules of persistent activity, the cellular correlate of working memory. Here, we investigate the mechanisms that underlie persistent activity emergence (ON) and termination (OFF) and search for the minimum network size required for expressing these states within physiological regimes. We show that (a) NMDA-mediated dendritic spikes gate the induction of persistent firing in the microcircuit. (b) The minimum network size required for persistent activity induction is inversely proportional to the synaptic drive of each excitatory neuron. (c) Relaxation of connectivity and synaptic delay constraints eliminates the gating effect of NMDA spikes, albeit at a cost of much larger networks. (d) Persistent activity termination by increased inhibition depends on the strength of the synaptic input and is negatively modulated by dADP. (e) Slow synaptic mechanisms and network activity contain predictive information regarding the ability of a given stimulus to turn ON and/or OFF persistent firing in the microcircuit model. Overall, this study zooms out from dendrites to cell assemblies and suggests a tight interaction between dendritic non-linearities and network properties (size/connectivity) that may facilitate the short-memory function of the PFC.

Working memory, the ability to retain information for a short period of time, is a fundamental cognitive function that shapes behavior. The cellular correlate of working memory is the prolonged spiking (persistent) activity of neurons in the prefrontal cortex. Impairments of prefrontal cortex functionalities and working memory have been associated with a variety of cognitive disorders, such as schizophrenia, the attention deficit hyperactivity disorder, and drug addiction. Hence, understanding how neurons embedded in the local circuitry support and maintain persistent activity is of outmost importance. Our work uses a multi-level integrative approach spanning from the dendritic, to the neuronal and network levels to identify the key biophysical and anatomical mechanisms contributing to persistent activity, leading to a number of high impact findings: it predicts a tradeoff between dendritic regenerative events and the size of a network expressing persistent activity. It also proposes when and how the persistent state can be stabilized, opening new avenues for pharmacological interventions. Finally, it describes decoding mechanisms for upcoming ON/OFF state transitions, furthering our understanding of information processing in the PFC and shedding new light on the emergence of anticipatory behaviors.

Modulatory effects of inhibition on persistent activity in a cortical microcircuit model

Neocortical network activity is generated through a dynamic balance between excitation, provided by pyramidal neurons, and inhibition, provided by interneurons. Imbalance of the excitation/inhibition ratio has been identified in several neuropsychiatric diseases, such as schizophrenia, autism and epilepsy, which also present with other cognitive deficits and symptoms associated with prefrontal cortical (PFC) dysfunction. We undertook a computational approach to study how changes in the excitation/inhibition balance in a PFC microcircuit model affect the properties of persistent activity, considered the cellular correlate of working memory function in PFC. To this end, we constructed a PFC microcircuit, consisting of pyramidal neuron models and all three different interneuron types: fast-spiking (FS), regular-spiking (RS), and irregular-spiking (IS) interneurons. Persistent activity was induced in the microcircuit model with a stimulus to the proximal apical dendrites of the pyramidal neuron models, and its properties were analyzed, such as the induction profile, the interspike intervals (ISIs) and neuronal synchronicity. Our simulations showed that (a) the induction but not the firing frequency or neuronal synchronicity is modulated by changes in the NMDA-to-AMPA ratio on FS interneuron model, (b) removing or decreasing the FS model input to the pyramidal neuron models greatly limited the biophysical modulation of persistent activity induction, decreased the ISIs and neuronal synchronicity during persistent activity, (c) the induction and firing properties could not be altered by the addition of other inhibitory inputs to the soma (from RS or IS models), and (d) the synchronicity change could be reversed by the addition of other inhibitory inputs to the soma, but beyond the levels of the control network. Thus, generic somatic inhibition acts as a pacemaker of persistent activity and FS specific inhibition modulates the output of the pacemaker.