Showing 1–50 of 212 results for author: Eduardo

Large-scale Epidemiological modeling: Scanning for Mosquito-Borne Diseases Spatio-temporal Patterns in Brazil

Authors: Eduardo C. Araujo, Claudia T. Codeço, Sandro Loch, Luã B. Vacaro, Laís P. Freitas, Raquel M. Lana, Leonardo S. Bastos, Iasmim F. de Almeida, Fernanda Valente, Luiz M. Carvalho, Flávio C. Coelho

Abstract: The influence of climate on mosquito-borne diseases like dengue and chikungunya is well-established, but comprehensively tracking long-term spatial and temporal trends across large areas has been hindered by fragmented data and limited analysis tools. This study presents an unprecedented analysis, in terms of breadth, estimating the SIR transmission parameters from incidence data in all 5,570 muni… ▽ More The influence of climate on mosquito-borne diseases like dengue and chikungunya is well-established, but comprehensively tracking long-term spatial and temporal trends across large areas has been hindered by fragmented data and limited analysis tools. This study presents an unprecedented analysis, in terms of breadth, estimating the SIR transmission parameters from incidence data in all 5,570 municipalities in Brazil over 14 years (2010-2023) for both dengue and chikungunya. We describe the Episcanner computational pipeline, developed to estimate these parameters, producing a reusable dataset describing all dengue and chikungunya epidemics that have taken place in this period, in Brazil. The analysis reveals new insights into the climate-epidemic nexus: We identify distinct geographical and temporal patterns of arbovirus disease incidence across Brazil, highlighting how climatic factors like temperature and precipitation influence the timing and intensity of dengue and chikungunya epidemics. The innovative Episcanner tool empowers researchers and public health officials to explore these patterns in detail, facilitating targeted interventions and risk assessments. This research offers a new perspective on the long-term dynamics of climate-driven mosquito-borne diseases and their geographical specificities linked to the effects of global temperature fluctuations such as those captured by the ENSO index. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: Submitted for peer review

A Scalable Tool For Analyzing Genomic Variants Of Humans Using Knowledge Graphs and Machine Learning

Authors: Shivika Prasanna, Ajay Kumar, Deepthi Rao, Eduardo Simoes, Praveen Rao

Abstract: The integration of knowledge graphs and graph machine learning (GML) in genomic data analysis offers several opportunities for understanding complex genetic relationships, especially at the RNA level. We present a comprehensive approach for leveraging these technologies to analyze genomic variants, specifically in the context of RNA sequencing (RNA-seq) data from COVID-19 patient samples. The prop… ▽ More The integration of knowledge graphs and graph machine learning (GML) in genomic data analysis offers several opportunities for understanding complex genetic relationships, especially at the RNA level. We present a comprehensive approach for leveraging these technologies to analyze genomic variants, specifically in the context of RNA sequencing (RNA-seq) data from COVID-19 patient samples. The proposed method involves extracting variant-level genetic information, annotating the data with additional metadata using SnpEff, and converting the enriched Variant Call Format (VCF) files into Resource Description Framework (RDF) triples. The resulting knowledge graph is further enhanced with patient metadata and stored in a graph database, facilitating efficient querying and indexing. We utilize the Deep Graph Library (DGL) to perform graph machine learning tasks, including node classification with GraphSAGE and Graph Convolutional Networks (GCNs). Our approach demonstrates significant utility using our proposed tool, VariantKG, in three key scenarios: enriching graphs with new VCF data, creating subgraphs based on user-defined features, and conducting graph machine learning for node classification. △ Less

Submitted 30 July, 2024; originally announced July 2024.

Comments: arXiv admin note: substantial text overlap with arXiv:2312.04423

Decoding the interaction mediators from landscape-induced spatial patterns

Authors: E. H. Colombo, L. Defaveri, C. Anteneodo

Abstract: Interactions between organisms are mediated by an intricate network of physico-chemical substances and other organisms. Understanding the dynamics of mediators and how they shape the population spatial distribution is key to predict ecological outcomes and how they would be transformed by changes in environmental constraints. However, due to the inherent complexity involved, this task is often unf… ▽ More Interactions between organisms are mediated by an intricate network of physico-chemical substances and other organisms. Understanding the dynamics of mediators and how they shape the population spatial distribution is key to predict ecological outcomes and how they would be transformed by changes in environmental constraints. However, due to the inherent complexity involved, this task is often unfeasible, from the empirical and theoretical perspectives. In this paper, we make progress in addressing this central issue, creating a bridge that provides a two-way connection between the features of the ensemble of underlying mediators and the wrinkles in the population density induced by a landscape defect (or spatial perturbation). The bridge is constructed by applying the Feynman-Vernon decomposition, which disentangles the influences among the focal population and the mediators in a compact way. This is achieved though an interaction kernel, which effectively incorporates the mediators' degrees of freedom, explaining the emergence of nonlocal influence between individuals, an ad hoc assumption in modeling population dynamics. Concrete examples are worked out and reveal the complexity behind a possible top-down inference procedure. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Mutation-Bias Learning in Games

Authors: Johann Bauer, Sheldon West, Eduardo Alonso, Mark Broom

Abstract: We present two variants of a multi-agent reinforcement learning algorithm based on evolutionary game theoretic considerations. The intentional simplicity of one variant enables us to prove results on its relationship to a system of ordinary differential equations of replicator-mutator dynamics type, allowing us to present proofs on the algorithm's convergence conditions in various settings via its… ▽ More We present two variants of a multi-agent reinforcement learning algorithm based on evolutionary game theoretic considerations. The intentional simplicity of one variant enables us to prove results on its relationship to a system of ordinary differential equations of replicator-mutator dynamics type, allowing us to present proofs on the algorithm's convergence conditions in various settings via its ODE counterpart. The more complicated variant enables comparisons to Q-learning based algorithms. We compare both variants experimentally to WoLF-PHC and frequency-adjusted Q-learning on a range of settings, illustrating cases of increasing dimensionality where our variants preserve convergence in contrast to more complicated algorithms. The availability of analytic results provides a degree of transferability of results as compared to purely empirical case studies, illustrating the general utility of a dynamical systems perspective on multi-agent reinforcement learning when addressing questions of convergence and reliable generalisation. △ Less

Submitted 28 May, 2024; originally announced May 2024.

MSC Class: 37N40 (Primary) 91A26 (Secondary)

Comparison of reaction networks of insulin signaling

Authors: Patrick Vincent N. Lubenia, Eduardo R. Mendoza, Angelyn R. Lao

Abstract: Understanding the insulin signaling cascade provides insights on the underlying mechanisms of biological phenomena such as insulin resistance, diabetes, Alzheimer's disease, and cancer. For this reason, previous studies utilized chemical reaction network theory to perform comparative analyses of reaction networks of insulin signaling in healthy (INSMS: INSulin Metabolic Signaling) and diabetic cel… ▽ More Understanding the insulin signaling cascade provides insights on the underlying mechanisms of biological phenomena such as insulin resistance, diabetes, Alzheimer's disease, and cancer. For this reason, previous studies utilized chemical reaction network theory to perform comparative analyses of reaction networks of insulin signaling in healthy (INSMS: INSulin Metabolic Signaling) and diabetic cells (INRES: INsulin RESistance). This study extends these analyses using various methods which give further insights regarding insulin signaling. Using embedded networks, we discuss evidence of the presence of a structural "bifurcation" in the signaling process between INSMS and INRES. Concordance profiles of INSMS and INRES show that both have a high propensity to remain monostationary. Moreover, the concordance properties allow us to present heuristic evidence that INRES has a higher level of stability beyond its monostationarity. Finally, we discuss a new way of analyzing reaction networks through network translation. This method gives rise to three new insights: (i) each stoichiometric class of INSMS and INRES contains a unique positive equilibrium; (ii) any positive equilibrium of INSMS is exponentially stable and is a global attractor in its stoichiometric class; and (iii) any positive equilibrium of INRES is locally asymptotically stable. These results open up opportunities for collaboration with experimental biologists to understand insulin signaling better. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 18 pages, 0 figure

Interplay between Contractivity and Monotonicity for Reaction Networks

Authors: Alon Duvall, M. Ali Al-Radhawi, Dhruv D. Jatkar, Eduardo Sontag

Abstract: This work studies relationships between monotonicity and contractivity, and applies the results to establish that many reaction networks are weakly contractive, and thus, under appropriate compactness conditions, globally convergent to equilibria. Verification of these properties is achieved through a novel algorithm that can be used to generate cones for monotone systems. The results given here a… ▽ More This work studies relationships between monotonicity and contractivity, and applies the results to establish that many reaction networks are weakly contractive, and thus, under appropriate compactness conditions, globally convergent to equilibria. Verification of these properties is achieved through a novel algorithm that can be used to generate cones for monotone systems. The results given here allow a unified proof of global convergence for several classes of networks that had been previously studied in the literature. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Embedding-based comparison of reaction networks of Wnt signaling

Authors: Bryan S. Hernandez, Patrick Vincent N. Lubenia, Eduardo R. Mendoza

