Multiscale representations of community structures in attractor neural networks.

Our cognition relies on the ability of the brain to segment hierarchically structured events on multiple scales. Recent evidence suggests that the brain performs this event segmentation based on the structure of state-transition graphs behind sequential experiences. However, the underlying circuit mechanisms are poorly understood. In this paper we propose an extended attractor network model for graph-based hierarchical computation which we call the Laplacian associative memory. This model generates multiscale representations for communities (clusters) of associative links between memory items, and the scale is regulated by the heterogenous modulation of inhibitory circuits. We analytically and numerically show that these representations correspond to graph Laplacian eigenvectors, a popular method for graph segmentation and dimensionality reduction. Finally, we demonstrate that our model exhibits chunked sequential activity patterns resembling hippocampal theta sequences. Our model connects graph theory and attractor dynamics to provide a biologically plausible mechanism for abstraction in the brain.

Extended Temporal Association Memory by Modulations of Inhibitory Circuits.

Hebbian learning of excitatory synapses plays a central role in storing activity patterns in associative memory models. Interstimulus Hebbian learning associates multiple items by converting temporal correlation to spatial correlation between attractors. Growing evidence suggests the importance of inhibitory plasticity in memory processing, but the consequence of such regulation in associative memory has not been understood. Noting that Hebbian learning of inhibitory synapses yields an anti-Hebbian effect, we show that the combination of Hebbian and anti-Hebbian learning can significantly increase the span of temporal association between correlated attractors as well as the sensitivity of these states to external input. Furthermore, these effects are regulated by changing the ratio of local and global recurrent inhibition after learning weights for excitation-inhibition balance. Our results suggest a nontrivial role of plasticity and modulation of inhibitory circuits in associative memory.

Structure of the human M2 muscarinic acetylcholine receptor bound to an antagonist.

The parasympathetic branch of the autonomic nervous system regulates the activity of multiple organ systems. Muscarinic receptors are G-protein-coupled receptors that mediate the response to acetylcholine released from parasympathetic nerves. Their role in the unconscious regulation of organ and central nervous system function makes them potential therapeutic targets for a broad spectrum of diseases. The M2 muscarinic acetylcholine receptor (M2 receptor) is essential for the physiological control of cardiovascular function through activation of G-protein-coupled inwardly rectifying potassium channels, and is of particular interest because of its extensive pharmacological characterization with both orthosteric and allosteric ligands. Here we report the structure of the antagonist-bound human M2 receptor, the first human acetylcholine receptor to be characterized structurally, to our knowledge. The antagonist 3-quinuclidinyl-benzilate binds in the middle of a long aqueous channel extending approximately two-thirds through the membrane. The orthosteric binding pocket is formed by amino acids that are identical in all five muscarinic receptor subtypes, and shares structural homology with other functionally unrelated acetylcholine binding proteins from different species. A layer of tyrosine residues forms an aromatic cap restricting dissociation of the bound ligand. A binding site for allosteric ligands has been mapped to residues at the entrance to the binding pocket near this aromatic cap. The structure of the M2 receptor provides insights into the challenges of developing subtype-selective ligands for muscarinic receptors and their propensity for allosteric regulation.

Temporal relation between neural activity and neurite pruning on a numerical model and a microchannel device with micro electrode array.

Synapse elimination and neurite pruning are essential processes for the formation of neuronal circuits. These regressive events depend on neural activity and occur in the early postnatal days known as the critical period, but what makes this temporal specificity is not well understood. One possibility is that the neural activities during the developmentally regulated shift of action of GABA inhibitory transmission lead to the critical period. Moreover, it has been reported that the shifting action of the inhibitory transmission on immature neurons overlaps with synapse elimination and neurite pruning and that increased inhibitory transmission by drug treatment could induce temporal shift of the critical period. However, the relationship among these phenomena remains unclear because it is difficult to experimentally show how the developmental shift of inhibitory transmission influences neural activities and whether the activities promote synapse elimination and neurite pruning. In this study, we modeled synapse elimination in neuronal circuits using the modified Izhikevich's model with functional shifting of GABAergic transmission. The simulation results show that synaptic pruning within a specified period like the critical period is spontaneously generated as a function of the developmentally shifting inhibitory transmission and that the specific firing rate and increasing synchronization of neural circuits are seen at the initial stage of the critical period. This temporal relationship was experimentally supported by an in vitro primary culture of rat cortical neurons in a microchannel on a multi-electrode array (MEA). The firing rate decreased remarkably between the 18-25 days in vitro (DIV), and following these changes in the firing rate, the neurite density was slightly reduced. Our simulation and experimental results suggest that decreasing neural activity due to developing inhibitory synaptic transmission could induce synapse elimination and neurite pruning at particular time such as the critical period. Additionally, these findings indicate that we can estimate the maturity level of inhibitory transmission and the critical period by measuring the firing rate and the degree of synchronization in engineered neural networks.

