Novel studies conducting cardiac safety assessment using human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are promising but might be limited by their specificity and predictivity. It is often challenging to correctly classify ion channel blockers or to sufficiently predict the risk for Torsade de Pointes (TdP). In this study, we developed a method combining in vitro and in silico experiments to improve machine learning approaches in delivering fast and reliable prediction of drug-induced ion-channel blockade and proarrhythmic behaviour. The algorithm is based on the construction of a dictionary and a greedy optimization, leading to the definition of optimal classifiers. Finally, we present a numerical tool that can accurately predict compound-induced pro-arrhythmic risk and involvement of sodium, calcium and potassium channels, based on hiPSC-CM field potential data.

Many research questions in visual perception involve determining whether stimulus properties are represented and processed independently. In visual neuroscience, there is great interest in determining whether important object dimensions are represented independently in the brain. For example, theories of face recognition have proposed either completely or partially independent processing of identity and emotional expression. Unfortunately, most previous research has only vaguely defined what is meant by "independence," which hinders its precise quantification and testing. This article develops a new quantitative framework that links signal detection theory from psychophysics and encoding models from computational neuroscience, focusing on a special form of independence defined in the psychophysics literature: perceptual separability. The new theory allowed us, for the first time, to precisely define separability of neural representations and to theoretically link behavioral and brain measures of separability. The framework formally specifies the relation between these different levels of perceptual and brain representation, providing the tools for a truly integrative research approach. In particular, the theory identifies exactly what valid inferences can be made about independent encoding of stimulus dimensions from the results of multivariate analyses of neuroimaging data and psychophysical studies. In addition, commonly used operational tests of independence are re-interpreted within this new theoretical framework, providing insights on their correct use and interpretation. Finally, we apply this new framework to the study of separability of brain representations of face identity and emotional expression (neutral/sad) in a human fMRI study with male and female participants.

For a chemical signal to propagate across a cell, it must navigate a tortuous environment involving a variety of organelle barriers. In this work we study mathematical models for a basic chemical signal, the arrival times at the nuclear membrane of proteins that are activated at the cell membrane and diffuse throughout the cytosol. Organelle surfaces within human B cells are reconstructed from soft X-ray tomographic images, and modeled as reflecting barriers to the molecules' diffusion. We show that signal inactivation sharpens signals, reducing variability in the arrival time at the nuclear membrane. Inactivation can also compensate for an observed slowdown in signal propagation induced by the presence of organelle barriers, leading to arrival times at the nuclear membrane that are comparable to models in which the cytosol is treated as an open, empty region. In the limit of strong signal inactivation this is achieved by filtering out molecules that traverse non-geodesic paths.

Resection of the bulk of a tumour often cannot eliminate all cancer cells, due to their infiltration into the surrounding healthy tissue. This may lead to recurrence of the tumour at a later time. We use a reaction-diffusion equation based model of tumour growth to investigate how the invasion front is delayed by resection, and how this depends on the density and behaviour of the remaining cancer cells. We show that the delay time is highly sensitive to qualitative details of the proliferation dynamics of the cancer cell population. The typically assumed logistic type proliferation leads to unrealistic results, predicting immediate recurrence. We find that in glioblastoma cell cultures the cell proliferation rate is an increasing function of the density at small cell densities. Our analysis suggests that cooperative behaviour of cancer cells, analogous to the Allee effect in ecology, can play a critical role in determining the time until tumour recurrence.

Morphogen gradients are crucial for the development of organisms. The biochemical properties of many morphogens prevent their extracellular free diffusion, indicating the need of an active mechanism for transport. The involvement of filopodial structures (cytonemes) has been proposed for morphogen signaling. Here, we describe an in silico model based on the main general features of cytoneme-meditated gradient formation and its implementation into Cytomorph, an open software tool. We have tested the spatial and temporal adaptability of our model quantifying Hedgehog (Hh) gradient formation in two Drosophila tissues. Cytomorph is able to reproduce the gradient and explain the different scaling between the two epithelia. After experimental validation, we studied the predicted impact of a range of features such as length, size, density, dynamics and contact behavior of cytonemes on Hh morphogen distribution. Our results illustrate Cytomorph as an adaptive tool to test different morphogen gradients and to generate hypotheses that are difficult to study experimentally.

