Paramecium is a large unicellular organism that swims in fresh water using cilia. When stimulated by various means (mechanically, chemically, optically, thermally), it often swims backward then turns and swims forward again in a new direction: this is called the avoiding reaction. This reaction is triggered by a calcium-based action potential. For this reason, several authors have called Paramecium the "swimming neuron". Here we present an empirically constrained model of its action potential based on electrophysiology experiments on live immobilized paramecia, together with simultaneous measurement of ciliary beating using particle image velocimetry. Using these measurements and additional behavioral measurements of free swimming, we extend the electrophysiological model by coupling calcium concentration to kinematic parameters, turning it into a swimming model. In this way, we obtain a model of autonomously behaving Paramecium. Finally, we demonstrate how the modeled organism interacts with an environment, can follow gradients and display collective behavior. This work provides a modeling basis for investigating the physiological basis of autonomous behavior of Paramecium in ecological environments.

Microbial specialised metabolism is full of valuable natural products that are applied clinically, agriculturally, and industrially. The genes that encode their biosynthesis are often physically clustered on the genome in biosynthetic gene clusters (BGCs). Many BGCs consist of multiple groups of co-evolving genes called sub-clusters that are responsible for the biosynthesis of a specific chemical moiety in a natural product. Sub-clusters therefore provide an important link between the structures of a natural product and its BGC, which can be leveraged for predicting natural product structures from sequence, as well as for linking chemical structures and metabolomics-derived mass features to BGCs. While some initial computational methodologies have been devised for sub-cluster detection, current approaches are not scalable, have only been run on small and outdated datasets, or produce an impractically large number of possible sub-clusters to mine through. Here, we constructed a scalable method for unsupervised sub-cluster detection, called iPRESTO, based on topic modelling and statistical analysis of co-occurrence patterns of enzyme-coding protein families. iPRESTO was used to mine sub-clusters across 150,000 prokaryotic BGCs from antiSMASH-DB. After annotating a fraction of the resulting sub-cluster families, we could predict a substructure for 16% of the antiSMASH-DB BGCs. Additionally, our method was able to confirm 83% of the experimentally characterised sub-clusters in MIBiG reference BGCs. Based on iPRESTO-detected sub-clusters, we could correctly identify the BGCs for xenorhabdin and salbostatin biosynthesis (which had not yet been annotated in BGC databases), as well as propose a candidate BGC for akashin biosynthesis. Additionally, we show for a collection of 145 actinobacteria how substructures can aid in linking BGCs to molecules by correlating iPRESTO-detected sub-clusters to MS/MS-derived Mass2Motifs substructure patterns. This work paves the way for deeper functional and structural annotation of microbial BGCs by improved linking of orphan molecules to their cognate gene clusters, thus facilitating accelerated natural product discovery.

Epidemiological and clinical evidence indicates that humans infected with the 1918 pandemic H1N1 influenza virus and highly pathogenic avian H5N1 influenza viruses often displayed severe lung pathology. High viral load and extensive infiltration of macrophages are the hallmarks of highly pathogenic (HP) influenza viral infections. However, it remains unclear what biological mechanisms primarily determine the observed difference in the kinetics of viral load and macrophages between HP and low pathogenic (LP) viral infections, and how the mechanistic differences are associated with viral pathogenicity. In this study, we develop a mathematical model of viral dynamics that includes the dynamics of different macrophage populations and interferon. We fit the model to in vivo kinetic data of viral load and macrophage level from BALB/c mice infected with an HP or LP strain of H1N1/H5N1 virus using Bayesian inference. Our primary finding is that HP viruses have a higher viral infection rate, a lower interferon production rate and a lower macrophage recruitment rate compared to LP viruses, which are strongly associated with more severe tissue damage (quantified by a higher percentage of epithelial cell loss). We also quantify the relative contribution of macrophages to viral clearance and find that macrophages do not play a dominant role in the direct clearance of free viruses although their role in mediating immune responses such as interferon production is crucial. Our work provides new insight into the mechanisms that convey the observed difference in viral and macrophage kinetics between HP and LP infections and establishes an improved model-fitting framework to enhance the analysis of new data on viral pathogenicity.

Humans constantly assess the appearance of materials to plan actions, such as stepping on icy roads without slipping. Visual inference of materials is important but challenging because a given material can appear dramatically different in various scenes. This problem especially stands out for translucent materials, whose appearance strongly depends on lighting, geometry, and viewpoint. Despite this, humans can still distinguish between different materials, and it remains unsolved how to systematically discover visual features pertinent to material inference from natural images. Here, we develop an unsupervised style-based image generation model to identify perceptually relevant dimensions for translucent material appearances from photographs. We find our model, with its layer-wise latent representation, can synthesize images of diverse and realistic materials. Importantly, without supervision, human-understandable scene attributes, including the object's shape, material, and body color, spontaneously emerge in the model's layer-wise latent space in a scale-specific manner. By embedding an image into the learned latent space, we can manipulate specific layers' latent code to modify the appearance of the object in the image. Specifically, we find that manipulation on the early-layers (coarse spatial scale) transforms the object's shape, while manipulation on the later-layers (fine spatial scale) modifies its body color. The middle-layers of the latent space selectively encode translucency features and manipulation of such layers coherently modifies the translucency appearance, without changing the object's shape or body color. Moreover, we find the middle-layers of the latent space can successfully predict human translucency ratings, suggesting that translucent impressions are established in mid-to-low spatial scale features. This layer-wise latent representation allows us to systematically discover perceptually relevant image features for human translucency perception. Together, our findings reveal that learning the scale-specific statistical structure of natural images might be crucial for humans to efficiently represent material properties across contexts.

The mouse gustatory cortex (GC) is involved in taste-guided decision-making in addition to sensory processing. Rodent GC exhibits metastable neural dynamics during ongoing and stimulus-evoked activity, but how these dynamics evolve in the context of a taste-based decision-making task remains unclear. Here we employ analytical and modeling approaches to i) extract metastable dynamics in ensemble spiking activity recorded from the GC of mice performing a perceptual decision-making task; ii) investigate the computational mechanisms underlying GC metastability in this task; and iii) establish a relationship between GC dynamics and behavioral performance. Our results show that activity in GC during perceptual decision-making is metastable and that this metastability may serve as a substrate for sequentially encoding sensory, abstract cue, and decision information over time. Perturbations of the model's metastable dynamics indicate that boosting inhibition in different coding epochs differentially impacts network performance, explaining a counterintuitive effect of GC optogenetic silencing on mouse behavior.