Transmembrane topology and oligomeric structure of the high-affinity choline transporter.

The high-affinity choline transporter CHT1 mediates choline uptake essential for acetylcholine synthesis in cholinergic nerve terminals. CHT1 belongs to the Na(+)/glucose cotransporter family (SLC5), which is postulated to have a common 13-transmembrane domain core; however, no direct experimental evidence for CHT1 transmembrane topology has yet been reported. We examined the transmembrane topology of human CHT1 using cysteine-scanning analysis. Single cysteine residues were introduced into the putative extra- and intracellular loops and probed for external accessibility for labeling with a membrane-impermeable, sulfhydryl-specific biotinylating reagent in intact cells expressing these mutants. The results provide experimental evidence for a topological model of a 13-transmembrane domain protein with an extracellular amino terminus and an intracellular carboxyl terminus. We also constructed a three-dimensional homology model of CHT1 based on the crystal structure of the bacterial Na(+)/galactose cotransporter, which supports our conclusion of CHT1 transmembrane topology. Furthermore, we examined whether CHT1 exists as a monomer or oligomer. Chemical cross-linking induces the formation of a higher molecular weight form of CHT1 on the cell surface in HEK293 cells. Two different epitope-tagged CHT1 proteins expressed in the same cells can be co-immunoprecipitated. Moreover, co-expression of an inactive mutant I89A with the wild type induces a dominant-negative effect on the overall choline uptake activity. These results indicate that CHT1 forms a homo-oligomer on the cell surface in cultured cells.

Molecular properties of muscarinic acetylcholine receptors.

Muscarinic acetylcholine receptors, which comprise five subtypes (M1-M5 receptors), are expressed in both the CNS and PNS (particularly the target organs of parasympathetic neurons). M1-M5 receptors are integral membrane proteins with seven transmembrane segments, bind with acetylcholine (ACh) in the extracellular phase, and thereafter interact with and activate GTP-binding regulatory proteins (G proteins) in the intracellular phase: M1, M3, and M5 receptors interact with Gq-type G proteins, and M2 and M4 receptors with Gi/Go-type G proteins. Activated G proteins initiate a number of intracellular signal transduction systems. Agonist-bound muscarinic receptors are phosphorylated by G protein-coupled receptor kinases, which initiate their desensitization through uncoupling from G proteins, receptor internalization, and receptor breakdown (down regulation). Recently the crystal structures of M2 and M3 receptors were determined and are expected to contribute to the development of drugs targeted to muscarinic receptors. This paper summarizes the molecular properties of muscarinic receptors with reference to the historical background and bias to studies performed in our laboratories.

Substrate-induced internalization of the high-affinity choline transporter.

Cholinergic neurons are endowed with a high-affinity choline uptake system for efficient synthesis of acetylcholine at the presynaptic terminals. The high-affinity choline transporter CHT1 is responsible for choline uptake, the rate-limiting step in acetylcholine synthesis. However, endogenous physiological factors that affect CHT1 expression or function and consequently regulate the acetylcholine synthesis rate are essentially unknown. Here we demonstrate that extracellular substrate decreases the cell-surface expression of CHT1 in rat brain synaptosomes, primary cultures from the basal forebrain, and mammalian cell lines transfected with CHT1. Extracellular choline rapidly decreases cell-surface CHT1 expression by accelerating its internalization, a process that is mediated by a dynamin-dependent endocytosis pathway in HEK293 cells. Specific inhibitor hemicholinium-3 decreases the constitutive internalization rate and thereby increases cell-surface CHT1 expression. We also demonstrate that the constitutive internalization of CHT1 depends on extracellular pH in cultured cells. Our results collectively suggest that the internalization of CHT1 is induced by extracellular substrate, providing a novel feedback mechanism for the regulation of acetylcholine synthesis at the cholinergic presynaptic terminals.

