Functionally Distinct Neuronal Ensembles within the Memory Engram.

Memories are believed to be encoded by sparse ensembles of neurons in the brain. However, it remains unclear whether there is functional heterogeneity within individual memory engrams, i.e., if separate neuronal subpopulations encode distinct aspects of the memory and drive memory expression differently. Here, we show that contextual fear memory engrams in the mouse dentate gyrus contain functionally distinct neuronal ensembles, genetically defined by the Fos- or Npas4-dependent transcriptional pathways. The Fos-dependent ensemble promotes memory generalization and receives enhanced excitatory synaptic inputs from the medial entorhinal cortex, which we find itself also mediates generalization. The Npas4-dependent ensemble promotes memory discrimination and receives enhanced inhibitory drive from local cholecystokinin-expressing interneurons, the activity of which is required for discrimination. Our study provides causal evidence for functional heterogeneity within the memory engram and reveals synaptic and circuit mechanisms used by each ensemble to regulate the memory discrimination-generalization balance.

Controlling Visually Guided Behavior by Holographic Recalling of Cortical Ensembles.

Neurons in cortical circuits are often coactivated as ensembles, yet it is unclear whether ensembles play a functional role in behavior. Some ensemble neurons have pattern completion properties, triggering the entire ensemble when activated. Using two-photon holographic optogenetics in mouse primary visual cortex, we tested whether recalling ensembles by activating pattern completion neurons alters behavioral performance in a visual task. Disruption of behaviorally relevant ensembles by activation of non-selective neurons decreased performance, whereas activation of only two pattern completion neurons from behaviorally relevant ensembles improved performance, by reliably recalling the whole ensemble. Also, inappropriate behavioral choices were evoked by the mistaken activation of behaviorally relevant ensembles. Finally, in absence of visual stimuli, optogenetic activation of two pattern completion neurons could trigger behaviorally relevant ensembles and correct behavioral responses. Our results demonstrate a causal role of neuronal ensembles in a visually guided behavior and suggest that ensembles implement internal representations of perceptual states.

Molecular determinants of protein evolvability.

The plethora of biological functions that sustain life is rooted in the remarkable evolvability of proteins. An emerging view highlights the importance of a protein's initial state in dictating evolutionary success. A deeper comprehension of the mechanisms that govern the evolvability of these initial states can provide invaluable insights into protein evolution. In this review, we describe several molecular determinants of protein evolvability, unveiled by experimental evolution and ancestral sequence reconstruction studies. We further discuss how genetic variation and epistasis can promote or constrain functional innovation and suggest putative underlying mechanisms. By establishing a clear framework for these determinants, we provide potential indicators enabling the forecast of suitable evolutionary starting points and delineate molecular mechanisms in need of deeper exploration.

Modeling Flexible Protein Structure With AlphaFold2 and Crosslinking Mass Spectrometry.

We propose a pipeline that combines AlphaFold2 (AF2) and crosslinking mass spectrometry (XL-MS) to model the structure of proteins with multiple conformations. The pipeline consists of two main steps: ensemble generation using AF2 and conformer selection using XL-MS data. For conformer selection, we developed two scores-the monolink probability score (MP) and the crosslink probability score (XLP)-both of which are based on residue depth from the protein surface. We benchmarked MP and XLP on a large dataset of decoy protein structures and showed that our scores outperform previously developed scores. We then tested our methodology on three proteins having an open and closed conformation in the Protein Data Bank: Complement component 3 (C3), luciferase, and glutamine-binding periplasmic protein, first generating ensembles using AF2, which were then screened for the open and closed conformations using experimental XL-MS data. In five out of six cases, the most accurate model within the AF2 ensembles-or a conformation within 1 Å of this model-was identified using crosslinks, as assessed through the XLP score. In the remaining case, only the monolinks (assessed through the MP score) successfully identified the open conformation of glutamine-binding periplasmic protein, and these results were further improved by including the "occupancy" of the monolinks. This serves as a compelling proof-of-concept for the effectiveness of monolinks. In contrast, the AF2 assessment score was only able to identify the most accurate conformation in two out of six cases. Our results highlight the complementarity of AF2 with experimental methods like XL-MS, with the MP and XLP scores providing reliable metrics to assess the quality of the predicted models. The MP and XLP scoring functions mentioned above are available at https://gitlab.com/topf-lab/xlms-tools.