The spread of seizures across brain networks is the main impairing factor, often leading to loss-of-consciousness, in people with epilepsy. Despite advances in recording and modeling brain activity, uncovering the nature of seizure spreading dynamics remains an important challenge to understanding and treating pharmacologically resistant epilepsy. To address this challenge, we introduce a new probabilistic model that captures the spreading dynamics in patient-specific complex networks. Network connectivity and interaction time delays between brain areas were estimated from white-matter tractography. The model's computational tractability allows it to play an important complementary role to more detailed models of seizure dynamics. We illustrate model fitting and predictive performance in the context of patient-specific Epileptor networks. We derive the phase diagram of spread size (order parameter) as a function of brain excitability and global connectivity strength, for different patient-specific networks. Phase diagrams allow the prediction of whether a seizure will spread depending on excitability and connectivity strength. In addition, model simulations predict the temporal order of seizure spread across network nodes. Furthermore, we show that the order parameter can exhibit both discontinuous and continuous (critical) phase transitions as neural excitability and connectivity strength are varied. Existence of a critical point, where response functions and fluctuations in spread size show power-law divergence with respect to control parameters, is supported by mean-field approximations and finite-size scaling analyses. Notably, the critical point separates two distinct regimes of spreading dynamics characterized by unimodal and bimodal spread-size distributions. Our study sheds new light on the nature of phase transitions and fluctuations in seizure spreading dynamics. We expect it to play an important role in the development of closed-loop stimulation approaches for preventing seizure spread in pharmacologically resistant epilepsy. Our findings may also be of interest to related models of spreading dynamics in epidemiology, biology, finance, and statistical physics.

Upper respiratory tract infections (URTIs) represent a large strain on primary health resources. To mitigate URTI transmission and public health burdens, it is important to pre-empt and provide forward guidance on URTI burden, while taking into account various facets which influence URTI transmission. This is so that appropriate public health measures can be taken to mitigate strain on primary care resources. This study describes a new approach to forecasting URTIs which can be used for national public health resource planning. Specifically, using environmental and disease data comprising more than 1000 dimensions, we developed sub-models which optimizes model explainability, in-sample model fit, predictive accuracy and combines many weaker predictors over a 2-month time horizon to generate direct, point forecasts over a 1-8 week ahead forecast horizon. Predictive performance was evaluated using rolling out-of-sample forecast assessment within both periods with/without structural breaks in transmission over the period of 2012-2022. We showed that forecast combinations of 5 other forecasting models had better and more consistent predictive performance than other modelling approaches, over periods with and without structural breaks in transmission dynamics. Furthermore, epidemiological analysis on high dimensional data was enabled using post-selection inference, to show the dynamic association between lower temperature, increases in past relative humidity and absolute humidity and increased URTIs attendance. The methods proposed can be used for outbreak preparedness and guide healthcare resource planning, in both stable periods of transmission and periods where structural breaks in data occur.

To adapt to their environments, animals learn associations between sensory stimuli and unconditioned stimuli. In invertebrates, olfactory associative learning primarily occurs in the mushroom body, which is segregated into separate compartments. Within each compartment, Kenyon cells (KCs) encoding sparse odor representations project onto mushroom body output neurons (MBONs) whose outputs guide behavior. Associated with each compartment is a dopamine neuron (DAN) that modulates plasticity of the KC-MBON synapses within the compartment. Interestingly, DAN-induced plasticity of the KC-MBON synapse is imbalanced in the sense that it only weakens the synapse and is temporally sparse. We propose a normative mechanistic model of the MBON as a linear discriminant analysis (LDA) classifier that predicts the presence of an unconditioned stimulus (class identity) given a KC odor representation (feature vector). Starting from a principled LDA objective function and under the assumption of temporally sparse DAN activity, we derive an online algorithm which maps onto the mushroom body compartment. Our model accounts for the imbalanced learning at the KC-MBON synapse and makes testable predictions that provide clear contrasts with existing models.

How cell specification can be controlled in a reproducible manner is a fundamental question in developmental biology. In ascidians, a group of invertebrate chordates, geometry plays a key role in achieving this control. Here, we use mathematical modeling to demonstrate that geometry dictates the neural-epidermal cell fate choice in the 32-cell stage ascidian embryo by a two-step process involving first the modulation of ERK signaling and second, the expression of the neural marker gene, Otx. The model describes signal transduction by the ERK pathway that is stimulated by FGF and attenuated by ephrin, and ERK-mediated control of Otx gene expression, which involves both an activator and a repressor of ETS-family transcription factors. Considering the measured area of cell surface contacts with FGF- or ephrin-expressing cells as inputs, the solutions of the model reproduce the experimental observations about ERK activation and Otx expression in the different cells under normal and perturbed conditions. Sensitivity analyses and computations of Hill coefficients allow us to quantify the robustness of the specification mechanism controlled by cell surface area and to identify the respective role played by each signaling input. Simulations also predict in which conditions the dual control of gene expression by an activator and a repressor that are both under the control of ERK can induce a robust ON/OFF control of neural fate induction.

