Components of immune systems face significant selective pressure to efficiently use organismal resources, mitigate infection, and resist parasitic manipulation. A theoretically optimal immune defense balances investment in constitutive and inducible immune components depending on the kinds of parasites encountered, but genetic and dynamic constraints can force deviation away from theoretical optima. One such potential constraint is pleiotropy, the phenomenon where a single gene affects multiple phenotypes. Although pleiotropy can prevent or dramatically slow adaptive evolution, it is prevalent in the signaling networks that compose metazoan immune systems. We hypothesized that pleiotropy is maintained in immune signaling networks despite slowed adaptive evolution because it provides some other advantage, such as forcing network evolution to compensate in ways that increase host fitness during infection. To study the effects of pleiotropy on the evolution of immune signaling networks, we used an agent-based modeling approach to evolve a population of host immune systems infected by simultaneously co-evolving parasites. Four kinds of pleiotropic restrictions on evolvability were incorporated into the networks, and their evolutionary outcomes were compared to, and competed against, non-pleiotropic networks. As the networks evolved, we tracked several metrics of immune network complexity, relative investment in inducible and constitutive defenses, and features associated with the winners and losers of competitive simulations. Our results suggest non-pleiotropic networks evolve to deploy highly constitutive immune responses regardless of parasite prevalence, but some implementations of pleiotropy favor the evolution of highly inducible immunity. These inducible pleiotropic networks are no less fit than non-pleiotropic networks and can out-compete non-pleiotropic networks in competitive simulations. These provide a theoretical explanation for the prevalence of pleiotropic genes in immune systems and highlight a mechanism that could facilitate the evolution of inducible immune responses.

Healthy skeletal muscle undergoes repair in response to mechanically localised strains during activities such as exercise. The ability of cells to transduce the external stimuli into a cascade of cell signalling responses is important to the process of muscle repair and regeneration. In chronic myopathies such as Duchenne muscular dystrophy and inflammatory myopathies, muscle is often subject to chronic necrosis and inflammation that perturbs tissue homeostasis and leads to non-localised, widespread damage across the tissue. Here we present an agent-based model that simulates muscle repair in response to both localised eccentric contractions similar to what would be experienced during exercise, and non-localised widespread inflammatory damage that is present in chronic disease. Computational modelling of muscle repair allows for in silico exploration of phenomena related to muscle disease. In our model, widespread inflammation led to delayed clearance of tissue damage, and delayed repair for recovery of initial fibril counts at all damage levels. Macrophage recruitment was delayed and significantly higher in widespread compared to localised damage. At higher damage percentages of 10%, widespread damage led to impaired muscle regeneration and changes in muscle geometry that represented alterations commonly observed in chronic myopathies, such as fibrosis. This computational work offers insight into the progression and aetiology of inflammatory muscle diseases, and suggests a focus on the muscle regeneration cascade in understanding the progression of muscle damage in inflammatory myopathies.

[This corrects the article DOI: 10.1371/journal.pcbi.1010651.].

Bacterial microcompartments (BMC) are complex macromolecular assemblies that participate in varied chemical processes in about one fourth of bacterial species. BMC-encapsulated enzymatic activities are segregated from other cell contents by means of semipermeable shells, justifying why BMC are viewed as prototype nano-reactors for biotechnological applications. Herein, we undertook a comparative study of bending propensities of BMC hexamers (BMC-H), the most abundant shell constituents. Published data show that some BMC-H, like β-carboxysomal CcmK, tend to assemble flat whereas other BMC-H often build curved objects. Inspection of available crystal structures presenting BMC-H in tiled arrangements permitted us to identify two major assembly modes with a striking connection with experimental trends. All-atom molecular dynamics (MD) supported that BMC-H bending is triggered robustly only from the arrangement adopted in crystals by BMC-H that experimentally form curved objects, leading to very similar arrangements to those found in structures of recomposed BMC shells. Simulations on triplets of planar-behaving hexamers, which were previously reconfigured to comply with such organization, confirmed that bending propensity is mostly defined by the precise lateral positioning of hexamers, rather than by BMC-H identity. Finally, an interfacial lysine was pinpointed as the most decisive residue in controlling PduA spontaneous curvature. Globally, results presented herein should contribute to improve our understanding of the variable mechanisms of biogenesis characterized for BMC, and of possible strategies to regulate BMC size and shape.

When a neuron is driven beyond its threshold, it spikes. The fact that it does not communicate its continuous membrane potential is usually seen as a computational liability. Here we show that this spiking mechanism allows neurons to produce an unbiased estimate of their causal influence, and a way of approximating gradient descent-based learning. Importantly, neither activity of upstream neurons, which act as confounders, nor downstream non-linearities bias the results. We show how spiking enables neurons to solve causal estimation problems and that local plasticity can approximate gradient descent using spike discontinuity learning.

The increase in microbial sequenced genomes from pure cultures and metagenomic samples reflects the current attainability of whole-genome and shotgun sequencing methods. However, software for genome visualization still lacks automation, integration of different analyses, and customizable options for non-experienced users. In this study, we introduce GenoVi, a Python command-line tool able to create custom circular genome representations for the analysis and visualization of microbial genomes and sequence elements. It is designed to work with complete or draft genomes, featuring customizable options including 25 different built-in color palettes (including 5 color-blind safe palettes), text formatting options, and automatic scaling for complete genomes or sequence elements with more than one replicon/sequence. Using a Genbank format file as the input file or multiple files within a directory, GenoVi (i) visualizes genomic features from the GenBank annotation file, (ii) integrates a Cluster of Orthologs Group (COG) categories analysis using DeepNOG, (iii) automatically scales the visualization of each replicon of complete genomes or multiple sequence elements, (iv) and generates COG histograms, COG frequency heatmaps and output tables including general stats of each replicon or contig processed. GenoVi's potential was assessed by analyzing single and multiple genomes of Bacteria and Archaea. Paraburkholderia genomes were analyzed to obtain a fast classification of replicons in large multipartite genomes. GenoVi works as an easy-to-use command-line tool and provides customizable options to automatically generate genomic maps for scientific publications, educational resources, and outreach activities. GenoVi is freely available and can be downloaded from https://github.com/robotoD/GenoVi.

Established prognostic tests based on limited numbers of transcripts can identify high-risk breast cancer patients, yet are approved only for individuals presenting with specific clinical features or disease characteristics. Deep learning algorithms could hold potential for stratifying patient cohorts based on full transcriptome data, yet the development of robust classifiers is hampered by the number of variables in omics datasets typically far exceeding the number of patients. To overcome this hurdle, we propose a classifier based on a data augmentation pipeline consisting of a Wasserstein generative adversarial network (GAN) with gradient penalty and an embedded auxiliary classifier to obtain a trained GAN discriminator (T-GAN-D). Applied to 1244 patients of the METABRIC breast cancer cohort, this classifier outperformed established breast cancer biomarkers in separating low- from high-risk patients (disease specific death, progression or relapse within 10 years from initial diagnosis). Importantly, the T-GAN-D also performed across independent, merged transcriptome datasets (METABRIC and TCGA-BRCA cohorts), and merging data improved overall patient stratification. In conclusion, the reiterative GAN-based training process allowed generating a robust classifier capable of stratifying low- vs high-risk patients based on full transcriptome data and across independent and heterogeneous breast cancer cohorts.