Enhancing the activity of Pd ensembles on graphene by manipulating coordination environment.

Atomic dispersion of metal catalysts on a substrate accounts for the increased atomic efficiency of single-atom catalysts (SACs) in various catalytic schemes compared to the nanoparticle counterparts. However, lacking neighboring metal sites has been shown to deteriorate the catalytic performance of SACs in a few industrially important reactions, such as dehalogenation, CO oxidation, and hydrogenation. Metal ensemble catalysts (Mn), an extended concept to SACs, have emerged as a promising alternative to overcome such limitation. Inspired by the fact that the performance of fully isolated SACs can be enhanced by tailoring their coordination environment (CE), we here evaluate whether the CE of Mn can also be manipulated in order to enhance their catalytic activity. We synthesized a set of Pd ensembles (Pdn) on doped graphene supports (Pdn/X-graphene where X = O, S, B, and N). We found that introducing S and N onto oxidized graphene modifies the first shell of Pdn converting Pd-O to Pd-S and Pd-N, respectively. We further found that the B dopant significantly affected the electronic structure of Pdn by serving as an electron donor in the second shell. We examined the performance of Pdn/X-graphene toward selective reductive catalysis, such as bromate reduction, brominated organic hydrogenation, and aqueous-phase CO2 reduction. We observed that Pdn/N-graphene exhibited superior performance by lowering the activation energy of the rate-limiting step, i.e., H2 dissociation into atomic hydrogen. The results collectively suggest controlling the CE of SACs in an ensemble configuration is a viable strategy to optimize and enhance their catalytic performance.

Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology.

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

The Neural Dynamics of Face Ensemble and Central Face Processing.

Extensive work has investigated the neural processing of single faces, including the role of shape and surface properties. However, much less is known about the neural basis of face ensemble perception (e.g., simultaneously viewing several faces in a crowd). Importantly, the contribution of shape and surface properties have not been elucidated in face ensemble processing. Furthermore, how single central faces are processed within the context of an ensemble remains unclear. Here, we probe the neural dynamics of ensemble representation using pattern analyses as applied to electrophysiology data in healthy adults (seven males, nine females). Our investigation relies on a unique set of stimuli, depicting different facial identities, which vary parametrically and independently along their shape and surface properties. These stimuli were organized into ensemble displays consisting of six surround faces arranged in a circle around one central face. Overall, our results indicate that both shape and surface properties play a significant role in face ensemble encoding, with the latter demonstrating a more pronounced contribution. Importantly, we find that the neural processing of the center face precedes that of the surround faces in an ensemble. Further, the temporal profile of center face decoding is similar to that of single faces, while those of single faces and face ensembles diverge extensively from each other. Thus, our work capitalizes on a new center-surround paradigm to elucidate the neural dynamics of ensemble processing and the information that underpins it. Critically, our results serve to bridge the study of single and ensemble face perception.

EnGens: a computational framework for generation and analysis of representative protein conformational ensembles.

Proteins are dynamic macromolecules that perform vital functions in cells. A protein structure determines its function, but this structure is not static, as proteins change their conformation to achieve various functions. Understanding the conformational landscapes of proteins is essential to understand their mechanism of action. Sets of carefully chosen conformations can summarize such complex landscapes and provide better insights into protein function than single conformations. We refer to these sets as representative conformational ensembles. Recent advances in computational methods have led to an increase in the number of available structural datasets spanning conformational landscapes. However, extracting representative conformational ensembles from such datasets is not an easy task and many methods have been developed to tackle it. Our new approach, EnGens (short for ensemble generation), collects these methods into a unified framework for generating and analyzing representative protein conformational ensembles. In this work, we: (1) provide an overview of existing methods and tools for representative protein structural ensemble generation and analysis; (2) unify existing approaches in an open-source Python package, and a portable Docker image, providing interactive visualizations within a Jupyter Notebook pipeline; (3) test our pipeline on a few canonical examples from the literature. Representative ensembles produced by EnGens can be used for many downstream tasks such as protein-ligand ensemble docking, Markov state modeling of protein dynamics and analysis of the effect of single-point mutations.

