Noise in biology

Probing the design principles of photosynthetic systems through fluorescence noise measurement

Elucidating the energetic processes which govern photosynthesis, the engine of life on earth, are an essential goal both for fundamental research and for cutting-edge biotechnological applications. Fluorescent signal of photosynthetic markers has long been utilised in this endeavour. In this research we demonstrate the use of fluorescent noise analysis to reveal further layers of intricacy in photosynthetic energy transfer. While noise is a common tool analysing dynamics in physics and engineering, its application in biology has thus far been limited. Here, a distinct behaviour in photosynthetic pigments across various chemical and biological environments is measured. These changes seem to elucidate quantum effects governing the generation of oxidative radicals. Although our method offers insights, it is important to note that the interpretation should be further validated expertly to support as conclusive theory. This innovative method is simple, non-invasive, and immediate, making it a promising tool to uncover further, more complex energetic events in photosynthesis, with potential uses in environmental monitoring, agriculture, and food-tech.

Responses of organ precursors to correct and incorrect inductive signals

During embryonic development, the inductive molecules produced by local origins normally arrive at their target tissues in a nondirectional, diffusion manner. The target organ precursor cells must correctly interpret these inductive signals to ensure proper specification/differentiation, which is dependent on two prerequisites: (i) obtaining cell-intrinsic competence; and (ii) receiving correct inductive signals while resisting incorrect ones. Gain of intrinsic competence could avoid a large number of misinductions because the incompetent cells are nonresponsive to inductive signals. However, in cases of different precursor cells with similar competence and located in close proximity, resistance to incorrect inductive signals is essential for accurate determination of cell fate. Here we outline the mechanisms of how organ precursors respond to correct and incorrect inductive signals.

FEP Augmentation as a Means to Solve Data Paucity Problems for Machine Learning in Chemical Biology

In the realm of medicinal chemistry, the primary objective is to swiftly optimize a multitude of chemical properties of a set of compounds to yield a clinical candidate poised for clinical trials. In recent years, two computational techniques, machine learning (ML) and physics-based methods, have evolved substantially and are now frequently incorporated into the medicinal chemist's toolbox to enhance the efficiency of both hit optimization and candidate design. Both computational methods come with their own set of limitations, and they are often used independently of each other. ML's capability to screen extensive compound libraries expediently is tempered by its reliance on quality data, which can be scarce especially during early-stage optimization. Contrarily, physics-based approaches like free energy perturbation (FEP) are frequently constrained by low throughput and high cost by comparison; however, physics-based methods are capable of making highly accurate binding affinity predictions. In this study, we harnessed the strength of FEP to overcome data paucity in ML by generating virtual activity data sets which then inform the training of algorithms. Here, we show that ML algorithms trained with an FEP-augmented data set could achieve comparable predictive accuracy to data sets trained on experimental data from biological assays. Throughout the paper, we emphasize key mechanistic considerations that must be taken into account when aiming to augment data sets and lay the groundwork for successful implementation. Ultimately, the study advocates for the synergy of physics-based methods and ML to expedite the lead optimization process. We believe that the physics-based augmentation of ML will significantly benefit drug discovery, as these techniques continue to evolve.

Biology in the 21st century: Natural selection is cognitive selection

Natural selection has a formal definition as the natural process that results in the survival and reproductive success of individuals or groups best adjusted to their environment, leading to the perpetuation of those genetic qualities best suited to that organism's environmental niche. Within conventional Neo-Darwinism, the largest source of those variations that can be selected is presumed to be secondary to random genetic mutations. As these arise, natural selection sustains adaptive traits in the context of a 'struggle for existence'. Consequently, in the 20th century, natural selection was generally portrayed as the primary evolutionary driver. The 21st century offers a comprehensive alternative to Neo-Darwinian dogma within Cognition-Based Evolution. The substantial differences between these respective evolutionary frameworks have been most recently articulated in a revision of Crick's Central Dogma, a former centerpiece of Neo-Darwinism. The argument is now advanced that the concept of natural selection should also be comprehensively reappraised. Cognitive selection is presented as a more precise term better suited to 21st century biology. Since cognition began with life's origin, natural selection represents cognitive selection.

Synergistic and antagonistic effects of deterministic and stochastic cell-cell variations

By diversifying, cells in a clonal population can together overcome the limits of individuals. Diversity in single-cell growth rates allows the population to survive environmental stresses, such as antibiotics, and grow faster than the undiversified population. These functional cell-cell variations can arise stochastically, from noise in biochemical reactions, or deterministically, by asymmetrically distributing damaged components. While each of the mechanisms is well understood, the effect of the combined mechanisms is unclear. To evaluate the contribution of the deterministic component we developed a mathematical model by mapping the growing population to the Ising model. To analyze the combined effects of stochastic and deterministic contributions we introduced the analytical results of the Ising-mapping into an Euler-Lotka framework. Model results, confirmed by simulations and experimental data, show that deterministic cell-cell variations increase near-linearly with stress. As a consequence, we predict that the gain in population doubling time from cell-cell variations is primarily stochastic at low stress but may cross over to deterministic at higher stresses. Furthermore, we find that while the deterministic component minimizes population damage, stochastic variations antagonize this effect. Together our results may help identifying stress-tolerant pathogenic cells and thus inspire novel antibiotic strategies.

Regulation of Chromatin Modifications through Coordination of Nucleus Size and Epithelial Cell Morphology Heterogeneity

Cell morphology heterogeneity within epithelial collectives is a pervasive phenomenon intertwined with tissue mechanical properties. Despite its widespread occurrence, the underlying mechanisms driving cell morphology heterogeneity and its consequential biological ramifications remain elusive. Here, we investigate the dynamic evolution of epithelial cell morphology and nucleus morphology during crowding, unveiling a consistent correlation between the two. Our investigation reveals a persistent log-normal probability distribution characterizing both cell and nucleus areas across diverse crowding stages and epithelial model systems. We showed that this morphological diversity arises from asymmetric partitioning during cell division and is perpetuated through actomyosin-mediated regulation of cell-nucleus size coordination. Moreover, we provide insights into the impact of nucleus morphology on chromatin dynamics, demonstrating that constraining nucleus area leads to downregulation of the euchromatic mark H3K9ac and upregulation of the heterochromatic mark H3K27me3 through modulation of histone demethylase UTX expression. These findings under-score the significance of cell morphology heterogeneity as a driver of chromatin state diversity, shaping functional variability within epithelial tissues.