Despite the considerable progress of in vivo neural recording techniques, inferring the biophysical mechanisms underlying large scale coordination of brain activity from neural data remains challenging. One obstacle is the difficulty to link high dimensional functional connectivity measures to mechanistic models of network activity. We address this issue by investigating spike-field coupling (SFC) measurements, which quantify the synchronization between, on the one hand, the action potentials produced by neurons, and on the other hand mesoscopic "field" signals, reflecting subthreshold activities at possibly multiple recording sites. As the number of recording sites gets large, the amount of pairwise SFC measurements becomes overwhelmingly challenging to interpret. We develop Generalized Phase Locking Analysis (GPLA) as an interpretable dimensionality reduction of this multivariate SFC. GPLA describes the dominant coupling between field activity and neural ensembles across space and frequencies. We show that GPLA features are biophysically interpretable when used in conjunction with appropriate network models, such that we can identify the influence of underlying circuit properties on these features. We demonstrate the statistical benefits and interpretability of this approach in various computational models and Utah array recordings. The results suggest that GPLA, used jointly with biophysical modeling, can help uncover the contribution of recurrent microcircuits to the spatio-temporal dynamics observed in multi-channel experimental recordings.

Human vision is still largely unexplained. Computer vision made impressive progress on this front, but it is still unclear to which extent artificial neural networks approximate human object vision at the behavioral and neural levels. Here, we investigated whether machine object vision mimics the representational hierarchy of human object vision with an experimental design that allows testing within-domain representations for animals and scenes, as well as across-domain representations reflecting their real-world contextual regularities such as animal-scene pairs that often co-occur in the visual environment. We found that DCNNs trained in object recognition acquire representations, in their late processing stage, that closely capture human conceptual judgements about the co-occurrence of animals and their typical scenes. Likewise, the DCNNs representational hierarchy shows surprising similarities with the representational transformations emerging in domain-specific ventrotemporal areas up to domain-general frontoparietal areas. Despite these remarkable similarities, the underlying information processing differs. The ability of neural networks to learn a human-like high-level conceptual representation of object-scene co-occurrence depends upon the amount of object-scene co-occurrence present in the image set thus highlighting the fundamental role of training history. Further, although mid/high-level DCNN layers represent the category division for animals and scenes as observed in VTC, its information content shows reduced domain-specific representational richness. To conclude, by testing within- and between-domain selectivity while manipulating contextual regularities we reveal unknown similarities and differences in the information processing strategies employed by human and artificial visual systems.

As infected and vaccinated population increases, some countries decided not to impose non-pharmaceutical intervention measures anymore and to coexist with COVID-19. However, we do not have a comprehensive understanding of its consequence, especially for China where most population has not been infected and most Omicron transmissions are silent. This paper aims to reveal the complete silent transmission dynamics of COVID-19 by agent-based simulations overlaying a big data of more than 0.7 million real individual mobility tracks without any intervention measures throughout a week in a Chinese city, with an extent of completeness and realism not attained in existing studies. Together with the empirically inferred transmission rate of COVID-19, we find surprisingly that with only 70 citizens to be infected initially, 0.33 million becomes infected silently at last. We also reveal a characteristic daily periodic pattern of the transmission dynamics, with peaks in mornings and afternoons. In addition, by inferring individual professions, visited locations and age group, we found that retailing, catering and hotel staff are more likely to get infected than other professions, and elderly and retirees are more likely to get infected at home than outside home.

Transcranial alternating current stimulation (tACS) represents a promising non-invasive treatment for an increasingly wide range of neurological and neuropsychiatric disorders. The ability to use periodically oscillating electric fields to non-invasively engage neural dynamics opens up the possibility of recruiting synaptic plasticity and to modulate brain function. However, despite consistent reports about tACS clinical effectiveness, strong state-dependence combined with the ubiquitous heterogeneity of cortical networks collectively results in high outcome variability. Introducing variations in intrinsic neuronal timescales, we explored how such heterogeneity influences stimulation-induced change in synaptic connectivity. We examined how spike timing dependent plasticity, at the level of cells, intra- and inter-laminar cortical networks, can be selectively and preferentially engaged by periodic stimulation. Using leaky integrate-and-fire neuron models, we analyzed cortical circuits comprised of multiple cell-types, alongside superficial multi-layered networks expressing distinct layer-specific timescales. Our results show that mismatch in neuronal timescales within and/or between cells-and the resulting variability in excitability, temporal integration properties and frequency tuning-enables selective and directional control on synaptic connectivity by tACS. Our work provides new vistas on how to recruit neural heterogeneity to guide brain plasticity using non-invasive stimulation paradigms.

The mosquito Aedes aegypti is the vector of a number of medically-important viruses, including dengue virus, yellow fever virus, chikungunya virus, and Zika virus, and as such vector control is a key approach to managing the diseases they cause. Understanding the impact of vector control on these diseases is aided by first understanding its impact on Ae. aegypti population dynamics. A number of detail-rich models have been developed to couple the dynamics of the immature and adult stages of Ae. aegypti. The numerous assumptions of these models enable them to realistically characterize impacts of mosquito control, but they also constrain the ability of such models to reproduce empirical patterns that do not conform to the models' behavior. In contrast, statistical models afford sufficient flexibility to extract nuanced signals from noisy data, yet they have limited ability to make predictions about impacts of mosquito control on disease caused by pathogens that the mosquitoes transmit without extensive data on mosquitoes and disease. Here, we demonstrate how the differing strengths of mechanistic realism and statistical flexibility can be fused into a single model. Our analysis utilizes data from 176,352 household-level Ae. aegypti aspirator collections conducted during 1999-2011 in Iquitos, Peru. The key step in our approach is to calibrate a single parameter of the model to spatio-temporal abundance patterns predicted by a generalized additive model (GAM). In effect, this calibrated parameter absorbs residual variation in the abundance time-series not captured by other features of the mechanistic model. We then used this calibrated parameter and the literature-derived parameters in the agent-based model to explore Ae. aegypti population dynamics and the impact of insecticide spraying to kill adult mosquitoes. The baseline abundance predicted by the agent-based model closely matched that predicted by the GAM. Following spraying, the agent-based model predicted that mosquito abundance rebounds within about two months, commensurate with recent experimental data from Iquitos. Our approach was able to accurately reproduce abundance patterns in Iquitos and produce a realistic response to adulticide spraying, while retaining sufficient flexibility to be applied across a range of settings.

[This corrects the article DOI: 10.1371/journal.pcbi.1010919.].

Motor symptoms of Parkinson's Disease (PD) are associated with dopamine deficits and pathological oscillation of basal ganglia (BG) neurons in the β range ([12-30] Hz). However, how dopamine depletion affects the oscillation dynamics of BG nuclei is still unclear. With a spiking neurons model, we here capture the features of BG nuclei interactions leading to oscillations in dopamine-depleted condition. We highlight that both the loop between subthalamic nucleus (STN) and Globus Pallidus pars externa (GPe) and the loop between striatal fast spiking and medium spiny neurons and GPe display resonances in the β range, and synchronize to a common β frequency through interaction. Crucially, the synchronization depends on dopamine depletion: the two loops are largely independent for high levels of dopamine, but progressively synchronize as dopamine is depleted due to the increased strength of the striatal loop. The model is validated against recent experimental reports on the role of cortical inputs, STN and GPe activity in the generation of β oscillations. Our results highlight the role of the interplay between the GPe-STN and the GPe-striatum loop in generating sustained β oscillations in PD subjects, and explain how this interplay depends on the level of dopamine. This paves the way to the design of therapies specifically addressing the onset of pathological β oscillations.

Controlled ovarian stimulation is tailored to the patient based on clinical parameters but estimating the number of retrieved metaphase II (MII) oocytes is a challenge. Here, we have developed a model that takes advantage of the patient's genetic and clinical characteristics simultaneously for predicting the stimulation outcome. Sequence variants in reproduction-related genes identified by next-generation sequencing were matched to groups of various MII oocyte counts using ranking, correspondence analysis, and self-organizing map methods. The gradient boosting machine technique was used to train models on a clinical dataset of 8,574 or a clinical-genetic dataset of 516 ovarian stimulations. The clinical-genetic model predicted the number of MII oocytes better than that based on clinical data. Anti-Müllerian hormone level and antral follicle count were the two most important predictors while a genetic feature consisting of sequence variants in the GDF9, LHCGR, FSHB, ESR1, and ESR2 genes was the third. The combined contribution of genetic features important for the prediction was over one-third of that revealed for anti-Müllerian hormone. Predictions of our clinical-genetic model accurately matched individuals' actual outcomes preventing over- or underestimation. The genetic data upgrades the personalized prediction of ovarian stimulation outcomes, thus improving the in vitro fertilization procedure.