A topic of growing interest in computational neuroscience is the discovery of fundamental principles underlying global dynamics and the self-organization of the brain. In particular, the notion that the brain operates near criticality has gained considerable support, and recent work has shown that the dynamics of different brain states may be modeled by pairwise maximum entropy Ising models at various distances from a phase transition, i.e., from criticality. Here we aim to characterize two brain states (psychedelics-induced and placebo) as captured by functional magnetic resonance imaging (fMRI), with features derived from the Ising spin model formalism (system temperature, critical point, susceptibility) and from algorithmic complexity. We hypothesized, along the lines of the entropic brain hypothesis, that psychedelics drive brain dynamics into a more disordered state at a higher Ising temperature and increased complexity. We analyze resting state blood-oxygen-level-dependent (BOLD) fMRI data collected in an earlier study from fifteen subjects in a control condition (placebo) and during ingestion of lysergic acid diethylamide (LSD). Working with the automated anatomical labeling (AAL) brain parcellation, we first create "archetype" Ising models representative of the entire dataset (global) and of the data in each condition. Remarkably, we find that such archetypes exhibit a strong correlation with an average structural connectome template obtained from dMRI (r = 0.6). We compare the archetypes from the two conditions and find that the Ising connectivity in the LSD condition is lower than the placebo one, especially in homotopic links (interhemispheric connectivity), reflecting a significant decrease of homotopic functional connectivity in the LSD condition. The global archetype is then personalized for each individual and condition by adjusting the system temperature. The resulting temperatures are all near but above the critical point of the model in the paramagnetic (disordered) phase. The individualized Ising temperatures are higher in the LSD condition than the placebo condition (p = 9 × 10-5). Next, we estimate the Lempel-Ziv-Welch (LZW) complexity of the binarized BOLD data and the synthetic data generated with the individualized model using the Metropolis algorithm for each participant and condition. The LZW complexity computed from experimental data reveals a weak statistical relationship with condition (p = 0.04 one-tailed Wilcoxon test) and none with Ising temperature (r(13) = 0.13, p = 0.65), presumably because of the limited length of the BOLD time series. Similarly, we explore complexity using the block decomposition method (BDM), a more advanced method for estimating algorithmic complexity. The BDM complexity of the experimental data displays a significant correlation with Ising temperature (r(13) = 0.56, p = 0.03) and a weak but significant correlation with condition (p = 0.04, one-tailed Wilcoxon test). This study suggests that the effects of LSD increase the complexity of brain dynamics by loosening interhemispheric connectivity-especially homotopic links. In agreement with earlier work using the Ising formalism with BOLD data, we find the brain state in the placebo condition is already above the critical point, with LSD resulting in a shift further away from criticality into a more disordered state.

Multiple coronaviruses including MERS-CoV causing Middle East Respiratory Syndrome, SARS-CoV causing SARS, and SARS-CoV-2 causing COVID-19, use a mechanism known as -1 programmed ribosomal frameshifting (-1 PRF) to replicate. SARS-CoV-2 possesses a unique RNA pseudoknotted structure that stimulates -1 PRF. Targeting -1 PRF in SARS-CoV-2 to impair viral replication can improve patients' prognoses. Crucial to developing these therapies is understanding the structure of the SARS-CoV-2 -1 PRF pseudoknot. Our goal is to expand knowledge of -1 PRF structural conformations. Following a structural alignment approach, we identify similarities in -1 PRF pseudoknots of SARS-CoV-2, SARS-CoV, and MERS-CoV. We provide in-depth analysis of the SARS-CoV-2 and MERS-CoV -1 PRF pseudoknots, including reference and noteworthy mutated sequences. To better understand the impact of mutations, we provide insight on -1 PRF pseudoknot sequence mutations and their effect on resulting structures. We introduce Shapify, a novel algorithm that given an RNA sequence incorporates structural reactivity (SHAPE) data and partial structure information to output an RNA secondary structure prediction within a biologically sound hierarchical folding approach. Shapify enhances our understanding of SARS-CoV-2 -1 PRF pseudoknot conformations by providing energetically favourable predictions that are relevant to structure-function and may correlate with -1 PRF efficiency. Applied to the SARS-CoV-2 -1 PRF pseudoknot, Shapify unveils previously unknown paths from initial stems to pseudoknotted structures. By contextualizing our work with available experimental data, our structure predictions motivate future RNA structure-function research and can aid 3-D modeling of pseudoknots.

Environmentally induced or epigenetic-related beta-cell dysfunction and insulin resistance play a critical role in the progression to diabetes. We developed a mathematical modeling framework capable of studying the progression to diabetes incorporating various diabetogenic factors. Considering the heightened risk of beta-cell defects induced by obesity, we focused on the obesity-diabetes model to further investigate the influence of obesity on beta-cell function and glucose regulation. The model characterizes individualized glucose and insulin dynamics over the span of a lifetime. We then fit the model to the longitudinal data of the Pima Indian population, which captures both the fluctuations and long-term trends of glucose levels. As predicted, controlling or eradicating the obesity-related factor can alleviate, postpone, or even reverse diabetes. Furthermore, our results reveal that distinct abnormalities of beta-cell function and levels of insulin resistance among individuals contribute to different risks of diabetes. This study may encourage precise interventions to prevent diabetes and facilitate individualized patient treatment.

Influenza pandemics typically occur in multiple waves of infection, often associated with initial emergence of a novel virus, followed (in temperate regions) by a resurgence accompanying the onset of the annual influenza season. Here, we examined whether data collected from an initial pandemic wave could be informative, for the need to implement non-pharmaceutical measures in any resurgent wave. Drawing from the 2009 H1N1 pandemic in 10 states in the USA, we calibrated simple mathematical models of influenza transmission dynamics to data for laboratory confirmed hospitalisations during the initial 'spring' wave. We then projected pandemic outcomes (cumulative hospitalisations) during the fall wave, and compared these projections with data. Model results showed reasonable agreement for all states that reported a substantial number of cases in the spring wave. Using this model we propose a probabilistic decision framework that can be used to determine the need for preemptive measures such as postponing school openings, in advance of a fall wave. This work illustrates how model-based evidence synthesis, in real-time during an early pandemic wave, could be used to inform timely decisions for pandemic response.

Two decades of in vivo imaging have revealed how diverse T-cell motion patterns can be. Such recordings have sparked the notion of search "strategies": T cells may have evolved ways to search for antigen efficiently depending on the task at hand. Mathematical models have indeed confirmed that several observed T-cell migration patterns resemble a theoretical optimum; for example, frequent turning, stop-and-go motion, or alternating short and long motile runs have all been interpreted as deliberately tuned behaviours, optimising the cell's chance of finding antigen. But the same behaviours could also arise simply because T cells cannot follow a straight, regular path through the tight spaces they navigate. Even if T cells do follow a theoretically optimal pattern, the question remains: which parts of that pattern have truly been evolved for search, and which merely reflect constraints from the cell's migration machinery and surroundings? We here employ an approach from the field of evolutionary biology to examine how cells might evolve search strategies under realistic constraints. Using a cellular Potts model (CPM), where motion arises from intracellular dynamics interacting with cell shape and a constraining environment, we simulate evolutionary optimization of a simple task: explore as much area as possible. We find that our simulated cells indeed evolve their motility patterns. But the evolved behaviors are not shaped solely by what is functionally optimal; importantly, they also reflect mechanistic constraints. Cells in our model evolve several motility characteristics previously attributed to search optimisation-even though these features are not beneficial for the task given here. Our results stress that search patterns may evolve for other reasons than being "optimal". In part, they may be the inevitable side effects of interactions between cell shape, intracellular dynamics, and the diverse environments T cells face in vivo.