Variability in Neural Circuit Formation

The study of neural development is usually concerned with the question of how nervous systems get put together. Variation in these processes is usually of interest as a means of revealing these normative mechanisms. However, variation itself can be an object of study and is of interest from multiple angles. First, the nature of variation in both the processes and the outcomes of neural development is relevant to our understanding of how these processes and outcomes are encoded in the genome. Second, variation in the wiring of the brain in humans may underlie variation in all kinds of psychological and behavioral traits, as well as neurodevelopmental disorders. And third, genetic variation that affects circuit development provides the raw material for evolutionary change. Here, I examine these different aspects of variation in circuit development and consider what they may tell us about these larger questions.

Design patterns of biological cells

Design patterns are generalized solutions to frequently recurring problems. They were initially developed by architects and computer scientists to create a higher level of abstraction for their designs. Here, we extend these concepts to cell biology to lend a new perspective on the evolved designs of cells' underlying reaction networks. We present a catalog of 21 design patterns divided into three categories: creational patterns describe processes that build the cell, structural patterns describe the layouts of reaction networks, and behavioral patterns describe reaction network function. Applying this pattern language to the E. coli central metabolic reaction network, the yeast pheromone response signaling network, and other examples lends new insights into these systems.

Life Is Uncertain: Inherent Variability Exhibited by Organisms, and at Higher Levels of Biological Organization

Organisms act stochastically. A not uncommon view in the ecological literature is that this is mainly due to the observer having insufficient information or a stochastic environment-and not partly because organisms themselves respond with inherent unpredictability. In this study, I compile the evidence that contradicts that view. Organisms generate uncertainty internally, which results in irreducible stochastic responses. I consider why: for instance, stochastic responses are associated with greater adaptability to changing environments and resource availability. Over longer timescales, biologically generated uncertainty influences behavior, evolution, and macroecological processes. Indeed, it could be stated that organisms are systems defined by the internal generation, magnification, and record-keeping of uncertainty as inputs to responses. Important practical implications arise if organisms can indeed be defined by an association with specific classes of inherent uncertainty: not least that isolating those signatures then provides a potential means for detecting life, for considering the forms that life could theoretically take, and for exploring the wider limits to how life might become distributed. These are all fundamental goals in astrobiology.

Accounting for network noise in graph-guided Bayesian modeling of structured high-dimensional data

There is a growing body of literature on knowledge-guided statistical learning methods for analysis of structured high-dimensional data (such as genomic and transcriptomic data) that can incorporate knowledge of underlying networks derived from functional genomics and functional proteomics. These methods have been shown to improve variable selection and prediction accuracy and yield more interpretable results. However, these methods typically use graphs extracted from existing databases or rely on subject matter expertise, which are known to be incomplete and may contain false edges. To address this gap, we propose a graph-guided Bayesian modeling framework to account for network noise in regression models involving structured high-dimensional predictors. Specifically, we use 2 sources of network information, including the noisy graph extracted from existing databases and the estimated graph from observed predictors in the dataset at hand, to inform the model for the true underlying network via a latent scale modeling framework. This model is coupled with the Bayesian regression model with structured high-dimensional predictors involving an adaptive structured shrinkage prior. We develop an efficient Markov chain Monte Carlo algorithm for posterior sampling. We demonstrate the advantages of our method over existing methods in simulations, and through analyses of a genomics dataset and another proteomics dataset for Alzheimer's disease.

Harnessing synthetic active particles for physical reservoir computing

The processing of information is an indispensable property of living systems realized by networks of active processes with enormous complexity. They have inspired many variants of modern machine learning, one of them being reservoir computing, in which stimulating a network of nodes with fading memory enables computations and complex predictions. Reservoirs are implemented on computer hardware, but also on unconventional physical substrates such as mechanical oscillators, spins, or bacteria often summarized as physical reservoir computing. Here we demonstrate physical reservoir computing with a synthetic active microparticle system that self-organizes from an active and passive component into inherently noisy nonlinear dynamical units. The self-organization and dynamical response of the unit are the results of a delayed propulsion of the microswimmer to a passive target. A reservoir of such units with a self-coupling via the delayed response can perform predictive tasks despite the strong noise resulting from the Brownian motion of the microswimmers. To achieve efficient noise suppression, we introduce a special architecture that uses historical reservoir states for output. Our results pave the way for the study of information processing in synthetic self-organized active particle systems.

Flagellar dynamics reveal fluctuations and kinetic limit in the Escherichia coli chemotaxis network

The Escherichia coli chemotaxis network, by which bacteria modulate their random run/tumble swimming pattern to navigate their environment, must cope with unavoidable number fluctuations ("noise") in its molecular constituents like other signaling networks. The probability of clockwise (CW) flagellar rotation, or CW bias, is a measure of the chemotaxis network's output, and its temporal fluctuations provide a proxy for network noise. Here we quantify fluctuations in the chemotaxis signaling network from the switching statistics of flagella, observed using time-resolved fluorescence microscopy of individual optically trapped E. coli cells. This approach allows noise to be quantified across the dynamic range of the network. Large CW bias fluctuations are revealed at steady state, which may play a critical role in driving flagellar switching and cell tumbling. When the network is stimulated chemically to higher activity, fluctuations dramatically decrease. A stochastic theoretical model, inspired by work on gene expression noise, points to CheY activation occurring in bursts, driving CW bias fluctuations. This model also shows that an intrinsic kinetic ceiling on network activity places an upper limit on activated CheY and CW bias, which when encountered suppresses network fluctuations. This limit may also prevent cells from tumbling unproductively in steep gradients.

Deep learning in bioinformatics

Deep learning is a powerful machine learning technique that can learn from large amounts of data using multiple layers of artificial neural networks. This paper reviews some applications of deep learning in bioinformatics, a field that deals with analyzing and interpreting biological data. We first introduce the basic concepts of deep learning and then survey the recent advances and challenges of applying deep learning to various bioinformatics problems, such as genome sequencing, gene expression analysis, protein structure prediction, drug discovery, and disease diagnosis. We also discuss future directions and opportunities for deep learning in bioinformatics. We aim to provide an overview of deep learning so that bioinformaticians applying deep learning models can consider all critical technical and ethical aspects. Thus, our target audience is biomedical informatics researchers who use deep learning models for inference. This review will inspire more bioinformatics researchers to adopt deep-learning methods for their research questions while considering fairness, potential biases, explainability, and accountability.