Identification of physiologically active substances as novel ligands for MRGPRD.

Mas-related G-protein coupled receptor member D (MRGPRD) is a G protein-coupled receptor (GPCR) which belongs to the Mas-related GPCRs expressed in the dorsal root ganglia (DRG). In this study, we investigated two novel ligands in addition to beta-alanine: (1) beta-aminoisobutyric acid, a physiologically active substance, with which possible relation to tumors has been seen together with beta-alanine; (2) diethylstilbestrol, a synthetic estrogen hormone. In addition to the novel ligands, we found that transfection of MRGPRD leads fibroblast cells to form spheroids, which would be related to oncogenicity. To understand the MRGPRD novel character, oncogenicity, a large chemical library was screened in order to obtain MRGPRD antagonists to utilize in exploring the character. The antagonist in turn inhibited the spheroid proliferation that is dependent on MRGPRD signaling as well as MRGPRD signals activated by beta-alanine. The antagonist, a small-molecule compound we found in this study, is a potential anticancer agent.

Platform for the rapid construction and evaluation of GPCRs for crystallography in Saccharomyces cerevisiae.

BACKGROUND: Recent successes in the determination of G-protein coupled receptor (GPCR) structures have relied on the ability of receptor variants to overcome difficulties in expression and purification. Therefore, the quick screening of functionally expressed stable receptor variants is vital. RESULTS: We developed a platform using Saccharomyces cerevisiae for the rapid construction and evaluation of functional GPCR variants for structural studies. This platform enables us to perform a screening cycle from construction to evaluation of variants within 6-7 days. We firstly confirmed the functional expression of 25 full-length class A GPCRs in this platform. Then, in order to improve the expression level and stability, we generated and evaluated the variants of the four GPCRs (hADRB2, hCHRM2, hHRH1 and hNTSR1). These stabilized receptor variants improved both functional activity and monodispersity. Finally, the expression level of the stabilized hHRH1 in Pichia pastoris was improved up to 65 pmol/mg from negligible expression of the functional full-length receptor in S. cerevisiae at first screening. The stabilized hHRH1 was able to be purified for use in crystallization trials. CONCLUSIONS: We demonstrated that the S. cerevisiae system should serve as an easy-to-handle and rapid platform for the construction and evaluation of GPCR variants. This platform can be a powerful prescreening method to identify a suitable GPCR variant for crystallography.

Evaluation of the Pichia pastoris expression system for the production of GPCRs for structural analysis.

BACKGROUND: Various protein expression systems, such as Escherichia coli (E. coli), Saccharomyces cerevisiae (S. cerevisiae), Pichia pastoris (P. pastoris), insect cells and mammalian cell lines, have been developed for the synthesis of G protein-coupled receptors (GPCRs) for structural studies. Recently, the crystal structures of four recombinant human GPCRs, namely ß2 adrenergic receptor, adenosine A2a receptor, CXCR4 and dopamine D3 receptor, were successfully determined using an insect cell expression system. GPCRs expressed in insect cells are believed to undergo mammalian-like posttranscriptional modifications and have similar functional properties than in mammals. Crystal structures of GPCRs have not yet been solved using yeast expression systems. In the present study, P. pastoris and insect cell expression systems for the human muscarinic acetylcholine receptor M2 subtype (CHRM2) were developed and the quantity and quality of CHRM2 synthesized by both expression systems were compared for the application in structural studies. RESULTS: The ideal conditions for the expression of CHRM2 in P. pastoris were 60 hr at 20°C in a buffer of pH 7.0. The specific activity of the expressed CHRM2 was 28.9 pmol/mg of membrane protein as determined by binding assays using [3H]-quinuclidinyl benzilate (QNB). Although the specific activity of the protein produced by P. pastoris was lower than that of Sf9 insect cells, CHRM2 yield in P. pastoris was 2-fold higher than in Sf9 insect cells because P. pastoris was cultured at high cell density. The dissociation constant (Kd) for QNB in P. pastoris was 101.14 ± 15.07 pM, which was similar to that in Sf9 insect cells (86.23 ± 8.57 pM). There were no differences in the binding affinity of CHRM2 for QNB between P. pastoris and Sf9 insect cells. CONCLUSION: Compared to insect cells, P. pastoris is easier to handle, can be grown at lower cost, and can be expressed quicker at a large scale. Yeast, P. pastoris, and insect cells are all effective expression systems for GPCRs. The results of the present study strongly suggested that protein expression in P. pastoris can be applied to the structural and biochemical studies of GPCRs.

