A Novel Milli-fluidic Liver Tissue Chip with Continuous Recirculation for Predictive Pharmacokinetics Applications.

A crucial step in lead selection during drug development is accurate estimation and optimization of hepatic clearance using in vitro methods. However, current methods are limited by factors such as lack of physiological relevance, short culture/incubation times that are not consistent with drug exposure patterns in patients, use of drug absorbing materials, and evaporation during long-term incubation. To address these technological needs, we developed a novel milli-fluidic human liver tissue chip (LTC) that was designed with continuous media recirculation and optimized for hepatic cultures using human primary hepatocytes. Here, we characterized the LTC using a series of physiologically relevant metrics and test compounds to demonstrate that we could accurately predict the PK of both low- and high-clearance compounds. The non-biological characterization indicated that the cyclic olefin copolymer (COC)-based LTC exhibited negligible evaporation and minimal non-specific binding of drugs of varying ionic states and lipophilicity. Biologically, the LTC exhibited functional and polarized hepatic culture with sustained metabolic CYP activity for at least 15 days. This long-term culture was then used for drug clearance studies for low- and high-clearance compounds for at least 12 days, and clearance was estimated for a range of compounds with high in vitro-in vivo correlation (IVIVC). We also demonstrated that LTC can be induced by rifampicin, and the culture age had insignificant effect on depletion kinetic and predicted clearance value. Thus, we used advances in bioengineering to develop a novel purpose-built platform with high reproducibility and minimal variability to address unmet needs for PK applications.

Emotional Memory Processing during REM Sleep with Implications for Post-Traumatic Stress Disorder.

REM sleep is important for the processing of emotional memories, including fear memories. Rhythmic interactions, especially in the theta band, between the medial prefrontal cortex (mPFC) and limbic structures are thought to play an important role, but the ways in which memory processing occurs at a mechanistic and circuits level are largely unknown. To investigate how rhythmic interactions lead to fear extinction during REM sleep, we used a biophysically based model that included the infralimbic cortex (IL), a part of the mPFC with a critical role in suppressing fear memories. Theta frequency (4-12 Hz) inputs to a given cell assembly in IL, representing an emotional memory, resulted in the strengthening of connections from the IL to the amygdala and the weakening of connections from the amygdala to the IL, resulting in the suppression of the activity of fear expression cells for the associated memory. Lower frequency (4 Hz) theta inputs effected these changes over a wider range of input strengths. In contrast, inputs at other frequencies were ineffective at causing these synaptic changes and did not suppress fear memories. Under post-traumatic stress disorder (PTSD) REM sleep conditions, rhythmic activity dissipated, and 4 Hz theta inputs to IL were ineffective, but higher-frequency (10 Hz) theta inputs to IL induced changes similar to those seen with 4 Hz inputs under normal REM sleep conditions, resulting in the suppression of fear expression cells. These results suggest why PTSD patients may repeatedly experience the same emotionally charged dreams and suggest potential neuromodulatory therapies for the amelioration of PTSD symptoms.SIGNIFICANCE STATEMENT Rhythmic interactions in the theta band between the mPFC and limbic structures are thought to play an important role in processing emotional memories, including fear memories, during REM sleep. The infralimbic cortex (IL) in the mPFC is thought to play a critical role in suppressing fear memories. We show that theta inputs to the IL, unlike other frequency inputs, are effective in producing synaptic changes that suppress the activity of fear expression cells associated with a given memory. Under PTSD REM sleep conditions, lower-frequency (4 Hz) theta inputs to the IL do not suppress the activity of fear expression cells associated with the given memory but, surprisingly, 10 Hz inputs do. These results suggest potential neuromodulatory therapies for PTSD.

Constraints on persistent activity in a biologically detailed network model of the prefrontal cortex with heterogeneities.

Persistent activity, the maintenance of neural activation over short periods of time in cortical networks, is widely thought to underlie the cognitive function of working memory. A large body of modeling studies has reproduced this kind of activity using cell assemblies with strengthened synaptic connections. However, almost all of these studies have considered persistent activity within networks with homogeneous neurons and synapses, making it difficult to judge the validity of such model results for cortical dynamics, which is based on highly heterogeneous neurons. Here, we consider persistent activity in a detailed, strongly data-driven network model of the prefrontal cortex with heterogeneous neuron and synapse parameters. Surprisingly, persistent activity could not be reproduced in this model without incorporating further constraints. We identified three factors that prevent successful persistent activity: heterogeneity in the cell parameters of interneurons, heterogeneity in the parameters of short-term synaptic plasticity and heterogeneity in the synaptic weights. We also discovered a general dynamic mechanism that prevents persistent activity in the presence of heterogeneities, namely a gradual drop-out of cell assembly neurons out of a bistable regime as input variability increases. Based on this mechanism, we found that persistent activity is recovered if heterogeneity is compensated, e.g., by a homeostatic plasticity mechanism. Cell assemblies shaped in this way may be potentially targeted by distinct inputs or become more responsive to specific tuning or spectral properties. Finally, we show that persistent activity in the model is robust against external noise, but the compensation of heterogeneities may prevent the dynamic generation of intrinsic in vivo-like irregular activity. These results may help informing the ongoing debate about the neural basis of working memory.