Transmission of many communicable diseases depends on proximity contacts among humans. Modeling the dynamics of proximity contacts can help determine whether an outbreak is likely to trigger an epidemic. While the advent of commodity mobile devices has eased the collection of proximity contact data, battery capacity and associated costs impose tradeoffs between the observation frequency and scanning duration used for contact detection. The choice of observation frequency should depend on the characteristics of a particular pathogen and accompanying disease. We downsampled data from five contact network studies, each measuring participant-participant contact every 5 minutes for durations of four or more weeks. These studies included a total of 284 participants and exhibited different community structures. We found that for epidemiological models employing high-resolution proximity data, both the observation method and observation frequency configured to collect proximity data impact the simulation results. This impact is subject to the population's characteristics as well as pathogen infectiousness. By comparing the performance of two observation methods, we found that in most cases, half-hourly Bluetooth discovery for one minute can collect proximity data that allows agent-based transmission models to produce a reasonable estimation of the attack rate, but more frequent Bluetooth discovery is preferred to model individual infection risks or for highly transmissible pathogens. Our findings inform the empirical basis for guidelines to inform data collection that is both efficient and effective.

Models of evolution of simple languages have typically assumed full alignment of the speaker and listeners interests, with perfect understanding representing the optimal outcome for both parties. In more realistic settings, communicating individuals will often desire different outcomes from one another. Previous work has shown that misalignment of speaker-listener interests reduces the maximum informativeness among Nash-equilibrium languages, and that multiple equilibrium languages (with different degrees of informativeness) are supported. We study the stochastic evolutionary dynamics of signaling games in which the alignment of speaker-listener interests can vary. We find that increased misalignment of speaker-listener interests is associated with a decrease in information transmission. Moreover, the most common languages to evolve are typically the most informative languages supportable as static Nash equilibria, suggesting a solution to the 'equilibrium selection problem'. In addition, our dynamics reveal the mechanism by which less informative languages evolve: words that previously signaled intense states come to be used hyperbolically for less intense states, with listeners' interpretation of these newly-ambiguous words evolving downward in response. We ground our results in linguistic data on intensifiers such as so and very, words which have unique dynamics-with constant recycling and innovation that match our theoretical results well.

From insects to mammals, oocytes and sperm develop within germline cysts comprising cells connected by intercellular bridges (ICBs). In numerous insects, formation of the cyst is accompanied by growth of the fusome-a membranous organelle that permeates the cyst. Fusome composition and function are best understood in Drosophila melanogaster: during oogenesis, the fusome dictates cyst topology and size and facilitates oocyte selection, while during spermatogenesis, the fusome synchronizes the cyst's response to DNA damage. Despite its distinct and sex-specific roles during insect gametogenesis, elucidating fusome growth and inheritance in females and its structure and connectivity in males has remained challenging. Here, we take advantage of advances in three-dimensional (3D) confocal microscopy and computational image processing tools to reconstruct the topology, growth, and distribution of the fusome in both sexes. In females, our experimental findings inform a theoretical model for fusome assembly and inheritance and suggest that oocyte selection proceeds through an 'equivalency with a bias' mechanism. In males, we find that cell divisions can deviate from the maximally branched pattern observed in females, leading to greater topological variability. Our work consolidates existing disjointed experimental observations and contributes a readily generalizable computational approach for quantitative studies of gametogenesis within and across species.

Sickle cell disease (SCD) is a common genetic disorder in Africa. Some ongoing work in SCD research includes the analysis and comparisons of variation in phenotypic presentations and disease outcomes with the genotypic signatures. This has contributed to the observed growth of molecular and genetic data in SCD. However, while the "omics" data continues to pile, the capacity to interpret and turn the genetic findings into clinical practice is still underdeveloped, especially in the developing region. Building bioinformatics infrastructure and capacity in the region is key to bridging the gap. This paper seeks to illustrate how the Sickle Cell Programme (SCP) at the Muhimbili University of Health and Allied Sciences (MUHAS) in Tanzania, modeled the integration of infrastructure for bioinformatics and clinical research while running day-to-day clinical care for SCD in Tanzania.

The structure of social networks strongly affects how different phenomena spread in human society, from the transmission of information to the propagation of contagious diseases. It is well-known that heterogeneous connectivity strongly favors spread, but a precise characterization of the redundancy present in social networks and its effect on the robustness of transmission is still lacking. This gap is addressed by the metric backbone, a weight- and connectivity-preserving subgraph that is sufficient to compute all shortest paths of weighted graphs. This subgraph is obtained via algebraically-principled axioms and does not require statistical sampling based on null-models. We show that the metric backbones of nine contact networks obtained from proximity sensors in a variety of social contexts are generally very small, 49% of the original graph for one and ranging from about 6% to 20% for the others. This reflects a surprising amount of redundancy and reveals that shortest paths on these networks are very robust to random attacks and failures. We also show that the metric backbone preserves the full distribution of shortest paths of the original contact networks-which must include the shortest inter- and intra-community distances that define any community structure-and is a primary subgraph for epidemic transmission based on pure diffusion processes. This suggests that the organization of social contact networks is based on large amounts of shortest-path redundancy which shapes epidemic spread in human populations. Thus, the metric backbone is an important subgraph with regard to epidemic spread, the robustness of social networks, and any communication dynamics that depend on complex network shortest paths.

The optomotor response (OMR) is central to the locomotory behavior in diverse animal species including insects, fish and mammals. Furthermore, the study of the OMR in larval zebrafish has become a key model system for investigating the neural basis of sensorimotor control. However, a comprehensive understanding of the underlying control algorithms is still outstanding. In fish it is often assumed that the OMR, by reducing average optic flow across the retina, serves to stabilize position with respect to the ground. Yet the degree to which this is achieved, and how it could emerge from the intermittent burst dynamics of larval zebrafish swimming, are unclear. Here, we combine detailed computational modeling with a new approach to free-swimming experiments in which we control the amount of visual feedback produced by a given motor effort by varying the height of the larva above a moving grid stimulus. We develop an account of underlying feedback control mechanisms that describes both the bout initiation process and the control of swim speed during bouts. We observe that the degree to which fish stabilize their position is only partial and height-dependent, raising questions about its function. We find the relative speed profile during bouts follows a fixed temporal pattern independent of absolute bout speed, suggesting that bout speed and bout termination are not separately controlled. We also find that the reverse optic flow, experienced when the fish is swimming faster than the stimulus, plays a minimal role in control of the OMR despite carrying most of the sensory information about self-movement. These results shed new light on the underlying dynamics of the OMR in larval zebrafish and will be crucial for future work aimed at identifying the neural basis of this behavior.