Bayesian inference for phylogenetics is a gold standard for computing distributions of phylogenies. However, Bayesian phylogenetics faces the challenging computational problem of moving throughout the high-dimensional space of trees. Fortunately, hyperbolic space offers a low dimensional representation of tree-like data. In this paper, we embed genomic sequences as points in hyperbolic space and perform hyperbolic Markov Chain Monte Carlo for Bayesian inference in this space. The posterior probability of an embedding is computed by decoding a neighbour-joining tree from the embedding locations of the sequences. We empirically demonstrate the fidelity of this method on eight data sets. We systematically investigated the effect of embedding dimension and hyperbolic curvature on the performance in these data sets. The sampled posterior distribution recovers the splits and branch lengths to a high degree over a range of curvatures and dimensions. We systematically investigated the effects of the embedding space's curvature and dimension on the Markov Chain's performance, demonstrating the suitability of hyperbolic space for phylogenetic inference.

Mathematical models are often used to explore network-driven cellular processes from a systems perspective. However, a dearth of quantitative data suitable for model calibration leads to models with parameter unidentifiability and questionable predictive power. Here we introduce a combined Bayesian and Machine Learning Measurement Model approach to explore how quantitative and non-quantitative data constrain models of apoptosis execution within a missing data context. We find model prediction accuracy and certainty strongly depend on rigorous data-driven formulations of the measurement, and the size and make-up of the datasets. For instance, two orders of magnitude more ordinal (e.g., immunoblot) data are necessary to achieve accuracy comparable to quantitative (e.g., fluorescence) data for calibration of an apoptosis execution model. Notably, ordinal and nominal (e.g., cell fate observations) non-quantitative data synergize to reduce model uncertainty and improve accuracy. Finally, we demonstrate the potential of a data-driven Measurement Model approach to identify model features that could lead to informative experimental measurements and improve model predictive power.

Clostridioides difficile pathogenesis is mediated through its two toxin proteins, TcdA and TcdB, which induce intestinal epithelial cell death and inflammation. It is possible to alter C. difficile toxin production by changing various metabolite concentrations within the extracellular environment. However, it is unknown which intracellular metabolic pathways are involved and how they regulate toxin production. To investigate the response of intracellular metabolic pathways to diverse nutritional environments and toxin production states, we use previously published genome-scale metabolic models of C. difficile strains CD630 and CDR20291 (iCdG709 and iCdR703). We integrated publicly available transcriptomic data with the models using the RIPTiDe algorithm to create 16 unique contextualized C. difficile models representing a range of nutritional environments and toxin states. We used Random Forest with flux sampling and shadow pricing analyses to identify metabolic patterns correlated with toxin states and environment. Specifically, we found that arginine and ornithine uptake is particularly active in low toxin states. Additionally, uptake of arginine and ornithine is highly dependent on intracellular fatty acid and large polymer metabolite pools. We also applied the metabolic transformation algorithm (MTA) to identify model perturbations that shift metabolism from a high toxin state to a low toxin state. This analysis expands our understanding of toxin production in C. difficile and identifies metabolic dependencies that could be leveraged to mitigate disease severity.

Rhodotorula toruloides is a non-conventional, oleaginous yeast able to naturally accumulate high amounts of microbial lipids. Constraint-based modeling of R. toruloides has been mainly focused on the comparison of experimentally measured and model predicted growth rates, while the intracellular flux patterns have been analyzed on a rather general level. Hence, the intrinsic metabolic properties of R. toruloides that make lipid synthesis possible are not thoroughly understood. At the same time, the lack of diverse physiological data sets has often been the bottleneck to predict accurate fluxes. In this study, we collected detailed physiology data sets of R. toruloides while growing on glucose, xylose, and acetate as the sole carbon source in chemically defined medium. Regardless of the carbon source, the growth was divided into two phases from which proteomic and lipidomic data were collected. Complemental physiological parameters were collected in these two phases and altogether implemented into metabolic models. Simulated intracellular flux patterns demonstrated the role of phosphoketolase in the generation of acetyl-CoA, one of the main precursors during lipid biosynthesis, while the role of ATP citrate lyase was not confirmed. Metabolic modeling on xylose as a carbon substrate was greatly improved by the detection of chirality of D-arabinitol, which together with D-ribulose were involved in an alternative xylose assimilation pathway. Further, flux patterns pointed to metabolic trade-offs associated with NADPH allocation between nitrogen assimilation and lipid biosynthetic pathways, which was linked to large-scale differences in protein and lipid content. This work includes the first extensive multi-condition analysis of R. toruloides using enzyme-constrained models and quantitative proteomics. Further, more precise kcat values should extend the application of the newly developed enzyme-constrained models that are publicly available for future studies.

Biochemical processes in cells, including enzyme-catalyzed reactions, occur in crowded conditions with various background macromolecules occupying up to 40% of cytoplasm's volume. Viral enzymes in the host cell also encounter such crowded conditions as they often function at the endoplasmic reticulum membranes. We focus on an enzyme encoded by the hepatitis C virus, the NS3/4A protease, which is crucial for viral replication. We have previously found experimentally that synthetic crowders, polyethylene glycol (PEG) and branched polysucrose (Ficoll), differently affect the kinetic parameters of peptide hydrolysis catalyzed by NS3/4A. To gain understanding of the reasons for such behavior, we perform atomistic molecular dynamics simulations of NS3/4A in the presence of either PEG or Ficoll crowders and with and without the peptide substrates. We find that both crowder types make nanosecond long contacts with the protease and slow down its diffusion. However, they also affect the enzyme structural dynamics; crowders induce functionally relevant helical structures in the disordered parts of the protease cofactor, NS4A, with the PEG effect being more pronounced. Overall, PEG interactions with NS3/4A are slightly stronger but Ficoll forms more hydrogen bonds with NS3. The crowders also interact with substrates; we find that the substrate diffusion is reduced much more in the presence of PEG than Ficoll. However, contrary to NS3, the substrate interacts more strongly with Ficoll than with PEG crowders, with the substrate diffusion being similar to crowder diffusion. Importantly, crowders also affect the substrate-enzyme interactions. We observe that both PEG and Ficoll enhance the presence of substrates near the active site, especially near catalytic H57 but Ficoll crowders increase substrate binding more than PEG molecules.

Neural system identification aims at learning the response function of neurons to arbitrary stimuli using experimentally recorded data, but typically does not leverage normative principles such as efficient coding of natural environments. Visual systems, however, have evolved to efficiently process input from the natural environment. Here, we present a normative network regularization for system identification models by incorporating, as a regularizer, the efficient coding hypothesis, which states that neural response properties of sensory representations are strongly shaped by the need to preserve most of the stimulus information with limited resources. Using this approach, we explored if a system identification model can be improved by sharing its convolutional filters with those of an autoencoder which aims to efficiently encode natural stimuli. To this end, we built a hybrid model to predict the responses of retinal neurons to noise stimuli. This approach did not only yield a higher performance than the "stand-alone" system identification model, it also produced more biologically plausible filters, meaning that they more closely resembled neural representation in early visual systems. We found these results applied to retinal responses to different artificial stimuli and across model architectures. Moreover, our normatively regularized model performed particularly well in predicting responses of direction-of-motion sensitive retinal neurons. The benefit of natural scene statistics became marginal, however, for predicting the responses to natural movies. In summary, our results indicate that efficiently encoding environmental inputs can improve system identification models, at least for noise stimuli, and point to the benefit of probing the visual system with naturalistic stimuli.