Neural mechanisms for learning hierarchical structures of information.

Spatial and temporal information from the environment is often hierarchically organized, so is our knowledge formed about the environment. Identifying the meaningful segments embedded in hierarchically structured information is crucial for cognitive functions, including visual, auditory, motor, memory, and language processing. Segmentation enables the grasping of the links between isolated entities, offering the basis for reasoning and thinking. Importantly, the brain learns such segmentation without external instructions. Here, we review the underlying computational mechanisms implemented at the single-cell and network levels. The network-level mechanism has an interesting similarity to machine-learning methods for graph segmentation. The brain possibly implements methods for the analysis of the hierarchical structures of the environment at multiple levels of its processing hierarchy.

Differential rate constants of racemization of aspartyl and asparaginyl residues in human alpha A-crystallin mutants.

Asp58 and Asp151 in alpha A-crystallin of human eye lenses become highly inverted and isomerized to d-beta-Asp residues with age. Racemization was previously shown to proceed rapidly when the residue on the carboxyl side of the Asp residue is small. Asn was also demonstrated to be more susceptible to racemization than Asp in protein. In this study, the changes of rate constants for racemization at Asp58 and Asp151 and at Asn58 and Asn151 were investigated using D58N, S59T, D151N and A152V mutants obtained through site-directed mutagenesis. The rate constant of racemization at Asn151 in D151N was found to be 1.5 times more rapid than Asp151 in the wild-type. For A152V, the rate constant at Asp151 was 1/4 that of the wild-type. There were no significant differences in the rate constants of racemization for both Asp58 and Asn58 residues. The aggregate size of D58N, S59T and D151N mutants increased or increased in polydispersity and their chaperone activities decreased. The size and chaperone activity of A152V was unchanged. These results suggest that structures close to Asp58 and Asp151 residues in the protein affect the rate constant of Asp racemization and the size and chaperone function of alpha A-crystallin.

A highly conserved tryptophan residue in the fourth transmembrane domain of the A adenosine receptor is essential for ligand binding but not receptor homodimerization.

Dimerization between G protein-coupled receptors (GPCRs) is a clearly established phenomenon. However, limited information is currently available on the interface essential for this process. Based on structural comparisons and sequence homology between rhodopsin and A(1) adenosine receptor (A(1)R), we initially hypothesized that four residues in transmembrane (TM) 4 and TM5 are involved in A(1)R homodimerization. Accordingly, these residues were substituted with Ala by site-directed mutagenesis. Interestingly, the mutant protein displayed no significant decrease in homodimer formation compared with wild-type A(1)R, as evident from coimmunoprecipitation and BRET(2) analyses (improved bioluminescence resonance energy transfer system offered by Perkin-Elmer Life Sciences), but lost ligand binding activity almost completely. Further studies disclosed that this effect was derived from the mutation of one particular residue, Trp132, which is highly conserved among many GPCRs. Confocal immunofluorescence and cell-surface biotinylation studies revealed that the mutant receptors localized normally at transfected cell membranes, signifying that loss of ligand binding was not because of defective cellular trafficking. Molecular modeling of the A(1)R-ligand complex disclosed that Trp132 interacted with several residues located in TM3 and TM5 that stabilized agonist binding. Thus, loss of interactions of Trp with these residues may, in turn, disrupt binding to agonists. Our study provides strong evidence of the essential role of the highly conserved Trp132 in TM4 of adenosine receptors.

Comparison of functional non-glycosylated GPCRs expression in Pichia pastoris.