A key challenge in mobilising growing numbers of digitised biological specimens for scientific research is finding high-throughput methods to extract phenotypic measurements on these datasets. In this paper, we test a pose estimation approach based on Deep Learning capable of accurately placing point labels to identify key locations on specimen images. We then apply the approach to two distinct challenges that each requires identification of key features in a 2D image: (i) identifying body region-specific plumage colouration on avian specimens and (ii) measuring morphometric shape variation in Littorina snail shells. For the avian dataset, 95% of images are correctly labelled and colour measurements derived from these predicted points are highly correlated with human-based measurements. For the Littorina dataset, more than 95% of landmarks were accurately placed relative to expert-labelled landmarks and predicted landmarks reliably captured shape variation between two distinct shell ecotypes ('crab' vs 'wave'). Overall, our study shows that pose estimation based on Deep Learning can generate high-quality and high-throughput point-based measurements for digitised image-based biodiversity datasets and could mark a step change in the mobilisation of such data. We also provide general guidelines for using pose estimation methods on large-scale biological datasets.

The impacts of disease on host vital rates can be demonstrated using longitudinal studies, but these studies can be expensive and logistically challenging. We examined the utility of hidden variable models to infer the individual effects of infectious disease from population-level measurements of survival when longitudinal studies are not possible. Our approach seeks to explain temporal deviations in population-level survival after introducing a disease causative agent when disease prevalence cannot be directly measured by coupling survival and epidemiological models. We tested this approach using an experimental host system (Drosophila melanogaster) with multiple distinct pathogens to validate the ability of the hidden variable model to infer per-capita disease rates. We then applied the approach to a disease outbreak in harbor seals (Phoca vituline) that had data on observed strandings but no epidemiological data. We found that our hidden variable modeling approach could successfully detect the per-capita effects of disease from monitored survival rates in both the experimental and wild populations. Our approach may prove useful for detecting epidemics from public health data in regions where standard surveillance techniques are not available and in the study of epidemics in wildlife populations, where longitudinal studies can be especially difficult to implement.

Large networks of interconnected components, such as genes or machines, can coordinate complex behavioral dynamics. One outstanding question has been to identify the design principles that allow such networks to learn new behaviors. Here, we use Boolean networks as prototypes to demonstrate how periodic activation of network hubs provides a network-level advantage in evolutionary learning. Surprisingly, we find that a network can simultaneously learn distinct target functions upon distinct hub oscillations. We term this emergent property resonant learning, as the new selected dynamical behaviors depend on the choice of the period of the hub oscillations. Furthermore, this procedure accelerates the learning of new behaviors by an order of magnitude faster than without oscillations. While it is well-established that modular network architecture can be selected through evolutionary learning to produce different network behaviors, forced hub oscillations emerge as an alternative evolutionary learning strategy for which network modularity is not necessarily required.

Epigenetic regulatory mechanisms allow multicellular organisms to develop distinct specialized cell identities despite having the same total genome. Cell-fate choices are based on gene expression programs and environmental cues that cells experience during embryonic development, and are usually maintained throughout the life of the organism despite new environmental cues. The evolutionarily conserved Polycomb group (PcG) proteins form Polycomb Repressive Complexes that help orchestrate these developmental choices. Post-development, these complexes actively maintain the resulting cell fate, even in the face of environmental perturbations. Given the crucial role of these polycomb mechanisms in providing phenotypic fidelity (i.e. maintenance of cell fate), we hypothesize that their dysregulation after development will lead to decreased phenotypic fidelity allowing dysregulated cells to sustainably switch their phenotype in response to environmental changes. We call this abnormal phenotypic switching phenotypic pliancy. We introduce a general computational evolutionary model that allows us to test our systems-level phenotypic pliancy hypothesis in-silico and in a context-independent manner. We find that 1) phenotypic fidelity is an emergent systems-level property of PcG-like mechanism evolution, and 2) phenotypic pliancy is an emergent systems-level property resulting from this mechanism's dysregulation. Since there is evidence that metastatic cells behave in a phenotypically pliant manner, we hypothesize that progression to metastasis is driven by the emergence of phenotypic pliancy in cancer cells as a result of PcG mechanism dysregulation. We corroborate our hypothesis using single-cell RNA-sequencing data from metastatic cancers. We find that metastatic cancer cells are phenotypically pliant in the same manner as predicted by our model.

Causal interactions and correlations between clinically-relevant biomarkers are important to understand, both for informing potential medical interventions as well as predicting the likely health trajectory of any individual as they age. These interactions and correlations can be hard to establish in humans, due to the difficulties of routine sampling and controlling for individual differences (e.g., diet, socio-economic status, medication). Because bottlenose dolphins are long-lived mammals that exhibit several age-related phenomena similar to humans, we analyzed data from a well controlled 25-year longitudinal cohort of 144 dolphins. The data from this study has been reported on earlier, and consists of 44 clinically relevant biomarkers. This time-series data exhibits three starkly different influences: (A) directed interactions between biomarkers, (B) sources of biological variation that can either correlate or decorrelate different biomarkers, and (C) random observation-noise which combines measurement error and very rapid fluctuations in the dolphin's biomarkers. Importantly, the sources of biological variation (type-B) are large in magnitude, often comparable to the observation errors (type-C) and larger than the effect of the directed interactions (type-A). Attempting to recover the type-A interactions without accounting for the type-B and type-C variation can result in an abundance of false-positives and false-negatives. Using a generalized regression which fits the longitudinal data with a linear model accounting for all three influences, we demonstrate that the dolphins exhibit many significant directed interactions (type-A), as well as strong correlated variation (type-B), between several pairs of biomarkers. Moreover, many of these interactions are associated with advanced age, suggesting that these interactions can be monitored and/or targeted to predict and potentially affect aging.