Embracing Heterogeneity: Challenging the Paradigm of Replisomes as Deterministic Machines

The paradigm of cellular systems as deterministic machines has long guided our understanding of biology. Advancements in technology and methodology, however, have revealed a world of stochasticity, challenging the notion of determinism. Here, we explore the stochastic behavior of multi-protein complexes, using the DNA replication system (replisome) as a prime example. The faithful and timely copying of DNA depends on the simultaneous action of a large set of enzymes and scaffolding factors. This fundamental cellular process is underpinned by dynamic protein-nucleic acid assemblies that must transition between distinct conformations and compositional states. Traditionally viewed as a well-orchestrated molecular machine, recent experimental evidence has unveiled significant variability and heterogeneity in the replication process. In this review, we discuss recent advances in single-molecule approaches and single-particle cryo-EM, which have provided insights into the dynamic processes of DNA replication. We comment on the new challenges faced by structural biologists and biophysicists as they attempt to describe the dynamic cascade of events leading to replisome assembly, activation, and progression. The fundamental principles uncovered and yet to be discovered through the study of DNA replication will inform on similar operating principles for other multi-protein complexes.

Actomyosin pulsing rescues embryonic tissue folding from disruption by myosin fluctuations

During early development, myosin II mechanically reshapes and folds embryo tissue. A muchstudied example is ventral furrow formation in Drosophila, marking the onset of gastrulation. Furrowing is driven by contraction of actomyosin networks on apical cell surfaces, but how the myosin patterning encodes tissue shape is unclear, and elastic models failed to reproduce essential features of experimental cell contraction profiles. The myosin patterning exhibits substantial cell-to-cell fluctuations with pulsatile time-dependence, a striking but unexplained feature of morphogenesis in many organisms. Here, using biophysical modeling we find viscous forces offer the principal resistance to actomyosin-driven apical constriction. In consequence, tissue shape is encoded in the direction-dependent curvature of the myosin patterning which orients an anterior-posterior furrow. Tissue contraction is highly sensitive to cell-to-cell myosin fluctuations, explaining furrowing failure in genetically perturbed embryos whose fluctuations are temporally persistent. In wild-type embryos this disastrous outcome is averted by pulsatile myosin time-dependence, which rescues furrowing by eliminating high frequencies in the fluctuation power spectrum. This low pass filter mechanism may underlie the usage of actomyosin pulsing in diverse morphogenetic processes across many organisms.

A bioelectrical phase transition patterns the first vertebrate heartbeats

A regular heartbeat is essential to vertebrate life. In the mature heart, this function is driven by an anatomically localized pacemaker. By contrast, pacemaking capability is broadly distributed in the early embryonic heart1-3, raising the question of how tissue-scale activity is first established and then maintained during embryonic development. The initial transition of the heart from silent to beating has never been characterized at the timescale of individual electrical events, and the structure in space and time of the early heartbeats remains poorly understood. Using all-optical electrophysiology, we captured the very first heartbeat of a zebrafish and analysed the development of cardiac excitability and conduction around this singular event. The first few beats appeared suddenly, had irregular interbeat intervals, propagated coherently across the primordial heart and emanated from loci that varied between animals and over time. The bioelectrical dynamics were well described by a noisy saddle-node on invariant circle bifurcation with action potential upstroke driven by CaV1.2. Our work shows how gradual and largely asynchronous development of single-cell bioelectrical properties produces a stereotyped and robust tissue-scale transition from quiescence to coordinated beating.

Mechanistic modelling of within-mosquito viral dynamics: Insights into infection and dissemination patterns

Vector or host competence can be defined as the ability of an individual to become infected and subsequently transmit a pathogen. Assays to measure competence play a key part in the assessment of the factors affecting mosquito-borne virus transmission and of potential pathogen-blocking control tools for these viruses. For mosquitoes, competence for arboviruses can be measured experimentally and results are usually analysed using standard statistical approaches. Here we develop a mechanistic approach to studying within-mosquito virus dynamics that occur during vector competence experiments. We begin by developing a deterministic model of virus replication in the mosquito midgut and subsequent escape and replication in the hemocoel. We then extend this to a stochastic model to capture the between-individual variation observed in vector competence experiments. We show that the dose-response of the probability of mosquito midgut infection and variation in the dissemination rate can be explained by stochastic processes generated from a small founding population of virions, caused by a relatively low rate of virion infection of susceptible cells. We also show that comparing treatments or species in competence experiments by fitting mechanistic models could provide further insight into potential differences. Generally, our work adds to the growing body of literature emphasizing the importance of intrinsic stochasticity in biological systems.

Theoretical and computational tools to model multistable gene regulatory networks

The last decade has witnessed a surge of theoretical and computational models to describe the dynamics of complex gene regulatory networks, and how these interactions can give rise to multistable and heterogeneous cell populations. As the use of theoretical modeling to describe genetic and biochemical circuits becomes more widespread, theoreticians with mathematical and physical backgrounds routinely apply concepts from statistical physics, non-linear dynamics, and network theory to biological systems. This review aims at providing a clear overview of the most important methodologies applied in the field while highlighting current and future challenges. It also includes hands-on tutorials to solve and simulate some of the archetypical biological system models used in the field. Furthermore, we provide concrete examples from the existing literature for theoreticians that wish to explore this fast-developing field. Whenever possible, we highlight the similarities and differences between biochemical and regulatory networks and 'classical' systems typically studied in non-equilibrium statistical and quantum mechanics.

Using Constrained-Disorder Principle-Based Systems to Improve the Performance of Digital Twins in Biological Systems

Digital twins are computer programs that use real-world data to create simulations that predict the performance of processes, products, and systems. Digital twins may integrate artificial intelligence to improve their outputs. Models for dealing with uncertainties and noise are used to improve the accuracy of digital twins. Most currently used systems aim to reduce noise to improve their outputs. Nevertheless, biological systems are characterized by inherent variability, which is necessary for their proper function. The constrained-disorder principle defines living systems as having a disorder as part of their existence and proper operation while kept within dynamic boundaries. In the present paper, we review the role of noise in complex systems and its use in bioengineering. We describe the use of digital twins for medical applications and current methods for dealing with noise and uncertainties in modeling. The paper presents methods to improve the accuracy and effectiveness of digital twin systems by continuously implementing variability signatures while simultaneously reducing unwanted noise in their inputs and outputs. Accounting for the noisy internal and external environments of complex biological systems is necessary for the future design of improved, more accurate digital twins.