Epithelial tissues form folded structures during embryonic development and organogenesis. Whereas substantial efforts have been devoted to identifying mechanical and biochemical mechanisms that induce folding, whether and how their interplay synergistically shapes epithelial folds remains poorly understood. Here we propose a mechano-biochemical model for dorsal fold formation in the early Drosophila embryo, an epithelial folding event induced by shifts of cell polarity. Based on experimentally observed apical domain homeostasis, we couple cell mechanics to polarity and find that mechanical changes following the initial polarity shifts alter cell geometry, which in turn influences the reaction-diffusion of polarity proteins, thus forming a feedback loop between cell mechanics and polarity. This model can induce spontaneous fold formation in silico, recapitulate polarity and shape changes observed in vivo, and confer robustness to tissue shape change against small fluctuations in mechanics and polarity. These findings reveal emergent properties of a developing epithelium under control of intracellular mechano-polarity coupling.

Genome-wide somatic mutation profiles of tumours can now be assessed efficiently and promise to move precision medicine forward. Statistical analysis of mutation profiles is however challenging due to the low frequency of most mutations, the varying mutation rates across tumours, and the presence of a majority of passenger events that hide the contribution of driver events. Here we propose a method, NetNorM, to represent whole-exome somatic mutation data in a form that enhances cancer-relevant information using a gene network as background knowledge. We evaluate its relevance for two tasks: survival prediction and unsupervised patient stratification. Using data from 8 cancer types from The Cancer Genome Atlas (TCGA), we show that it improves over the raw binary mutation data and network diffusion for these two tasks. In doing so, we also provide a thorough assessment of somatic mutations prognostic power which has been overlooked by previous studies because of the sparse and binary nature of mutations.

Adaptive cognitive-control involves a hierarchical cortico-striatal gating system that supports selective updating, maintenance, and retrieval of useful cognitive and motor information. Here, we developed a task that independently manipulates selective gating operations into working-memory (input gating), from working-memory (output gating), and of responses (motor gating) and tested the neural dynamics and computational principles that support them. Increases in gating demands, captured by gate switches, were expressed by distinct EEG correlates at each gating level that evolved dynamically in partially overlapping time windows. Further, categorical representations of specific maintained items and of motor responses could be decoded from EEG when the corresponding gate was switching, thereby linking gating operations to prioritization. Finally, gate switching at all levels was related to increases in the motor decision threshold as quantified by the drift diffusion model. Together these results support the notion that cognitive gating operations scaffold on top of mechanisms involved in motor gating.

In Escherichia coli, translocation of RNA polymerase (RNAP) during transcription introduces supercoiling to DNA, which influences the initiation and elongation behaviors of RNAP. To quantify the role of supercoiling in transcription regulation, we developed a spatially resolved supercoiling model of transcription. The integrated model describes how RNAP activity feeds back with the local DNA supercoiling and how this mechanochemical feedback controls transcription, subject to topoisomerase activities and stochastic topological domain formation. This model establishes that transcription-induced supercoiling mediates the cooperation of co-transcribing RNAP molecules in highly expressed genes, and this cooperation is achieved under moderate supercoiling diffusion and high topoisomerase unbinding rates. It predicts that a topological domain could serve as a transcription regulator, generating substantial transcriptional noise. It also shows the relative orientation of two closely arranged genes plays an important role in regulating their transcription. The model provides a quantitative platform for investigating how genome organization impacts transcription.

Formation of functionally adequate vascular networks by angiogenesis presents a problem in biological patterning. Generated without predetermined spatial patterns, networks must develop hierarchical tree-like structures for efficient convective transport over large distances, combined with dense space-filling meshes for short diffusion distances to every point in the tissue. Moreover, networks must be capable of restructuring in response to changing functional demands without interruption of blood flow. Here, theoretical simulations based on experimental data are used to demonstrate that this patterning problem can be solved through over-abundant stochastic generation of vessels in response to a growth factor generated in hypoxic tissue regions, in parallel with refinement by structural adaptation and pruning. Essential biological mechanisms for generation of adequate and efficient vascular patterns are identified and impairments in vascular properties resulting from defects in these mechanisms are predicted. The results provide a framework for understanding vascular network formation in normal or pathological conditions and for predicting effects of therapies targeting angiogenesis.

We often need to rapidly change our mind about perceptual decisions in order to account for new information and correct mistakes. One fundamental, unresolved question is whether information processed prior to a decision being made ('pre-decisional information') has any influence on the likelihood and speed with which that decision is reversed. We investigated this using a luminance discrimination task in which participants indicated which of two flickering greyscale squares was brightest. Following an initial decision, the stimuli briefly remained on screen, and participants could change their response. Using psychophysical reverse correlation, we examined how moment-to-moment fluctuations in stimulus luminance affected participants' decisions. This revealed that the strength of even the very earliest (pre-decisional) evidence was associated with the likelihood and speed of later changes of mind. To account for this effect, we propose an extended diffusion model in which an initial 'snapshot' of sensory information biases ongoing evidence accumulation.