In order to treat degenerative diseases, the importance of advanced therapy medicinal products has increased in recent years. The newly developed treatment strategies require a rethinking of the appropriate analytical methods. Current standards are missing the complete and sterile analysis of the product of interest to make the drug manufacturing effort worthwhile. They only consider partial areas of the sample or product while also irreversibly damaging the investigated specimen. Two-dimensional T1 / T2 MR relaxometry meets these requirements and is therefore a promising in-process control during the manufacturing and classification process of cell-based treatments. In this study a tabletop MR scanner was used to perform two-dimensional MR relaxometry. Throughput was increased by developing an automation platform based on a low-cost robotic arm, resulting in the acquisition of a large dataset of cell-based measurements. Two-dimensional inverse Laplace transformation was used for post-processing, followed by data classification performed with support vector machines (SVM) as well as optimized artificial neural networks (ANN). The trained networks were able to distinguish non-differentiated from differentiated MSCs with a prediction accuracy of 85%. To increase versatility, an ANN was trained on 354 independent, biological replicates distributed across ten different cell lines, resulting in a prediction accuracy of up to 98% depending on data composition. The present study provides a proof of principle for the application of T1 / T2 relaxometry as a non-destructive cell classification method. It does not require labeling of cells and can perform whole mount analysis of each sample. Since all measurements can be performed under sterile conditions, it can be used as an in-process control for cellular differentiation. This distinguishes it from other characterization techniques, as most are destructive or require some type of cell labeling. These advantages highlight the technique's potential for preclinical screening of patient-specific cell-based transplants and drugs.

Protein kinases play a vital role in a wide range of cellular processes, and compounds that inhibit kinase activity emerging as a primary focus for targeted therapy development, especially in cancer. Consequently, efforts to characterize the behavior of kinases in response to inhibitor treatment, as well as downstream cellular responses, have been performed at increasingly large scales. Previous work with smaller datasets have used baseline profiling of cell lines and limited kinome profiling data to attempt to predict small molecule effects on cell viability, but these efforts did not use multi-dose kinase profiles and achieved low accuracy with very limited external validation. This work focuses on two large-scale primary data types, kinase inhibitor profiles and gene expression, to predict the results of cell viability screening. We describe the process by which we combined these data sets, examined their properties in relation to cell viability and finally developed a set of computational models that achieve a reasonably high prediction accuracy (R2 of 0.78 and RMSE of 0.154). Using these models, we identified a set of kinases, several of which are understudied, that are strongly influential in the cell viability prediction models. In addition, we also tested to see if a wider range of multiomics data sets could improve the model results and found that proteomic kinase inhibitor profiles were the single most informative data type. Finally, we validated a small subset of the model predictions in several triple-negative and HER2 positive breast cancer cell lines demonstrating that the model performs well with compounds and cell lines that were not included in the training data set. Overall, this result demonstrates that generic knowledge of the kinome is predictive of very specific cell phenotypes, and has the potential to be integrated into targeted therapy development pipelines.

Infectious diseases of plants present an ongoing and increasing threat to international biosecurity, with wide-ranging implications. An important challenge in plant disease management is achieving early detection of invading pathogens, which requires effective surveillance through the implementation of appropriate monitoring programmes. However, when monitoring relies on visual inspection as a means of detection, surveillance is often hindered by a long incubation period (delay from infection to symptom onset) during which plants may be infectious but not displaying visible symptoms. 'Sentinel' plants-alternative susceptible host species that display visible symptoms of infection more rapidly-could be introduced to at-risk populations and included in monitoring programmes to act as early warning beacons for infection. However, while sentinel hosts exhibit faster disease progression and so allow pathogens to be detected earlier, this often comes at a cost: faster disease progression typically promotes earlier onward transmission. Here, we construct a computational model of pathogen transmission to explore this trade-off and investigate how including sentinel plants in monitoring programmes could facilitate earlier detection of invasive plant pathogens. Using Xylella fastidiosa infection in Olea europaea (European olive) as a current high profile case study, for which Catharanthus roseus (Madagascan periwinkle) is a candidate sentinel host, we apply a Bayesian optimisation algorithm to determine the optimal number of sentinel hosts to introduce for a given sampling effort, as well as the optimal division of limited surveillance resources between crop and sentinel plants. Our results demonstrate that including sentinel plants in monitoring programmes can reduce the expected prevalence of infection upon outbreak detection substantially, increasing the feasibility of local outbreak containment.

Design of peptide binders is an attractive strategy for targeting "undruggable" protein-protein interfaces. Current design protocols rely on the extraction of an initial sequence from one known protein interactor of the target protein, followed by in-silico or in-vitro mutagenesis-based optimization of its binding affinity. Wet lab protocols can explore only a minor portion of the vast sequence space and cannot efficiently screen for other desirable properties such as high specificity and low toxicity, while in-silico design requires intensive computational resources and often relies on simplified binding models. Yet, for a multivalent protein target, dozens to hundreds of natural protein partners already exist in the cellular environment. Here, we describe a peptide design protocol that harnesses this diversity via a machine learning generative model. After identifying putative natural binding fragments by literature and homology search, a compositional Restricted Boltzmann Machine is trained and sampled to yield hundreds of diverse candidate peptides. The latter are further filtered via flexible molecular docking and an in-vitro microchip-based binding assay. We validate and test our protocol on calcineurin, a calcium-dependent protein phosphatase involved in various cellular pathways in health and disease. In a single screening round, we identified multiple 16-length peptides with up to six mutations from their closest natural sequence that successfully interfere with the binding of calcineurin to its substrates. In summary, integrating protein interaction and sequence databases, generative modeling, molecular docking and interaction assays enables the discovery of novel protein-protein interaction modulators.

Numerous models have been developed to account for the complex properties of the random walks of biomolecules. However, when analysing experimental data, conditions are rarely met to ensure model identification. The dynamics may simultaneously be influenced by spatial and temporal heterogeneities of the environment, out-of-equilibrium fluxes and conformal changes of the tracked molecules. Recorded trajectories are often too short to reliably discern such multi-scale dynamics, which precludes unambiguous assessment of the type of random walk and its parameters. Furthermore, the motion of biomolecules may not be well described by a single, canonical random walk model. Here, we develop a two-step statistical testing scheme for comparing biomolecule dynamics observed in different experimental conditions without having to identify or make strong prior assumptions about the model generating the recorded random walks. We first train a graph neural network to perform simulation-based inference and thus learn a rich summary statistics vector describing individual trajectories. We then compare trajectories obtained in different biological conditions using a non-parametric maximum mean discrepancy (MMD) statistical test on their so-obtained summary statistics. This procedure allows us to characterise sets of random walks regardless of their generating models, without resorting to model-specific physical quantities or estimators. We first validate the relevance of our approach on numerically simulated trajectories. This demonstrates both the statistical power of the MMD test and the descriptive power of the learnt summary statistics compared to estimates of physical quantities. We then illustrate the ability of our framework to detect changes in α-synuclein dynamics at synapses in cultured cortical neurons, in response to membrane depolarisation, and show that detected differences are largely driven by increased protein mobility in the depolarised state, in agreement with previous findings. The method provides a means of interpreting the differences it detects in terms of single trajectory characteristics. Finally, we emphasise the interest of performing various comparisons to probe the heterogeneity of experimentally acquired datasets at different levels of granularity (e.g., biological replicates, fields of view, and organelles).