Prefrontal oscillations modulate the propagation of neuronal activity required for working memory.

Cognition involves using attended information, maintained in working memory (WM), to guide action. During a cognitive task, a correct response requires flexible, selective gating so that only the appropriate information flows from WM to downstream effectors that carry out the response. In this work, we used biophysically-detailed modeling to explore the hypothesis that network oscillations in prefrontal cortex (PFC), leveraging local inhibition, can independently gate responses to items in WM. The key role of local inhibition was to control the period between spike bursts in the outputs, and to produce an oscillatory response no matter whether the WM item was maintained in an asynchronous or oscillatory state. We found that the WM item that induced an oscillatory population response in the PFC output layer with the shortest period between spike bursts was most reliably propagated. The network resonant frequency (i.e., the input frequency that produces the largest response) of the output layer can be flexibly tuned by varying the excitability of deep layer principal cells. Our model suggests that experimentally-observed modulation of PFC beta-frequency (15-30 Hz) and gamma-frequency (30-80 Hz) oscillations could leverage network resonance and local inhibition to govern the flexible routing of signals in service to cognitive processes like gating outputs from working memory and the selection of rule-based actions. Importantly, we show for the first time that nonspecific changes in deep layer excitability can tune the output gate's resonant frequency, enabling the specific selection of signals encoded by populations in asynchronous or fast oscillatory states. More generally, this represents a dynamic mechanism by which adjusting network excitability can govern the propagation of asynchronous and oscillatory signals throughout neocortex.

Biased competition in the absence of input bias revealed through corticostriatal computation.

Classical accounts of biased competition require an input bias to resolve the competition between neuronal ensembles driving downstream processing. However, flexible and reliable selection of behaviorally relevant ensembles can occur with unbiased stimulation: striatal D1 and D2 spiny projection neurons (SPNs) receive balanced cortical input, yet their activity determines the choice between GO and NO-GO pathways in the basal ganglia. We here present a corticostriatal model identifying three mechanisms that rely on physiological asymmetries to effect rate- and time-coded biased competition in the presence of balanced inputs. First, tonic input strength determines which one of the two SPN phenotypes exhibits a higher mean firing rate. Second, low-strength oscillatory inputs induce higher firing rate in D2 SPNs but higher coherence between D1 SPNs. Third, high-strength inputs oscillating at distinct frequencies can preferentially activate D1 or D2 SPN populations. Of these mechanisms, only the latter accommodates observed rhythmic activity supporting rule-based decision making in prefrontal cortex.

Flexible resonance in prefrontal networks with strong feedback inhibition.

Oscillations are ubiquitous features of brain dynamics that undergo task-related changes in synchrony, power, and frequency. The impact of those changes on target networks is poorly understood. In this work, we used a biophysically detailed model of prefrontal cortex (PFC) to explore the effects of varying the spike rate, synchrony, and waveform of strong oscillatory inputs on the behavior of cortical networks driven by them. Interacting populations of excitatory and inhibitory neurons with strong feedback inhibition are inhibition-based network oscillators that exhibit resonance (i.e., larger responses to preferred input frequencies). We quantified network responses in terms of mean firing rates and the population frequency of network oscillation; and characterized their behavior in terms of the natural response to asynchronous input and the resonant response to oscillatory inputs. We show that strong feedback inhibition causes the PFC to generate internal (natural) oscillations in the beta/gamma frequency range (>15 Hz) and to maximize principal cell spiking in response to external oscillations at slightly higher frequencies. Importantly, we found that the fastest oscillation frequency that can be relayed by the network maximizes local inhibition and is equal to a frequency even higher than that which maximizes the firing rate of excitatory cells; we call this phenomenon population frequency resonance. This form of resonance is shown to determine the optimal driving frequency for suppressing responses to asynchronous activity. Lastly, we demonstrate that the natural and resonant frequencies can be tuned by changes in neuronal excitability, the duration of feedback inhibition, and dynamic properties of the input. Our results predict that PFC networks are tuned for generating and selectively responding to beta- and gamma-rhythmic signals due to the natural and resonant properties of inhibition-based oscillators. They also suggest strategies for optimizing transcranial stimulation and using oscillatory networks in neuromorphic engineering.