Eukaryotic cells partition a wide variety of important materials and processes into biomolecular condensates-phase-separated droplets that lack a membrane. In addition to nonspecific electrostatic or hydrophobic interactions, phase separation also depends on specific binding motifs that link together constituent molecules. Nevertheless, few rules have been established for how these ubiquitous specific, saturating, motif-motif interactions drive phase separation. By integrating Monte Carlo simulations of lattice-polymers with mean-field theory, we show that the sequence of heterotypic binding motifs strongly affects a polymer's ability to phase separate, influencing both phase boundaries and condensate properties (e.g. viscosity and polymer diffusion). We find that sequences with large blocks of single motifs typically form more inter-polymer bonds, which promotes phase separation. Notably, the sequence of binding motifs influences phase separation primarily by determining the conformational entropy of self-bonding by single polymers. This contrasts with systems where the molecular architecture primarily affects the energy of the dense phase, providing a new entropy-based mechanism for the biological control of phase separation.

The bacterial flagellar motor is a highly efficient rotary machine used by many bacteria to propel themselves. It has recently been shown that at low speeds its rotation proceeds in steps. Here we propose a simple physical model, based on the storage of energy in protein springs, that accounts for this stepping behavior as a random walk in a tilted corrugated potential that combines torque and contact forces. We argue that the absolute angular position of the rotor is crucial for understanding step properties and show this hypothesis to be consistent with the available data, in particular the observation that backward steps are smaller on average than forward steps. We also predict a sublinear speed versus torque relationship for fixed load at low torque, and a peak in rotor diffusion as a function of torque. Our model provides a comprehensive framework for understanding and analyzing stepping behavior in the bacterial flagellar motor and proposes novel, testable predictions. More broadly, the storage of energy in protein springs by the flagellar motor may provide useful general insights into the design of highly efficient molecular machines.

Perceptual decisions are thought to be mediated by a mechanism of sequential sampling and integration of noisy evidence whose temporal weighting profile affects the decision quality. To examine temporal weighting, participants were presented with two brightness-fluctuating disks for 1, 2 or 3 seconds and were requested to choose the overall brighter disk at the end of each trial. By employing a signal-perturbation method, which deploys across trials a set of systematically controlled temporal dispersions of the same overall signal, we were able to quantify the participants' temporal weighting profile. Results indicate that, for intervals of 1 or 2 sec, participants exhibit a primacy-bias. However, for longer stimuli (3-sec) the temporal weighting profile is non-monotonic, with concurrent primacy and recency, which is inconsistent with the predictions of previously suggested computational models of perceptual decision-making (drift-diffusion and Ornstein-Uhlenbeck processes). We propose a novel, dynamic variant of the leaky-competing accumulator model as a potential account for this finding, and we discuss potential neural mechanisms.

Many large-scale functional connectivity studies have emphasized the importance of communication through increased inter-region correlations during task states. In contrast, local circuit studies have demonstrated that task states primarily reduce correlations among pairs of neurons, likely enhancing their information coding by suppressing shared spontaneous activity. Here we sought to adjudicate between these conflicting perspectives, assessing whether co-active brain regions during task states tend to increase or decrease their correlations. We found that variability and correlations primarily decrease across a variety of cortical regions in two highly distinct data sets: non-human primate spiking data and human functional magnetic resonance imaging data. Moreover, this observed variability and correlation reduction was accompanied by an overall increase in dimensionality (reflecting less information redundancy) during task states, suggesting that decreased correlations increased information coding capacity. We further found in both spiking and neural mass computational models that task-evoked activity increased the stability around a stable attractor, globally quenching neural variability and correlations. Together, our results provide an integrative mechanistic account that encompasses measures of large-scale neural activity, variability, and correlations during resting and task states.

Rhythmic voltage oscillations resulting from the summed activity of neuronal populations occur in many nervous systems. Contemporary observations suggest that coexistent oscillations interact and, in time, may switch in dominance. We recently reported an example of these interactions recorded from in vitro preparations of rat somatosensory cortex. We found that following an initial interval of coexistent gamma ( approximately 25 ms period) and beta2 ( approximately 40 ms period) rhythms in the superficial and deep cortical layers, respectively, a transition to a synchronous beta1 ( approximately 65 ms period) rhythm in all cortical layers occurred. We proposed that the switch to beta1 activity resulted from the novel mechanism of period concatenation of the faster rhythms: gamma period (25 ms)+beta2 period (40 ms) = beta1 period (65 ms). In this article, we investigate in greater detail the fundamental mechanisms of the beta1 rhythm. To do so we describe additional in vitro experiments that constrain a biologically realistic, yet simplified, computational model of the activity. We use the model to suggest that the dynamic building blocks (or motifs) of the gamma and beta2 rhythms combine to produce a beta1 oscillation that exhibits cross-frequency interactions. Through the combined approach of in vitro experiments and mathematical modeling we isolate the specific components that promote or destroy each rhythm. We propose that mechanisms vital to establishing the beta1 oscillation include strengthened connections between a population of deep layer intrinsically bursting cells and a transition from antidromic to orthodromic spike generation in these cells. We conclude that neural activity in the superficial and deep cortical layers may temporally combine to generate a slower oscillation.