With the evolution of multicellularity, communication among cells in different tissues and organs became pivotal to life. Molecular basis of such communication has long been studied, but genome-wide screens for genes and other biomolecules mediating tissue-tissue signaling are lacking. To systematically identify inter-tissue mediators, we present a novel computational approach MultiCens (Multilayer/Multi-tissue network Centrality measures). Unlike single-layer network methods, MultiCens can distinguish within- vs. across-layer connectivity to quantify the "influence" of any gene in a tissue on a query set of genes of interest in another tissue. MultiCens enjoys theoretical guarantees on convergence and decomposability, and performs well on synthetic benchmarks. On human multi-tissue datasets, MultiCens predicts known and novel genes linked to hormones. MultiCens further reveals shifts in gene network architecture among four brain regions in Alzheimer's disease. MultiCens-prioritized hypotheses from these two diverse applications, and potential future ones like "Multi-tissue-expanded Gene Ontology" analysis, can enable whole-body yet molecular-level systems investigations in humans.

Computational fluid dynamics (CFD) can be used to simulate vascular haemodynamics and analyse potential treatment options. CFD has shown to be beneficial in improving patient outcomes. However, the implementation of CFD for routine clinical use is yet to be realised. Barriers for CFD include high computational resources, specialist experience needed for designing simulation set-ups, and long processing times. The aim of this study was to explore the use of machine learning (ML) to replicate conventional aortic CFD with automatic and fast regression models. Data used to train/test the model consisted of 3,000 CFD simulations performed on synthetically generated 3D aortic shapes. These subjects were generated from a statistical shape model (SSM) built on real patient-specific aortas (N = 67). Inference performed on 200 test shapes resulted in average errors of 6.01% ±3.12 SD and 3.99% ±0.93 SD for pressure and velocity, respectively. Our ML-based models performed CFD in ∼0.075 seconds (4,000x faster than the solver). This proof-of-concept study shows that results from conventional vascular CFD can be reproduced using ML at a much faster rate, in an automatic process, and with reasonable accuracy.

Flies detect and track moving targets among visual clutter, and this process mainly relies on visual motion. Visual motion is analyzed or computed with the pathway from the retina to T4/T5 cells. The computation of local directional motion was formulated as an elementary movement detector (EMD) model more than half a century ago. Solving target detection or figure-ground discrimination problems can be equivalent to extracting boundaries between a target and the background based on the motion discontinuities in the output of a retinotopic array of EMDs. Individual EMDs cannot measure true velocities, however, due to their sensitivity to pattern properties such as luminance contrast and spatial frequency content. It remains unclear how local directional motion signals are further integrated to enable figure-ground discrimination. Here, we present a computational model inspired by fly motion vision. Simulations suggest that the heavily fluctuating output of an EMD array is naturally surmounted by a lobula network, which is hypothesized to be downstream of the local motion detectors and have parallel pathways with distinct directional selectivity. The lobula network carries out a spatiotemporal smoothing operation for visual motion, especially across time, enabling the segmentation of moving figures from the background. The model qualitatively reproduces experimental observations in the visually evoked response characteristics of one type of lobula columnar (LC) cell. The model is further shown to be robust to natural scene variability. Our results suggest that the lobula is involved in local motion-based target detection.

Agent-based models (ABMs) have enabled great advances in the study of tumor development and therapeutic response, allowing researchers to explore the spatiotemporal evolution of the tumor and its microenvironment. However, these models face serious drawbacks in the realm of parameterization - ABM parameters are typically set individually based on various data and literature sources, rather than through a rigorous parameter estimation approach. While ABMs can be fit to simple time-course data (such as tumor volume), that type of data loses the spatial information that is a defining feature of ABMs. While tumor images provide spatial information, it is exceedingly difficult to compare tumor images to ABM simulations beyond a qualitative visual comparison. Without a quantitative method of comparing the similarity of tumor images to ABM simulations, a rigorous parameter fitting is not possible. Here, we present a novel approach that applies neural networks to represent both tumor images and ABM simulations as low dimensional points, with the distance between points acting as a quantitative measure of difference between the two. This enables a quantitative comparison of tumor images and ABM simulations, where the distance between simulated and experimental images can be minimized using standard parameter-fitting algorithms. Here, we describe this method and present two examples to demonstrate the application of the approach to estimate parameters for two distinct ABMs. Overall, we provide a novel method to robustly estimate ABM parameters.

Predator-prey theory is commonly used to describe tumor growth in the presence of selective pressure from the adaptive immune system. These interactions are mediated by the tumor immunopeptidome (what the tumor "shows" the body) and the T-cell receptor (TCR) repertoire (how well the body "sees" cancer cells). The tumor immunopeptidome comprises neoantigens which can be gained and lost throughout tumorigenesis and treatment. Heterogeneity in the immunopeptidome is predictive of poor response to immunotherapy in some tumor types, suggesting that the TCR repertoire is unable to support a fully polyclonal response against every neoantigen. Importantly, while tumor and T-cell populations are known to compete with each other for intratumoral resources, whether between-lineage competition among peripheral T cells influences the TCR repertoire is unknown and difficult to interrogate experimentally. Computational models may offer a way to investigate these phenomena and deepen our understanding of the tumor-immune axis. Here, we construct a predator-prey-like model and calibrate it to preclinical and clinical data to describe tumor growth and immunopeptidome diversification. Simultaneously, we model the expansion of antigen-specific T-cell lineages and their consumption of both lineage-specific antigenic resources and lineage-agnostic, shared resources. This predator-prey-like framework accurately described clinically observed immunopeptidomes; recapitulated response-associated effects of immunotherapy, including immunoediting; and allowed exploration of treatment of tumors with varying growth and mutation rates.

Mitochondria form a network in the cell that rapidly changes through fission, fusion, and motility. Dysregulation of this four-dimensional (4D: x,y,z,time) network is implicated in numerous diseases ranging from cancer to neurodegeneration. While lattice light-sheet microscopy has recently made it possible to image mitochondria in 4D, quantitative analysis methods for the resulting datasets have been lacking. Here we present MitoTNT, the first-in-class software for Mitochondrial Temporal Network Tracking in 4D live-cell fluorescence microscopy data. MitoTNT uses spatial proximity and network topology to compute an optimal tracking assignment. To validate the accuracy of tracking, we created a reaction-diffusion simulation to model mitochondrial network motion and remodeling events. We found that our tracking is >90% accurate for ground-truth simulations and agrees well with published motility results for experimental data. We used MitoTNT to quantify 4D mitochondrial networks from human induced pluripotent stem cells. First, we characterized sub-fragment motility and analyzed network branch motion patterns. We revealed that the skeleton node motion is correlated along branch nodes and is uncorrelated in time. Second, we identified fission and fusion events with high spatiotemporal resolution. We found that mitochondrial skeleton nodes near the fission/fusion sites move nearly twice as fast as random skeleton nodes and that microtubules play a role in mediating selective fission/fusion. Finally, we developed graph-based transport simulations that model how material would distribute on experimentally measured mitochondrial temporal networks. We showed that pharmacological perturbations increase network reachability but decrease network resilience through a combination of altered mitochondrial fission/fusion dynamics and motility. MitoTNT's easy-to-use tracking module, interactive 4D visualization capability, and powerful post-tracking analyses aim at making temporal network tracking accessible to the wider mitochondria research community.

