Hundreds of millions of structured proteins sustain life through chemical interactions and catalytic reactions1. Though dynamic, these proteins are assumed to be built upon fixed scaffolds of secondary structure, α-helices and β-sheets. Experimentally determined structures of over >58,000 non-redundant proteins support this assumption, though it has recently been challenged by ∼100 fold-switching proteins2. These “metamorphic3” proteins, though ostensibly rare, raise the question of how many uncharacterized proteins have shapeshifting–rather than fixed–secondary structures. To address this question, we developed a comparative sequence-based approach that predicts fold-switching proteins from differences in secondary structure propensity. We applied this approach to the universally conserved NusG transcription factor family of ∼15,000 proteins, one of which has a 50-residue regulatory subunit experimentally shown to switch between α-helical and β-sheet folds4. Our approach predicted that 25% of the sequences in this family undergo similar α-helix ⇌ β-sheet transitions, a frequency two orders of magnitude larger than previously observed. Our predictions evade state-of-the-art computational methods but were confirmed experimentally by circular dichroism and nuclear magnetic resonance spectroscopy for all 10 assiduously chosen dissimilar variants. These results suggest that fold switching is a pervasive mechanism of transcriptional regulation in all kingdoms of life and imply that numerous uncharacterized proteins may also switch folds.

Wistar-Kyoto (WKY) rats as an endogenous depression model partially lack a response to classic selective serotonin reuptake inhibitors (SSRIs). Thus, this strain has the potential to be established as a model of treatment-resistant depression (TRD). However, the SSRI resistance in WKY rats is still not fully understood. In this study, WKY and control rats were subjected to a series of tests, namely, a forced swim test (FST), a sucrose preference test (SPT), and an open field test (OFT), and were scanned in a 7.0-T MRI scanner before and after three-week citalopram or saline administration. Behavioral results demonstrated that WKY rats had increased immobility in the FST and decreased sucrose preference in the SPT and central time spent in the OFT. However, citalopram did not improve immobility in the FST. The amplitude of low-frequency fluctuation (ALFF) analysis showed regional changes in the striatum and hippocampus of WKY rats. However, citalopram partially reversed the ALFF value in the dorsal part of the two regions. Functional connectivity (FC) analysis showed that FC strengths were decreased in WKY rats compared with controls. Nevertheless, citalopram partially increased FC strengths in WKY rats. Based on FC, global graph analysis demonstrated decreased network efficiency in WKY + saline group compared with control + saline group, but citalopram showed weak network efficiency improvement. In conclusion, resting-state fMRI results implied widely affected brain function at both regional and global levels in WKY rats. Citalopram had only partial effects on these functional changes, indicating a potential treatment resistance mechanism.

The PV2 (Celio 1990), a cluster of parvalbumin-positive neurons located in the ventromedial region of the distal periaqueductal gray (PAG) has not been previously described as its own entity, leading us to study its extent, connections, and gene expression. It is an oval, bilateral, elongated cluster composed of approximately 475 parvalbumin-expressing neurons in a single mouse hemisphere. In its anterior portion it impinges upon the paratrochlear nucleus (Par4) and in its distal portion it is harbored in the posterodorsal raphe nucleus (PDR). It is known to receive inputs from the orbitofrontal cortex and from the parvafox nucleus in the ventrolateral hypothalamus. Using anterograde tracing methods in parvalbumin-Cre mice, the main projections of the PV2 cluster innervate the supraoculomotor periaqueductal gray (Su3) of the PAG, the parvafox nucleus of the lateral hypothalamus, the gemini nuclei of the posterior hypothalamus, the septal regions, and the diagonal band in the forebrain, as well as various nuclei within the reticular formation in the midbrain and brainstem. Within the brainstem, projections were discrete, but involved areas implicated in autonomic control. The PV2 cluster expressed various peptides and receptors, including the receptor for Adcyap1, a peptide secreted by one of its main afferences, namely, the parvafox nucleus. The expression of GAD1 and GAD2 in the region of the PV2, the presence of Vgat-1 in a subpopulation of PV2-neurons as well as the coexistence of GAD67 immunoreactivity with parvalbumin in terminal endings indicates the inhibitory nature of a subpopulation of PV2-neurons. The PV2 cluster may be part of a feedback controlling the activity of the hypothalamic parvafox and the Su3 nuclei in the periaqueductal gray.

Neural computation in biological and artificial networks relies on nonlinear synaptic integration. The structural connectivity matrix of synaptic weights between neurons is a critical determinant of overall network function. However, quantitative links between neural network structure and function are complex and subtle. For example, many networks can give rise to similar functional responses, and the same network can function differently depending on context. Whether certain patterns of synaptic connectivity are required to generate specific network-level computations is largely unknown. Here we introduce a geometric framework for identifying synaptic connections required by steady-state responses in recurrent networks of rectified-linear neurons. Assuming that the number of specified response patterns does not exceed the number of input synapses, we analytically calculate all feedforward and recurrent connectivity matrices that can generate the specified responses from the network inputs. We then use this analytical characterization to rigorously analyze the solution space geometry and derive certainty conditions guaranteeing a non-zero synapse between neurons. Numerical simulations of feedforward and recurrent networks verify our analytical results. Our theoretical framework could be applied to neural activity data to make anatomical predictions that follow generally from the model architecture. It thus provides novel opportunities for discerning what model features are required to accurately relate neural network structure and function.