The σB alternative sigma factor circuit modulates noise to generate different types of pulsing dynamics

Single-cell approaches are revealing a high degree of heterogeneity, or noise, in gene expression in isogenic bacteria. How gene circuits modulate this noise in gene expression to generate robust output dynamics is unclear. Here we use the Bacillus subtilis alternative sigma factor σB as a model system for understanding the role of noise in generating circuit output dynamics. σB controls the general stress response in B. subtilis and is activated by a range of energy and environmental stresses. Recent single-cell studies have revealed that the circuit can generate two distinct outputs, stochastic pulsing and a single pulse response, but the conditions under which each response is generated are under debate. We implement a stochastic mathematical model of the σB circuit to investigate this and find that the system's core circuit can generate both response types. This is despite one response (stochastic pulsing) being stochastic in nature, and the other (single response pulse) being deterministic. We demonstrate that the main determinant for whichever response is generated is the degree with which the input pathway activates the core circuit, although the noise properties of the input pathway also biases the system towards one or the other type of output. Thus, our work shows how stochastic modelling can reveal the mechanisms behind non-intuitive gene circuit output dynamics.

Increased gene expression variability hinders the formation of regional mechanical conflicts leading to reduced organ shape robustness

To relate gene networks and organ shape, one needs to address two wicked problems: i) Gene expression is often variable locally, and shape is reproducible globally; ii) gene expression can have cascading effects on tissue mechanics, with possibly counterintuitive consequences for the final organ shape. Here, we address such wicked problems, taking advantage of simpler plant organ development where shape only emerges from cell division and elongation. We confirm that mutation in VERNALIZATION INDEPENDENCE 3 (VIP3), a subunit of the conserved polymerase-associated factor 1 complex (Paf1C), increases gene expression variability in Arabidopsis. Then, we focused on the Arabidopsis sepal, which exhibits a reproducible shape and stereotypical regional growth patterns. In vip3 sepals, we measured higher growth heterogeneity between adjacent cells. This even culminated in the presence of negatively growing cells in specific growth conditions. Interestingly, such increased local noise interfered with the stereotypical regional pattern of growth. We previously showed that regional differential growth at the wild-type sepal tip triggers a mechanical conflict, to which cells resist by reinforcing their walls, leading to growth arrest. In vip3, the disturbed regional growth pattern delayed organ growth arrest and increased final organ shape variability. Altogether, we propose that gene expression variability is managed by Paf1C to ensure organ robustness by building up mechanical conflicts at the regional scale, instead of the local scale.

Systematic analysis and optimization of early warning signals for critical transitions using distribution data

Abrupt shifts between alternative regimes occur in complex systems, from cell regulation to brain functions to ecosystems. Several model-free early warning signals (EWS) have been proposed to detect impending transitions, but failure or poor performance in some systems have called for better investigation of their generic applicability. Notably, there are still ongoing debates whether such signals can be successfully extracted from data in particular from biological experiments. In this work, we systematically investigate properties and performance of dynamical EWS in different deteriorating conditions, and we propose an optimized combination to trigger warnings as early as possible, eventually verified on experimental data from microbiological populations. Our results explain discrepancies observed in the literature between warning signs extracted from simulated models and from real data, provide guidance for EWS selection based on desired systems and suggest an optimized composite indicator to alert for impending critical transitions using distribution data.

Deep neural network analysis models for complex random telegraph signals

Time-fluctuating signals are ubiquitous and diverse in many physical, chemical, and biological systems, among which random telegraph signals (RTSs) refer to a series of instantaneous switching events between two discrete levels from single-particle movements. A reliable RTS analysis is a crucial prerequisite to identify underlying mechanisms related to device performance and sensitivity. When numerous levels are involved, complex patterns of multilevel RTSs occur and make their quantitative analysis exponentially difficult, hereby systematic approaches are often elusive. In this work, we present a three-step analysis protocol via progressive knowledge-transfer, where the outputs of the early step are passed onto a subsequent step. Especially, to quantify complex RTSs, we resort to three deep neural network architectures whose trained models can process raw temporal data directly. We furthermore demonstrate the model accuracy extensively with a large dataset of different RTS types in terms of additional background noise types and amplitude size. Our protocol offers structured schemes to extract the parameter values of complex RTSs as imperative information with which researchers can draw meaningful and relevant interpretations and inferences of given devices and systems.

Self-organized collective cell behaviors as design principles for synthetic developmental biology

Over the past two decades, molecular cell biology has graduated from a mostly analytic science to one with substantial synthetic capability. This success is built on a deep understanding of the structure and function of biomolecules and molecular mechanisms. For synthetic biology to achieve similar success at the scale of tissues and organs, an equally deep understanding of the principles of development is required. Here, we review some of the central concepts and recent progress in tissue patterning, morphogenesis and collective cell migration and discuss their value for synthetic developmental biology, emphasizing in particular the power of (guided) self-organization and the role of theoretical advances in making developmental insights applicable in synthesis.

Noise and vestibular perception of passive self-motion

Noise defined as random disturbances is ubiquitous in both the external environment and the nervous system. Depending on the context, noise can degrade or improve information processing and performance. In all cases, it contributes to neural systems dynamics. We review some effects of various sources of noise on the neural processing of self-motion signals at different stages of the vestibular pathways and the resulting perceptual responses. Hair cells in the inner ear reduce the impact of noise by means of mechanical and neural filtering. Hair cells synapse on regular and irregular afferents. Variability of discharge (noise) is low in regular afferents and high in irregular units. The high variability of irregular units provides information about the envelope of naturalistic head motion stimuli. A subset of neurons in the vestibular nuclei and thalamus are optimally tuned to noisy motion stimuli that reproduce the statistics of naturalistic head movements. In the thalamus, variability of neural discharge increases with increasing motion amplitude but saturates at high amplitudes, accounting for behavioral violation of Weber's law. In general, the precision of individual vestibular neurons in encoding head motion is worse than the perceptual precision measured behaviorally. However, the global precision predicted by neural population codes matches the high behavioral precision. The latter is estimated by means of psychometric functions for detection or discrimination of whole-body displacements. Vestibular motion thresholds (inverse of precision) reflect the contribution of intrinsic and extrinsic noise to perception. Vestibular motion thresholds tend to deteriorate progressively after the age of 40 years, possibly due to oxidative stress resulting from high discharge rates and metabolic loads of vestibular afferents. In the elderly, vestibular thresholds correlate with postural stability: the higher the threshold, the greater is the postural imbalance and risk of falling. Experimental application of optimal levels of either galvanic noise or whole-body oscillations can ameliorate vestibular function with a mechanism reminiscent of stochastic resonance. Assessment of vestibular thresholds is diagnostic in several types of vestibulopathies, and vestibular stimulation might be useful in vestibular rehabilitation.