N-linked glycosylation is the most common post-translational modification of G-protein-coupled receptors (GPCRs) and is correlated to the localization and function of the receptors depending on each receptor. However, heterogeneity of glycosylation can interfere with protein crystallization. The removal of N-linked glycosylation from membrane proteins improves the ability to crystallize these proteins. We screened 25 non-glycosylated GPCRs for functional receptor production in the methylotrophic yeast Pichia pastoris using specific ligand-receptor binding assays. We found that five clones were expressed at greater than 10 pmol/mg, 9 clones at 1-10 pmol/mg and 11 clones at less than 1 pmol/mg of membrane protein. Further optimization of culture parameters including culture scale, induction time, pH and temperature enabled us to achieve expression of a functional human muscarinic acetylcholine receptor subtype 2 (CHRM2) with a B(max) value of 51.2 pmol/mg of membrane protein. Approximately 1.9 mg of the human CHRM2 was produced from a 1-L culture.

Acetylcholine synthesis and release in NIH3T3 cells coexpressing the high-affinity choline transporter and choline acetyltransferase.

Acetylcholine (ACh) is known to be a key neurotransmitter in the central and peripheral nervous systems, but it is also produced in a variety of non-neuronal tissues and cells, including lymphocytes, placenta, amniotic membrane, vascular endothelial cells, keratinocytes, and epithelial cells in the digestive and respiratory tracts. To investigate contribution made by the high-affinity choline transporter (CHT1) to ACh synthesis in both cholinergic neurons and nonneuronal cells, we transfected rat CHT1 cDNA into NIH3T3ChAT cells, a mouse fibroblast line expressing mouse choline acetyltransferase (ChAT), to establish the NIH3T3ChAT 112-1 cell line, which stably expresses both CHT1 and ChAT. NIH3T3ChAT 112-1 cells showed increased binding of the CHT1 inhibitor [(3)H]hemicholinium-3 (HC-3) and greater [(3)H]choline uptake and ACh synthesis than NIH3T3ChAT 103-1 cells, a CHT1-negative control cell line. HC-3 significantly inhibited ACh synthesis in NIH3T3ChAT 112-1 cells but did not affect synthesis in NIH3T3ChAT 103-1 cells. ACh synthesis in NIH3T3ChAT 112-1 cells was also reduced by amiloride, an inhibitor of organic cation transporters (OCTs) involved in low-affinity choline uptake, and by procaine and lidocaine, two local anesthetics that inhibit plasma membrane phospholipid metabolism. These results suggest that CHT1 plays a key role in ACh synthesis in NIH3T3ChAT 112-1 cells and that choline taken up by OCTs or derived from the plasma membrane is also utilized for ACh synthesis in both cholinergic neurons and nonneuronal cholinergic cells, such as lymphocytes.

Unsupervised Detection of Cell-Assembly Sequences by Similarity-Based Clustering.

Neurons which fire in a fixed temporal pattern (i.e., "cell assemblies") are hypothesized to be a fundamental unit of neural information processing. Several methods are available for the detection of cell assemblies without a time structure. However, the systematic detection of cell assemblies with time structure has been challenging, especially in large datasets, due to the lack of efficient methods for handling the time structure. Here, we show a method to detect a variety of cell-assembly activity patterns, recurring in noisy neural population activities at multiple timescales. The key innovation is the use of a computer science method to comparing strings ("edit similarity"), to group spikes into assemblies. We validated the method using artificial data and experimental data, which were previously recorded from the hippocampus of male Long-Evans rats and the prefrontal cortex of male Brown Norway/Fisher hybrid rats. From the hippocampus, we could simultaneously extract place-cell sequences occurring on different timescales during navigation and awake replay. From the prefrontal cortex, we could discover multiple spike sequences of neurons encoding different segments of a goal-directed task. Unlike conventional event-driven statistical approaches, our method detects cell assemblies without creating event-locked averages. Thus, the method offers a novel analytical tool for deciphering the neural code during arbitrary behavioral and mental processes.

Muscarinic M4 receptor recycling requires a motif in the third intracellular loop.