Task-related neural activity is widespread across populations of neurons during goal-directed behaviors. However, little is known about the synaptic reorganization and circuit mechanisms that lead to broad activity changes. Here we trained a limited subset of neurons in a spiking network with strong synaptic interactions to reproduce the activity of neurons in the motor cortex during a decision-making task. We found that task-related activity, resembling the neural data, emerged across the network, even in the untrained neurons. Analysis of trained networks showed that strong untrained synapses, which were independent of the task and determined the dynamical state of the network, mediated the spread of task-related activity. Optogenetic perturbations suggest that the motor cortex is strongly-coupled, supporting the applicability of the mechanism to cortical networks. Our results reveal a cortical mechanism that facilitates distributed representations of task-variables by spreading the activity from a subset of plastic neurons to the entire network through task-independent strong synapses.

Small-molecule fluorescent stains enable the imaging of cellular structures without the need for genetic manipulation. Here, we introduce 2,7-diaminobenzopyrylium (DAB) dyes as live-cell mitochondrial stains excited with violet light. This amalgam of the coumarin and rhodamine fluorophore structures yields dyes with high photostability and tunable spectral properties.

Insects like Drosophila produce a second brain adapted to the form and behavior of a larva. Neurons for both larval and adult brains are produced by the same stem cells (neuroblasts) but the larva possesses only the earliest born neurons produced from each. To understand how a functional larval brain is made from this reduced set of neurons, we examined the origins and metamorphic fates of the neurons of the larval and adult mushroom body circuits. The adult mushroom body core is built sequentially of γ Kenyon cells, that form a medial lobe, followed by α’β’, and αβ Kenyon cells that form additional medial lobes and two vertical lobes. Extrinsic input (MBINs) and output (MBONs) neurons divide this core into computational compartments. The larval mushroom body contains only γ neurons. Its medial lobe compartments are roughly homologous to those of the adult and same MBONs are used for both. The larval vertical lobe, however, is an analogous “facsimile” that uses a larval-specific branch on the γ neurons to make up for the missing α’β’, and αβ neurons. The extrinsic cells for the facsimile are early-born neurons that trans-differentiate to serve a mushroom body function in the larva and then shift to other brain circuits in the adult. These findings are discussed in the context of the evolution of a larval brain in insects with complete metamorphosis.

Signaling through the TNF-family receptor Fas/CD95 can trigger apoptosis or non-apoptotic cellular responses and is essential for protection from autoimmunity. Receptor clustering has been observed following interaction with Fas ligand (FasL), but the stoichiometry of Fas, particularly when triggered by membrane-bound FasL, the only form of FasL competent at inducing programmed cell death, is not known. Here we used super-resolution microscopy to study the behavior of single molecules of Fas/CD95 on the plasma membrane after interaction of Fas with FasL on planar lipid bilayers. We observed rapid formation of Fas protein superclusters containing more than 20 receptors after interactions with membrane-bound FasL. Fluorescence correlation imaging demonstrated recruitment of FADD dependent on an intact Fas death domain, with lipid raft association playing a secondary role. Flow-cytometric FRET analysis confirmed these results, and also showed that some Fas clustering can occur in the absence of FADD and caspase-8. Point mutations in the Fas death domain associated with autoimmune lymphoproliferative syndrome (ALPS) completely disrupted Fas reorganization and FADD recruitment, confirming structure-based predictions of the critical role that these residues play in Fas-Fas and Fas-FADD interactions. Finally, we showed that induction of apoptosis correlated with the ability to form superclusters and recruit FADD.

Summary Many embryonic organs undergo epithelial morphogenesis to form tree-like hierarchical structures. However, it remains unclear what drives the budding and branching of stratified epithelia, such as in the embryonic salivary gland and pancreas. Here, we performed live-organ imaging of mouse embryonic salivary glands at single-cell resolution to reveal that budding morphogenesis is driven by expansion and folding of a distinct epithelial surface cell sheet characterized by strong cell-matrix adhesions and weak cell-cell adhesions. Profiling of single-cell transcriptomes of this epithelium revealed spatial patterns of transcription underlying these cell adhesion differences. We then synthetically reconstituted budding morphogenesis by experimentally suppressing E-cadherin expression and inducing basement membrane formation in 3D spheroid cultures of engineered cells, which required β1-integrin-mediated cell-matrix adhesion for successful budding. Thus, stratified epithelial budding, the key first step of branching morphogenesis, is driven by an overall combination of strong cell-matrix adhesion and weak cell-cell adhesion by peripheral epithelial cells.

Dysfunctional sociability is a core symptom in autism spectrum disorder (ASD) that may arise from neural-network dysconnectivity between multiple brain regions. However, pathogenic neural-network mechanisms underlying social dysfunction are largely unknown. Here, we demonstrate that circuit-selective mutation (ctMUT) of ASD-risk Shank3 gene within a unidirectional projection from the prefrontal cortex to the basolateral amygdala alters spine morphology and excitatory-inhibitory balance of the circuit. Shank3 ctMUT mice show reduced sociability as well as elevated neural activity and its amplitude variability, which is consistent with the neuroimaging results from human ASD patients. Moreover, the circuit hyper-activity disrupts the temporal correlation of socially tuned neurons to the events of social interactions. Finally, optogenetic circuit activation in wild-type mice partially recapitulates the reduced sociability of Shank3 ctMUT mice, while circuit inhibition in Shank3 ctMUT mice partially rescues social behavior. Collectively, these results highlight a circuit-level pathogenic mechanism of Shank3 mutation that drives social dysfunction.