Dynamic fluctuations in a bacterial metabolic network

The operation of the central metabolism is typically assumed to be deterministic, but dynamics and high connectivity of the metabolic network make it potentially prone to generating fluctuations. However, time-resolved measurements of metabolite levels in individual cells that are required to characterize such fluctuations remained a challenge, particularly in small bacterial cells. Here we use single-cell metabolite measurements based on Förster resonance energy transfer, combined with computer simulations, to explore the real-time dynamics of the metabolic network of Escherichia coli. We observe that steplike exposure of starved E. coli to glycolytic carbon sources elicits large periodic fluctuations in the intracellular concentration of pyruvate in individual cells. These fluctuations are consistent with predicted oscillatory dynamics of E. coli metabolic network, and they are primarily controlled by biochemical reactions around the pyruvate node. Our results further indicate that fluctuations in glycolysis propagate to other cellular processes, possibly leading to temporal heterogeneity of cellular states within a population.

Autonomous optimization of neuroprosthetic stimulation parameters that drive the motor cortex and spinal cord outputs in rats and monkeys

Neural stimulation can alleviate paralysis and sensory deficits. Novel high-density neural interfaces can enable refined and multipronged neurostimulation interventions. To achieve this, it is essential to develop algorithmic frameworks capable of handling optimization in large parameter spaces. Here, we leveraged an algorithmic class, Gaussian-process (GP)-based Bayesian optimization (BO), to solve this problem. We show that GP-BO efficiently explores the neurostimulation space, outperforming other search strategies after testing only a fraction of the possible combinations. Through a series of real-time multi-dimensional neurostimulation experiments, we demonstrate optimization across diverse biological targets (brain, spinal cord), animal models (rats, non-human primates), in healthy subjects, and in neuroprosthetic intervention after injury, for both immediate and continual learning over multiple sessions. GP-BO can embed and improve "prior" expert/clinical knowledge to dramatically enhance its performance. These results advocate for broader establishment of learning agents as structural elements of neuroprosthetic design, enabling personalization and maximization of therapeutic effectiveness.

Actomyosin pulsing rescues embryonic tissue folding from disruption by myosin fluctuations

During early development, myosin II mechanically reshapes and folds embryo tissue. A much-studied example is ventral furrow formation in Drosophila , marking the onset of gastrulation. Furrowing is driven by contraction of actomyosin networks on apical cell surfaces, but how the myosin patterning encodes tissue shape is unclear, and elastic models failed to reproduce essential features of experimental cell contraction profiles. The myosin patterning exhibits substantial cell-to-cell fluctuations with pulsatile time-dependence, a striking but unexplained feature of morphogenesis in many organisms. Here, using biophysical modeling we find viscous forces offer the principle resistance to actomyosin-driven apical constriction. In consequence, tissue shape is encoded in the direction-dependent curvature of the myosin patterning which orients an anterior-posterior furrow. Tissue contraction is highly sensitive to cell-to-cell myosin fluctuations, explaining furrowing failure in genetically perturbed embryos whose fluctuations are temporally persistent. In wild-type embryos, this catastrophic outcome is averted by pulsatile myosin time-dependence, a time-averaging effect that rescues furrowing. This low pass filter mechanism may underlie the usage of actomyosin pulsing in diverse morphogenetic processes across many organisms.

Cell state transitions: catch them if you can

The Company of Biologists' 2022 workshop on 'Cell State Transitions: Approaches, Experimental Systems and Models' brought together an international and interdisciplinary team of investigators spanning the fields of cell and developmental biology, stem cell biology, physics, mathematics and engineering to tackle the question of how cells precisely navigate between distinct identities and do so in a dynamic manner. This second edition of the workshop was organized after a successful virtual workshop on the same topic that took place in 2021.

On multistability and constitutive relations of cell motion on fibronectin lanes

Cell motility on flat substrates exhibits coexisting steady and oscillatory morphodynamics, the biphasic adhesion-velocity relation, and the universal correlation between speed and persistence (UCSP) as simultaneous observations common to many cell types. Their universality and concurrency suggest a unifying mechanism causing all three of them. Stick-slip models for cells on one-dimensional lanes suggest multistability to arise from the nonlinear friction of retrograde flow. This study suggests a mechanical mechanism controlled by integrin signaling on the basis of a biophysical model and analysis of trajectories of MDA-MB-231 cells on fibronectin lanes, which additionally explains the constitutive relations. The experiments exhibit cells with steady or oscillatory morphodynamics and either spread or moving with spontaneous transitions between the dynamic regimes, spread and moving, and spontaneous direction reversals. Our biophysical model is based on the force balance at the protrusion edge, the noisy clutch of retrograde flow, and a response function of friction and membrane drag to integrin signaling. The theory reproduces the experimentally observed cell states, characteristics of oscillations, and state probabilities. Analysis of experiments with the biophysical model establishes a stick-slip oscillation mechanism, and explains multistability of cell states and the statistics of state transitions. It suggests protrusion competition to cause direction reversal events, the statistics of which explain the UCSP. The effect of integrin signaling on drag and friction explains the adhesion-velocity relation and cell behavior at fibronectin density steps. The dynamics of our mechanism are nonlinear flow mechanics driven by F-actin polymerization and shaped by the noisy clutch of retrograde flow friction, protrusion competition via membrane tension, and drag forces. Integrin signaling controls the parameters of the mechanical system.

Mechanistic modeling of viral particle production

Viral systems such as wild-type viruses, viral vectors, and virus-like particles are essential components of modern biotechnology and medicine. Despite their importance, the commercial-scale production of viral systems remains highly inefficient for multiple reasons. Computational strategies are a promising avenue for improving process development, optimization, and control, but require a mathematical description of the system. This article reviews mechanistic modeling strategies for the production of viral particles, both at the cellular and bioreactor scales. In many cases, techniques and models from adjacent fields such as epidemiology and wild-type viral infection kinetics can be adapted to construct a suitable process model. These process models can then be employed for various purposes such as in-silico testing of novel process operating strategies and/or advanced process control.