The present study was performed to identify sequence(s) in the third intracellular loop (i3) of the muscarinic acetylcholine receptor M4 subtype (M4 receptor) involved in its internalization and recycling. In transiently transfected human embryonic kidney 293-tsA201 cells, 40 to 50% of cell-surface M4 receptors are internalized in an agonist-dependent manner, and approximately 65% of internalized receptors are recycled back to the cell surface after removal of the agonist. We examined the internalization and recycling of M4 receptor mutants with partial deletion in i3 and found that various mutants (M4del-K(235)-K(240), M4del-T(241)-K(271), and M4del-W(339)-N(372)) showed internalization and cell-surface recycling in a similar manner to the M4 receptor. We also found that the mutant M4del-L(272)-R(338) was internalized to only half the extent of the M4 receptor and was recycled after agonist removal, and the mutant M4del-V(373)-A(393) was also internalized to half the extent of the wild type but was not recycled back to the cell surface after agonist removal. When the sequence corresponding to Val(373)-Ala(393) was grafted onto the i3 portion of a recycling-negative mutant of muscarinic M2 receptor with deletion of almost the whole of the i3 sequence, approximately 40% of the chimeric receptor on the cell surface was internalized, and more than 65% of the internalized receptors were recycled back to the cell surface. These results indicate that the regions including Leu(272)-Arg(338) and Val(373)-Ala(393) are involved in internalization of the M4 receptor, and the region including Val(373)-Ala(393) is indispensable for its recycling, whereas the other regions of i3 are dispensable for internalization and recycling.

Dendritic processing of spontaneous neuronal sequences for single-trial learning.

Spontaneous firing sequences are ubiquitous in cortical networks, but their roles in cellular and network-level computations remain unexplored. In the hippocampus, such sequences, conventionally called preplay, have been hypothesized to participate in learning and memory. Here, we present a computational model for encoding input sequence patterns into internal network states based on the propagation of preplay sequences in recurrent neuronal networks. The model instantiates two synaptic pathways in cortical neurons, one for proximal dendrite-somatic interactions to generate intrinsic preplay sequences and the other for distal dendritic processing of extrinsic signals. The core dendritic computation is the maximization of matching between patterned activities in the two compartments through nonlinear spike generation. The model performs robust single-trial learning with long-term stability and independence that are modulated by the plasticity of dendrite-targeted inhibition. Our results demonstrate that dendritic computation enables somatic spontaneous firing sequences to act as templates for rapid and stable memory formation.

Recurrent network model for learning goal-directed sequences through reverse replay.

Reverse replay of hippocampal place cells occurs frequently at rewarded locations, suggesting its contribution to goal-directed path learning. Symmetric spike-timing dependent plasticity (STDP) in CA3 likely potentiates recurrent synapses for both forward (start to goal) and reverse (goal to start) replays during sequential activation of place cells. However, how reverse replay selectively strengthens forward synaptic pathway is unclear. Here, we show computationally that firing sequences bias synaptic transmissions to the opposite direction of propagation under symmetric STDP in the co-presence of short-term synaptic depression or afterdepolarization. We demonstrate that significant biases are created in biologically realistic simulation settings, and this bias enables reverse replay to enhance goal-directed spatial memory on a W-maze. Further, we show that essentially the same mechanism works in a two-dimensional open field. Our model for the first time provides the mechanistic account for the way reverse replay contributes to hippocampal sequence learning for reward-seeking spatial navigation.

Neuronal stability in medial frontal cortex sets individual variability in decision-making.

In the brain, decision making is instantiated in dedicated neural circuits. However, there is considerable individual variability in decision-making behavior, particularly under uncertainty. The origins of decision variability within these conserved neural circuits are not known. Here we demonstrate in the rat medial frontal cortex (MFC) that individual variability is a consequence of altered stability in neuronal populations. In a sensory-guided choice task, rats trained on familiar stimuli were exposed to unfamiliar stimuli, resulting in variable choice responses across individuals. We created a recurrent network model to examine the source of variability in MFC neurons, and found that the landscape of neural population trajectories explained choice variability across different unfamiliar stimuli. We experimentally confirmed model predictions showing that trial-by-trial variability in neuronal activity indexes the landscape and predicts individual variation. These results show that neural stability is a critical component of the MFC neural dynamics that underpins individual variation in decision-making.