DynaSim: A MATLAB Toolbox for Neural Modeling and Simulation.

DynaSim is an open-source MATLAB/GNU Octave toolbox for rapid prototyping of neural models and batch simulation management. It is designed to speed up and simplify the process of generating, sharing, and exploring network models of neurons with one or more compartments. Models can be specified by equations directly (similar to XPP or the Brian simulator) or by lists of predefined or custom model components. The higher-level specification supports arbitrarily complex population models and networks of interconnected populations. DynaSim also includes a large set of features that simplify exploring model dynamics over parameter spaces, running simulations in parallel using both multicore processors and high-performance computer clusters, and analyzing and plotting large numbers of simulated data sets in parallel. It also includes a graphical user interface (DynaSim GUI) that supports full functionality without requiring user programming. The software has been implemented in MATLAB to enable advanced neural modeling using MATLAB, given its popularity and a growing interest in modeling neural systems. The design of DynaSim incorporates a novel schema for model specification to facilitate future interoperability with other specifications (e.g., NeuroML, SBML), simulators (e.g., NEURON, Brian, NEST), and web-based applications (e.g., Geppetto) outside MATLAB. DynaSim is freely available at http://dynasimtoolbox.org. This tool promises to reduce barriers for investigating dynamics in large neural models, facilitate collaborative modeling, and complement other tools being developed in the neuroinformatics community.

Thalamocortical control of propofol phase-amplitude coupling.

The anesthetic propofol elicits many different spectral properties on the EEG, including alpha oscillations (8-12 Hz), Slow Wave Oscillations (SWO, 0.1-1.5 Hz), and dose-dependent phase-amplitude coupling (PAC) between alpha and SWO. Propofol is known to increase GABAA inhibition and decrease H-current strength, but how it generates these rhythms and their interactions is still unknown. To investigate both generation of the alpha rhythm and its PAC to SWO, we simulate a Hodgkin-Huxley network model of a hyperpolarized thalamus and corticothalamic inputs. We find, for the first time, that the model thalamic network is capable of independently generating the sustained alpha seen in propofol, which may then be relayed to cortex and expressed on the EEG. This dose-dependent sustained alpha critically relies on propofol GABAA potentiation to alter the intrinsic spindling mechanisms of the thalamus. Furthermore, the H-current conductance and background excitation of these thalamic cells must be within specific ranges to exhibit any intrinsic oscillations, including sustained alpha. We also find that, under corticothalamic SWO UP and DOWN states, thalamocortical output can exhibit maximum alpha power at either the peak or trough of this SWO; this implies the thalamus may be the source of propofol-induced PAC. Hyperpolarization level is the main determinant of whether the thalamus exhibits trough-max PAC, which is associated with lower propofol dose, or peak-max PAC, associated with higher dose. These findings suggest: the thalamus generates a novel rhythm under GABAA potentiation such as under propofol, its hyperpolarization may determine whether a patient experiences trough-max or peak-max PAC, and the thalamus is a critical component of propofol-induced cortical spectral phenomena. Changes to the thalamus may be a critical part of how propofol accomplishes its effects, including unconsciousness.

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 GABAAergic 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.

Sequential then interactive processing of letters and words in the left fusiform gyrus.

Despite decades of cognitive, neuropsychological and neuroimaging studies, it is unclear if letters are identified before word-form encoding during reading, or if letters and their combinations are encoded simultaneously and interactively. Here using functional magnetic resonance imaging, we show that a 'letter-form' area (responding more to consonant strings than false fonts) can be distinguished from an immediately anterior 'visual word-form area' in ventral occipito-temporal cortex (responding more to words than consonant strings). Letter-selective magnetoencephalographic responses begin in the letter-form area ∼60 ms earlier than word-selective responses in the word-form area. Local field potentials confirm the latency and location of letter-selective responses. This area shows increased high-gamma power for ∼400 ms, and strong phase-locking with more anterior areas supporting lexico-semantic processing. These findings suggest that during reading, visual stimuli are first encoded as letters before their combinations are encoded as words. Activity then rapidly spreads anteriorly, and the entire network is engaged in sustained integrative processing.

Theta oscillations are sensitive to both early and late conflict processing stages: effects of alcohol intoxication.