Making use of noise in biological systems

Disorder and noise are inherent in biological systems. They are required to provide systems with the advantages required for proper functioning. Noise is a part of the flexibility and plasticity of biological systems. It provides systems with increased routes, improves information transfer, and assists in response triggers. This paper reviews recent studies on noise at the genome, cellular, and whole organ levels. We focus on the need to use noise in system engineering. We present some of the challenges faced in studying noise. Optimizing the efficiency of complex systems requires a degree of variability in their functions within certain limits. Constrained noise can be considered a method for improving system robustness by regulating noise levels in continuously dynamic settings. The digital pill-based artificial intelligence (AI)-based platform is the first to implement second-generation AI comprising variability-based signatures. This platform enhances the efficacy of the therapeutic regimens. Systems requiring variability and mechanisms regulating noise are mandatory for understanding biological functions.

Precision of morphogen-driven tissue patterning during development is enhanced through contact-mediated cellular interactions

Cells in developing embryos reliably differentiate to attain location-specific fates, despite fluctuations in morphogen concentrations that provide positional information and in molecular processes that interpret it. We show that local contact-mediated cell-cell interactions utilize inherent asymmetry in the response of patterning genes to the global morphogen signal yielding a bimodal response. This results in robust developmental outcomes with a consistent identity for the dominant gene at each cell, substantially reducing the uncertainty in the location of boundaries between distinct fates.

Learning in Transcriptional Network Models: Computational Discovery of Pathway-Level Memory and Effective Interventions

Trainability, in any substrate, refers to the ability to change future behavior based on past experiences. An understanding of such capacity within biological cells and tissues would enable a particularly powerful set of methods for prediction and control of their behavior through specific patterns of stimuli. This top-down mode of control (as an alternative to bottom-up modification of hardware) has been extensively exploited by computer science and the behavioral sciences; in biology however, it is usually reserved for organism-level behavior in animals with brains, such as training animals towards a desired response. Exciting work in the field of basal cognition has begun to reveal degrees and forms of unconventional memory in non-neural tissues and even in subcellular biochemical dynamics. Here, we characterize biological gene regulatory circuit models and protein pathways and find them capable of several different kinds of memory. We extend prior results on learning in binary transcriptional networks to continuous models and identify specific interventions (regimes of stimulation, as opposed to network rewiring) that abolish undesirable network behavior such as drug pharmacoresistance and drug sensitization. We also explore the stability of created memories by assessing their long-term behavior and find that most memories do not decay over long time periods. Additionally, we find that the memory properties are quite robust to noise; surprisingly, in many cases noise actually increases memory potential. We examine various network properties associated with these behaviors and find that no one network property is indicative of memory. Random networks do not show similar memory behavior as models of biological processes, indicating that generic network dynamics are not solely responsible for trainability. Rational control of dynamic pathway function using stimuli derived from computational models opens the door to empirical studies of proto-cognitive capacities in unconventional embodiments and suggests numerous possible applications in biomedicine, where behavior shaping of pathway responses stand as a potential alternative to gene therapy.

Bistability and noise-induced transient behaviour of steady states in a cancer network with the regulation of microRNA

MicroRNAs (miRs) regulatory network models are highly non-linear due to the negative regulation of gene expression at the post-transcriptional level by miRs and can produce interesting dynamics of the system such as bistability for miR-transcriptional factors interactions. In this article, we focus on the miR-17-92 cluster and its interaction with transcriptional factors (proteins) E2F and Myc. Environmental fluctuations (noise) and randomness in the bio-chemical reactions can be very important and change the dynamical role of miR-17-92 in the regulatory network. We have formulated a stochastically forced miR-17-92 and E2F-Myc interaction model and study the phenomena of intrinsic and extrinsic noise which can induce random switching between steady states or the destruction of the bistability. Using a method that employs stochastic sensitivity functions we have constructed confidence ellipses to determine the configurational arrangements of equilibrium and spatial arrangements of random states near stable equilibria. Simulations are carried out to numerically show the flow of the solution trajectories under noise. Finally, we summarize the simulation results and the impact of noise on the dual non-linear role of miR-17-92 cluster to act as an oncogene or as a tumour suppressor gene.

Computationally efficient mechanism discovery for cell invasion with uncertainty quantification

Parameter estimation for mathematical models of biological processes is often difficult and depends significantly on the quality and quantity of available data. We introduce an efficient framework using Gaussian processes to discover mechanisms underlying delay, migration, and proliferation in a cell invasion experiment. Gaussian processes are leveraged with bootstrapping to provide uncertainty quantification for the mechanisms that drive the invasion process. Our framework is efficient, parallelisable, and can be applied to other biological problems. We illustrate our methods using a canonical scratch assay experiment, demonstrating how simply we can explore different functional forms and develop and test hypotheses about underlying mechanisms, such as whether delay is present. All code and data to reproduce this work are available at https://github.com/DanielVandH/EquationLearning.jl.

3D Tomographic Analysis of the Order-Disorder Interplay in the Pachyrhynchus congestus mirabilis Weevil

The bright colors of Pachyrhynchus weevils originate from complex dielectric nanostructures within their elytral scales. In contrast to previous work exhibiting highly ordered single-network diamond-type photonic crystals, here, it is shown by combining optical microscopy and spectroscopy measurements with 3D focused ion beam (FIB) tomography that the blue scales of P. congestus mirabilis differ from that of an ordered diamond structure. Through the use of FIB tomography on elytral scales filled with platinum (Pt) by electron beam-assisted deposition, it is revealed that the red scales of this weevil possess a periodic diamond structure, while the network morphology of the blue scales exhibit diamond morphology only on the single scattering unit level with disorder on longer length scales. Full wave simulations performed on the reconstructed volumes indicate that this local order is sufficient to open a partial photonic bandgap even at low dielectric constant contrast between chitin and air in the absence of long-range or translational order. The observation of disordered and ordered photonic crystals within a single organism opens up interesting questions on the cellular origin of coloration and studies on bio-inspired replication of angle-independent colors.

Confinement-induced accumulation and de-mixing of microscopic active-passive mixtures