Recent advances in molecular transduction of odorants in the Olfactory Sensory Neurons (OSNs) of the Drosophila Antenna have shown that the odorant object identity is multiplicatively coupled with the odorant concentration waveform. The resulting combinatorial neural code is a confounding representation of odorant semantic information (identity) and syntactic information (concentration). To distill the functional logic of odor information processing in the Antennal Lobe (AL) a number of challenges need to be addressed including 1) how is the odorant semantic information decoupled from the syntactic information at the level of the AL, 2) how are these two information streams processed by the diverse AL Local Neurons (LNs) and 3) what is the end-to-end functional logic of the AL? By analyzing single-channel physiology recordings at the output of the AL, we found that the Projection Neuron responses can be decomposed into a concentration-invariant component, and two transient components boosting the positive/negative concentration contrast that indicate onset/offset timing information of the odorant object. We hypothesized that the concentration-invariant component, in the multi-channel context, is the recovered odorant identity vector presented between onset/offset timing events. We developed a model of LN pathways in the Antennal Lobe termed the differential Divisive Normalization Processors (DNPs), which robustly extract the semantics (the identity of the odorant object) and the ON/OFF semantic timing events indicating the presence/absence of an odorant object. For real-time processing with spiking PN models, we showed that the phase-space of the biological spike generator of the PN offers an intuit perspective for the representation of recovered odorant semantics and examined the dynamics induced by the odorant semantic timing events. Finally, we provided theoretical and computational evidence for the functional logic of the AL as a robust ON-OFF odorant object identity recovery processor across odorant identities, concentration amplitudes and waveform profiles.

Multistep protein-protein interactions underlie most biological processes, but their characterization through methods such as isothermal titration calorimetry (ITC) is largely confined to simple models that provide little information on the intermediate, individual steps. In this study, we primarily examine the essential hub protein LC8, a small dimer that binds disordered regions of 100+ client proteins in two symmetrical grooves at the dimer interface. Mechanistic details of LC8 binding have remained elusive, hampered in part by ITC data analyses employing simple models that treat bivalent binding as a single event with a single binding affinity. We build on existing Bayesian ITC approaches to quantify thermodynamic parameters for multi-site binding interactions impacted by significant uncertainty in protein concentration. Using a two-site binding model, we identify positive cooperativity with high confidence for LC8 binding to multiple client peptides. In contrast, application of an identical model to the two-site binding between the coiled-coil NudE dimer and the intermediate chain of dynein reveals little evidence of cooperativity. We propose that cooperativity in the LC8 system drives the formation of saturated induced-dimer structures, the functional units of most LC8 complexes. In addition to these system-specific findings, our work advances general ITC analysis in two ways. First, we describe a previously unrecognized mathematical ambiguity in concentrations in standard binding models and clarify how it impacts the precision with which binding parameters are determinable in cases of high uncertainty in analyte concentrations. Second, building on observations in the LC8 system, we develop a system-agnostic heat map of practical parameter identifiability calculated from synthetic data which demonstrates that the ability to determine microscopic binding parameters is strongly dependent on both the parameters themselves and experimental conditions. The work serves as a foundation for determination of multi-step binding interactions, and we outline best practices for Bayesian analysis of ITC experiments.

Multicellularity has evolved several independent times over the past hundreds of millions of years and given rise to a wide diversity of complex life. Recent studies have found that large differences in the fundamental structure of early multicellular life cycles can affect fitness and influence multicellular adaptation. Yet, there is an underlying assumption that at some scale or categorization multicellular life cycles are similar in terms of their adaptive potential. Here, we consider this possibility by exploring adaptation in a class of simple multicellular life cycles of filamentous organisms that only differ in one respect, how many daughter filaments are produced. We use mathematical models and evolutionary simulations to show that despite the similarities, qualitatively different mutations fix. In particular, we find that mutations with a tradeoff between cell growth and group survival, i.e. "selfish" or "altruistic" traits, spread differently. Specifically, altruistic mutations more readily spread in life cycles that produce few daughters while in life cycles producing many daughters either type of mutation can spread depending on the environment. Our results show that subtle changes in multicellular life cycles can fundamentally alter adaptation.

Mechanistic models have been used for centuries to describe complex interconnected processes, including biological ones. As the scope of these models has widened, so have their computational demands. This complexity can limit their suitability when running many simulations or when real-time results are required. Surrogate machine learning (ML) models can be used to approximate the behaviour of complex mechanistic models, and once built, their computational demands are several orders of magnitude lower. This paper provides an overview of the relevant literature, both from an applicability and a theoretical perspective. For the latter, the paper focuses on the design and training of the underlying ML models. Application-wise, we show how ML surrogates have been used to approximate different mechanistic models. We present a perspective on how these approaches can be applied to models representing biological processes with potential industrial applications (e.g., metabolism and whole-cell modelling) and show why surrogate ML models may hold the key to making the simulation of complex biological systems possible using a typical desktop computer.

Expression noise, the variability of the amount of gene product among isogenic cells grown in identical conditions, originates from the inherent stochasticity of diffusion and binding of the molecular players involved in transcription and translation. It has been shown that expression noise is an evolvable trait and that central genes exhibit less noise than peripheral genes in gene networks. A possible explanation for this pattern is increased selective pressure on central genes since they propagate their noise to downstream targets, leading to noise amplification. To test this hypothesis, we developed a new gene regulatory network model with inheritable stochastic gene expression and simulated the evolution of gene-specific expression noise under constraint at the network level. Stabilizing selection was imposed on the expression level of all genes in the network and rounds of mutation, selection, replication and recombination were performed. We observed that local network features affect both the probability to respond to selection, and the strength of the selective pressure acting on individual genes. In particular, the reduction of gene-specific expression noise as a response to stabilizing selection on the gene expression level is higher in genes with higher centrality metrics. Furthermore, global topological structures such as network diameter, centralization and average degree affect the average expression variance and average selective pressure acting on constituent genes. Our results demonstrate that selection at the network level leads to differential selective pressure at the gene level, and local and global network characteristics are an essential component of gene-specific expression noise evolution.

Addressing many of the major outstanding questions in the fields of microbial evolution and pathogenesis will require analyses of populations of microbial genomes. Although population genomic studies provide the analytical resolution to investigate evolutionary and mechanistic processes at fine spatial and temporal scales-precisely the scales at which these processes occur-microbial population genomic research is currently hindered by the practicalities of obtaining sufficient quantities of the relatively pure microbial genomic DNA necessary for next-generation sequencing. Here we present swga2.0, an optimized and parallelized pipeline to design selective whole genome amplification (SWGA) primer sets. Unlike previous methods, swga2.0 incorporates active and machine learning methods to evaluate the amplification efficacy of individual primers and primer sets. Additionally, swga2.0 optimizes primer set search and evaluation strategies, including parallelization at each stage of the pipeline, to dramatically decrease program runtime. Here we describe the swga2.0 pipeline, including the empirical data used to identify primer and primer set characteristics, that improve amplification performance. Additionally, we evaluate the novel swga2.0 pipeline by designing primer sets that successfully amplify Prevotella melaninogenica, an important component of the lung microbiome in cystic fibrosis patients, from samples dominated by human DNA.

The genetic code refers to a rule that maps 64 codons to 20 amino acids. Nearly all organisms, with few exceptions, share the same genetic code, the standard genetic code (SGC). While it remains unclear why this universal code has arisen and been maintained during evolution, it may have been preserved under selection pressure. Theoretical studies comparing the SGC and numerically created hypothetical random genetic codes have suggested that the SGC has been subject to strong selection pressure for being robust against translation errors. However, these prior studies have searched for random genetic codes in only a small subspace of the possible code space due to limitations in computation time. Thus, how the genetic code has evolved, and the characteristics of the genetic code fitness landscape, remain unclear. By applying multicanonical Monte Carlo, an efficient rare-event sampling method, we efficiently sampled random codes from a much broader random ensemble of genetic codes than in previous studies, estimating that only one out of every 1020 random codes is more robust than the SGC. This estimate is significantly smaller than the previous estimate, one in a million. We also characterized the fitness landscape of the genetic code that has four major fitness peaks, one of which includes the SGC. Furthermore, genetic algorithm analysis revealed that evolution under such a multi-peaked fitness landscape could be strongly biased toward a narrow peak, in an evolutionary path-dependent manner.