Twenty-first century hydroclimate: A continually changing baseline, with more frequent extremes.

SignificanceTwenty-first century trends in hydroclimate are so large that future average conditions will, in most cases, fall into the range of what we would today consider extreme drought or pluvial states. Using large climate model ensembles, we remove the background trend and find that the risk of droughts and pluvials relative to that (changing) baseline is fairly similar to the 20th century risk. By continually adapting to long-term background changes, these risks could therefore perhaps be minimized. However, increases in the frequency of extremely wet and dry years are still present even after removing the trend, indicating that sustainably managing hydroclimate-driven risks in a warmer world will face increasingly difficult challenges.

Differential Effects of Astrocyte Manipulations on Learned Motor Behavior and Neuronal Ensembles in the Motor Cortex.

Although motor cortex is crucial for learning precise and reliable movements, whether and how astrocytes contribute to its plasticity and function during motor learning is unknown. Here, we report that astrocyte-specific manipulations in primary motor cortex (M1) during a lever push task alter motor learning and execution, as well as the underlying neuronal population coding. Mice that express decreased levels of the astrocyte glutamate transporter 1 (GLT1) show impaired and variable movement trajectories, whereas mice with increased astrocyte Gq signaling show decreased performance rates, delayed response times, and impaired trajectories. In both groups, which include male and female mice, M1 neurons have altered interneuronal correlations and impaired population representations of task parameters, including response time and movement trajectories. RNA sequencing further supports a role for M1 astrocytes in motor learning and shows changes in astrocytic expression of glutamate transporter genes, GABA transporter genes, and extracellular matrix protein genes in mice that have acquired this learned behavior. Thus, astrocytes coordinate M1 neuronal activity during motor learning, and our results suggest that this contributes to learned movement execution and dexterity through mechanisms that include regulation of neurotransmitter transport and calcium signaling.SIGNIFICANCE STATEMENT We demonstrate for the first time that in the M1 of mice, astrocyte function is critical for coordinating neuronal population activity during motor learning. We demonstrate that knockdown of astrocyte glutamate transporter GLT1 affects specific components of learning, such as smooth trajectory formation. Altering astrocyte calcium signaling by activation of Gq-DREADD upregulates GLT1 and affects other components of learning, such as response rates and reaction times as well as trajectory smoothness. In both manipulations, neuronal activity in motor cortex is dysregulated, but in different ways. Thus, astrocytes have a crucial role in motor learning via their influence on motor cortex neurons, and they do so by mechanisms that include regulation of glutamate transport and calcium signals.

Creation of Neuronal Ensembles and Cell-Specific Homeostatic Plasticity through Chronic Sparse Optogenetic Stimulation.

Cortical computations emerge from the dynamics of neurons embedded in complex cortical circuits. Within these circuits, neuronal ensembles, which represent subnetworks with shared functional connectivity, emerge in an experience-dependent manner. Here we induced ensembles in ex vivo cortical circuits from mice of either sex by differentially activating subpopulations through chronic optogenetic stimulation. We observed a decrease in voltage correlation, and importantly a synaptic decoupling between the stimulated and nonstimulated populations. We also observed a decrease in firing rate during Up-states in the stimulated population. These ensemble-specific changes were accompanied by decreases in intrinsic excitability in the stimulated population, and a decrease in connectivity between stimulated and nonstimulated pyramidal neurons. By incorporating the empirically observed changes in intrinsic excitability and connectivity into a spiking neural network model, we were able to demonstrate that changes in both intrinsic excitability and connectivity accounted for the decreased firing rate, but only changes in connectivity accounted for the observed decorrelation. Our findings help ascertain the mechanisms underlying the ability of chronic patterned stimulation to create ensembles within cortical circuits and, importantly, show that while Up-states are a global network-wide phenomenon, functionally distinct ensembles can preserve their identity during Up-states through differential firing rates and correlations.SIGNIFICANCE STATEMENT The connectivity and activity patterns of local cortical circuits are shaped by experience. This experience-dependent reorganization of cortical circuits is driven by complex interactions between different local learning rules, external input, and reciprocal feedback between many distinct brain areas. Here we used an ex vivo approach to demonstrate how simple forms of chronic external stimulation can shape local cortical circuits in terms of their correlated activity and functional connectivity. The absence of feedback between different brain areas and full control of external input allowed for a tractable system to study the underlying mechanisms and development of a computational model. Results show that differential stimulation of subpopulations of neurons significantly reshapes cortical circuits and forms subnetworks referred to as neuronal ensembles.