Prior neuroimaging evidence indicates that decision conflict activates medial and lateral prefrontal and parietal cortices. Theoretical accounts of cognitive control highlight anterior cingulate cortex (ACC) as a central node in this network. However, a better understanding of the relative primacy and functional contributions of these areas to decision conflict requires insight into the neural dynamics of successive processing stages including conflict detection, response selection and execution. Moderate alcohol intoxication impairs cognitive control as it interferes with the ability to inhibit dominant, prepotent responses when they are no longer correct. To examine the effects of moderate intoxication on successive processing stages during cognitive control, spatio-temporal changes in total event-related theta power were measured during Stroop-induced conflict. Healthy social drinkers served as their own controls by participating in both alcohol (0.6 g/kg ethanol for men, 0.55 g/kg women) and placebo conditions in a counterbalanced design. Anatomically-constrained magnetoencephalography (aMEG) approach was applied to complex power spectra for theta (4-7 Hz) frequencies. The principal generator of event-related theta power to conflict was estimated to ACC, with contributions from fronto-parietal areas. The ACC was uniquely sensitive to conflict during both early conflict detection, and later response selection and execution stages. Alcohol attenuated theta power to conflict across successive processing stages, suggesting that alcohol-induced deficits in cognitive control may result from theta suppression in the executive network. Slower RTs were associated with attenuated theta power estimated to ACC, indicating that alcohol impairs motor preparation and execution subserved by the ACC. In addition to their relevance for the currently prevailing accounts of cognitive control, our results suggest that alcohol-induced impairment of top-down strategic processing underlies poor self-control and inability to refrain from drinking.

Sequential temporo-fronto-temporal activation during monitoring of the auditory environment for temporal patterns.

Subjects detected rarely occurring shifts between two simple tone-patterns, in a paradigm that dissociated the effects of rarity from those of pitch, habituation, and attention. Whole-head magnetoencephalography suggested that rare attended pattern-shifts evoked activity first in the superior temporal plane (sTp, peak ~100 ms), then superior temporal sulcus (sTs, peak ~130 ms), then posteroventral prefrontal (pvpF, peak ~230 ms), and anterior temporal cortices (aT, peak ~370 ms). Activity was more prominent in the right hemisphere. After subtracting the effects of nonshift tones (balanced for pitch and habituation status), weak but consistent differential effects of pattern-shifts began in aT at 90-130 ms, spread to sTs and sTp at ∼130 ms, then pvpF, and finally returned to aT. Cingulate activity resembled prefrontal. Responses to pattern shifts were greatly attenuated when the same stimuli were ignored, suggesting that the initial superior temporal activity reflected an attention-related mismatch negativity. The prefrontal activity at ~230 ms corresponded in latency and task correlates with simultaneously recorded event-related potential components N2b and P3a; the subsequent temporal activity corresponded to the P3b. These results were confirmed in sensors specific for frontal or temporal cortex, and thus are independent of the inverse method used. Overall, these results suggest that auditory working memory for temporal patterns begins with detection of the pattern change by an interaction of anterior and superior temporal structures, followed by identification of the event and its consequences led by posteroventral prefrontal and cingulate cortices, and finally, definitive encoding of the event in anterior temporal areas.

Multimodal imaging of repetition priming: Using fMRI, MEG, and intracranial EEG to reveal spatiotemporal profiles of word processing.

Repetition priming is a core feature of memory processing whose anatomical correlates remain poorly understood. In this study, we use advanced multimodal imaging (functional magnetic resonance imaging (fMRI) and magnetoencephalography; MEG) to investigate the spatiotemporal profile of repetition priming. We use intracranial electroencephalography (iEEG) to validate our fMRI/MEG measurements. Twelve controls completed a semantic judgment task with fMRI and MEG that included words presented once (new, 'N') and words that repeated (old, 'O'). Six patients with epilepsy completed the same task during iEEG recordings. Blood-oxygen level dependent (BOLD) responses for N vs. O words were examined across the cortical surface and within regions of interest. MEG waveforms for N vs. O words were estimated using a noise-normalized minimum norm solution, and used to interpret the timecourse of fMRI. Spatial concordance was observed between fMRI and MEG repetition effects from 350 to 450 ms within bilateral occipitotemporal and medial temporal, left prefrontal, and left posterior temporal cortex. Additionally, MEG revealed widespread sources within left temporoparietal regions, whereas fMRI revealed bilateral reductions in occipitotemporal and left superior frontal, and increases in inferior parietal, precuneus, and dorsolateral prefrontal activity. BOLD suppression in left posterior temporal, left inferior prefrontal, and right occipitotemporal cortex correlated with MEG repetition-related reductions. IEEG responses from all three regions supported the timecourse of MEG and localization of fMRI. Furthermore, iEEG decreases to repeated words were associated with decreased gamma power in several regions, providing evidence that gamma oscillations are tightly coupled to cognitive phenomena and reflect regional activations seen in the BOLD signal.