Multi-view data can be generated from diverse sources, by different technologies, and in multiple modalities. In various fields, integrating information from multi-view data has pushed the frontier of discovery. In this paper, we develop a new approach for multi-view clustering, which overcomes the limitations of existing methods such as the need of pooling data across views, restrictions on the clustering algorithms allowed within each view, and the disregard for complementary information between views. Our new method, called CPS-merge analysis, merges clusters formed by the Cartesian product of single-view cluster labels, guided by the principle of maximizing clustering stability as evaluated by CPS analysis. In addition, we introduce measures to quantify the contribution of each view to the formation of any cluster. CPS-merge analysis can be easily incorporated into an existing clustering pipeline because it only requires single-view cluster labels instead of the original data. We can thus readily apply advanced single-view clustering algorithms. Importantly, our approach accounts for both consensus and complementary effects between different views, whereas existing ensemble methods focus on finding a consensus for multiple clustering results, implying that results from different views are variations of one clustering structure. Through experiments on single-cell datasets, we demonstrate that our approach frequently outperforms other state-of-the-art methods.

Neurons integrate from thousands of synapses whose strengths span an order of magnitude. Intriguingly, in mouse neocortex, the few 'strong' synapses are formed between similarly tuned cells, suggesting they determine spiking output. This raises the question of how other computational primitives, including 'background' activity from the many 'weak' synapses, short-term plasticity, and temporal factors contribute to spiking. We used paired recordings and extracellular stimulation experiments to map excitatory postsynaptic potential (EPSP) amplitudes and paired-pulse ratios of synaptic connections formed between pyramidal neurons in layer 2/3 (L2/3) of barrel cortex. While net short-term plasticity was weak, strong synaptic connections were exclusively depressing. Importantly, we found no evidence for clustering of synaptic properties on individual neurons. Instead, EPSPs and paired-pulse ratios of connections converging onto the same cells spanned the full range observed across L2/3, which critically constrains theoretical models of cortical filtering. To investigate how different computational primitives of synaptic information processing interact to shape spiking, we developed a computational model of a pyramidal neuron in the excitatory L2/3 circuitry, which was constrained by our experiments and published in vivo data. We found that strong synapses were substantially depressed during ongoing activation and their ability to evoke correlated spiking primarily depended on their high temporal synchrony and high firing rates observed in vivo. However, despite this depression, their larger EPSP amplitudes strongly amplified information transfer and responsiveness. Thus, our results contribute to a nuanced framework of how cortical neurons exploit synergies between temporal coding, synaptic properties, and noise to transform synaptic inputs into spikes.

Making no use of physical laws or co-evolutionary information, de novo deep learning (DL) models for RNA secondary structure prediction have achieved far superior performances than traditional algorithms. However, their statistical underpinning raises the crucial question of generalizability. We present a quantitative study of the performance and generalizability of a series of de novo DL models, with a minimal two-module architecture and no post-processing, under varied similarities between seen and unseen sequences. Our models demonstrate excellent expressive capacities and outperform existing methods on common benchmark datasets. However, model generalizability, i.e., the performance gap between the seen and unseen sets, degrades rapidly as the sequence similarity decreases. The same trends are observed from several recent DL and machine learning models. And an inverse correlation between performance and generalizability is revealed collectively across all learning-based models with wide-ranging architectures and sizes. We further quantitate how generalizability depends on sequence and structure identity scores via pairwise alignment, providing unique quantitative insights into the limitations of statistical learning. Generalizability thus poses a major hurdle for deploying de novo DL models in practice and various pathways for future advances are discussed.

Our understanding of how speed and persistence of cell migration affects the growth rate and size of tumors remains incomplete. To address this, we developed a mathematical model wherein cells migrate in two-dimensional space, divide, die or intravasate into the vasculature. Exploring a wide range of speed and persistence combinations, we find that tumor growth positively correlates with increasing speed and higher persistence. As a biologically relevant example, we focused on Golgi fragmentation, a phenomenon often linked to alterations of cell migration. Golgi fragmentation was induced by depletion of Giantin, a Golgi matrix protein, the downregulation of which correlates with poor patient survival. Applying the experimentally obtained migration and invasion traits of Giantin depleted breast cancer cells to our mathematical model, we predict that loss of Giantin increases the number of intravasating cells. This prediction was validated, by showing that circulating tumor cells express significantly less Giantin than primary tumor cells. Altogether, our computational model identifies cell migration traits that regulate tumor progression and uncovers a role of Giantin in breast cancer progression.

Dorsal horn of the spinal cord is an important crossroad of pain neuraxis, especially for the neuronal plasticity mechanisms that can lead to chronic pain states. Windup is a well-known spinal pain facilitation process initially described several decades ago, but its exact mechanism is still not fully understood. Here, we combine both ex vivo and in vivo electrophysiological recordings of rat spinal neurons with computational modeling to demonstrate a role for ASIC1a-containing channels in the windup process. Spinal application of the ASIC1a inhibitory venom peptides mambalgin-1 and psalmotoxin-1 (PcTx1) significantly reduces the ability of deep wide dynamic range (WDR) neurons to develop windup in vivo. All deep WDR-like neurons recorded from spinal slices exhibit an ASIC current with biophysical and pharmacological characteristics consistent with functional expression of ASIC1a homomeric channels. A computational model of WDR neuron supplemented with different ASIC1a channel parameters accurately reproduces the experimental data, further supporting a positive contribution of these channels to windup. It also predicts a calcium-dependent windup decrease for elevated ASIC conductances, a phenomenon that was experimentally validated using the Texas coral snake ASIC-activating toxin (MitTx) and calcium-activated potassium channel inhibitory peptides (apamin and iberiotoxin). This study supports a dual contribution to windup of calcium permeable ASIC1a channels in deep laminae projecting neurons, promoting it upon moderate channel activity, but ultimately leading to calcium-dependent windup inhibition associated to potassium channels when activity increases.

The computational principles adopted by the hippocampus in associative memory (AM) tasks have been one of the most studied topics in computational and theoretical neuroscience. Recent theories suggested that AM and the predictive activities of the hippocampus could be described within a unitary account, and that predictive coding underlies the computations supporting AM in the hippocampus. Following this theory, a computational model based on classical hierarchical predictive networks was proposed and was shown to perform well in various AM tasks. However, this fully hierarchical model did not incorporate recurrent connections, an architectural component of the CA3 region of the hippocampus that is crucial for AM. This makes the structure of the model inconsistent with the known connectivity of CA3 and classical recurrent models such as Hopfield Networks, which learn the covariance of inputs through their recurrent connections to perform AM. Earlier PC models that learn the covariance information of inputs explicitly via recurrent connections seem to be a solution to these issues. Here, we show that although these models can perform AM, they do it in an implausible and numerically unstable way. Instead, we propose alternatives to these earlier covariance-learning predictive coding networks, which learn the covariance information implicitly and plausibly, and can use dendritic structures to encode prediction errors. We show analytically that our proposed models are perfectly equivalent to the earlier predictive coding model learning covariance explicitly, and encounter no numerical issues when performing AM tasks in practice. We further show that our models can be combined with hierarchical predictive coding networks to model the hippocampo-neocortical interactions. Our models provide a biologically plausible approach to modelling the hippocampal network, pointing to a potential computational mechanism during hippocampal memory formation and recall, which employs both predictive coding and covariance learning based on the recurrent network structure of the hippocampus.