Brain networks store new memories using functional and structural synaptic plasticity. Memory formation is generally attributed to Hebbian plasticity, while homeostatic plasticity is thought to have an ancillary role in stabilizing network dynamics. Here we report that homeostatic plasticity alone can also lead to the formation of stable memories. We analyze this phenomenon using a new theory of network remodeling, combined with numerical simulations of recurrent spiking neural networks that exhibit structural plasticity based on firing rate homeostasis. These networks are able to store repeatedly presented patterns and recall them upon the presentation of incomplete cues. Storage is fast, governed by the homeostatic drift. In contrast, forgetting is slow, driven by a diffusion process. Joint stimulation of neurons induces the growth of associative connections between them, leading to the formation of memory engrams. These memories are stored in a distributed fashion throughout connectivity matrix, and individual synaptic connections have only a small influence. Although memory-specific connections are increased in number, the total number of inputs and outputs of neurons undergo only small changes during stimulation. We find that homeostatic structural plasticity induces a specific type of "silent memories", different from conventional attractor states.

As nascent proteins are synthesized by the ribosome, they depart via an exit tunnel running through the center of the large subunit. The exit tunnel likely plays an important part in various aspects of translation. Although water plays a key role in many bio-molecular processes, the nature of water confined to the exit tunnel has remained unknown. Furthermore, solvent in biological cavities has traditionally been characterized as either a continuous dielectric fluid, or a discrete tightly bound molecule. Using atomistic molecular dynamics simulations, we predict that the thermodynamic and kinetic properties of water confined within the ribosome exit tunnel are quite different from this simple two-state model. We find that the tunnel creates a complex microenvironment for the solvent resulting in perturbed rotational dynamics and heterogenous dielectric behavior. This gives rise to a very rugged solvation landscape and significantly retarded solvent diffusion. We discuss how this non-bulk-like solvent is likely to affect important biophysical processes such as sequence dependent stalling, co-translational folding, and antibiotic binding. We conclude with a discussion of the general applicability of these results to other biological cavities.

Neurodegenerative diseases such as Alzheimer's or Parkinson's are associated with the prion-like propagation and aggregation of toxic proteins. A long standing hypothesis that amyloid-beta drives Alzheimer's disease has proven the subject of contemporary controversy; leading to new research in both the role of tau protein and its interaction with amyloid-beta. Conversely, recent work in mathematical modeling has demonstrated the relevance of nonlinear reaction-diffusion type equations to capture essential features of the disease. Such approaches have been further simplified, to network-based models, and offer researchers a powerful set of computationally tractable tools with which to investigate neurodegenerative disease dynamics. Here, we propose a novel, coupled network-based model for a two-protein system that includes an enzymatic interaction term alongside a simple model of aggregate transneuronal damage. We apply this theoretical model to test the possible interactions between tau proteins and amyloid-beta and study the resulting coupled behavior between toxic protein clearance and proteopathic phenomenology. Our analysis reveals ways in which amyloid-beta and tau proteins may conspire with each other to enhance the nucleation and propagation of different diseases, thus shedding new light on the importance of protein clearance and protein interaction mechanisms in prion-like models of neurodegenerative disease.

How splicing factors are recruited to nascent transcripts in the nucleus in order to assemble spliceosomes on newly synthesised pre-mRNAs is unknown. To address this question, we compared the intranuclear trafficking kinetics of small nuclear ribonucleoprotein particles (snRNP) and non-snRNP proteins in the presence and absence of splicing activity. Photobleaching experiments clearly show that spliceosomal proteins move continuously throughout the entire nucleus independently of ongoing transcription or splicing. Using quantitative experimental data, a mathematical model was applied for spliceosome assembly and recycling in the nucleus. The model assumes that splicing proteins move by Brownian diffusion and interact stochastically with binding sites located at different subnuclear compartments. Inhibition of splicing, which reduces the number of pre-mRNA binding sites available for spliceosome assembly, was modeled as a decrease in the on-rate binding constant in the nucleoplasm. Simulation of microscopy experiments before and after splicing inhibition yielded results consistent with the experimental observations. Taken together, our data argue against the view that spliceosomal components are stored in nuclear speckles until a signal triggers their recruitment to nascent transcripts. Rather, the results suggest that splicing proteins are constantly diffusing throughout the entire nucleus and collide randomly and transiently with pre-mRNAs.