Neuronal ensembles in memory processes.

A neuronal ensemble represents the concomitant activity of a specific group of neurons that could encompass a broad repertoire of brain functions such as motor, perceptual, memory or cognitive states. On the other hand, a memory engram portrays the physical manifestation of memory or the changes that enable learning and retrieval. Engram studies focused for many years on finding where memories are stored as in, which cells or brain regions represent a memory trace, and disregarded the investigation of how neuronal activity patterns give rise to such memories. Recent experiments suggest that the association and reactivation of specific neuronal groups could be the main mechanism underlying the brain's ability to remember past experiences and envision future actions. Thus, the growing consensus is that the interaction between neuronal ensembles could allow sequential activity patterns to become memories and recurrent memories to compose complex behaviors. The goal of this review is to propose how the neuronal ensemble framework could be translated and useful to understand memory processes.

Extended N-centered ensemble density functional theory of double electronic excitations.

A recent work (arXiv:2401.04685) has merged N-centered ensembles of neutral and charged electronic ground states with ensembles of neutral ground and excited states, thus providing a general and in-principle exact (so-called extended N-centered) ensemble density functional theory of neutral and charged electronic excitations. This formalism made it possible to revisit the concept of density-functional derivative discontinuity, in the particular case of single excitations from the highest occupied Kohn-Sham (KS) molecular orbital, without invoking the usual "asymptotic behavior of the density" argument. In this work, we address a broader class of excitations, with a particular focus on double excitations. An exact implementation of the theory is presented for the two-electron Hubbard dimer model. A thorough comparison of the true physical ground- and excited-state electronic structures with that of the fictitious ensemble density-functional KS system is also presented. Depending on the choice of the density-functional ensemble as well as the asymmetry of the dimer and the correlation strength, an inversion of states can be observed. In some other cases, the strong mixture of KS states within the true physical system makes the assignment "single excitation" or "double excitation" irrelevant.

MutateX: an automated pipeline for in silico saturation mutagenesis of protein structures and structural ensembles.

Mutations, which result in amino acid substitutions, influence the stability of proteins and their binding to biomolecules. A molecular understanding of the effects of protein mutations is both of biotechnological and medical relevance. Empirical free energy functions that quickly estimate the free energy change upon mutation (ΔΔG) can be exploited for systematic screenings of proteins and protein complexes. In silico saturation mutagenesis can guide the design of new experiments or rationalize the consequences of known mutations. Often software such as FoldX, while fast and reliable, lack the necessary automation features to apply them in a high-throughput manner. We introduce MutateX, a software to automate the prediction of ΔΔGs associated with the systematic mutation of each residue within a protein, or protein complex to all other possible residue types, using the FoldX energy function. MutateX also supports ΔΔG calculations over protein ensembles, upon post-translational modifications and in multimeric assemblies. At the heart of MutateX lies an automated pipeline engine that handles input preparation, parallelization and outputs publication-ready figures. We illustrate the MutateX protocol applied to different case studies. The results of the high-throughput scan provided by our tools can help in different applications, such as the analysis of disease-associated mutations, to complement experimental deep mutational scans, or assist the design of variants for industrial applications. MutateX is a collection of Python tools that relies on open-source libraries. It is available free of charge under the GNU General Public License from https://github.com/ELELAB/mutatex.