[This corrects the article DOI: 10.1371/journal.pcbi.1010383.].

Anovulation refers to a menstrual cycle characterized by the absence of ovulation. Exogenous hormones such as synthetic progesterone and estrogen have been used to attain this state to achieve contraception. However, large doses are associated with adverse effects such as increased risk for thrombosis and myocardial infarction. This study utilizes optimal control theory on a modified menstrual cycle model to determine the minimum total exogenous estrogen/progesterone dose, and timing of administration to induce anovulation. The mathematical model correctly predicts the mean daily levels of pituitary hormones LH and FSH, and ovarian hormones E2, P4, and Inh throughout a normal menstrual cycle and reflects the reduction in these hormone levels caused by exogenous estrogen and/or progesterone. Results show that it is possible to reduce the total dose by 92% in estrogen monotherapy, 43% in progesterone monotherapy, and that it is most effective to deliver the estrogen contraceptive in the mid follicular phase. Finally, we show that by combining estrogen and progesterone the dose can be lowered even more. These results may give clinicians insights into optimal formulations and schedule of therapy that can suppress ovulation.

The long-term behaviors of biochemical systems are often described by their steady states. Deriving these states directly for complex networks arising from real-world applications, however, is often challenging. Recent work has consequently focused on network-based approaches. Specifically, biochemical reaction networks are transformed into weakly reversible and deficiency zero generalized networks, which allows the derivation of their analytic steady states. Identifying this transformation, however, can be challenging for large and complex networks. In this paper, we address this difficulty by breaking the complex network into smaller independent subnetworks and then transforming the subnetworks to derive the analytic steady states of each subnetwork. We show that stitching these solutions together leads to the analytic steady states of the original network. To facilitate this process, we develop a user-friendly and publicly available package, COMPILES (COMPutIng anaLytic stEady States). With COMPILES, we can easily test the presence of bistability of a CRISPRi toggle switch model, which was previously investigated via tremendous number of numerical simulations and within a limited range of parameters. Furthermore, COMPILES can be used to identify absolute concentration robustness (ACR), the property of a system that maintains the concentration of particular species at a steady state regardless of any initial concentrations. Specifically, our approach completely identifies all the species with and without ACR in a complex insulin model. Our method provides an effective approach to analyzing and understanding complex biochemical systems.

BACKGROUND: Antibiotic treatments are often associated with a late slowdown in bacterial killing. This separates the killing of bacteria into at least two distinct phases: a quick phase followed by a slower phase, the latter of which is linked to treatment success. Current mechanistic explanations for the in vitro slowdown are either antibiotic persistence or heteroresistance. Persistence is defined as the switching back and forth between susceptible and non-susceptible states, while heteroresistance is defined as the coexistence of bacteria with heterogeneous susceptibilities. Both are also thought to cause a slowdown in the decline of bacterial populations in patients and therefore complicate and prolong antibiotic treatments. Reduced bacterial death rates over time are also observed within tuberculosis patients, yet the mechanistic reasons for this are unknown and therefore the strategies to mitigate them are also unknown. METHODS AND FINDINGS: We analyse a dose ranging trial for rifampicin in tuberculosis patients and show that there is a slowdown in the decline of bacteria. We show that the late phase of bacterial killing depends more on the peak drug concentrations than the total drug exposure. We compare these to pharmacokinetic-pharmacodynamic models of rifampicin heteroresistance and persistence. We find that the observation on the slow phase's dependence on pharmacokinetic measures, specifically peak concentrations are only compatible with models of heteroresistance and incompatible with models of persistence. The quantitative agreement between heteroresistance models and observations is very good ([Formula: see text]). To corroborate the importance of the slowdown, we validate our results by estimating the time to sputum culture conversion and compare the results to a different dose ranging trial. CONCLUSIONS: Our findings indicate that higher doses, specifically higher peak concentrations may be used to optimize rifampicin treatments by accelerating bacterial killing in the slow phase. It adds to the growing body of literature supporting higher rifampicin doses for shortening tuberculosis treatments.

Despite the accumulation of data and studies, deciphering animal vocal communication remains challenging. In most cases, researchers must deal with the sparse recordings composing Small, Unbalanced, Noisy, but Genuine (SUNG) datasets. SUNG datasets are characterized by a limited number of recordings, most often noisy, and unbalanced in number between the individuals or categories of vocalizations. SUNG datasets therefore offer a valuable but inevitably distorted vision of communication systems. Adopting the best practices in their analysis is essential to effectively extract the available information and draw reliable conclusions. Here we show that the most recent advances in machine learning applied to a SUNG dataset succeed in unraveling the complex vocal repertoire of the bonobo, and we propose a workflow that can be effective with other animal species. We implement acoustic parameterization in three feature spaces and run a Supervised Uniform Manifold Approximation and Projection (S-UMAP) to evaluate how call types and individual signatures cluster in the bonobo acoustic space. We then implement three classification algorithms (Support Vector Machine, xgboost, neural networks) and their combination to explore the structure and variability of bonobo calls, as well as the robustness of the individual signature they encode. We underscore how classification performance is affected by the feature set and identify the most informative features. In addition, we highlight the need to address data leakage in the evaluation of classification performance to avoid misleading interpretations. Our results lead to identifying several practical approaches that are generalizable to any other animal communication system. To improve the reliability and replicability of vocal communication studies with SUNG datasets, we thus recommend: i) comparing several acoustic parameterizations; ii) visualizing the dataset with supervised UMAP to examine the species acoustic space; iii) adopting Support Vector Machines as the baseline classification approach; iv) explicitly evaluating data leakage and possibly implementing a mitigation strategy.

Protein design is a technique to engineer proteins by permuting amino acids in the sequence to obtain novel functionalities. However, exploring all possible combinations of amino acids is generally impossible due to the exponential growth of possibilities with the number of designable sites. The present work introduces circuits implementing a pure quantum approach, Grover's algorithm, to solve protein design problems. Our algorithms can adjust to implement any custom pair-wise energy tables and protein structure models. Moreover, the algorithm's oracle is designed to consist of only adder functions. Quantum computer simulators validate the practicality of our circuits, containing up to 234 qubits. However, a smaller circuit is implemented on real quantum devices. Our results show that using [Formula: see text] iterations, the circuits find the correct results among all N possibilities, providing the expected quadratic speed up of Grover's algorithm over classical methods (i.e., [Formula: see text]).

Memories may be encoded in the brain via strongly interconnected groups of neurons, called assemblies. The concept of Hebbian plasticity suggests that these assemblies are generated through synaptic plasticity, strengthening the recurrent connections within select groups of neurons that receive correlated stimulation. To remain stable in absence of such stimulation, the assemblies need to be self-reinforcing under the plasticity rule. Previous models of such assembly maintenance require additional mechanisms of fast homeostatic plasticity often with biologically implausible timescales. Here we provide a model of neuronal assembly generation and maintenance purely based on spike-timing-dependent plasticity (STDP) between excitatory neurons. It uses irregularly and stochastically spiking neurons and STDP that depresses connections of uncorrelated neurons. We find that assemblies do not grow beyond a certain size, because temporally imprecisely correlated spikes dominate the plasticity in large assemblies. Assemblies in the model can be learned or spontaneously emerge. The model allows for prominent, stable overlap structures between static assemblies. Further, assemblies can drift, particularly according to a novel, transient overlap-based mechanism. Finally the model indicates that assemblies grow in the aging brain, where connectivity decreases.