Abstract: This work introduces a new method for comparing two reaction networks of the same or closely related systems through their embedded networks in terms of the shared set of species. Hence, we call this method the Common Species Embedded Networks (CSEN) analysis. Using this approach, we conduct a comparison of existing reaction networks associated with Wnt signaling models (Lee, Schmitz, MacLean, and… ▽ More This work introduces a new method for comparing two reaction networks of the same or closely related systems through their embedded networks in terms of the shared set of species. Hence, we call this method the Common Species Embedded Networks (CSEN) analysis. Using this approach, we conduct a comparison of existing reaction networks associated with Wnt signaling models (Lee, Schmitz, MacLean, and Feinberg) that we have identified. The analysis yields three important results for these Wnt models. First, the CSEN analysis of the Lee (mono-stationary) and Feinberg (multi-stationary) shows a strong similarity, justifying the study of the Feinberg model, which was a modified Lee model constructed to study an important network property called "concordance". It also challenge the absoluteness of discrimination of the models into mono-stationarity versus multi-stationarity, which is a main result of Maclean et al. (PNAS USA 2015). Second, the CSEN analysis provides evidence supporting a strong similarity between the Schmitz and MacLean models, as indicated by the "proximate equivalence" that we have identified. Third, the analysis underscores the absence of a comparable relationship between the Feinberg and MacLean models, highlighting distinctive differences between the two. Thus, our approach could be a useful tool to compare mathematical models of the same or closely related systems. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Competition for binding targets results in paradoxical effects for simultaneous activator and repressor action -- Extended Version

Authors: M. Ali Al-Radhawi, Krishna Manoj, Dhruv D. Jatkar, Alon Duvall, Domitilla Del Vecchio, Eduardo D. Sontag

Abstract: In the context of epigenetic transformations in cancer metastasis, a puzzling effect was recently discovered, in which the elimination (knock-out) of an activating regulatory element leads to increased (rather than decreased) activity of the element being regulated. It has been postulated that this paradoxical behavior can be explained by activating and repressing transcription factors competing f… ▽ More In the context of epigenetic transformations in cancer metastasis, a puzzling effect was recently discovered, in which the elimination (knock-out) of an activating regulatory element leads to increased (rather than decreased) activity of the element being regulated. It has been postulated that this paradoxical behavior can be explained by activating and repressing transcription factors competing for binding to other possible targets. It is very difficult to prove this hypothesis in mammalian cells, due to the large number of potential players and the complexity of endogenous intracellular regulatory networks. Instead, this paper analyzes this issue through an analogous synthetic biology construct which aims to reproduce the paradoxical behavior using standard bacterial gene expression networks. The paper first reviews the motivating cancer biology work, and then describes a proposed synthetic construct. A mathematical model is formulated, and basic properties of uniqueness of steady states and convergence to equilibria are established, as well as an identification of parameter regimes which should lead to observing such paradoxical phenomena (more activator leads to less activity at steady state). A proof is also given to show that this is a steady-state property, and for initial transients the phenomenon will not be observed. This work adds to the general line of work of resource competition in synthetic circuits. △ Less

Submitted 21 March, 2024; originally announced March 2024.

Comments: 14 pages, 8 figures

A necessary condition for non-monotonic dose response, with an application to a kinetic proofreading model -- Extended version

Authors: Polly Y. Yu, Eduardo D. Sontag

Abstract: Steady state non-monotonic ("biphasic") dose responses are often observed in experimental biology, which raises the control-theoretic question of identifying which possible mechanisms might underlie such behaviors. It is well known that the presence of an incoherent feedforward loop (IFFL) in a network may give rise to a non-monotonic response. It has been conjectured that this condition is also n… ▽ More Steady state non-monotonic ("biphasic") dose responses are often observed in experimental biology, which raises the control-theoretic question of identifying which possible mechanisms might underlie such behaviors. It is well known that the presence of an incoherent feedforward loop (IFFL) in a network may give rise to a non-monotonic response. It has been conjectured that this condition is also necessary, i.e. that a non-monotonic response implies the existence of an IFFL. In this paper, we show that this conjecture is false, and in the process prove a weaker version: that either an IFFL must exist or both a positive loop and a negative feedback loop must exist. Towards this aim, we give necessary and sufficient conditions for when minors of a symbolic matrix have mixed signs. Finally, we study in full generality when a model of immune T-cell activation could exhibit a steady state non-monotonic dose response. △ Less

Submitted 18 April, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: Appendix included

Beyond Sight: Probing Alignment Between Image Models and Blind V1

Authors: Jacob Granley, Galen Pogoncheff, Alfonso Rodil, Leili Soo, Lily Marie Turkstra, Lucas Gil Nadolskis, Arantxa Alfaro Saez, Cristina Soto Sanchez, Eduardo Fernandez Jover, Michael Beyeler

Abstract: Neural activity in the visual cortex of blind humans persists in the absence of visual stimuli. However, little is known about the preservation of visual representation capacity in these cortical regions, which could have significant implications for neural interfaces such as visual prostheses. In this work, we present a series of analyses on the shared representations between evoked neural activi… ▽ More Neural activity in the visual cortex of blind humans persists in the absence of visual stimuli. However, little is known about the preservation of visual representation capacity in these cortical regions, which could have significant implications for neural interfaces such as visual prostheses. In this work, we present a series of analyses on the shared representations between evoked neural activity in the primary visual cortex (V1) of a blind human with an intracortical visual prosthesis, and latent visual representations computed in deep neural networks (DNNs). In the absence of natural visual input, we examine two alternative forms of inducing neural activity: electrical stimulation and mental imagery. We first quantitatively demonstrate that latent DNN activations are aligned with neural activity measured in blind V1. On average, DNNs with higher ImageNet accuracy or higher sighted primate neural predictivity are more predictive of blind V1 activity. We further probe blind V1 alignment in ResNet-50 and propose a proof-of-concept approach towards interpretability of blind V1 neurons. The results of these studies suggest the presence of some form of natural visual processing in blind V1 during electrically evoked visual perception and present unique directions in mechanistically understanding and interfacing with blind V1. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted preprint version

Leveraging The Finite States of Emotion Processing to Study Late-Life Mental Health

Authors: Yuanzhe Huang, Saurab Faruque, Minjie Wu, Akiko Mizuno, Eduardo Diniz, Shaolin Yang, George Dewitt Stetten, Noah Schweitzer, Hecheng Jin, Linghai Wang, Howard J. Aizenstein

Abstract: Traditional approaches in mental health research apply General Linear Models (GLM) to describe the longitudinal dynamics of observed psycho-behavioral measurements (questionnaire summary scores). Similarly, GLMs are also applied to characterize relationships between neurobiological measurements (regional fMRI signals) and perceptual stimuli or other regional signals. While these methods are useful… ▽ More Traditional approaches in mental health research apply General Linear Models (GLM) to describe the longitudinal dynamics of observed psycho-behavioral measurements (questionnaire summary scores). Similarly, GLMs are also applied to characterize relationships between neurobiological measurements (regional fMRI signals) and perceptual stimuli or other regional signals. While these methods are useful for exploring linear correlations among the isolated signals of those constructs (i.e., summary scores or fMRI signals), these classical frameworks fall short in providing insights into the comprehensive system-level dynamics underlying observable changes. Hidden Markov Models (HMM) are a statistical model that enable us to describe the sequential relations among multiple observable constructs, and when applied through the lens of Finite State Automata (FSA), can provide a more integrated and intuitive framework for modeling and understanding the underlying controller (the prescription for how to respond to inputs) that fundamentally defines any system, as opposed to linearly correlating output signals produced by the controller. We present a simple and intuitive HMM processing pipeline vcHMM (See Preliminary Data) that highlights FSA theory and is applicable for both behavioral analysis of questionnaire data and fMRI data. HMMs offer theoretic promise as they are computationally equivalent to the FSA, the control processor of a Turing Machine (TM) The dynamic programming Viterbi algorithm is used to leverage the HMM model. It efficiently identifies the most likely sequence of hidden states. The vcHMM pipeline leverages this grammar to understand how behavior and neural activity relate to depression. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Few-shot genes selection: subset of PAM50 genes for breast cancer subtypes classification

Authors: Leandro Y. S. Okimoto, Rayol Mendonca-Neto, Fabíola G. Nakamura, Eduardo F. Nakamura, David Fenyö, Claudio T. Silva

Abstract: Background: In recent years, researchers have made significant strides in understanding the heterogeneity of breast cancer and its various subtypes. However, the wealth of genomic and proteomic data available today necessitates efficient frameworks, instruments, and computational tools for meaningful analysis. Despite its success as a prognostic tool, the PAM50 gene signature's reliance on many ge… ▽ More Background: In recent years, researchers have made significant strides in understanding the heterogeneity of breast cancer and its various subtypes. However, the wealth of genomic and proteomic data available today necessitates efficient frameworks, instruments, and computational tools for meaningful analysis. Despite its success as a prognostic tool, the PAM50 gene signature's reliance on many genes presents challenges in terms of cost and complexity. Consequently, there is a need for more efficient methods to classify breast cancer subtypes using a reduced gene set accurately. Results: This study explores the potential of achieving precise breast cancer subtype categorization using a reduced gene set derived from the PAM50 gene signature. By employing a "Few-Shot Genes Selection" method, we randomly select smaller subsets from PAM50 and evaluate their performance using metrics and a linear model, specifically the Support Vector Machine (SVM) classifier. In addition, we aim to assess whether a more compact gene set can maintain performance while simplifying the classification process. Our findings demonstrate that certain reduced gene subsets can perform comparable or superior to the full PAM50 gene signature. Conclusions: The identified gene subsets, with 36 genes, have the potential to contribute to the development of more cost-effective and streamlined diagnostic tools in breast cancer research and clinical settings. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Journal ref: BMC Bioinformatics 25, 92 (2024)

The impact of high frequency-based stability on the onset of action potentials in neuron models

Authors: Eduardo Cerpa, Nathaly Corrales, Matías Courdurier, Leonel E. Medina, Esteban Paduro

Abstract: This paper studies the phenomenon of conduction block in model neurons using high-frequency biphasic stimulation (HFBS). The focus is investigating the triggering of undesired onset action potentials when the HFBS is turned on. The approach analyzes the transient behavior of an averaged system corresponding to the FitzHugh-Nagumo neuron model using Lyapunov and quasi-static methods. The first resu… ▽ More This paper studies the phenomenon of conduction block in model neurons using high-frequency biphasic stimulation (HFBS). The focus is investigating the triggering of undesired onset action potentials when the HFBS is turned on. The approach analyzes the transient behavior of an averaged system corresponding to the FitzHugh-Nagumo neuron model using Lyapunov and quasi-static methods. The first result provides a more comprehensive understanding of the onset activation through a mathematical proof of how to avoid it using a ramp in the amplitude of the oscillatory source. The second result tests the response of the blocked system to a piecewise linear stimulus, providing a quantitative description of how the HFBS strength translates into conduction block robustness. The results of this work can provide insights for the design of electrical neurostimulation therapies. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 32 pages, 3 figures

MSC Class: 37N25; 92-10

Targeting of chondrocyte plasticity via connexin43 modulation attenuates cellular senescence and fosters a pro-regenerative environment in osteoarthritis

Authors: Marta Varela-Eirin, Adrian Varela-Vazquez, Amanda Guitian-Caamano, Carlos Luis Paino, Virginia Mato, Raquel Largo, Trond Aasen, Arantxa Tabernero, Eduardo Fonseca, Mustapha Kandouz, Jose Ramon Caeiro, Alfonso Blanco, Maria D. Mayan

Abstract: Osteoarthritis (OA), a chronic disease characterized by articular cartilage degeneration, is a leading cause of disability and pain worldwide. In OA, chondrocytes in cartilage undergo phenotypic changes and senescence, restricting cartilage regeneration and favouring disease progression. Similar to other wound-healing disorders, chondrocytes from OA patients show a chronic increase in the gap junc… ▽ More Osteoarthritis (OA), a chronic disease characterized by articular cartilage degeneration, is a leading cause of disability and pain worldwide. In OA, chondrocytes in cartilage undergo phenotypic changes and senescence, restricting cartilage regeneration and favouring disease progression. Similar to other wound-healing disorders, chondrocytes from OA patients show a chronic increase in the gap junction channel protein connexin43 (Cx43), which regulates signal transduction through the exchange of elements or recruitment/release of signalling factors. Although immature or stem-like cells are present in cartilage from OA patients, their origin and role in disease progression are unknown. In this study, we found that Cx43 acts as a positive regulator of chondrocyte-mesenchymal transition. Downregulation of either Cx43 by CRISPR/Cas9 or Cx43-mediated gap junctional intercellular communication (GJIC) by carbenoxolone treatment triggered rediferentiation of osteoarthritic chondrocytes into a more differentiated state, associated with decreased synthesis of MMPs and proinflammatory factors, and reduced senescence. We have identified causal Cx43-sensitive circuit in chondrocytes that regulates dedifferentiation, redifferentiation and senescence. We propose that chondrocytes undergo chondrocyte-mesenchymal transition where increased Cx43-mediated GJIC during OA facilitates Twist-1 nuclear translocation as a novel mechanism involved in OA progression. These findings support the use of Cx43 as an appropriate therapeutic target to halt OA progression and to promote cartilage regeneration. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 16 pages, 5 figures

Forecasting insect abundance using time series embedding and machine learning

Authors: Gabriel R. Palma, Rodrigo F. Mello, Wesley A. C. Godoy, Eduardo Engel, Douglas Lau, Charles Markham, Rafael A. Moral

Abstract: Implementing insect monitoring systems provides an excellent opportunity to create accurate interventions for insect control. However, selecting the appropriate time for an intervention is still an open question due to the inherent difficulty of implementing on-site monitoring in real-time. This decision is even more critical with insect species that can abruptly increase population size. A possib… ▽ More Implementing insect monitoring systems provides an excellent opportunity to create accurate interventions for insect control. However, selecting the appropriate time for an intervention is still an open question due to the inherent difficulty of implementing on-site monitoring in real-time. This decision is even more critical with insect species that can abruptly increase population size. A possible solution to enhance decision-making is to apply forecasting methods to predict insect abundance. However, another layer of complexity is added when other covariates are considered in the forecasting, such as climate time series collected along the monitoring system. Multiple possible combinations of climate time series and their lags can be used to build a forecasting method. Therefore, this research paper proposes a new approach to address this problem by combining statistics, machine learning, and time series embedding. We used two datasets containing a time series of aphids and climate data collected weekly in Coxilha and Passo Fundo municipalities in Southern Brazil for eight years. We conduct a simulation study based on a probabilistic autoregressive model with exogenous time series based on Poisson and negative binomial distributions to check the influence of incorporating climate time series on the performance of our approach. We pre-processed the data using our newly proposed approach and more straightforward approaches commonly used to train machine learning algorithms in time series problems. We evaluate the performance of the selected machine algorithms by looking at the Root Mean Squared Error obtained using one-step-ahead forecasting. Based on Random Forests, Lasso-regularised linear regression, and LightGBM regression algorithms, our novel approach yields competitive forecasts while automatically selecting insect abundances, climate time series and their lags to aid forecasting. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 27 pages

Exploring the Relationship between Fractional Hill Coefficient, Intermediate Processes and Mesoscopic Fluctuations

Authors: Manuel Eduardo Hernández-García, Jorge Velázquez-Castro

Abstract: The Hill function is relevant for describing enzyme binding and other processes in gene regulatory networks. Despite its theoretical foundation, it is often used empirically as a useful fitting function. Theoretical predictions suggest that the Hill coefficient should be an integer; however, it is often assigned a fractional value. This study aims to show that the use of fractional Hill coefficien… ▽ More The Hill function is relevant for describing enzyme binding and other processes in gene regulatory networks. Despite its theoretical foundation, it is often used empirically as a useful fitting function. Theoretical predictions suggest that the Hill coefficient should be an integer; however, it is often assigned a fractional value. This study aims to show that the use of fractional Hill coefficients actually indicates the presence of intermediate processes or, in some cases, the effect of noise on these systems. The deterministic approximation of enzymatic processes leads to the derivation of the Hill function, which can be expanded around the noise magnitude to derive mesoscopic corrections. This study establishes relationships between the intermediate processes and the decimal Hill coefficient, both with and without fluctuations. This outcome contributes to a deeper understanding of the underlying processes associated with the fractional Hill coefficient, while also enabling the prediction of an effective value of the Hill coefficient from the underlying mechanism, allowing us to have a simplified description of complex systems. △ Less

Submitted 25 December, 2023; originally announced December 2023.

Comments: 12 pages, 9 figures

MSC Class: 92-11 ACM Class: G.3; I.6

Scalable Knowledge Graph Construction and Inference on Human Genome Variants

Authors: Shivika Prasanna, Deepthi Rao, Eduardo Simoes, Praveen Rao

Abstract: Real-world knowledge can be represented as a graph consisting of entities and relationships between the entities. The need for efficient and scalable solutions arises when dealing with vast genomic data, like RNA-sequencing. Knowledge graphs offer a powerful approach for various tasks in such large-scale genomic data, such as analysis and inference. In this work, variant-level information extracte… ▽ More Real-world knowledge can be represented as a graph consisting of entities and relationships between the entities. The need for efficient and scalable solutions arises when dealing with vast genomic data, like RNA-sequencing. Knowledge graphs offer a powerful approach for various tasks in such large-scale genomic data, such as analysis and inference. In this work, variant-level information extracted from the RNA-sequences of vaccine-naïve COVID-19 patients have been represented as a unified, large knowledge graph. Variant call format (VCF) files containing the variant-level information were annotated to include further information for each variant. The data records in the annotated files were then converted to Resource Description Framework (RDF) triples. Each VCF file obtained had an associated CADD scores file that contained the raw and Phred-scaled scores for each variant. An ontology was defined for the VCF and CADD scores files. Using this ontology and the extracted information, a large, scalable knowledge graph was created. Available graph storage was then leveraged to query and create datasets for further downstream tasks. We also present a case study using the knowledge graph and perform a classification task using graph machine learning. We also draw comparisons between different Graph Neural Networks (GNNs) for the case study. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Pulsed interactions unify reaction-diffusion and spatial nonlocal models for biological pattern formation

Authors: Eduardo H. Colombo, Ricardo Martinez-Garcia, Justin M. Calabrese, Cristóbal López, Emilio Hernández-García

Abstract: The emergence of a spatially-organized population distribution depends on the dynamics of the population and mediators of interaction (activators and inhibitors). Two broad classes of models have been used to investigate when and how self-organization is triggered, namely, reaction-diffusion and spatially nonlocal models. Nevertheless, these models implicitly assume smooth propagation scenarios, n… ▽ More The emergence of a spatially-organized population distribution depends on the dynamics of the population and mediators of interaction (activators and inhibitors). Two broad classes of models have been used to investigate when and how self-organization is triggered, namely, reaction-diffusion and spatially nonlocal models. Nevertheless, these models implicitly assume smooth propagation scenarios, neglecting that individuals many times interact by exchanging short and abrupt pulses of the mediating substance. A recently proposed framework advances in the direction of properly accounting for these short-scale fluctuations by applying a coarse-graining procedure on the pulse dynamics. In this paper, we generalize the coarse-graining procedure and apply the extended formalism to new scenarios in which mediators influence individuals' reproductive success or their motility. We show that, in the slow- and fast-mediator limits, pulsed interactions recover, respectively, the reaction-diffusion and nonlocal models, providing a mechanistic connection between them. Furthermore, at each limit, the spatial stability condition is qualitatively different, leading to a timescale-induced transition where spatial patterns emerge as mediator dynamics becomes sufficiently fast. △ Less

Submitted 21 March, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Journal ref: Journal of Statistical Mechanics: Theory and Experiment 2024, 034001 (2024)

Quantifying Hierarchical Selection

Authors: Hardik Rajpal, Clem von Stengel, Pedro A. M. Mediano, Fernando E. Rosas, Eduardo Viegas, Pablo A. Marquet, Henrik J. Jensen

Abstract: At what level does selective pressure effectively act? When considering the reproductive dynamics of interacting and mutating agents, it has long been debated whether selection is better understood by focusing on the individual or if hierarchical selection emerges as a consequence of joint adaptation. Despite longstanding efforts in theoretical ecology there is still no consensus on this fundament… ▽ More At what level does selective pressure effectively act? When considering the reproductive dynamics of interacting and mutating agents, it has long been debated whether selection is better understood by focusing on the individual or if hierarchical selection emerges as a consequence of joint adaptation. Despite longstanding efforts in theoretical ecology there is still no consensus on this fundamental issue, most likely due to the difficulty in obtaining adequate data spanning sufficient number of generations and the lack of adequate tools to quantify the effect of hierarchical selection. Here we capitalise on recent advances in information-theoretic data analysis to advance this state of affairs by investigating the emergence of high-order structures -- such as groups of species -- in the collective dynamics of the Tangled Nature model of evolutionary ecology. Our results show that evolutionary dynamics can lead to clusters of species that act as a selective group, that acquire information-theoretic agency. Overall, our findings provide quantitative evidence supporting the relevance of high-order structures in evolutionary ecology, which can emerge even from relatively simple processes of adaptation and selection. △ Less

Submitted 1 November, 2023; v1 submitted 31 October, 2023; originally announced October 2023.

Improving Molecular Properties Prediction Through Latent Space Fusion

Authors: Eduardo Soares, Akihiro Kishimoto, Emilio Vital Brazil, Seiji Takeda, Hiroshi Kajino, Renato Cerqueira

Abstract: Pre-trained Language Models have emerged as promising tools for predicting molecular properties, yet their development is in its early stages, necessitating further research to enhance their efficacy and address challenges such as generalization and sample efficiency. In this paper, we present a multi-view approach that combines latent spaces derived from state-of-the-art chemical models. Our appr… ▽ More Pre-trained Language Models have emerged as promising tools for predicting molecular properties, yet their development is in its early stages, necessitating further research to enhance their efficacy and address challenges such as generalization and sample efficiency. In this paper, we present a multi-view approach that combines latent spaces derived from state-of-the-art chemical models. Our approach relies on two pivotal elements: the embeddings derived from MHG-GNN, which represent molecular structures as graphs, and MoLFormer embeddings rooted in chemical language. The attention mechanism of MoLFormer is able to identify relations between two atoms even when their distance is far apart, while the GNN of MHG-GNN can more precisely capture relations among multiple atoms closely located. In this work, we demonstrate the superior performance of our proposed multi-view approach compared to existing state-of-the-art methods, including MoLFormer-XL, which was trained on 1.1 billion molecules, particularly in intricate tasks such as predicting clinical trial drug toxicity and inhibiting HIV replication. We assessed our approach using six benchmark datasets from MoleculeNet, where it outperformed competitors in five of them. Our study highlights the potential of latent space fusion and feature integration for advancing molecular property prediction. In this work, we use small versions of MHG-GNN and MoLFormer, which opens up an opportunity for further improvement when our approach uses a larger-scale dataset. △ Less

Submitted 20 October, 2023; originally announced October 2023.

Comments: 8 Pages, 4 Figures - Submited to the AI4Science Workshop - Neurips 2023

Network transformation-based analysis of biochemical systems

Authors: Dylan Antonio Talabis, Eduardo Mendoza

Abstract: A dynamical system obtains a wide variety of kinetic realizations, which is advantageous for the analysis of biochemical systems. A reaction network, derived from a dynamical system, may or may not possess some properties needed for a thorough analysis. We improve and extend the work of M. Johnston \cite{JOHN2014} and Hong et al. \cite{HONG2023} on network translations to network transformations,… ▽ More A dynamical system obtains a wide variety of kinetic realizations, which is advantageous for the analysis of biochemical systems. A reaction network, derived from a dynamical system, may or may not possess some properties needed for a thorough analysis. We improve and extend the work of M. Johnston \cite{JOHN2014} and Hong et al. \cite{HONG2023} on network translations to network transformations, where the network is modified while preserving the dynamical system. These transformations can shrink, extend, or retain the stoichiometric subspace. Here, we show that positive dependent network can be translated to a weakly reversible network. Using the kinetic realizations of (1) calcium signaling in the olfactory system and (2) metabolic insulin signaling, we demonstrate the benefits of transformed systems with positive deficiency for analyzing biochemical systems. Furthermore, we present an algorithm for a network transformation of a weakly reversible non-complex factorizable kinetic (NFK) system to a weakly reversible complex factorizable kinetic (CFK) system, thereby enhancing the Subspace Coincidence Theorem for NFK systems of Nazareno et al. \cite{NAZA2019}. Finally, using the transformed kinetic realization of monolignol biosynthesis in \textit{Populus xylem}, we study the structural and kinetic properties of transformed systems,including the invariance of concordance and variation of injectivity and mono-/multi-stationarity under network transformation. △ Less

Submitted 30 April, 2024; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 40 pages. arXiv admin note: text overlap with arXiv:1109.2923 by other authors

MSC Class: 92C45; 65H10; 80A30

What it takes to solve the Origin(s) of Life: An integrated review of techniques

Authors: OoLEN, Silke Asche, Carla Bautista, David Boulesteix, Alexandre Champagne-Ruel, Cole Mathis, Omer Markovitch, Zhen Peng, Alyssa Adams, Avinash Vicholous Dass, Arnaud Buch, Eloi Camprubi, Enrico Sandro Colizzi, Stephanie Colón-Santos, Hannah Dromiack, Valentina Erastova, Amanda Garcia, Ghjuvan Grimaud, Aaron Halpern, Stuart A Harrison, Seán F. Jordan, Tony Z Jia, Amit Kahana, Artemy Kolchinsky, Odin Moron-Garcia , et al. (13 additional authors not shown)

Abstract: Understanding the origin(s) of life (OoL) is a fundamental challenge for science in the 21st century. Research on OoL spans many disciplines, including chemistry, physics, biology, planetary sciences, computer science, mathematics and philosophy. The sheer number of different scientific perspectives relevant to the problem has resulted in the coexistence of diverse tools, techniques, data, and sof… ▽ More Understanding the origin(s) of life (OoL) is a fundamental challenge for science in the 21st century. Research on OoL spans many disciplines, including chemistry, physics, biology, planetary sciences, computer science, mathematics and philosophy. The sheer number of different scientific perspectives relevant to the problem has resulted in the coexistence of diverse tools, techniques, data, and software in OoL studies. This has made communication between the disciplines relevant to the OoL extremely difficult because the interpretation of data, analyses, or standards of evidence can vary dramatically. Here, we hope to bridge this wide field of study by providing common ground via the consolidation of tools and techniques rather than positing a unifying view on how life emerges. We review the common tools and techniques that have been used significantly in OoL studies in recent years. In particular, we aim to identify which information is most relevant for comparing and integrating the results of experimental analyses into mathematical and computational models. This review aims to provide a baseline expectation and understanding of technical aspects of origins research, rather than being a primer on any particular topic. As such, it spans broadly -- from analytical chemistry to mathematical models -- and highlights areas of future work that will benefit from a multidisciplinary approach to tackling the mystery of life's origin. Ultimately, we hope to empower a new generation of OoL scientists by reviewing how they can investigate life's origin, rather than dictating how to think about the problem. △ Less

Submitted 24 August, 2023; v1 submitted 22 August, 2023; originally announced August 2023.

TrajPy: empowering feature engineering for trajectory analysis across domains

Authors: Maurício Moreira-Soares, Eduardo Mossmann, Rui D. M. Travasso, José Rafael Bordin

Abstract: Trajectories, sequentially measured quantities that form a path, are an important presence in many different fields, from hadronic beams in physics to electrocardiograms in medicine. Trajectory anal-ysis requires the quantification and classification of curves either using statistical descriptors or physics-based features. To date, there is no extensive and user-friendly package for trajectory ana… ▽ More Trajectories, sequentially measured quantities that form a path, are an important presence in many different fields, from hadronic beams in physics to electrocardiograms in medicine. Trajectory anal-ysis requires the quantification and classification of curves either using statistical descriptors or physics-based features. To date, there is no extensive and user-friendly package for trajectory anal-ysis available, despite its importance and potential application across domains. We developed a free open-source python package named TrajPy as a complementary tool to empower trajectory analysis. The package showcases a friendly graphic user interface and provides a set of physical descriptors that help characterizing these intricate structures. In combina-tion with image analysis, it was already successfully applied to the study of mitochondrial motility in neuroblastoma cell lines and to the analysis of in silico models for cell migration. The TrajPy package was developed in Python 3 and released under the GNU GPL-3 license. Easy installation is available through PyPi and the development source code can be found in the repository https://github.com/ocbe-uio/TrajPy/. The package release is automatically archived under the DOI 10.5281/zenodo.3656044. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 4 pages, 1 figure

To reverse engineer an entire nervous system

Authors: Gal Haspel, Edward S Boyden, Jeffrey Brown, George Church, Netta Cohen, Christopher Fang-Yen, Steven Flavell, Miriam B Goodman, Anne C Hart, Oliver Hobert, Eduardo J Izquierdo, Konstantinos Kagias, Shawn Lockery, Yangning Lu, Adam Marblestone, Jordan Matelsky, Hanspeter Pfister, Horacio G Rotstein, Monika Scholz, Eli Shlizerman, Quilee Simeon, Michael A Skuhersky, Vineet Tiruvadi, Vivek Venkatachalam, Guangyu Robert Yang , et al. (3 additional authors not shown)

Abstract: A primary goal of neuroscience is to understand how nervous systems, or assemblies of neural circuits, generate and control behavior. Testing and refining our theories of neural control would be greatly facilitated if we could reliably simulate an entire nervous system so we could replicate the brain dynamics in response to any stimuli and different contexts. More fundamentally, reconstructing or… ▽ More A primary goal of neuroscience is to understand how nervous systems, or assemblies of neural circuits, generate and control behavior. Testing and refining our theories of neural control would be greatly facilitated if we could reliably simulate an entire nervous system so we could replicate the brain dynamics in response to any stimuli and different contexts. More fundamentally, reconstructing or modeling a system is an important milestone in understanding it, and so, simulating an entire nervous system is in itself one of the goals, indeed dreams, of systems neuroscience. To do so requires us to identify how each neuron's output depends on its inputs, within some nervous system. This deconstruction, understanding function from input-output pairs, falls into the realm of reverse engineering. Current efforts at reverse engineering the brain focus on the mammalian nervous system, but these brains are complex, allowing only recordings of tiny subsystems. Here we argue that the time is ripe to embark on a concerted effort to reverse engineer a smaller system and that the nematode C. elegans is the ideal candidate system. In particular, the established and growing toolkit of optophysiology techniques can non-invasively capture and control each neuron's activity and scale to hundreds of thousands of experiments, across a large population of animals. Data across populations and behaviors can be combined because across individuals neuronal identities are largely conserved in form and function. Modern machine-learning-based model training should then enable a simulation of C. elegans' impressive breadth of brain states and behaviors. The ability to reverse engineer an entire nervous system will benefit systems neuroscience as well as the design of artificial intelligence systems, enabling fundamental insights as well as new approaches for investigations of progressively larger nervous systems. △ Less

Submitted 9 December, 2023; v1 submitted 12 August, 2023; originally announced August 2023.

Comments: 23 pages, 2 figures, opinion paper

On structural contraction of biological interaction networks

Authors: M. Ali Al-Radhawi, David Angeli, Eduardo Sontag

Abstract: In previous work, we have developed an approach for characterizing the long-term dynamics of classes of Biological Interaction Networks (BINs), based on "rate-dependent Lyapunov functions". In this work, we show that stronger notions of convergence can be established by proving structural contractivity with respect to non-standard norms. We illustrate our theory with examples from signaling pathwa… ▽ More In previous work, we have developed an approach for characterizing the long-term dynamics of classes of Biological Interaction Networks (BINs), based on "rate-dependent Lyapunov functions". In this work, we show that stronger notions of convergence can be established by proving structural contractivity with respect to non-standard norms. We illustrate our theory with examples from signaling pathways. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Sparse electrophysiological source imaging predicts aging-related gait speed slowing

Authors: Vega-Hernandez, Mayrim, Galan-Garcia, Lidice, Perez-Hidalgo-Gato, Jhoanna, Ontivero-Ortega, Marlis, Garcia-Agustin, Daysi, Garcia-Reyes, Ronaldo, Bosch-Bayard, Jorge, Marinazzo, Daniele, Martinez-Montes, Eduardo, Valdes-Sosa, Pedro A

Abstract: Objective: We seek stable Electrophysiological Source Imaging (ESI) biomarkers associated with Gait Speed (GS) as a measure of functional decline. Towards this end we determine the predictive value of ESI activation and connectivity patterns of resting-state EEG Theta rhythm on physical performance decline measured by a slowing GS in aging individuals. Methods: As potential biomarkers related to G… ▽ More Objective: We seek stable Electrophysiological Source Imaging (ESI) biomarkers associated with Gait Speed (GS) as a measure of functional decline. Towards this end we determine the predictive value of ESI activation and connectivity patterns of resting-state EEG Theta rhythm on physical performance decline measured by a slowing GS in aging individuals. Methods: As potential biomarkers related to GS changes, we estimate ESI using flexible sparse/smooth/non-negative models (NN-SLASSO), from which activation ESI (aESI) and connectivity ESI (cESI) features are selected using the Stable Sparse Classifier method. Results and Conclusions: Novel sparse aESI models outperformed traditional methods such as the LORETA family. The models combining aESI and cESI features improved the predictability of GS changes. Selected biomarkers from activation/connectivity patterns were localized to orbitofrontal and temporal cortical regions. Significance: The proposed methodology contributes to understanding the activation and connectivity of ESI complex patterns related to GS, providing potential biomarker features for GS slowing. Given the known relationship between GS decline and cognitive impairment, this preliminary work suggests it might be applied to other, more complex measures of healthy and pathological aging. Importantly, it might allow an ESI-based evaluation of rehabilitation programs. △ Less

Submitted 22 August, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comparative Analysis of Kinetic Realizations of Insulin Signaling

Authors: Patrick Vincent N. Lubenia, Eduardo R. Mendoza, Angelyn R. Lao

Abstract: Several studies have developed dynamical models to understand the underlying mechanisms of insulin signaling, a signaling cascade that leads to the translocation of glucose, the human body's main source of energy. Fortunately, reaction network analysis allows us to extract properties of dynamical systems without depending on their model parameter values. This study focuses on the comparison of ins… ▽ More Several studies have developed dynamical models to understand the underlying mechanisms of insulin signaling, a signaling cascade that leads to the translocation of glucose, the human body's main source of energy. Fortunately, reaction network analysis allows us to extract properties of dynamical systems without depending on their model parameter values. This study focuses on the comparison of insulin signaling in healthy state (INSMS or INSulin Metabolic Signaling) and in type 2 diabetes (INRES or INsulin RESistance) using reaction network analysis. The analysis uses network decomposition to identify the different subsystems involved in insulin signaling (e.g., insulin receptor binding and recycling, GLUT4 translocation, and ERK signaling pathway, among others). Furthermore, results show that INSMS and INRES are similar with respect to some network, structo-kinetic, and kinetic properties. Their differences, however, provide insights into what happens when insulin resistance occurs. First, the variation in the number of species involved in INSMS and INRES suggests that when irregularities occur in the insulin signaling pathway, other complexes (and, hence, other processes) get involved, characterizing insulin resistance. Second, the loss of concordance exhibited by INRES suggests less restrictive interplay between the species involved in insulin signaling, leading to unusual activities in the signaling cascade. Lastly, GLUT4 losing its absolute concentration robustness in INRES may signify that the transporter has lost its reliability in shuttling glucose to the cell, inhibiting efficient cellular energy production. This study also suggests possible applications of the equilibria parametrization and network decomposition, resulting from the analysis, to potentially establish absolute concentration robustness in a species. △ Less

Submitted 23 October, 2023; v1 submitted 7 July, 2023; originally announced July 2023.

Comments: 32 pages, 1 figure

Corrected Hill Function in Stochastic Gene Regulatory Networks

Authors: Manuel Eduardo Hernández-García, Jorge Velázquez-Castro

Abstract: Describing reaction rates in stochastic bio-circuits is commonly done by directly introducing the deterministically deduced Hill function into the master equation. However, when fluctuations in enzymatic reaction rates are not neglectable, the Hill function must be derived, considering all the involved stochastic reactions. In this work, we derived the stochastic version of the Hill function from… ▽ More Describing reaction rates in stochastic bio-circuits is commonly done by directly introducing the deterministically deduced Hill function into the master equation. However, when fluctuations in enzymatic reaction rates are not neglectable, the Hill function must be derived, considering all the involved stochastic reactions. In this work, we derived the stochastic version of the Hill function from the master equation of the complete set of reactions that, in the macroscopic limit, lead to the Hill function reaction rate. We performed a series expansion around the average values of the concentrations, which allowed us to find corrections for the deterministic Hill function. This process allowed us to quantify the fluctuations of enzymatic reactions. We found that the underlying variability in propensity rates of gene regulatory networks has an important non-linear effect that reduces the intrinsic fluctuations of the mRNA and protein concentrations. △ Less

Submitted 3 November, 2023; v1 submitted 6 July, 2023; originally announced July 2023.

Comments: 23 pages, 14 figures

MSC Class: 92-11 ACM Class: I.6

Beyond Chemical Language: A Multimodal Approach to Enhance Molecular Property Prediction

Authors: Eduardo Soares, Emilio Vital Brazil, Karen Fiorela Aquino Gutierrez, Renato Cerqueira, Dan Sanders, Kristin Schmidt, Dmitry Zubarev

Abstract: We present a novel multimodal language model approach for predicting molecular properties by combining chemical language representation with physicochemical features. Our approach, MULTIMODAL-MOLFORMER, utilizes a causal multistage feature selection method that identifies physicochemical features based on their direct causal effect on a specific target property. These causal features are then inte… ▽ More We present a novel multimodal language model approach for predicting molecular properties by combining chemical language representation with physicochemical features. Our approach, MULTIMODAL-MOLFORMER, utilizes a causal multistage feature selection method that identifies physicochemical features based on their direct causal effect on a specific target property. These causal features are then integrated with the vector space generated by molecular embeddings from MOLFORMER. In particular, we employ Mordred descriptors as physicochemical features and identify the Markov blanket of the target property, which theoretically contains the most relevant features for accurate prediction. Our results demonstrate a superior performance of our proposed approach compared to existing state-of-the-art algorithms, including the chemical language-based MOLFORMER and graph neural networks, in predicting complex tasks such as biodegradability and PFAS toxicity estimation. Moreover, we demonstrate the effectiveness of our feature selection method in reducing the dimensionality of the Mordred feature space while maintaining or improving the model's performance. Our approach opens up promising avenues for future research in molecular property prediction by harnessing the synergistic potential of both chemical language and physicochemical features, leading to enhanced performance and advancements in the field. △ Less

Submitted 22 June, 2023; originally announced June 2023.

Comments: 14 pages, 6 Figures, 5 tables. Submited to NEURIPS 2023, Under review

ACM Class: J.2; I.2.1

Estimating the severity of COVID-19 Omicron Variant in USA, India, Brazil, France, Germany and UK

Authors: Tulio Eduardo Rodrigues, Otaviano Helene

Abstract: This work evaluates the severity of COVID-19 Omicron variant in terms of the case-fatality-rates (CFR) with respective uncertainty intervals via a simultaneous fitting of confirmed cases and deaths in the USA, India, Brazil, France, Germany and United Kingdom. The CFRs were calculated combining Monte Carlo simulations and analytical methods based on Gompertz functions under the framework of the Le… ▽ More This work evaluates the severity of COVID-19 Omicron variant in terms of the case-fatality-rates (CFR) with respective uncertainty intervals via a simultaneous fitting of confirmed cases and deaths in the USA, India, Brazil, France, Germany and United Kingdom. The CFRs were calculated combining Monte Carlo simulations and analytical methods based on Gompertz functions under the framework of the Least Square Method and assuming that the deaths can be described by a convolution of the confirmed cases with a common gamma function to describe the case to death period. Linear backgrounds both for cases and deaths were included in the fitting to account for the contributions from other strains within the Omicron peaks. The fitting included 125 and 113 epidemiological weeks, for cases and deaths, respectively, and 64 parameters, resulting in a chi^2 of 176.5 for 174 degrees of freedom (p = 0.434). The CFRs with 95% confidence intervals for USA, India, Brazil, France, Germany and United Kingdom were 0.295 (0.154-0.436)%, 0.232 (0.134-0.331)%, 0.49 (0.27-0.71)%, 0.056 (0.028-0.084)%, 0.129 (0.074-0.184)% and 0.168 (0.107-0.229)%, respectively. The case to death period was satisfactorily described by a common gamma function with mean of 15.71 +/- 0.55 days and coefficient of variation of 0.354 +/- 0.070. The proposed calculations provided accurate and reliable information about the respective CFRs and the case to death period, the latter being consistent with previous estimates for the symptom onset to death made in the early stages of the pandemic. The CFRs thus obtained are considerably lower than previous measurements available in the literature, suggesting that the latter may have been overestimated, as the probability of deaths from other strains of the virus under the generally prominent Omicron peak, here accounted for in terms of linear backgrounds, was not considered. △ Less

Submitted 20 February, 2023; originally announced February 2023.

Comments: 13 pages, 4 figures and 2 tables

Community-developed checklists for publishing images and image analysis

Authors: Christopher Schmied, Michael Nelson, Sergiy Avilov, Gert-Jan Bakker, Cristina Bertocchi, Johanna Bischof, Ulrike Boehm, Jan Brocher, Mariana Carvalho, Catalin Chiritescu, Jana Christopher, Beth Cimini, Eduardo Conde-Sousa, Michael Ebner, Rupert Ecker, Kevin Eliceiri, Julia Fernandez-Rodriguez, Nathalie Gaudreault, Laurent Gelman, David Grunwald, Tingting Gu, Nadia Halidi, Mathias Hammer, Matthew Hartley, Marie Held , et al. (29 additional authors not shown)

Abstract: Images document scientific discoveries and are prevalent in modern biomedical research. Microscopy imaging in particular is currently undergoing rapid technological advancements. However for scientists wishing to publish the obtained images and image analyses results, there are to date no unified guidelines. Consequently, microscopy images and image data in publications may be unclear or difficult… ▽ More Images document scientific discoveries and are prevalent in modern biomedical research. Microscopy imaging in particular is currently undergoing rapid technological advancements. However for scientists wishing to publish the obtained images and image analyses results, there are to date no unified guidelines. Consequently, microscopy images and image data in publications may be unclear or difficult to interpret. Here we present community-developed checklists for preparing light microscopy images and image analysis for publications. These checklists offer authors, readers, and publishers key recommendations for image formatting and annotation, color selection, data availability, and for reporting image analysis workflows. The goal of our guidelines is to increase the clarity and reproducibility of image figures and thereby heighten the quality of microscopy data is in publications. △ Less

Submitted 14 September, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 28 pages, 8 Figures, 3 Supplmentary Figures, Manuscript, Essential recommendations for publication of microscopy image data

Subharmonic solutions for a class of predator-prey models with degenerate weights in periodic environments

Authors: Julián López-Gómez, Eduardo Muñoz-Hernández, Fabio Zanolin

Abstract: This paper deals with the existence, multiplicity, minimal complexity and global structure of the subharmonic solutions to a class of planar Hamiltonian systems with periodic coefficients, being the classical predator-prey model of V. Volterra its most paradigmatic example. By means of a topological approach based on techniques from global bifurcation theory, the first part of the paper ascertains… ▽ More This paper deals with the existence, multiplicity, minimal complexity and global structure of the subharmonic solutions to a class of planar Hamiltonian systems with periodic coefficients, being the classical predator-prey model of V. Volterra its most paradigmatic example. By means of a topological approach based on techniques from global bifurcation theory, the first part of the paper ascertains their nature, multiplicity and minimal complexity, as well as their global minimal structure, in terms of the configuration of the function coefficients in the setting of the model. The second part of the paper introduces a dynamical system approach based on the theory of topological horseshoes that permits to detect, besides subharmonic solutions, ``chaotic-type'' solutions. As a byproduct of our analysis, the simplest predator-prey prototype models in periodic environments can provoke chaotic dynamics. This cannot occur in cooperative and quasi-cooperative dynamics, as a consequence of the ordering imposed by the maximum principle. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: 72 pages, 21 figures

MSC Class: 34C25; 37B55; 37E40; 37J12

On the major role played by the curvature of intracranial aneurysms walls in determining their mechanical response, local hemodynamics, and rupture likelihood

Authors: Iago Oliveira, Philip Cardiff, Carlos Eduardo Baccin, Rafael Tatit, José Luiz Gasche

Abstract: The properties of intracranial aneurysms (IAs) walls are known to be driven by the underlying hemodynamics adjacent to the IA sac. Different pathways exist explaining the connections between hemodynamics and local tissue properties. The emergence of such theories is essential if one wishes to compute the mechanical response of a patient-specific IA wall and predict its rupture. Apart from the hemo… ▽ More The properties of intracranial aneurysms (IAs) walls are known to be driven by the underlying hemodynamics adjacent to the IA sac. Different pathways exist explaining the connections between hemodynamics and local tissue properties. The emergence of such theories is essential if one wishes to compute the mechanical response of a patient-specific IA wall and predict its rupture. Apart from the hemodynamics and tissue properties, one could assume that the mechanical response also depends on the local morphology, more specifically, the wall curvature, with larger values at highly-curved wall portions. Nonetheless, this contradicts observations of IA rupture sites more often found at the dome, where the curvature is lower. This seeming contradiction indicates a complex interaction between local hemodynamics, wall morphology, and mechanical response, which warrants further investigation. This was the main goal of this work. We accomplished this by analysing the stress and stretch fields in different regions of the wall for a sample of IAs, which have been classified based on particular local hemodynamics and local curvature. Pulsatile numerical simulations were performed using the one-way fluid-solid interaction strategy implemented in OpenFOAM (solids4foam toolbox). We found that the variable best correlated with regions of high stress and stretch was the wall curvature. Additionally, our data suggest a connection between the local curvature and local hemodynamics, indicating that the curvature is a property that could be used to assess both mechanical response and hemodynamic conditions, and, moreover, to suggest new metrics based on the curvature to predict the likelihood of rupture. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: Preprint submitted to Acta Biomaterialia, with 27 pages and 11 figures

A compartmental model for Xylella fastidiosa diseases with explicit vector seasonal dynamics

Authors: Àlex Giménez-Romero, Eduardo Moralejo, Manuel A. Matías

Abstract: The bacterium Xylella fastidiosa (Xf) is mainly transmitted by the spittlebug, Philaenus spumarius, in Europe, where it has caused significant economic damage to olive and almond trees. Understanding the factors that determine disease dynamics in pathosystems that share similarities can help design control strategies focused on minimizing transmission chains. Here we introduce a compartmental mode… ▽ More The bacterium Xylella fastidiosa (Xf) is mainly transmitted by the spittlebug, Philaenus spumarius, in Europe, where it has caused significant economic damage to olive and almond trees. Understanding the factors that determine disease dynamics in pathosystems that share similarities can help design control strategies focused on minimizing transmission chains. Here we introduce a compartmental model for Xf-caused diseases in Europe that accounts for the main relevant epidemiological processes, including the seasonal dynamics of P. spumarius. The model was confronted with epidemiological data from the two major outbreaks of Xf in Europe, the olive quick disease syndrome (OQDS) in Apulia, Italy, caused by the subspecies pauca, and the almond leaf scorch disease (ALSD) in Majorca, Spain, caused by subspecies multiplex and fastidiosa. Using a Bayesian inference framework, we show how the model successfully reproduces the general field data in both diseases. In a global sensitivity analysis, the vector-plant and plant-vector transmission rates, together with the vector removal rate, were the most influential parameters in determining the time of the infected host population peak, the incidence peak and the final number of dead hosts. We also used our model to check different vector-based control strategies, showing that a joint strategy focused on increasing the rate of vector removal while lowering the number of annual newborn vectors is optimal for disease control. △ Less

Submitted 10 December, 2022; originally announced December 2022.

Comments: 18 pages, 8 figures

Estimation of fibre architecture and scar in myocardial tissue using electrograms: an in-silico study

Authors: Konstantinos Ntagiantas, Eduardo Pignatelli, Nicholas S. Peters, Chris D. Cantwell, Rasheda A. Chowdhury, Anil A. Bharath

Abstract: Atrial Fibrillation (AF) is characterized by disorganised electrical activity in the atria and is known to be sustained by the presence of regions of fibrosis (scars) or functional cellular remodeling, both of which may lead to areas of slow conduction. Estimating the effective conductivity of the myocardium and identifying regions of abnormal propagation is therefore crucial for the effective tre… ▽ More Atrial Fibrillation (AF) is characterized by disorganised electrical activity in the atria and is known to be sustained by the presence of regions of fibrosis (scars) or functional cellular remodeling, both of which may lead to areas of slow conduction. Estimating the effective conductivity of the myocardium and identifying regions of abnormal propagation is therefore crucial for the effective treatment of AF. We hypothesise that the spatial distribution of tissue conductivity can be directly inferred from an array of concurrently acquired contact electrograms (EGMs). We generate a dataset of simulated cardiac AP propagation using randomised scar distributions and a phenomenological cardiac model and calculate contact EGMs at various positions on the field. EGMs are enriched with noise extracted from biological data acquired in the lab. A deep neural network, based on a modified U-net architecture, is trained to estimate the location of the scar and quantify conductivity of the tissue with a Jaccard index of 91%. We adapt a wavelet-based surrogate testing analysis to confirm that the inferred conductivity distribution is an accurate representation of the ground truth input to the model. We find that the root mean square error (RMSE) between the ground truth and our predictions is significantly smaller ($p_{val}<0.01$) than the RMSE between the ground truth and surrogate samples. △ Less

Submitted 21 February, 2023; v1 submitted 6 December, 2022; originally announced December 2022.

Comments: 16 pages, 15 figures

Epigenetic factor competition reshapes the EMT landscape

Authors: M. Ali Al-Radhawi, Shubham Tripathi, Yun Zhang, Eduardo D. Sontag, Herbert Levine

Abstract: The emergence of and transitions between distinct phenotypes in isogenic cells can be attributed to the intricate interplay of epigenetic marks, external signals, and gene regulatory elements. These elements include chromatin remodelers, histone modifiers, transcription factors, and regulatory RNAs. Mathematical models known as Gene Regulatory Networks (GRNs) are an increasingly important tool to… ▽ More The emergence of and transitions between distinct phenotypes in isogenic cells can be attributed to the intricate interplay of epigenetic marks, external signals, and gene regulatory elements. These elements include chromatin remodelers, histone modifiers, transcription factors, and regulatory RNAs. Mathematical models known as Gene Regulatory Networks (GRNs) are an increasingly important tool to unravel the workings of such complex networks. In such models, epigenetic factors are usually proposed to act on the chromatin regions directly involved in the expression of relevant genes. However, it has been well-established that these factors operate globally and compete with each other for targets genome-wide. Therefore, a perturbation of the activity of a regulator can redistribute epigenetic marks across the genome and modulate the levels of competing regulators. In this paper, we propose a conceptual and mathematical modeling framework that incorporates both local and global competition effects between antagonistic epigenetic regulators in addition to local transcription factors, and show the counter-intuitive consequences of such interactions. We apply our approach to recent experimental findings on the Epithelial-Mesenchymal Transition (EMT). We show that it can explain the puzzling experimental data as well provide new verifiable predictions. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Journal ref: Proc Natl Acad Sci USA, 119:e2210844119, 2022

Pattern-Based Prediction of Population Outbreaks

Authors: Gabriel R. Palma, Wesley A. C. Godoy, Eduardo Engel, Douglas Lau, Edgar Galvan, Oliver Mason, Charles Markham, Rafael A. Moral

Abstract: Insect outbreaks are biotic disturbances in forests and agroecosystems that cause economic and ecological damage. This phenomenon depends on a variety of biological and physical factors. The complexity and practical importance of the issue have made the problem of predicting outbreaks a focus of recent research. Here, we propose the Pattern-Based Prediction (PBP) method for predicting population o… ▽ More Insect outbreaks are biotic disturbances in forests and agroecosystems that cause economic and ecological damage. This phenomenon depends on a variety of biological and physical factors. The complexity and practical importance of the issue have made the problem of predicting outbreaks a focus of recent research. Here, we propose the Pattern-Based Prediction (PBP) method for predicting population outbreaks. It is based on the Alert Zone Procedure, combined with elements from machine learning. It uses information on previous time series values that precede an outbreak event as predictors of future outbreaks, which can be useful when monitoring pest species. We illustrate the methodology using simulated datasets and real time series data obtained by monitoring aphids in wheat crops in Southern Brazil. We obtained an average test accuracy of $84.6\%$ in the simulation studies implemented with stochastic models, and $95.0\%$ for predicting outbreaks using the real dataset. This shows the feasibility of the PBP method in predicting outbreaks in population dynamics. We benchmarked our results against established state-of-the-art machine learning methods, namely Support Vector Machines, Deep Neural Networks, Long Short Term Memory and Random Forests. The PBP method yielded a competitive performance, associated with higher true-positive rates in most comparisons, while being able to provide interpretability rather than being a black-box method. This is an improvement over current state-of-the-art machine learning tools, especially when being used by non-specialists, such as ecologists aiming to use a quantitative approach for pest monitoring. We provide open-source code to implement the PBP method in Python, through the \texttt{pypbp} package, which may be directly downloaded from the Python Package Index server or accessed through \url{https://pypbp-documentation.readthedocs.io} △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 40 pages

Pulsed interaction signals as a route to biological pattern formation

Authors: Eduardo H. Colombo, Cristóbal López, Emilio Hernández-García

Abstract: We identify a mechanism for biological spatial pattern formation arising when the signals that mediate interactions between individuals in a population have pulsed character. Our general population-signal framework shows that while for a slow signal-dynamics limit no pattern formation is observed for any values of the model parameters, for a fast limit, on the contrary, pattern formation can occur… ▽ More We identify a mechanism for biological spatial pattern formation arising when the signals that mediate interactions between individuals in a population have pulsed character. Our general population-signal framework shows that while for a slow signal-dynamics limit no pattern formation is observed for any values of the model parameters, for a fast limit, on the contrary, pattern formation can occur. Furthermore, at these limits, our framework reduces, respectively, to reaction-diffusion and spatially nonlocal models, thus bridging these approaches. △ Less

Submitted 13 August, 2023; v1 submitted 11 July, 2022; originally announced July 2022.

Journal ref: Phys. Rev. Lett. 130, 058401 (2023)

Quantitative CT texture-based method to predict diagnosis and prognosis of fibrosing interstitial lung disease patterns

Authors: Babak Haghighi, Warren B. Gefter, Lauren Pantalone, Despina Kontos, Eduardo Mortani Barbosa Jr

Abstract: Purpose: To utilize high-resolution quantitative CT (QCT) imaging features for prediction of diagnosis and prognosis in fibrosing interstitial lung diseases (ILD). Approach: 40 ILD patients (20 usual interstitial pneumonia (UIP), 20 non-UIP pattern ILD) were classified by expert consensus of 2 radiologists and followed for 7 years. Clinical variables were recorded. Following segmentation of the lu… ▽ More Purpose: To utilize high-resolution quantitative CT (QCT) imaging features for prediction of diagnosis and prognosis in fibrosing interstitial lung diseases (ILD). Approach: 40 ILD patients (20 usual interstitial pneumonia (UIP), 20 non-UIP pattern ILD) were classified by expert consensus of 2 radiologists and followed for 7 years. Clinical variables were recorded. Following segmentation of the lung field, a total of 26 texture features were extracted using a lattice-based approach (TM model). The TM model was compared with previously histogram-based model (HM) for their abilities to classify UIP vs non-UIP. For prognostic assessment, survival analysis was performed comparing the expert diagnostic labels versus TM metrics. Results: In the classification analysis, the TM model outperformed the HM method with AUC of 0.70. While survival curves of UIP vs non-UIP expert labels in Cox regression analysis were not statistically different, TM QCT features allowed statistically significant partition of the cohort. Conclusions: TM model outperformed HM model in distinguishing UIP from non-UIP patterns. Most importantly, TM allows for partitioning of the cohort into distinct survival groups, whereas expert UIP vs non-UIP labeling does not. QCT TM models may improve diagnosis of ILD and offer more accurate prognostication, better guiding patient management. △ Less

Submitted 20 June, 2022; originally announced June 2022.

Model for prognostic of symptomatic, asymptomatic and hospitalized COVID-19 cases with correct demography evolution

Authors: Antonio Rafael Selva Castañeda, Erick Eduardo Ramirez-Torres, Luis Eugenio Valdés-García, Hilda María Morandeira-Padrón, Diana Sedal Yanez, Juan I. Montijano, Luis Enrique Bergues Cabrales

Abstract: The aim of this study is to propose a modified Susceptible-Exposed-Infectious-Removed (SEIR) model that describes the behaviour of symptomatic, asymptomatic and hospitalized patients of COVID-19 epidemic, including the effect of demographic variation of population. It is shown that considering a population growth proportional to the total population leads to solutions with a qualitative behaviour… ▽ More The aim of this study is to propose a modified Susceptible-Exposed-Infectious-Removed (SEIR) model that describes the behaviour of symptomatic, asymptomatic and hospitalized patients of COVID-19 epidemic, including the effect of demographic variation of population. It is shown that considering a population growth proportional to the total population leads to solutions with a qualitative behaviour different from the behaviour obtained in many studies, where constant growth ratio is assumed. An exhaustive theoretical study is carried out and the basic reproduction number $R_0$ is computed from the model equations. It is proved that if $R_0<1$ then the disease-free manifold is globally asymptotically stable, that is, the epidemics remits. Global and local stability of the equilibrium points is also studied. Numerical simulations are used to show the agreement between numerical results and theoretical properties. The model is fitted to experimental data corresponding to the pandemic evolution in the República de Cuba, showing a proper behaviour of infected cases which let us think that can provide a correct estimation of asymptomatic cases. In conclusion, the model seems to be an adequate tool for the study and control of infectious diseases in particular the COVID-19 disease transmission. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Reaction Network Analysis of Metabolic Insulin Signaling

Authors: Patrick Vincent N. Lubenia, Eduardo R. Mendoza, Angelyn R. Lao

Abstract: Absolute concentration robustness (ACR) and concordance are novel concepts in the theory of robustness and stability within Chemical Reaction Network Theory. In this paper, we have extended Shinar and Feinberg's reaction network analysis approach to the insulin signaling system based on recent advances in decomposing reaction networks. We have shown that the network with 20 species, 35 complexes,… ▽ More Absolute concentration robustness (ACR) and concordance are novel concepts in the theory of robustness and stability within Chemical Reaction Network Theory. In this paper, we have extended Shinar and Feinberg's reaction network analysis approach to the insulin signaling system based on recent advances in decomposing reaction networks. We have shown that the network with 20 species, 35 complexes, and 35 reactions is concordant, implying at most one positive equilibrium in each of its stoichiometric compatibility class. We have obtained the system's finest independent decomposition consisting of 10 subnetworks, a coarsening of which reveals three subnetworks which are not only functionally but also structurally important. Utilizing the network's deficiency-oriented coarsening, we have developed a method to determine positive equilibria for the entire network. Our analysis has also shown that the system has ACR in 8 species all coming from a deficiency zero subnetwork. Interestingly, we have shown that, for a set of rate constants, the insulin-regulated glucose transporter GLUT4 (important in glucose energy metabolism), has stable ACR. △ Less

Submitted 12 September, 2022; v1 submitted 10 April, 2022; originally announced April 2022.

Comments: 34 pages, 1 figure

Ketamine-Medetomidine General Anesthesia Occurs With Alternation of Cortical Electrophysiological Activity Between High and Low Complex States

Authors: Eduardo Cocca Padovani

Abstract: Anesthetic agents are known to induce a range of alterations in cortical electrophysiological activity, such as the rise of signature patterns, changes in statistical properties, and altered dynamic behavior of neural records. Plenty of methods can be used to monitor these changes, among them complexity metrics demonstrated to have the power to discriminate states involving distinct levels of awar… ▽ More Anesthetic agents are known to induce a range of alterations in cortical electrophysiological activity, such as the rise of signature patterns, changes in statistical properties, and altered dynamic behavior of neural records. Plenty of methods can be used to monitor these changes, among them complexity metrics demonstrated to have the power to discriminate states involving distinct levels of awareness. There is a consensus that anesthetic drugs can interfere with neural activities at different levels and time scales, being able to induce alterations both locally and in the spatiotemporal patterns established throughout the whole cortex. However, it is still unclear how such changes in the complexity of cortical activity are supposed to occur, and experimental evidence is still needed. For this purpose, we have analyzed an ECoG records database of a Ketamine-Medetomidine anesthetic induction in a non-human primate subject. The MDR-ECoG technique provided records of cortical activity with both high temporal and spatial resolution allied with extensive coverage of the cortical surface. The Permutation Entropy and the Fractal Dimension were employed to evaluate the complexity of the neural time series. It was found that the complexity of cortical activity was relatively constant during awakened conditions. The transition to unconsciousness occurred relatively fast; it required about 30 to 40 seconds for the first remarkable changes to take place. During anesthesia, the complexity assumed considerable variation at local levels and fluctuated apparently without a defined period. The cortex dynamically alternated between high and low complex states on a global scale. This study provides novel evidence of the effects of anesthetics on neural activity and cortical dynamics, contributing to the elucidation of anesthetic mechanisms in terms of circuits, pathways, and global brain functioning. △ Less

Submitted 5 September, 2023; v1 submitted 9 February, 2022; originally announced February 2022.

Comments: 22 pages, 11 figures

Predicting radiotherapy patient outcomes with real-time clinical data using mathematical modelling

Authors: Alexander P. Browning, Thomas D. Lewin, Ruth E. Baker, Philip K. Maini, Eduardo G. Moros, Jimmy Caudell, Helen M. Byrne, Heiko Enderling

Abstract: Longitudinal tumour volume data from head-and-neck cancer patients show that tumours of comparable pre-treatment size and stage may respond very differently to the same radiotherapy fractionation protocol. Mathematical models are often proposed to predict treatment outcome in this context, and have the potential to guide clinical decision-making and inform personalised fractionation protocols. Hin… ▽ More Longitudinal tumour volume data from head-and-neck cancer patients show that tumours of comparable pre-treatment size and stage may respond very differently to the same radiotherapy fractionation protocol. Mathematical models are often proposed to predict treatment outcome in this context, and have the potential to guide clinical decision-making and inform personalised fractionation protocols. Hindering effective use of models in this context is the sparsity of clinical measurements juxtaposed with the model complexity required to produce the full range of possible patient responses. In this work, we present a compartment model of tumour volume and tumour composition, which, despite relative simplicity, is capable of producing a wide range of patient responses. We then develop novel statistical methodology and leverage a cohort of existing clinical data to produce a predictive model of both tumour volume progression and the associated level of uncertainty that evolves throughout a patient's course of treatment. To capture inter-patient variability, all model parameters are patient specific, with a bootstrap particle filter-like Bayesian approach developed to model a set of training data as prior knowledge. We validate our approach against a subset of unseen data, and demonstrate both the predictive ability of our trained model and its limitations. △ Less

Submitted 13 December, 2023; v1 submitted 6 January, 2022; originally announced January 2022.

An approach to interfacing the brain with quantum computers: practical steps and caveats

Authors: Eduardo Reck Miranda, Satvik Venkatesh, Jose D. Martın-Guerrero, Carlos Hernani-Morales, Lucas Lamata, Enrique Solano

Abstract: We report on the first proof-of-concept system demonstrating how one can control a qubit with mental activity. We developed a method to encode neural correlates of mental activity as instructions for a quantum computer. Brain signals are detected utilizing electrodes placed on the scalp of a person, who learns how to produce the required mental activity to issue instructions to rotate and measure… ▽ More We report on the first proof-of-concept system demonstrating how one can control a qubit with mental activity. We developed a method to encode neural correlates of mental activity as instructions for a quantum computer. Brain signals are detected utilizing electrodes placed on the scalp of a person, who learns how to produce the required mental activity to issue instructions to rotate and measure a qubit. Currently, our proof-of-concept runs on a software simulation of a quantum computer. At the time of writing, available quantum computing hardware and brain activity sensing technology are not sufficiently developed for real-time control of quantum states with the brain. But we are one step closer to interfacing the brain with real quantum machines, as improvements in hardware technology at both fronts become available in time to come. The paper ends with a discussion on some of the challenging problems that need to be addressed before we can interface the brain with quantum hardware. △ Less

Submitted 4 January, 2022; originally announced January 2022.

Comments: Journal pre-submission draft. Revision pending

From single to collective motion of social amoebae: a computational study of interacting cells

Authors: Eduardo Moreno, Robert Großmann, Carsten Beta, Sergio Alonso

Abstract: The coupling of the internal mechanisms of cell polarization to cell shape deformations and subsequent cell crawling poses many interdisciplinary scientific challenges. Several mathematical approaches have been proposed to model the coupling of both processes, where one of the most successful methods relies on a phase field that encodes the morphology of the cell, together with the integration of… ▽ More The coupling of the internal mechanisms of cell polarization to cell shape deformations and subsequent cell crawling poses many interdisciplinary scientific challenges. Several mathematical approaches have been proposed to model the coupling of both processes, where one of the most successful methods relies on a phase field that encodes the morphology of the cell, together with the integration of partial differential equations that account for the polarization mechanism inside the cell domain as defined by the phase field. This approach has been previously employed to model the motion of single cells of the social amoeba Dictyostelium discoideum, a widely used model organism to study actin-driven motility and chemotaxis of eukaryotic cells. Besides single cell motility, Dictyostelium discoideum is also well-known for its collective behavior. Here, we extend the previously introduced model for single cell motility to describe the collective motion of large populations of interacting amoebae by including repulsive interactions between the cells. We performed numerical simulations of this model, first characterizing the motion of single cells in terms of their polarity and velocity vectors. We then systematically studied the collisions between two cells that provided the basic interaction scenarios also observed in larger ensembles of interacting amoebae. Finally, the relevance of the cell density was analyzed, revealing a systematic decrease of the motility with density, associated with the formation of transient cell clusters that emerge in this system. This model is a prototypical active matter system for the investigation of the emergent collective dynamics of deformable, self-driven cells with a highly complex, nonlinear coupling of cell shape deformations, self-propulsion and repulsive cell-cell interactions. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: 25 pages, 12 Figures and 1 Table

Journal ref: Front. Phys. 9 (2022) 750187

Transient frailty induced by cell division. Observation, reasons and implications

Authors: Zengru Di, Eduardo M. Garcia-Roger, Peter Richmond, Bertrand M. Roehner, Stephane Tronche

Abstract: We know that stress-factors, e.g. X-rays, have an effect on cells that is more lethal in rapid exponential growth than in stationary phase. It is this effect which makes radiotherapy effective in cancer treatment. This stress effect can be explained in two ways: (a) more vulnerability in the growth phase, (b) improved protection capacity and repair mechanisms in the stationary phase. Although the… ▽ More We know that stress-factors, e.g. X-rays, have an effect on cells that is more lethal in rapid exponential growth than in stationary phase. It is this effect which makes radiotherapy effective in cancer treatment. This stress effect can be explained in two ways: (a) more vulnerability in the growth phase, (b) improved protection capacity and repair mechanisms in the stationary phase. Although the two explanations do not exclude each other, they are very different in the sense that (a) is a general mechanism whereas (b) is strain and stress-factor dependent. In this paper we explore major facets of (a). Firstly, we emphasize that (a) can account for known experimental stress-factor evidence. Secondly, we observe that (a) rightly predicts that slow exponential growth (meaning with a doubling time of several hours) results in a lower death rate than fast exponential growth (doubling time of a fraction of one hour), an effect that cannot be explained in the (b) framework because both organisms are in the same phase. Thirdly, we conjecture that the stress-factor effect can be extended to situations of chronic stress due to non-optimal environmental conditions. If correct, this conjecture would imply that even in normal culture conditions the natural death rate is lower in the stationary phase than in the growth phase. Finally, the paper closes with the description of several open questions and of appropriate test-experiments meant to address them. △ Less

Submitted 11 December, 2021; originally announced December 2021.

Comments: 21 pages, 4 figures

HEROHE Challenge: assessing HER2 status in breast cancer without immunohistochemistry or in situ hybridization

Authors: Eduardo Conde-Sousa, João Vale, Ming Feng, Kele Xu, Yin Wang, Vincenzo Della Mea, David La Barbera, Ehsan Montahaei, Mahdieh Soleymani Baghshah, Andreas Turzynski, Jacob Gildenblat, Eldad Klaiman, Yiyu Hong, Guilherme Aresta, Teresa Araújo, Paulo Aguiar, Catarina Eloy, António Polónia

Abstract: Breast cancer is the most common malignancy in women, being responsible for more than half a million deaths every year. As such, early and accurate diagnosis is of paramount importance. Human expertise is required to diagnose and correctly classify breast cancer and define appropriate therapy, which depends on the evaluation of the expression of different biomarkers such as the transmembrane prote… ▽ More Breast cancer is the most common malignancy in women, being responsible for more than half a million deaths every year. As such, early and accurate diagnosis is of paramount importance. Human expertise is required to diagnose and correctly classify breast cancer and define appropriate therapy, which depends on the evaluation of the expression of different biomarkers such as the transmembrane protein receptor HER2. This evaluation requires several steps, including special techniques such as immunohistochemistry or in situ hybridization to assess HER2 status. With the goal of reducing the number of steps and human bias in diagnosis, the HEROHE Challenge was organized, as a parallel event of the 16th European Congress on Digital Pathology, aiming to automate the assessment of the HER2 status based only on hematoxylin and eosin stained tissue sample of invasive breast cancer. Methods to assess HER2 status were presented by 21 teams worldwide and the results achieved by some of the proposed methods open potential perspectives to advance the state-of-the-art. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Common Complexes of Decompositions and Complex Balanced Equilibria of Chemical Reaction Networks

Authors: Lauro L. Fontanil, Eduardo R. Mendoza

Abstract: A decomposition of a chemical reaction network (CRN) is produced by partitioning its set of reactions. The partition induces networks, called subnetworks, that are "smaller" than the given CRN which, at this point, can be called parent network. A complex is called a common complex if it occurs in at least two subnetworks in a decomposition. A decomposition is said to be incidence independent if th… ▽ More A decomposition of a chemical reaction network (CRN) is produced by partitioning its set of reactions. The partition induces networks, called subnetworks, that are "smaller" than the given CRN which, at this point, can be called parent network. A complex is called a common complex if it occurs in at least two subnetworks in a decomposition. A decomposition is said to be incidence independent if the image of the incidence map of the parent network is the direct sum of the images of the subnetworks' incidence maps. It has been recently discovered that the complex balanced equilibria of the parent network and its subnetworks are fundamentally connected in an incidence independent decomposition. In this paper, we utilized the set of common complexes and a developed criterion to investigate decomposition's incidence independence properties. A framework was also developed to analyze decomposition classes with similar structure and incidence independence properties. We identified decomposition classes that can be characterized by their sets of common complexes and studied their incidence independence. Some of these decomposition classes occur in some biological and chemical models. Finally, a sufficient condition was obtained for the complex balancing of some power law kinetic (PLK) systems with incidence independent and complex balanced decompositions. This condition led to a generalization of the Defficiency Zero Theorem for some PLK systems. △ Less

Submitted 12 September, 2021; originally announced September 2021.

Macaque's Cortical Functional Connectivity Dynamics at the Onset of Propofol-Induced Anesthesia

Authors: Eduardo C. Padovani

Abstract: Propofol, when administered for general anesthesia, induces oscillatory dynamic brain states that are thought to underlie the drug's pharmacological effects. Despite the elucidation of propofol's mechanisms of action at the molecular level, its effects on neural circuits and overall cortical functioning, which eventually lead to unconsciousness, are still unclear. To identify possible mechanisms,… ▽ More Propofol, when administered for general anesthesia, induces oscillatory dynamic brain states that are thought to underlie the drug's pharmacological effects. Despite the elucidation of propofol's mechanisms of action at the molecular level, its effects on neural circuits and overall cortical functioning, which eventually lead to unconsciousness, are still unclear. To identify possible mechanisms, the spatial-temporal patterns of functional connectivity established among specialized cortical areas in anesthetized subjects need to be described. Within this context, the present research involved the analysis of dense sub-dural ECoG electrode array recordings from macaques under propofol anesthetic induction. Granger causality methodology was used to infer functional connectivity interactions in five physiological frequency bands serially over time, every five seconds throughout the experiment. The time-resolved networks obtained permitted us to observe the unfolding of the anesthetic induction and compare the networks obtained under different experimental conditions. About one minute after administering propofol, functional connectivity started to increase for 4-5 minutes, then decreased until the LOC was achieved. A predominant Granger causality flow from the occipital and temporal areas towards the frontal and parietal regions was also evidenced during the transition. During general anesthesia, the local connectivity of the occipital lobe increased, as did the interactions between the occipital and temporal lobes. Conversely, the functional connectivity from the frontal and parietal lobes toward the temporal and occipital regions was mainly impaired. The research is one of the first studies to describe the dynamics of the functional connectivity during the transitional state that precedes the LOC. △ Less

Submitted 5 September, 2023; v1 submitted 1 August, 2021; originally announced August 2021.

Comments: 20 pages, 5 figures

Is there an infant mortality in bacteria?

Authors: Eduardo M. Garcia-Roger, Peter Richmond, Bertrand M. Roehner

Abstract: This manuscript proposes a significant step in our long-run investigation of infant mortality across species. Since 2016 (Berrut et al. 2016) a succession of studies (Bois et al. 2019) have traced infant mortality from organisms of high complexity (e.g. mammals) down to unicellular organisms. Infant mortality may be considered as a filtering process through which organisms with potentially lethal… ▽ More This manuscript proposes a significant step in our long-run investigation of infant mortality across species. Since 2016 (Berrut et al. 2016) a succession of studies (Bois et al. 2019) have traced infant mortality from organisms of high complexity (e.g. mammals) down to unicellular organisms. Infant mortality may be considered as a filtering process through which organisms with potentially lethal congenital defects are eliminated. Such defects may have many causes but here we focus particularly on mishaps resulting from non-optimal conditions in the production of proteins, enzymes and other crucial macromolecules. The statistical signature of infant mortality consists in a falling age-specific death rate. The question we address here is whether infant mortality episodes take place in bacteria in the minutes precededing or following cell division. It will be shown that while experiments carried out in the 20th century tried but failed to detect such an effect (mostly because of limited sample size), more recent observations provided consistent evidence of a sizeable mortality, with a rate of the order of 0.7 per 1,000 per hour, in the exponential growth phase of E. coli. A further crucial test will be to measure the age-specific, post-division death rate. An experiment is outlined for that purpose. It is based on the selection of stained cells through flow cytometry and the derivation of their ages at death from their sizes. If an infant mortality effect can be identified in E. coli it can be conjectured that a similar effect also exists in other unicellular organisms, both prokaryote and eukaryote. △ Less

Submitted 18 March, 2021; originally announced March 2021.

Comments: 16 p., 5 figures