Understanding the out-of-equilibrium properties of noisy microscale systems and the extent to which they can be modulated externally, is a crucial scientific and technological challenge. It holds the promise to unlock disruptive new technologies ranging from targeted delivery of chemicals within the body to directed assembly of new materials. Here we focus on how active matter can be harnessed to transport passive microscopic systems in a statistically predictable way. Using a minimal active-passive system of weakly Brownian particles and swimming microalgae, we show that spatial confinement leads to a complex non-monotonic steady-state distribution of colloids, with a pronounced peak at the boundary. The particles’ emergent active dynamics is well captured by a space-dependent Poisson process resulting from the space-dependent motion of the algae. Based on our findings, we then realise experimentally the de-mixing of the active-passive suspension, opening the way for manipulating colloidal objects via controlled activity fields.

Understanding how order emerges in active matter may facilitate macroscopic control of microscopic objects. Here, Williams et al. show how to control the transport of passive microscopic particles in presence of motile algae in conjunction with boundary-induced accumulation of microswimmers.

Persistent Correlation in Cellular Noise Determines Longevity of Viral Infections

The slowly decaying viral dynamics, even after 2-3 weeks from diagnosis, is one of the characteristics of COVID-19 infection that is still unexplored in theoretical and experimental studies. This long-lived characteristic of viral infections in the framework of inherent variations or noise present at the cellular level is often overlooked. Therefore, in this work, we aim to understand the effect of these variations by proposing a stochastic non-Markovian model that not only captures the coupled dynamics between the immune cells and the virus but also enables the study of the effect of fluctuations. Numerical simulations of our model reveal that the long-range temporal correlations in fluctuations dictate the long-lived dynamics of a viral infection and, in turn, also affect the rates of immune response. Furthermore, predictions of our model system are in agreement with the experimental viral load data of COVID-19 patients from various countries.

A general theory for temperature dependence in biology

At present, there is no simple, first principles-based, and general model for quantitatively describing the full range of observed biological temperature responses. Here we derive a general theory for temperature dependence in biology based on Eyring-Evans-Polanyi's theory for chemical reaction rates. Assuming only that the conformational entropy of molecules changes with temperature, we derive a theory for the temperature dependence of enzyme reaction rates which takes the form of an exponential function modified by a power law and that describes the characteristic asymmetric curved temperature response. Based on a few additional principles, our model can be used to predict the temperature response above the enzyme level, thus spanning quantum to classical scales. Our theory provides an analytical description for the shape of temperature response curves and demonstrates its generality by showing the convergence of all temperature dependence responses onto universal relationships-a universal data collapse-under appropriate normalization and by identifying a general optimal temperature, around 25 ^∘C, characterizing all temperature response curves. The model provides a good fit to empirical data for a wide variety of biological rates, times, and steady-state quantities, from molecular to ecological scales and across multiple taxonomic groups (from viruses to mammals). This theory provides a simple framework to understand and predict the impact of temperature on biological quantities based on the first principles of thermodynamics, bridging quantum to classical scales.

Stochastic fluctuations promote ordered pattern formation of cells in the Notch-Delta signaling pathway

The Notch-Delta signaling pathway mediates cell differentiation implicated in many regulatory processes including spatiotemporal patterning in tissues by promoting alternate cell fates between neighboring cells. At the multicellular level, this "lateral inhibition” principle leads to checkerboard patterns with alternation of Sender and Receiver cells. While it is well known that stochasticity modulates cell fate specification, little is known about how stochastic fluctuations at the cellular level propagate during multicell pattern formation. Here, we model stochastic fluctuations in the Notch-Delta pathway in the presence of two different noise types–shot and white–for a multicell system. Our results show that intermediate fluctuations reduce disorder and guide the multicell lattice toward checkerboard-like patterns. By further analyzing cell fate transition events, we demonstrate that intermediate noise amplitudes provide enough perturbation to facilitate “proofreading” of disordered patterns and cause cells to switch to the correct ordered state (Sender surrounded by Receivers, and vice versa). Conversely, high noise can override environmental signals coming from neighboring cells and lead to switching between ordered and disordered patterns. Therefore, in analogy with spin glass systems, intermediate noise levels allow the multicell Notch system to escape frustrated patterns and relax towards the lower energy checkerboard pattern while at large noise levels the system is unable to find this ordered base of attraction.

The Notch pathway is involved in many biological processes and is known to form precise spatial patterns alternating Sender and Receiver cell states. Quantifying the implications of stochastic fluctuations provided insight that patterns formed in Notch-mediated pathways must follow a predetermined path towards checkerboard or exist in a noisy environment which promotes order through error correction.

We model Notch pattern formation stochastically and analyze the spatiotemporal dynamics. Our results show multicellular systems equilibrate towards ordered systems, but mistakes in the initial lattice propagate causing the systems to relax into frustrated systems. Only through existing in a noisy environment are the systems able to relax into the checkerboard pattern. Analyzing the temporal dynamics confirms, in environments with intermediate noise, the “incorrect” cells (Sender in a Sender environment, and vice versa) can be flipped to the correct state (Sender in a Receiver environment, and vice versa). Comparing with the spin glass energy landscape, we suggest the multicellular model follows a rugged landscape to form patterns with stochastic fluctuations required to enforce order throughout the system.

Hypertension: Constraining the Expression of ACE-II by Adopting Optimal Macronutrients Diet Predicted via Support Vector Machine

The recent elevation of cases infected from novel COVID-19 has placed the human life in trepidation mode, especially for those suffering from comorbidities. Most of the studies in the last few months have undeniably raised concerns for hypertensive patients that face greater risk of fatality from COVID-19. Furthermore, one of the recent WHO reports has estimated a total of 1.13 billion people are at a risk of hypertension of which two-thirds live in low and middle income countries. The gradual escalation of the hypertension problem andthe sudden rise of COVID-19 cases have placed an increasingly higher number of human lives at risk in low and middle income countries. To lower the risk of hypertension, most physicians recommend drugs that have angiotensin-converting enzyme (ACE) inhibitors. However, prolonged use of such drugs is not recommended due to metabolic risks and the increase in the expression of ACE-II which could facilitate COVID-19 infection. In contrast, the intake of optimal macronutrients is one of the possible alternatives to naturally control hypertension. In the present study, a nontrivial feature selection and machine learning algorithm is adopted to intelligently predict the food-derived antihypertensive peptide. The proposed idea of the paper lies in reducing the computational power while retaining the performance of the support vector machine (SVM) by estimating the dominant pattern in the features space through feature filtering. The proposed feature filtering algorithm has reported a trade-off performance by reducing the chances of Type I error, which is desirable when recommending a dietary food to patients suffering from hypertension. The maximum achievable accuracy of the best performing SVM models through feature selection are 86.17% and 85.61%, respectively.