Critical phenomena are wildly observed in living systems. If the system is at criticality, it can quickly transfer information and achieve optimal response to external stimuli. Especially, animal collective behavior has numerous critical properties, which are related to other research regions, such as the brain system. Although the critical phenomena influencing collective behavior have been extensively studied, two important aspects require clarification. First, these critical phenomena never occur on a single scale but are instead nested from the micro- to macro-levels (e.g., from a Lévy walk to scale-free correlation). Second, the functional role of group criticality is unclear. To elucidate these aspects, the ambiguous interaction model is constructed in this study; this model has a common framework and is a natural extension of previous representative models (such as the Boids and Vicsek models). We demonstrate that our model can explain the nested criticality of collective behavior across several scales (considering scale-free correlation, super diffusion, Lévy walks, and 1/f fluctuation for relative velocities). Our model can also explain the relationship between scale-free correlation and group turns. To examine this relation, we propose a new method, applying partial information decomposition (PID) to two scale-free induced subgroups. Using PID, we construct information flows between two scale-free induced subgroups and find that coupling of the group morphology (i.e., the velocity distributions) and its fluctuation power (i.e., the fluctuation distributions) likely enable rapid group turning. Thus, the flock morphology may help its internal fluctuation convert to dynamic behavior. Our result sheds new light on the role of group morphology, which is relatively unheeded, retaining the importance of fluctuation dynamics in group criticality.

The basic building blocks of the electrophysiology of cardiomyocytes are ion channels integrated in the cell membranes. Close to the ion channels there are very strong electrical and chemical gradients. However, these gradients extend for only a few nano-meters and are therefore commonly ignored in mathematical models. The full complexity of the dynamics is modelled by the Poisson-Nernst-Planck (PNP) equations but these equations must be solved using temporal and spatial scales of nano-seconds and nano-meters. Here we report solutions of the PNP equations in a fraction of two abuttal cells separated by a tiny extracellular space. We show that when only the potassium channels of the two cells are open, a stationary solution is reached with the well-known Debye layer close to the membranes. When the sodium channels of one of the cells are opened, a very strong and brief electrochemical wave emanates from the channels. If the extracellular space is sufficiently small and the number of sodium channels is sufficiently high, the wave extends all the way over to the neighboring cell and may therefore explain cardiac conduction even at very low levels of gap junctional coupling.

Accurate treatment adjustment to physical activity (PA) remains a challenging problem in type 1 diabetes (T1D) management. Exercise-driven effects on glucose metabolism depend strongly on duration and intensity of the activity, and are highly variable between patients. In-silico evaluation can support the development of improved treatment strategies, and can facilitate personalized treatment optimization. This requires models of the glucose-insulin system that capture relevant exercise-related processes. We developed a model of glucose-insulin regulation that describes changes in glucose metabolism for aerobic moderate- to high-intensity PA of short and prolonged duration. In particular, we incorporated the insulin-independent increase in glucose uptake and production, including glycogen depletion, and the prolonged rise in insulin sensitivity. The model further includes meal absorption and insulin kinetics, allowing simulation of everyday scenarios. The model accurately predicts glucose dynamics for varying PA scenarios in a range of independent validation data sets, and full-day simulations with PA of different timing, duration and intensity agree with clinical observations. We personalized the model on data from a multi-day free-living study of children with T1D by adjusting a small number of model parameters to each child. To assess the use of the personalized models for individual treatment evaluation, we compared subject-specific treatment options for PA management in replay simulations of the recorded data with altered meal, insulin and PA inputs.

Theories of efficient coding propose that the auditory system is optimized for the statistical structure of natural sounds, yet the transformations underlying optimal acoustic representations are not well understood. Using a database of natural sounds including human speech and a physiologically-inspired auditory model, we explore the consequences of peripheral (cochlear) and mid-level (auditory midbrain) filter tuning transformations on the representation of natural sound spectra and modulation statistics. Whereas Fourier-based sound decompositions have constant time-frequency resolution at all frequencies, cochlear and auditory midbrain filters bandwidths increase proportional to the filter center frequency. This form of bandwidth scaling produces a systematic decrease in spectral resolution and increase in temporal resolution with increasing frequency. Here we demonstrate that cochlear bandwidth scaling produces a frequency-dependent gain that counteracts the tendency of natural sound power to decrease with frequency, resulting in a whitened output representation. Similarly, bandwidth scaling in mid-level auditory filters further enhances the representation of natural sounds by producing a whitened modulation power spectrum (MPS) with higher modulation entropy than both the cochlear outputs and the conventional Fourier MPS. These findings suggest that the tuning characteristics of the peripheral and mid-level auditory system together produce a whitened output representation in three dimensions (frequency, temporal and spectral modulation) that reduces redundancies and allows for a more efficient use of neural resources. This hierarchical multi-stage tuning strategy is thus likely optimized to extract available information and may underlies perceptual sensitivity to natural sounds.

Modeling biological mechanisms is a key for disease understanding and drug-target identification. However, formulating quantitative models in the field of Alzheimer's Disease is challenged by a lack of detailed knowledge of relevant biochemical processes. Additionally, fitting differential equation systems usually requires time resolved data and the possibility to perform intervention experiments, which is difficult in neurological disorders. This work addresses these challenges by employing the recently published Variational Autoencoder Modular Bayesian Networks (VAMBN) method, which we here trained on combined clinical and patient level gene expression data while incorporating a disease focused knowledge graph. Our approach, called iVAMBN, resulted in a quantitative model that allowed us to simulate a down-expression of the putative drug target CD33, including potential impact on cognitive impairment and brain pathophysiology. Experimental validation demonstrated a high overlap of molecular mechanism predicted to be altered by CD33 perturbation with cell line data. Altogether, our modeling approach may help to select promising drug targets.