Spatiotemporal oscillations underlie all cognitive brain functions. Large-scale brain models, constrained by neuroimaging data, aim to trace the principles underlying such macroscopic neural activity from the intricate and multi-scale structure of the brain. Despite substantial progress in the field, many aspects about the mechanisms behind the onset of spatiotemporal neural dynamics are still unknown. In this work we establish a simple framework for the emergence of complex brain dynamics, including high-dimensional chaos and travelling waves. The model consists of a complex network of 90 brain regions, whose structural connectivity is obtained from tractography data. The activity of each brain area is governed by a Jansen neural mass model and we normalize the total input received by each node so it amounts the same across all brain areas. This assumption allows for the existence of an homogeneous invariant manifold, i.e., a set of different stationary and oscillatory states in which all nodes behave identically. Stability analysis of these homogeneous solutions unveils a transverse instability of the synchronized state, which gives rise to different types of spatiotemporal dynamics, such as chaotic alpha activity. Additionally, we illustrate the ubiquity of this route towards complex spatiotemporal activity in a network of next generation neural mass models. Altogehter, our results unveil the bifurcation landscape that underlies the emergence of function from structure in the brain.

Transcranial direct current stimulation (tDCS) can noninvasively modulate behavior, cognition, and physiologic brain functions depending on polarity and dose of stimulation as well as montage of electrodes. Concurrent tDCS-fMRI presents a novel way to explore the parameter space of non-invasive brain stimulation and to inform the experimenter as well as the participant if a targeted brain region or a network of spatially separate brain regions has been engaged and modulated. We compared a multi-electrode (ME) with a single electrode (SE) montage and both active conditions with a no-stimulation (NS) control condition to assess the engagement of a brain network and the ability of different electrode montages to modulate network activity. The multi-electrode montage targeted nodal regions of the right Arcuate Fasciculus Network (AFN) with anodal electrodes placed over the skull position of the posterior superior temporal/middle temporal gyrus (STG/MTG), supramarginal gyrus (SMG), posterior inferior frontal gyrus (IFG) and a return cathodal electrode over the left supraorbital region. In comparison, the single electrode montage used only one anodal electrode over a nodal brain region of the AFN, but varied the location between STG/MTG, SMG, and posterior IFG for different participants. Whole-brain rs-fMRI was obtained approximately every three seconds. The tDCS-stimulator was turned on at 3 minutes after the scanning started. A 4D rs-fMRI data set was converted to dynamic functional connectivity (DFC) matrices using a set of ROI pairs belonging to the AFN as well as other unrelated brain networks. In this study, we evaluated the performance of five algorithms to classify the DFC matrices from the three conditions (ME, SE, NS) into three different categories. The highest accuracy of 0.92 was obtained for the classification of the ME condition using the K Nearest Neighbor (KNN) algorithm. In other words, applying the classification algorithm allowed us to identify the engagement of the AFN and the ME condition was the best montage to achieve such an engagement. The top 5 ROI pairs that made a major contribution to the classification of participant's rs-fMRI data were identified using model performance parameters; ROI pairs were mainly located within the right AFN. This proof-of-concept study using a classification algorithm approach can be expanded to create a near real-time feedback system at a participant level to detect the engagement and modulation of a brain network that spans multiple brain lobes.

Adenosine-to-inosine RNA editing is essential to prevent undesired immune activation. This diverse process alters the genetic content of the RNA and may recode proteins, change splice sites and miRNA targets, and mimic genomic mutations. Recent studies have associated or implicated aberrant editing with pathological conditions, including cancer, autoimmune diseases, and neurological and psychiatric conditions. RNA editing patterns in cardiovascular tissues have not been investigated systematically so far, and little is known about its potential role in cardiac diseases. Some hints suggest robust editing in this system, including the fact that ADARB1 (ADAR2), the main coding-sequence editor, is most highly expressed in these tissues. Here we characterized RNA editing in the heart and arteries and examined a contributory role to the development of atherosclerosis and two structural heart diseases -Ischemic and Dilated Cardiomyopathies. Analyzing hundreds of RNA-seq samples taken from the heart and arteries of cardiac patients and controls, we find that global editing, alongside inflammatory gene expression, is increased in patients with atherosclerosis, cardiomyopathies, and heart failure. We describe a single recoding editing site and suggest it as a target for focused research. This recoding editing site in the IGFBP7 gene is one of the only evolutionary conserved sites between mammals, and we found it exhibits consistently increased levels of editing in these patients. Our findings reveal that RNA editing is abundant in arteries and is elevated in several key cardiovascular conditions. They thus provide a roadmap for basic and translational research of RNA as a mediator of atherosclerosis and non-genetic cardiomyopathies.

Neurons, represented as a tree structure of morphology, have various distinguished branches of dendrites. Different types of synaptic receptors distributed over dendrites are responsible for receiving inputs from other neurons. NMDA receptors (NMDARs) are expressed as excitatory units, and play a key physiological role in synaptic function. Although NMDARs are widely expressed in most types of neurons, they play a different role in the cerebellar Purkinje cells (PCs). Utilizing a computational PC model with detailed dendritic morphology, we explored the role of NMDARs at different parts of dendritic branches and regions. We found somatic responses can switch from silent, to simple spikes and complex spikes, depending on specific dendritic branches. Detailed examination of the dendrites regarding their diameters and distance to soma revealed diverse response patterns, yet explain two firing modes, simple and complex spike. Taken together, these results suggest that NMDARs play an important role in controlling excitability sensitivity while taking into account the factor of dendritic properties. Given the complexity of neural morphology varying in cell types, our work suggests that the functional role of NMDARs is not stereotyped but highly interwoven with local properties of neuronal structure.

Reduced cortical inhibition by somatostatin-expressing (SST) interneurons has been strongly associated with treatment-resistant depression. However, due to technical limitations it is impossible to establish experimentally in humans whether the effects of reduced SST interneuron inhibition on microcircuit activity have signatures detectable in clinically-relevant brain signals such as electroencephalography (EEG). To overcome these limitations, we simulated resting-state activity and EEG using detailed models of human cortical microcircuits with normal (healthy) or reduced SST interneuron inhibition (depression), and found that depression microcircuits exhibited increased theta, alpha and low beta power (4-16 Hz). The changes in depression involved a combination of an aperiodic broadband and periodic theta components. We then demonstrated the specificity of the EEG signatures of reduced SST interneuron inhibition by showing they were distinct from those corresponding to reduced parvalbumin-expressing (PV) interneuron inhibition. Our study thus links SST interneuron inhibition level to distinct features in EEG simulated from detailed human microcircuits, which can serve to better identify mechanistic subtypes of depression using EEG, and non-invasively monitor modulation of cortical inhibition.

Cellular distributions of the sphingolipid ceramide-1-phosphate (C1P) impact essential biological processes. C1P levels are spatiotemporally regulated by ceramide-1-phosphate transfer protein (CPTP), which efficiently shuttles C1P between organelle membranes. Yet, how CPTP rapidly extracts and inserts C1P into a membrane remains unknown. Here, we devise a multiscale simulation approach to elucidate biophysical details of CPTP-mediated C1P transport. We find that CPTP binds a membrane poised to extract and insert C1P and that membrane binding promotes conformational changes in CPTP that facilitate C1P uptake and release. By significantly disrupting a lipid's local hydrophobic environment in the membrane, CPTP lowers the activation free energy barrier for passive C1P desorption and enhances C1P extraction from the membrane. Upon uptake of C1P, further conformational changes may aid membrane unbinding in a manner reminiscent of the electrostatic switching mechanism used by other lipid transfer proteins. Insertion of C1P into an acceptor membrane, eased by a decrease in membrane order by CPTP, restarts the transfer cycle. Most notably, we provide molecular evidence for CPTP's ability to catalyze C1P extraction by breaking hydrophobic C1P-membrane contacts with compensatory hydrophobic lipid-protein contacts. Our work, thus, provides biophysical insights into how CPTP efficiently traffics C1P between membranes to maintain sphingolipid homeostasis and, additionally, presents a simulation method aptly suited for uncovering the catalytic mechanisms of other lipid transfer proteins.