Predictive uncertainty in deep learning-based MR image reconstruction using deep ensembles: Evaluation on the fastMRI data set.

PURPOSE: To estimate pixel-wise predictive uncertainty for deep learning-based MR image reconstruction and to examine the impact of domain shifts and architecture robustness. METHODS: Uncertainty prediction could provide a measure for robustness of deep learning (DL)-based MR image reconstruction from undersampled data. DL methods bear the risk of inducing reconstruction errors like in-painting of unrealistic structures or missing pathologies. These errors may be obscured by visual realism of DL reconstruction and thus remain undiscovered. Furthermore, most methods are task-agnostic and not well calibrated to domain shifts. We propose a strategy that estimates aleatoric (data) and epistemic (model) uncertainty, which entails training a deep ensemble (epistemic) with nonnegative log-likelihood (aleatoric) loss in addition to the conventional applied losses terms. The proposed procedure can be paired with any DL reconstruction, enabling investigations of their predictive uncertainties on a pixel level. Five different architectures were investigated on the fastMRI database. The impact on the examined uncertainty of in-distributional and out-of-distributional data with changes to undersampling pattern, imaging contrast, imaging orientation, anatomy, and pathology were explored. RESULTS: Predictive uncertainty could be captured and showed good correlation to normalized mean squared error. Uncertainty was primarily focused along the aliased anatomies and on hyperintense and hypointense regions. The proposed uncertainty measure was able to detect disease prevalence shifts. Distinct predictive uncertainty patterns were observed for changing network architectures. CONCLUSION: The proposed approach enables aleatoric and epistemic uncertainty prediction for DL-based MR reconstruction with an interpretable examination on a pixel level.

In vivo ephaptic coupling allows memory network formation.

It is increasingly clear that memories are distributed across multiple brain areas. Such "engram complexes" are important features of memory formation and consolidation. Here, we test the hypothesis that engram complexes are formed in part by bioelectric fields that sculpt and guide the neural activity and tie together the areas that participate in engram complexes. Like the conductor of an orchestra, the fields influence each musician or neuron and orchestrate the output, the symphony. Our results use the theory of synergetics, machine learning, and data from a spatial delayed saccade task and provide evidence for in vivo ephaptic coupling in memory representations.

Kinetic description of site ensembles on catalytic surfaces.

We demonstrate that the Langmuir-Hinshelwood formalism is an incomplete kinetic description and, in particular, that the Hinshelwood assumption (i.e., that adsorbates are randomly distributed on the surface) is inappropriate even in catalytic reactions as simple as A + A → A2 The Hinshelwood assumption results in miscounting of site pairs (e.g., A*-A*) and, consequently, in erroneous rates, reaction orders, and identification of rate-determining steps. The clustering and isolation of surface species unnoticed by the Langmuir-Hinshelwood model is rigorously accounted for by derivation of higher-order rate terms containing statistical factors specific to each site ensemble. Ensemble-specific statistical rate terms arise irrespective of and couple with lateral adsorbate interactions, are distinct for each elementary step including surface diffusion events (e.g., A* + * → * + A*), and provide physical insight obscured by the nonanalytical nature of the kinetic Monte Carlo (kMC) method-with which the higher-order formalism quantitatively agrees. The limitations of the Langmuir-Hinshelwood model are attributed to the incorrect assertion that the rate of an elementary step is the same with respect to each site ensemble. In actuality, each elementary step-including adsorbate diffusion-traverses through each ensemble with unique rate, reversibility, and kinetic-relevance to the overall reaction rate. Explicit kinetic description of ensemble-specific paths is key to the improvements of the higher-order formalism; enables quantification of ensemble-specific rate, reversibility, and degree of rate control of surface diffusion; and reveals that a single elementary step can, counter intuitively, be both equilibrated and rate determining.

Machine learning ensembles, neural network, hybrid and sparse regression approaches for weather based rainfed cotton yield forecast.