Fluctuating asymmetry as an indicator of stress

Fluctuating asymmetry as a special kind of asymmetry can be defined as deviations from a known predetermined ratio of the parts of morphological structure under study. As a special type of phenotypic variability fluctuating asymmetry is a manifestation of ontogenetic noise or developmental variability. This type of variability is ubiquitous and plays a significant role in the observed phenotypic diversity. The level of fluctuating asymmetry turns out to be an indicator of optimal developmental conditions and genetic coadaptation. It is also considered as a parameter of fitness. Thus, fluctuating asymmetry acts as a measure of developmental stability in developmental biology and as a measure of population condition in population biology.

Mutual information in changing environments: Nonlinear interactions, out-of-equilibrium systems, and continuously varying diffusivities

Biochemistry, ecology, and neuroscience are examples of prominent fields aiming at describing interacting systems that exhibit nontrivial couplings to complex, ever-changing environments. We have recently shown that linear interactions and a switching environment are encoded separately in the mutual information of the overall system. Here we first generalize these findings to a broad class of nonlinear interacting models. We find that a new term in the mutual information appears, quantifying the interplay between nonlinear interactions and environmental changes, and leading to either constructive or destructive information interference. Furthermore, we show that a higher mutual information emerges in out-of-equilibrium environments with respect to an equilibrium scenario. Finally, we generalize our framework to the case of continuously varying environments. We find that environmental changes can be mapped exactly into an effective spatially varying diffusion coefficient, shedding light on modeling of biophysical systems in inhomogeneous media.

Evaluation of statistical approaches for association testing in noisy drug screening data

BACKGROUND

Identifying associations among biological variables is a major challenge in modern quantitative biological research, particularly given the systemic and statistical noise endemic to biological systems. Drug sensitivity data has proven to be a particularly challenging field for identifying associations to inform patient treatment.

RESULTS

To address this, we introduce two semi-parametric variations on the commonly used concordance index: the robust concordance index and the kernelized concordance index (rCI, kCI), which incorporate measurements about the noise distribution from the data. We demonstrate that common statistical tests applied to the concordance index and its variations fail to control for false positives, and introduce efficient implementations to compute p-values using adaptive permutation testing. We then evaluate the statistical power of these coefficients under simulation and compare with Pearson and Spearman correlation coefficients. Finally, we evaluate the various statistics in matching drugs across pharmacogenomic datasets.

CONCLUSIONS

We observe that the rCI and kCI are better powered than the concordance index in simulation and show some improvement on real data. Surprisingly, we observe that the Pearson correlation was the most robust to measurement noise among the different metrics.

Enhancing Synchronization by Optimal Correlated Noise

From the flashes of fireflies to Josephson junctions and power infrastructure, networks of coupled phase oscillators provide a powerful framework to describe synchronization phenomena in many natural and engineered systems. Most real-world networks are under the influence of noisy, random inputs, potentially inhibiting synchronization. While noise is unavoidable, here we show that there exist optimal noise patterns which minimize desynchronizing effects and even enhance order. Specifically, using analytical arguments we show that in the case of a two-oscillator model, there exists a sharp transition from a regime where the optimal synchrony-enhancing noise is perfectly anticorrelated, to one where the optimal noise is correlated. More generally, we then use numerical optimization methods to demonstrate that there exist anticorrelated noise patterns that optimally enhance synchronization in large complex oscillator networks. Our results may have implications in networks such as power grids and neuronal networks, which are subject to significant amounts of correlated input noise.

Pathway-resolved decomposition demonstrates correlation and noise dependencies of redundant information processing in recurrent feed-forward topologies

In a biochemical assay that converts fan-in networks into feed-forward loops (FFLs), we show that the inter-regulator redundant information about the output gene product can be decomposed into finer components, mediated by the constituent pathways. Variance-based information within the linear noise regime facilitates quantifying these submodular redundancies. Contrary to the conventional wisdom on information decomposition, we report that information redundancy depends nontrivially on inter-regulator correlation. For the type-1 coherent (C1) and incoherent (I1) FFLs, the direct regulatory path-mediated redundancy is certainly correlation independent. However, components induced by the indirect regulatory path and interpathway interference are correlation dependent in (non)linear fashion. The trade-off between information redundancy and similarly decomposable extrinsic noise from input to output node has been demonstrated for the pathways and full motifs. Our analyses suggest that the interpathway cross redundancy positively and negatively influences the superposition of elementary redundancies in the C1- and I1-FFLs, respectively. Their corresponding total extrinsic noise is produced by the weighted sum and difference of the pathway-specific components. We find that the I1-FFL is able to manufacture more varied redundancy and extrinsic noise responses compared to the C1-FFL. Underlying the differing characteristics of the composite metrics across FFL variants, there exist uniformly behaving pathway-dependent elements. The decomposition framework has been meticulously explored in biologically rational parametric realizations through analytical estimates and stochastic simulations.

Noise resistant synchronization and collective rhythm switching in a model of animal group locomotion

Biology is suffused with rhythmic behaviour, and interacting biological oscillators often synchronize their rhythms with one another. Colonies of some ant species are able to synchronize their activity to fall into coherent bursts, but models of this phenomenon have neglected the potential effects of intrinsic noise and interspecific differences in individual-level behaviour. We investigated the individual and collective activity patterns of two Leptothorax ant species. We show that in one species (Leptothorax sp. W), ants converge onto rhythmic cycles of synchronized collective activity with a period of about 20 min. A second species (Leptothorax crassipilis) exhibits more complex collective dynamics, where dominant collective cycle periods range from 16 min to 2.8 h. Recordings that last 35 h reveal that, in both species, the same colony can exhibit multiple oscillation frequencies. We observe that workers of both species can be stimulated by nest-mates to become active after a refractory resting period, but the durations of refractory periods differ between the species and can be highly variable. We model the emergence of synchronized rhythms using an agent-based model informed by our empirical data. This simple model successfully generates synchronized group oscillations despite the addition of noise to ants' refractory periods. We also find that adding noise reduces the likelihood that the model will spontaneously switch between distinct collective cycle frequencies.