Nowadays, in vitro three-dimensional (3D) neuronal networks are becoming a consolidated experimental model to overcome most of the intrinsic limitations of bi-dimensional (2D) assemblies. In the 3D environment, experimental evidence revealed a wider repertoire of activity patterns, characterized by a modulation of the bursting features, than the one observed in 2D cultures. However, it is not totally clear and understood what pushes the neuronal networks towards different dynamical regimes. One possible explanation could be the underlying connectivity, which could involve a larger number of neurons in a 3D rather than a 2D space and could organize following well-defined topological schemes. Driven by experimental findings, achieved by recording 3D cortical networks organized in multi-layered structures coupled to Micro-Electrode Arrays (MEAs), in the present work we developed a large-scale computational network model made up of leaky integrate-and-fire (LIF) neurons to investigate possible structural configurations able to sustain the emerging patterns of electrophysiological activity. In particular, we investigated the role of the number of layers defining a 3D assembly and the spatial distribution of the connections within and among the layers. These configurations give rise to different patterns of activity that could be compared to the ones emerging from real in vitro 3D neuronal populations. Our results suggest that the introduction of three-dimensionality induced a global reduction in both firing and bursting rates with respect to 2D models. In addition, we found that there is a minimum number of layers necessary to obtain a change in the dynamics of the network. However, the effects produced by a 3D organization of the cells is somewhat mitigated if a scale-free connectivity is implemented in either one or all the layers of the network. Finally, the best matching of the experimental data is achieved supposing a 3D connectivity organized in structured bundles of links located in different areas of the 2D network.

In Italian universities, bioinformatics courses are increasingly being incorporated into different study paths. However, the content of bioinformatics courses is usually selected by the professor teaching the course, in the absence of national guidelines that identify the minimum indispensable knowledge in bioinformatics that undergraduate students from different scientific fields should achieve. The Training&Teaching group of the Bioinformatics Italian Society (BITS) proposed to university professors a survey aimed at portraying the current situation of bioinformatics courses within undergraduate curricula in Italy (i.e., bioinformatics courses activated within both bachelor's and master's degrees). Furthermore, the Training&Teaching group took a cue from the survey outcomes to develop recommendations for the design and the inclusion of bioinformatics courses in academic curricula. Here, we present the outcomes of the survey, as well as the BITS recommendations, with the hope that they may support BITS members in identifying learning outcomes and selecting content for their bioinformatics courses. As we share our effort with the broader international community involved in teaching bioinformatics at academic level, we seek feedback and thoughts on our proposal and hope to start a fruitful debate on the topic, including how to better fulfill the real bioinformatics knowledge needs of the research and the labor market at both the national and international level.

Spiking model neurons can be set up to respond selectively to specific spatio-temporal spike patterns by optimization of their input weights. It is unknown, however, if existing synaptic plasticity mechanisms can achieve this temporal mode of neuronal coding and computation. Here it is shown that changes of synaptic efficacies which tend to balance excitatory and inhibitory synaptic inputs can make neurons sensitive to particular input spike patterns. Simulations demonstrate that a combination of Hebbian mechanisms, hetero-synaptic plasticity and synaptic scaling is sufficient for self-organizing sensitivity for spatio-temporal spike patterns that repeat in the input. In networks inclusion of hetero-synaptic plasticity that depends on the pre-synaptic neurons leads to specialization and faithful representation of pattern sequences by a group of target neurons. Pattern detection is robust against a range of distortions and noise. The proposed combination of Hebbian mechanisms, hetero-synaptic plasticity and synaptic scaling is found to protect the memories for specific patterns from being overwritten by ongoing learning during extended periods when the patterns are not present. This suggests a novel explanation for the long term robustness of memory traces despite ongoing activity with substantial synaptic plasticity. Taken together, our results promote the plausibility of precise temporal coding in the brain.

Neural mass models are used to simulate cortical dynamics and to explain the electrical and magnetic fields measured using electro- and magnetoencephalography. Simulations evince a complex phase-space structure for these kinds of models; including stationary points and limit cycles and the possibility for bifurcations and transitions among different modes of activity. This complexity allows neural mass models to describe the itinerant features of brain dynamics. However, expressive, nonlinear neural mass models are often difficult to fit to empirical data without additional simplifying assumptions: e.g., that the system can be modelled as linear perturbations around a fixed point. In this study we offer a mathematical analysis of neural mass models, specifically the canonical microcircuit model, providing analytical solutions describing slow changes in the type of cortical activity, i.e. dynamical itinerancy. We derive a perturbation analysis up to second order of the phase flow, together with adiabatic approximations. This allows us to describe amplitude modulations in a relatively simple mathematical format providing analytic proof-of-principle for the existence of semi-stable states of cortical dynamics at the scale of a cortical column. This work allows for model inversion of neural mass models, not only around fixed points, but over regions of phase space that encompass transitions among semi or multi-stable states of oscillatory activity. Crucially, these theoretical results speak to model inversion in the context of multiple semi-stable brain states, such as the transition between interictal, pre-ictal and ictal activity in epilepsy.

The synaptic organization of the brain is constantly modified by activity-dependent synaptic plasticity. In several neurological disorders, abnormal neuronal activity and pathological synaptic connectivity may significantly impair normal brain function. Reorganization of neuronal circuits by therapeutic stimulation has the potential to restore normal brain dynamics. Increasing evidence suggests that the temporal stimulation pattern crucially determines the long-lasting therapeutic effects of stimulation. Here, we tested whether a specific pattern of brain stimulation can enable the suppression of pathologically strong inter-population synaptic connectivity through spike-timing-dependent plasticity (STDP). More specifically, we tested how introducing a time shift between stimuli delivered to two interacting populations of neurons can effectively decouple them. To that end, we first used a tractable model, i.e., two bidirectionally coupled leaky integrate-and-fire (LIF) neurons, to theoretically analyze the optimal range of stimulation frequency and time shift for decoupling. We then extended our results to two reciprocally connected neuronal populations (modules) where inter-population delayed connections were modified by STDP. As predicted by the theoretical results, appropriately time-shifted stimulation causes a decoupling of the two-module system through STDP, i.e., by unlearning pathologically strong synaptic interactions between the two populations. Based on the overall topology of the connections, the decoupling of the two modules, in turn, causes a desynchronization of the populations that outlasts the cessation of stimulation. Decoupling effects of the time-shifted stimulation can be realized by time-shifted burst stimulation as well as time-shifted continuous simulation. Our results provide insight into the further optimization of a variety of multichannel stimulation protocols aiming at a therapeutic reshaping of diseased brain networks.