Cotton is a major economic crop predominantly cultivated under rainfed situations. The accurate prediction of cotton yield invariably helps farmers, industries, and policy makers. The final cotton yield is mostly determined by the weather patterns that prevail during the crop growing phase. Crop yield prediction with greater accuracy is possible due to the development of innovative technologies which analyses the bigdata with its high-performance computing abilities. Machine learning technologies can make yield prediction reasonable and faster and with greater flexibility than process based complex crop simulation models. The present study demonstrates the usability of ML algorithms for yield forecasting and facilitates the comparison of different models. The cotton yield was simulated by employing the weekly weather indices as inputs and the model performance was assessed by nRMSE, MAPE and EF values. Results show that stacked generalised ensemble model and artificial neural networks predicted the cotton yield with lower nRMSE, MAPE and higher efficiency compared to other models. Variable importance studies in LASSO and ENET model found minimum temperature and relative humidity as the main determinates of cotton yield in all districts. The models were ranked based these performance metrics in the order of Stacked generalised ensemble > ANN > PCA ANN > SMLR ANN > LASSO> ENET > SVM > PCA SMLR > SMLR SVM > SMLR. This study shows that stacked generalised ensembling and ANN method can be used for reliable yield forecasting at district or county level and helps stakeholders in timely decision-making.

Explainable AI: Machine Learning Interpretation in Blackcurrant Powders.

Recently, explainability in machine and deep learning has become an important area in the field of research as well as interest, both due to the increasing use of artificial intelligence (AI) methods and understanding of the decisions made by models. The explainability of artificial intelligence (XAI) is due to the increasing consciousness in, among other things, data mining, error elimination, and learning performance by various AI algorithms. Moreover, XAI will allow the decisions made by models in problems to be more transparent as well as effective. In this study, models from the 'glass box' group of Decision Tree, among others, and the 'black box' group of Random Forest, among others, were proposed to understand the identification of selected types of currant powders. The learning process of these models was carried out to determine accuracy indicators such as accuracy, precision, recall, and F1-score. It was visualized using Local Interpretable Model Agnostic Explanations (LIMEs) to predict the effectiveness of identifying specific types of blackcurrant powders based on texture descriptors such as entropy, contrast, correlation, dissimilarity, and homogeneity. Bagging (Bagging_100), Decision Tree (DT0), and Random Forest (RF7_gini) proved to be the most effective models in the framework of currant powder interpretability. The measures of classifier performance in terms of accuracy, precision, recall, and F1-score for Bagging_100, respectively, reached values of approximately 0.979. In comparison, DT0 reached values of 0.968, 0.972, 0.968, and 0.969, and RF7_gini reached values of 0.963, 0.964, 0.963, and 0.963. These models achieved classifier performance measures of greater than 96%. In the future, XAI using agnostic models can be an additional important tool to help analyze data, including food products, even online.

NodeFlow: Towards End-to-End Flexible Probabilistic Regression on Tabular Data.

We introduce NodeFlow, a flexible framework for probabilistic regression on tabular data that combines Neural Oblivious Decision Ensembles (NODEs) and Conditional Continuous Normalizing Flows (CNFs). It offers improved modeling capabilities for arbitrary probabilistic distributions, addressing the limitations of traditional parametric approaches. In NodeFlow, the NODE captures complex relationships in tabular data through a tree-like structure, while the conditional CNF utilizes the NODE's output space as a conditioning factor. The training process of NodeFlow employs standard gradient-based learning, facilitating the end-to-end optimization of the NODEs and CNF-based density estimation. This approach ensures outstanding performance, ease of implementation, and scalability, making NodeFlow an appealing choice for practitioners and researchers. Comprehensive assessments on benchmark datasets underscore NodeFlow's efficacy, revealing its achievement of state-of-the-art outcomes in multivariate probabilistic regression setup and its strong performance in univariate regression tasks. Furthermore, ablation studies are conducted to justify the design choices of NodeFlow. In conclusion, NodeFlow's end-to-end training process and strong performance make it a compelling solution for practitioners and researchers. Additionally, it opens new avenues for research and application in the field of probabilistic regression on tabular data.