  Prognostic signatures play an important role in clinical research, offering
insights into the potential health outcomes of patients and guiding therapeutic
decisions. Although single-gene prognostic biomarkers are valuable, multi-gene
prognostic signatures offer deeper insights into disease progression. In this
paper, we propose MulMarker, a comprehensive framework for identifying
multi-gene prognostic signatures across various diseases. MulMarker comprises
three core modules: a chatbot for addressing user queries, a module for
identifying multi-gene prognostic signatures, and a module for generating
tailored reports. Employing MulMarker, we identified a cell cycle-related
prognostic signature that consists of CCNA1/2, CCNB1/2/3, CCNC, CCND1/2/3,
CCNE1/2, CCNF, CCNG1/2, and CCNH. Based on the prognostic signature, we
successfully stratified patients into high-risk and low-risk groups. Notably,
our analysis revealed that patients in the low-risk group demonstrated a
significantly higher survival rate compared to those in the high-risk group.
Overall, MulMarker represents an efficient approach for the identification of
multi-gene prognostic signatures. We release the code of MulMarker at
https://github.com/Tina9/MulMarker.


  In biochemical networks, complex dynamical features such as superlinear
growth and oscillations are classically considered a consequence of
autocatalysis. For the large class of parameter-rich kinetic models, which
includes Generalized Mass Action kinetics and Michaelis-Menten kinetics, we
show that certain submatrices of the stoichiometric matrix, so-called unstable
cores, are sufficient for a reaction network to admit instability and
potentially give rise to such complex dynamical behavior. The determinant of
the submatrix distinguishes unstable-positive feedbacks, with a single
real-positive eigenvalue, and unstable-negative feedbacks without real-positive
eigenvalues. Autocatalytic cores turn out to be exactly the unstable-positive
feedbacks that are Metzler matrices. Thus there are sources of dynamical
instability in chemical networks that are unrelated to autocatalysis. We use
such intuition to design non-autocatalytic biochemical networks with
superlinear growth and oscillations.


  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.


  Networks of interconnected agents are essential to study complex networked
systems' state evolution, stability, resilience, and control. Nevertheless, the
high dimensionality and nonlinear dynamics are vital factors preventing us from
theoretically analyzing them. Recently, the dimension-reduction approaches
reduced the system's size by mapping the original system to a one-dimensional
system such that only one effective representative can capture its macroscopic
dynamics. However, the approaches dramatically fail as the network becomes
heterogeneous and has multiple community structures. Here, we bridge the gap by
developing a generalized dimension reduction approach, which enables us to map
the original system to a $m$-dimensional system that consists of $m$
interacting components. Notably, by validating it on various dynamical models,
this approach accurately predicts the original system state and the tipping
point, if any. Furthermore, the numerical results demonstrate that this
approach approximates the system evolution and identifies the critical points
for complex networks with time delay.


  Sequence modelling approaches for epigenetic profile prediction have recently
expanded in terms of sequence length, model size, and profile diversity.
However, current models cannot infer on many experimentally feasible tissue and
assay pairs due to poor usage of contextual information, limiting $\textit{in
silico}$ understanding of regulatory genomics. We demonstrate that strong
correlation can be achieved across a large range of experimental conditions by
integrating tissue and assay embeddings into a Contextualised Genomic Network
(CGN). In contrast to previous approaches, we enhance long-range sequence
embeddings with contextual information in the input space, rather than
expanding the output space. We exhibit the efficacy of our approach across a
broad set of epigenetic profiles and provide the first insights into the effect
of genetic variants on epigenetic sequence model training. Our general approach
to context integration exceeds state of the art in multiple settings while
employing a more rigorous validation procedure.


  Ligand-based drug design aims to identify novel drug candidates of similar
shapes with known active molecules. In this paper, we formulated an in silico
shape-conditioned molecule generation problem to generate 3D molecule
structures conditioned on the shape of a given molecule. To address this
problem, we developed a translation- and rotation-equivariant shape-guided
generative model ShapeMol. ShapeMol consists of an equivariant shape encoder
that maps molecular surface shapes into latent embeddings, and an equivariant
diffusion model that generates 3D molecules based on these embeddings.
Experimental results show that ShapeMol can generate novel, diverse, drug-like
molecules that retain 3D molecular shapes similar to the given shape condition.
These results demonstrate the potential of ShapeMol in designing drug
candidates of desired 3D shapes binding to protein target pockets.


  This work presents a general chemical reaction network theory for olfactory
sensing processes that employ G-protein-coupled receptors as olfactory
receptors (ORs). The theory is applicable to general mixtures of odorants and
an arbitrary number of ORs. Reactions of ORs with G-proteins, both in the
presence and the absence of odorants, are explicitly considered. A unique
feature of the theory is the definition of an odor activity vector consisting
of strengths of odorant-induced signals from ORs relative to those due to
background G-protein activity in the absence of odorants. It is demonstrated
that each component of the odor activity defined this way reduces to a
Michaelis-Menten form capable of accounting for cooperation or competition
effects between different odorants. The main features of the theory are
illustrated for a two-odorant mixture. Known and potential mixture effects,
such as suppression, shadowing, inhibition, and synergy are quantitatively
described. Effects of relative values of rate constants, basal activity, and
G-protein concentration are also demonstrated.


  Graph generation poses a significant challenge as it involves predicting a
complete graph with multiple nodes and edges based on simply a given label.
This task also carries fundamental importance to numerous real-world
applications, including de-novo drug and molecular design. In recent years,
several successful methods have emerged in the field of graph generation.
However, these approaches suffer from two significant shortcomings: (1) the
underlying Graph Neural Network (GNN) architectures used in these methods are
often underexplored; and (2) these methods are often evaluated on only a
limited number of metrics. To fill this gap, we investigate the expressiveness
of GNNs under the context of the molecular graph generation task, by replacing
the underlying GNNs of graph generative models with more expressive GNNs.
Specifically, we analyse the performance of six GNNs on six different molecular
generative objectives on the ZINC-250k dataset in two different generative
frameworks: autoregressive generation models, such as GCPN and GraphAF, and
one-shot generation models, such as GraphEBM. Through our extensive
experiments, we demonstrate that advanced GNNs can indeed improve the
performance of GCPN, GraphAF, and GraphEBM on molecular generation tasks, but
GNN expressiveness is not a necessary condition for a good GNN-based generative
model. Moreover, we show that GCPN and GraphAF with advanced GNNs can achieve
state-of-the-art results across 17 other non-GNN-based graph generative
approaches, such as variational autoencoders and Bayesian optimisation models,
on the proposed molecular generative objectives (DRD2, Median1, Median2), which
are important metrics for de-novo molecular design.


  Machine learning (ML) is a promising approach for predicting small molecule
properties in drug discovery. Here, we provide a comprehensive overview of
various ML methods introduced for this purpose in recent years. We review a
wide range of properties, including binding affinities, solubility, and ADMET
(Absorption, Distribution, Metabolism, Excretion, and Toxicity). We discuss
existing popular datasets and molecular descriptors and embeddings, such as
chemical fingerprints and graph-based neural networks. We highlight also
challenges of predicting and optimizing multiple properties during hit-to-lead
and lead optimization stages of drug discovery and explore briefly possible
multi-objective optimization techniques that can be used to balance diverse
properties while optimizing lead candidates. Finally, techniques to provide an
understanding of model predictions, especially for critical decision-making in
drug discovery are assessed. Overall, this review provides insights into the
landscape of ML models for small molecule property predictions in drug
discovery. So far, there are multiple diverse approaches, but their
performances are often comparable. Neural networks, while more flexible, do not
always outperform simpler models. This shows that the availability of
high-quality training data remains crucial for training accurate models and
there is a need for standardized benchmarks, additional performance metrics,
and best practices to enable richer comparisons between the different
techniques and models that can shed a better light on the differences between
the many techniques.


  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 present eight reaction network operations that can
modify a reaction network while preserving the dynamical system. In particular,
we study the structural and kinetic properties of these transformed systems,
such as stoichiometric alignment, concordance and existence of equilibria. To
illustrate our approach, we use the kinetic realizations of the following
biochemical systems: (1) calcium signaling in the olfactory system, (2)
metabolic insulin signaling and (3) monolignol biosynthesis in \textit{Populus
xylem}. Furthermore, we present a sufficient condition for the invariance of
concordance under a network transformation and introduce an algorithm for a
network transformation of a weakly reversible non-CFK system to a weakly
reversible CFK system. We review the work of Hong et al. on network
translation, discussing its relevance in this study. Finally, we end by
highlighting some benefits of positive deficiency kinetic realizations for the
analysis of biochemical systems.


  The comprehension of how local interactions arise in global collective
behavior is of utmost importance in both biological and physical research.
Traditional agent-based models often rely on static rules that fail to capture
the dynamic strategies of the biological world. Reinforcement learning has been
proposed as a solution, but most previous methods adopt handcrafted reward
functions that implicitly or explicitly encourage the emergence of swarming
behaviors. In this study, we propose a minimal predator-prey coevolution
framework based on mixed cooperative-competitive multiagent reinforcement
learning, and adopt a reward function that is solely based on the fundamental
survival pressure, that is, prey receive a reward of $-1$ if caught by
predators while predators receive a reward of $+1$. Surprisingly, our analysis
of this approach reveals an unexpectedly rich diversity of emergent behaviors
for both prey and predators, including flocking and swirling behaviors for
prey, as well as dispersion tactics, confusion, and marginal predation
phenomena for predators. Overall, our study provides novel insights into the
collective behavior of organisms and highlights the potential applications in
swarm robotics.


  Our main goal is to examine the role of communities in epidemic spread in a
random graph model. More precisely, we consider a random graph model which
consists of overlapping complete graphs, representing households, workplaces,
school classes, and which also has a simple geometric structure. We study the
model's sensitivity to infection parameters and other tunable parameters of the
model, which might be helpful in finding efficient social distancing
strategies. We also quantitatively compare different vaccination strategies to
see which order is the best to defend the most vulnerable groups or the
population in general, and how important it is to gather and use information on
the position of infected individuals in the network.


  Biological cells replicate their genomes in a well-planned manner. The DNA
replication program of an organism determines the timing at which different
genomic regions are replicated, with fundamental consequences for cell
homeostasis and genome stability. Qualitatively, in a growing cell culture, one
expects that genomic regions that are replicated early should be more abundant
than regions that are replicated late. This abundance pattern can be
experimentally measured using deep sequencing. However, a general quantitative
theory to explain these data is still lacking. In this paper, we predict the
abundance of DNA fragments in asynchronously growing cultures from any given
stochastic model of the DNA replication program. As key examples, we present
stochastic models of the DNA replication programs in Escherichia coli and in
budding yeast. In both cases, our approach leads to analytical predictions that
are in excellent agreement with experimental data and permit to infer key
information about the replication program. In particular, our method is able to
infer the locations of known replication origins in budding yeast with high
accuracy. These examples demonstrate that our method can provide insight into a
broad range of organisms, from bacteria to eukaryotes.


  Information on the structure of molecules, retrieved via biochemical
databases, plays a pivotal role in various disciplines, such as metabolomics,
systems biology, and drug discovery. However, no such database can be complete,
and the chemical structure for a given compound is not necessarily consistent
between databases. This paper presents StructRecon, a novel tool for resolving
unique and correct molecular structures from database identifiers. StructRecon
traverses the cross-links between database entries in different databases to
construct what we call an identifier graph, which offers a more complete view
of the total information available on a particular compound across all the
databases. In order to reconcile discrepancies between databases, we first
present an extensible model for chemical structure which supports multiple
independent levels of detail, allowing standardisation of the structure to be
applied iteratively. In some cases, our standardisation approach results in
multiple structures for a given compound, in which case a random walk-based
algorithm is used to select the most likely structure among incompatible
alternates. We applied StructRecon to the EColiCore2 model, resolving a unique
chemical structure for 85.11 % of identifiers. StructRecon is open-source and
modular, which enables the potential support for more databases in the future.


  An important application of Synthetic Biology is the engineering of the host
cell system to yield useful products. However, an increase in the scale of the
host system leads to huge design space and requires a large number of
validation trials with high experimental costs. A comprehensible machine
learning approach that efficiently explores the hypothesis space and guides
experimental design is urgently needed for the Design-Build-Test-Learn (DBTL)
cycle of the host cell system. We introduce a novel machine learning framework
ILP-iML1515 based on Inductive Logic Programming (ILP) that performs abductive
logical reasoning and actively learns from training examples. In contrast to
numerical models, ILP-iML1515 is built on comprehensible logical
representations of a genome-scale metabolic model and can update the model by
learning new logical structures from auxotrophic mutant trials. The ILP-iML1515
framework 1) allows high-throughput simulations and 2) actively selects
experiments that reduce the experimental cost of learning gene functions in
comparison to randomly selected experiments.


  Phenotype switching can be triggered by external stimuli and by intrinsic
stochasticity. Here, we focus on the motility-matrix production switch in
Bacillus subtilis. We use modeling to describe the SinR-SlrR bistable switch
its regulation by SinI, and to distinguish different sources of stochasticity.
Our simulations indicate that intrinsic fluctuations in the synthesis of SinI
are insufficient to drive spontaneous switching and suggest that switching is
triggered by upstream noise from the Spo0A phosphorelay.


  PHF5A is a member of the zinc-finger proteins. To advance knowledge on their
role in carcinogenesis, data from experimental studies, animal models and
clinical studies in different tumorigenesis have been reviewed. Furthermore,
PHF5A as an oncogenic function, is frequently expressed in tumor cells and a
potential prognostic marker for different cancers. PHF5A is implicated in the
regulation of cancer cell proliferation, invasion, migration and metastasis.
Knockdown of PHF5A prevented the invasion and metastasis of tumor cells. Here,
the role of PHF5A in different cancers and their possible mechanism in relation
to recent literature is reviewed and discussed. However, there is an open
promising perspective to their therapeutic management for different cancer
types.


  Optimizing the properties of molecules (materials or drugs) for stronger
toughness, lower toxicity, or better bioavailability has been a long-standing
challenge. In this context, we propose a molecular optimization framework
called Q-Drug (Quantum-inspired optimization algorithm for Drugs) that
leverages quantum-inspired algorithms to optimize molecules on discrete binary
domain variables. The framework begins by encoding the molecules into binary
embeddings using a discrete VAE. The binary embeddings are then used to
construct an Ising energy-like objective function, over which the
state-of-the-art quantum-inspired optimization algorithm is adopted to find the
optima. The binary embeddings corresponding to the optima are decoded to obtain
the optimized molecules. We have tested the framework for optimizing drug
molecule properties and have found that it outperforms other molecular
optimization methods, finding molecules with better properties in 1/20th to
1/10th of the time previously required. The framework can also be deployed
directly on various quantum computing equipment, such as laser pulses CIMs,
FPGA Ising Machines, and quantum computers based on quantum annealing, among
others. Our work demonstrates a new paradigm that leverages the advantages of
quantum computing and AI to solve practically useful problems.


  We study the effect of noisy infection (contact) and recovery rates on the
distribution of outbreak sizes in the stochastic SIR model. The rates are
modeled as Ornstein-Uhlenbeck processes with finite correlation time and
variance, which we illustrate using outbreak data from the RSV 2019-2020 season
in the US. In the limit of large populations, we find analytical solutions for
the outbreak-size distribution in the long-correlated (adiabatic) and
short-correlated (white) noise regimes, and demonstrate that the distribution
can be highly skewed with significant probabilities for large fluctuations away
from mean-field theory. Furthermore, we assess the relative contribution of
demographic and reaction-rate noise on the outbreak-size variance, and show
that demographic noise becomes irrelevant in the presence of slowly varying
reaction-rate noise but persists for large system sizes if the noise is fast.
Finally, we show that the crossover to the white-noise regime typically occurs
for correlation times that are on the same order as the characteristic recovery
time in the model.


  The theory of large random matrices has proved an invaluable tool for the
study of systems with disordered interactions in many quite disparate research
areas. Widely applicable results, such as the celebrated elliptic law for dense
random matrices, allow one to deduce the statistical properties of the
interactions in a complex dynamical system that permit stability. However, such
simple and universal results have so far proved difficult to come by in the
case of sparse random matrices. Here, we present a new approach, which maps the
hermitized resolvent of a random matrix onto the response functions of a linear
dynamical system. The response functions are then evaluated using a path
integral formalism, enabling one to construct Feynman diagrams and perform a
perturbative analysis. This approach provides simple closed-form expressions
for the eigenvalue spectrum, allowing one to derive modified versions of the
classic elliptic and semi-circle laws that take into account the sparse
correction. Additionally, in order to demonstrate the broad utility of the path
integral framework, we derive a non-Hermitian generalization of the
Marchenko-Pastur law, and we also show how one can handle non-negligible
higher-order statistics (i.e. non-Gaussian statistics) in dense ensembles.


  In the modern world, technology is at its peak. Different avenues in
programming and technology have been explored for data analysis, automation,
and robotics. Machine learning is key to optimize data analysis, make accurate
predictions, and hasten/improve existing functions. Thus, presently, the field
of machine learning in artificial intelligence is being developed and its uses
in varying fields are being explored. One field in which its uses stand out is
that of microbial biosynthesis. In this paper, a comprehensive overview of the
differing machine learning programs used in biosynthesis is provided, alongside
brief descriptions of the fields of machine learning and microbial biosynthesis
separately. This information includes past trends, modern developments, future
improvements, explanations of processes, and current problems they face. Thus,
this paper's main contribution is to distill developments in, and provide a
holistic explanation of, 2 key fields and their applicability to improve
industry/research. It also highlights challenges and research directions,
acting to instigate more research and development in the growing fields.
Finally, the paper aims to act as a reference for academics performing
research, industry professionals improving their processes, and students
looking to understand the concept of machine learning in biosynthesis.


  We investigate the mechanisms behind the quasi-periodic outbursts on the
Covid-19 epidemics. Data for France and Germany show that the patterns of
outbursts exhibit a qualitative change in early 2022, which appears in a change
in their average period, and which is confirmed by time-frequency analysis.
This provides a signal which can be used to discriminate among several
mechanisms. Two main ideas have been proposed to explain periodicity in
epidemics. One involves memory effects and another considers exchanges between
epidemic clusters and a reservoir of population. We test these two approaches
in the particular case of the Covid-19 epidemics and show that the "cluster
model" is the only one which appears to be able to explain the observed pattern
with realistic parameters. A last section discusses our results in the context
of early studies of epidemics, and we stress the importance to work with models
with a limited number of parameters, which moreover can be sufficiently well
estimated, to draw conclusions on the general mechanisms behind the
observations.


  In this paper, we develop a computational approach for computing most likely
trajectories describing rare events that correspond to the emergence of
non-dominant genotypes. This work is based on the large deviations approach for
discrete Markov chains describing the genetic evolution of large bacterial
populations. We demonstrate that a gradient descent algorithm developed in this
paper results in the fast and accurate computation of most-likely trajectories
for a large number of bacterial genotypes. We supplement our analysis with
extensive numerical simulations demonstrating the computational advantage of
the designed gradient descent algorithm over other, more simplified,
approaches.


  Current antibody language models are limited by their use of unpaired
antibody sequence data and the biases in publicly available antibody sequence
datasets, which are skewed toward antibodies against a relatively small number
of pathogens. A recently published dataset (by Jaffe, et al) of approximately
1.6 x 10^6 natively paired human antibody sequences from healthy donors
represents by far the largest dataset of its kind and offers a unique
opportunity to evaluate how antibody language models can be improved by
training with natively paired antibody sequence data. We trained two Baseline
Antibody Language Models (BALM), using natively paired (BALM-paired) or
unpaired (BALM-unpaired) sequences from the Jaffe dataset. We provide evidence
that training with natively paired sequences substantially improves model
performance and that this improvement results from the model learning
immunologically relevant features that span the light and heavy chains. We also
show that ESM-2, a state-of-the-art general protein language model, learns
similar cross-chain features when fine-tuned with natively paired antibody
sequence data.


  Seawater temperature rise is damaging coral reef ecosystems. There is growing
evidence for the negative impact of rising temperatures on the survival of
adult corals and their reproductive success. However, the effect of elevated
temperatures on gametes remains scarcely studied. Here we tested the effect of
the thermal priming of gametes on the fertilization success in experimentally
tested populations of Acropora cytherea corals in French Polynesia. As
expected, a temperature of 30 {\textdegree}C (ambient +3 {\textdegree}C)
reduces coral fertilization success. However, the thermal exposure of gametes
to 30 {\textdegree}C after their release in seawater prior to fertilization
limited fertilization failure, with a greater impact of oocytes in comparison
to sperm. This temperature is similar to temperatures observed in nature under
the changing climate. Our findings imply that the thermal priming of early life
stages, such as gametes may play a role in maintaining the coral fertilization
success in spite of increasing seawater temperature.


  In this paper, we first reviewed several biclustering methods that are used
to identify the most significant clusters in gene expression data. Here we
mainly focused on the SSVD(sparse SVD) method and tried a new sparse penalty
named "Prenet penalty" which has been used only in factor analysis to gain
sparsity. Then in the simulation study, we tried different types of generated
datasets (with different sparsity and dimension) and tried 1-layer
approximation then for k-layers which shows the mixed Prenet penalty is very
effective for non-overlapped data. Finally, we used some real gene expression
data to show the behavior of our methods.


  Our study is devoted to a four-compartment epidemic model of a constant
population of independent random walkers. Each walker is in one of four
compartments (S-susceptible, C-infected but not infectious (period of
incubation), I-infected and infectious, R-recovered and immune) characterizing
the states of health. The walkers navigate independently on a periodic 2D
lattice. Infections occur by collisions of susceptible and infectious walkers.
Once infected, a walker undergoes the delayed cyclic transition pathway S $\to$
C $\to$ I $\to$ R $\to$ S. The random delay times between the transitions
(sojourn times in the compartments) are drawn from independent probability
density functions (PDFs). We analyze the existence of the endemic equilibrium
and stability of the globally healthy state and derive a condition for the
spread of the epidemics which we connect with the basic reproduction number
$R_0>1$. We give quantitative numerical evidence that a simple approach based
on random walkers offers an appropriate microscopic picture of the dynamics for
this class of epidemics.


  The classical model for the genealogies of a neutrally evolving population in
a fixed environment is due to Kingman. Kingman's coalescent process, which
produces a binary tree, universally emerges from many microscopic models in
which the variance in the number of offspring is finite. It is understood that
power-law offspring distributions with infinite variance can result in a very
different type of coalescent structure with merging of more than two lineages.
Here we investigate the regime where the variance of the offspring distribution
is finite but comparable to the population size. This is achieved by studying a
model in which the log offspring sizes have a stretched exponential form. Such
offspring distributions are motivated by biology, where they emerge from a toy
model of growth in a heterogenous environment, but also mathematics and
statistical physics, where limit theorems and phase transitions for sums over
random exponentials have received considerable attention due to their
appearance in the partition function of Derrida's Random Energy Model (REM). We
find that the limit coalescent is a $\beta$-coalescent -- a previously studied
model emerging from evolutionary dynamics models with heavy-tailed offspring
distributions. We also discuss the connection to previous results on the REM.


  Recent works have shown the promise of learning pre-trained models for 3D
molecular representation. However, existing pre-training models focus
predominantly on equilibrium data and largely overlook off-equilibrium
conformations. It is challenging to extend these methods to off-equilibrium
data because their training objective relies on assumptions of conformations
being the local energy minima. We address this gap by proposing a force-centric
pretraining model for 3D molecular conformations covering both equilibrium and
off-equilibrium data. For off-equilibrium data, our model learns directly from
their atomic forces. For equilibrium data, we introduce zero-force
regularization and forced-based denoising techniques to approximate
near-equilibrium forces. We obtain a unified pre-trained model for 3D molecular
representation with over 15 million diverse conformations. Experiments show
that, with our pre-training objective, we increase forces accuracy by around 3
times compared to the un-pre-trained Equivariant Transformer model. By
incorporating regularizations on equilibrium data, we solved the problem of
unstable MD simulations in vanilla Equivariant Transformers, achieving
state-of-the-art simulation performance with 2.45 times faster inference time
than NequIP. As a powerful molecular encoder, our pre-trained model achieves
on-par performance with state-of-the-art property prediction tasks.


  It is estimated that approximately 15% of cancers worldwide can be linked to
viral infections. The viruses that can cause or increase the risk of cancer
include human papillomavirus, hepatitis B and C viruses, Epstein-Barr virus,
and human immunodeficiency virus, to name a few. The computational analysis of
the massive amounts of tumor DNA data, whose collection is enabled by the
recent advancements in sequencing technologies, have allowed studies of the
potential association between cancers and viral pathogens. However, the high
diversity of oncoviral families makes reliable detection of viral DNA difficult
and thus, renders such analysis challenging. In this paper, we introduce XVir,
a data pipeline that relies on a transformer-based deep learning architecture
to reliably identify viral DNA present in human tumors. In particular, XVir is
trained on genomic sequencing reads from viral and human genomes and may be
used with tumor sequence information to find evidence of viral DNA in human
cancers. Results on semi-experimental data demonstrate that XVir is capable of
achieving high detection accuracy, generally outperforming state-of-the-art
competing methods while being more compact and less computationally demanding.


  The genome-scale metabolic model with protein constraint (PC-model) has been
increasingly popular for microbial metabolic simulations. We present PROSO
Toolbox, a unified and simple-to-use PC-model toolbox that takes any
high-quality genome-scale metabolic reconstruction as the input. The toolbox
can construct a PC-model automatically, apply various algorithms for
computational strain design and simulation, and help unveil metabolism from
gene expression data through a state-of-the-art OVERLAY workflow. It also has
detailed tutorials and documentation for maximum accessibility to researchers
from diverse backgrounds. PROSO Toolbox, tutorials, and documentation are
freely available online: https://github.com/QCSB/PROSO-Toolbox.


  We present a stochastic framework to decipher fluctuations propagation in
classes of coherent feed-forward loops. The systematic contribution of the
direct (one-step) and indirect (two-step) pathways is considered to quantify
fluctuations of the output node. We also consider both additive and
multiplicative integration mechanisms of the two parallel pathways (one-step
and two-step). Analytical expression of the output node's coefficient of
variation shows contributions of intrinsic, one-step, two-step, and
cross-interaction in closed form. We observe a diverse range of degeneracy and
non-degeneracy in each of the decomposed fluctuations term and their
contribution to the overall output fluctuations of each coherent feed-forward
loop motif. Analysis of output fluctuations reveals a maximal level of
fluctuations of the coherent feed-forward loop motif of type 1.


  Molecular dynamics simulations have emerged as a fundamental instrument for
studying biomolecules. At the same time, it is desirable to perform simulations
of a collection of particles under various conditions in which the molecules
can fluctuate. In this paper, we explore and adapt the soft prompt-based
learning method to molecular dynamics tasks. Our model can remarkably
generalize to unseen and out-of-distribution scenarios with limited training
data. While our work focuses on temperature as a test case, the versatility of
our approach allows for efficient simulation through any continuous dynamic
conditions, such as pressure and volumes. Our framework has two stages: 1)
Pre-trains with data mixing technique, augments molecular structure data and
temperature prompts, then applies a curriculum learning method by increasing
the ratio of them smoothly. 2) Meta-learning-based fine-tuning framework
improves sample-efficiency of fine-tuning process and gives the soft
prompt-tuning better initialization points. Comprehensive experiments reveal
that our framework excels in accuracy for in-domain data and demonstrates
strong generalization capabilities for unseen and out-of-distribution samples.


  Non-equilibrium cluster-cluster aggregation of particles diffusing in or at
the cell membrane has been hypothesized to lead to domains of finite size in
different biological contexts such as lipid rafts, cell adhesion complexes, or
postsynaptic domains in neurons. In this scenario, the desorption of particles
balances a continuous flux to the membrane, imposing a cut-off on possible
aggregate sizes and giving rise to a stationary size distribution. Here, we
investigate the case of non-equilibrium cluster-cluster aggregation in two
dimensions where diffusing particles and/or clusters remain fixed in space at
specific anchoring sites, which should be particularly relevant for synapses
but may also be present in other biological or physical systems. Using an
effective mean-field description of the concentration field around anchored
clusters, we derive an expression for their average size as a function of
parameters such as the anchoring site density. We furthermore propose and solve
appropriate rate equations that allow us to predict the size distributions of
both diffusing and fixed clusters. We confirm our results with particle-based
simulations, and discuss potential implications for biological and physical
systems.


  In drug discovery, mapping interactions between genes within cellular systems
is a crucial early step. This helps formulate hypotheses regarding molecular
mechanisms that could potentially be targeted by future medicines. The
CausalBench Challenge was an initiative to invite the machine learning
community to advance the state of the art in constructing gene-gene interaction
networks. These networks, derived from large-scale, real-world datasets of
single cells under various perturbations, are crucial for understanding the
causal mechanisms underlying disease biology. Using the framework provided by
the CausalBench benchmark, participants were tasked with enhancing the capacity
of the state of the art methods to leverage large-scale genetic perturbation
data. This report provides an analysis and summary of the methods submitted
during the challenge to give a partial image of the state of the art at the
time of the challenge. The winning solutions significantly improved performance
compared to previous baselines, establishing a new state of the art for this
critical task in biology and medicine.


  The ability to precisely edit genomes by deleting or adding genetic
information enables the study of biological functions and the building of
efficient cell factories. In many unconventional yeasts, such as promising new
hosts for cell factory design but also human pathogenic yeasts and food
spoilers, this progress has been limited by the fact that most yeasts favor
non-homologous end joining (NHEJ) over homologous recombination (HR) as DNA
repair mechanism, impairing genetic access to these hosts. In mammalian cells,
small molecules that either inhibit proteins involved in NHEJ, enhance protein
function in HR, or molecules that arrest the cell cycle in HR-dominant phases
are regarded as promising agents for the simple and transient increase of
HR-mediated genome editing without the need for a priori host engineering. Only
a few of these chemicals have been applied to the engineering of yeast although
the targeted proteins are mostly conserved; making chemical agents a yet
underexplored area in enhancing yeast engineering. Here, we consolidate
knowledge of available small molecules that have been used to improve HR
efficiency in mammalian cells and the few ones that have been used in yeast. We
include available high throughput (HTP)-compatible NHEJ/HR quantification
assays that could be used to screen for and isolate yeast-specific inhibitors.


  Concentrations of pathogen genomes measured in wastewater have recently
become available as a new data source to use when modeling the spread of
infectious diseases. One promising use for this data source is inference of the
effective reproduction number, the average number of individuals a newly
infected person will infect. We propose a model where new infections arrive
according to a time-varying immigration rate which can be interpreted as a
compound parameter equal to the product of the proportion of susceptibles in
the population and the transmission rate. This model allows us to estimate the
effective reproduction number from concentrations of pathogen genomes while
avoiding difficult to verify assumptions about the dynamics of the susceptible
population. As a byproduct of our primary goal, we also produce a new model for
estimating the effective reproduction number from case data using the same
framework. We test this modeling framework in an agent-based simulation study
with a realistic data generating mechanism which accounts for the time-varying
dynamics of pathogen shedding. Finally, we apply our new model to estimating
the effective reproduction number of SARS-CoV-2 in Los Angeles, California,
using pathogen RNA concentrations collected from a large wastewater treatment
facility.


  Infectious disease forecasting has been a key focus and proved to be crucial
in controlling epidemic. A recent trend is to develop forecast-ing models based
on graph neural networks (GNNs). However, existing GNN-based methods suffer
from two key limitations: (1) Current models broaden receptive fields by
scaling the depth of GNNs, which is insuffi-cient to preserve the semantics of
long-range connectivity between distant but epidemic related areas. (2)
Previous approaches model epidemics within single spatial scale, while ignoring
the multi-scale epidemic pat-terns derived from different scales. To address
these deficiencies, we devise the Multi-scale Spatio-temporal Graph Neural
Network (MSGNN) based on an innovative multi-scale view. To be specific, in the
proposed MSGNN model, we first devise a novel graph learning module, which
directly captures long-range connectivity from trans-regional epidemic signals
and integrates them into a multi-scale graph. Based on the learned multi-scale
graph, we utilize a newly designed graph convolution module to exploit
multi-scale epidemic patterns. This module allows us to facilitate multi-scale
epidemic modeling by mining both scale-shared and scale-specific pat-terns.
Experimental results on forecasting new cases of COVID-19 in United State
demonstrate the superiority of our method over state-of-arts. Further analyses
and visualization also show that MSGNN offers not only accurate, but also
robust and interpretable forecasting result.


  Polycyclic aromatic hydrocarbons (PAHs) are highly toxic, carcinogenic
substances. On soils contaminated with PAHs, crop cultivation, animal husbandry
and even the survival of microflora in the soil are greatly perturbed,
depending on the degree of contamination. Most microorganisms cannot tolerate
PAH-contaminated soils, however, some microbial strains can adapt to these
harsh conditions and survive on contaminated soils. Analysis of the metagenomes
of contaminated environmental samples may lead to discovery of PAH-degrading
enzymes suitable for green biotechnology methodologies ranging from
biocatalysis to pollution control. In the present study, our goal was to apply
a metagenomic data search to identify efficient novel enzymes in remediation of
PAH-contaminated soils. The metagenomic hits were further analyzed using a set
of bioinformatics tools to select protein sequences predicted to encode
well-folded soluble enzymes. Three novel enzymes (two dioxygenases and one
peroxidase) were cloned and used in soil remediation microcosms experiments.
The novel enzymes were found to be efficient for degradation of naphthalene and
phenanthrene. Adding the inorganic oxidant CaO2 further increased the degrading
potential of the novel enzymes for anthracene and pyrene. We conclude that
metagenome mining paired with bioinformatic predictions, structural modelling
and functional assays constitutes a powerful approach towards novel enzymes for
soil remediation.


  Throughout the history of electron microscopy, ribosomes have served as an
ideal subject for imaging and technological development, which in turn has
driven our understanding of ribosomal biology. Here, we provide a historical
perspective at the intersection of electron microscopy technology development
and ribosome biology and reflect on how this technique has shed light on each
stage of the life cycle of this dynamic macromolecular machine. With an
emphasis on prokaryotic systems, we specifically describe how pairing cryo-EM
with clever experimental design, time-resolved techniques, and next-generation
heterogeneous structural analysis has afforded insights into the modular nature
of assembly, the roles of the many transient biogenesis and translation
co-factors, and the subtle variations in structure and function between strains
and species. The work concludes with a prospective outlook on the field,
highlighting the pivotal role cryogenic electron tomography is playing in
adding cellular context to our understanding of ribosomal life cycles, and
noting how this exciting technology promises to bridge the gap between cellular
and structural biology.


  Biofilms inhabit a range of environments, such as dental plaques or soil
micropores, often characterized by intricate, non-even surfaces. However, the
impact of surface irregularities on the population dynamics of biofilms remains
elusive as most biofilm experiments are conducted on flat surfaces. Here, we
show that the shape of the surface on which a biofilm grows influences genetic
drift and selection within the biofilm. We culture E. coli biofilms in
micro-wells with an undulating bottom surface and observe the emergence of
clonal sectors whose size corresponds to that of the undulations, despite no
physical barrier separating different areas of the biofilm. The sectors are
remarkably stable over time and do not invade each other; we attribute this
stability to the characteristics of the velocity field within the growing
biofilm, which hinders mixing and clonal expansion. A microscopically-detailed
computer model fully reproduces these findings and highlights the role of
mechanical (physical) interactions such as adhesion and friction in microbial
evolution. The model also predicts clonal expansion to be severely limited even
for clones with a significant growth advantage - a finding which we
subsequently confirm experimentally using a mixture of antibiotic-sensitive and
antibiotic-resistant mutants in the presence of sub-lethal concentrations of
the antibiotic rifampicin. The strong suppression of selection contrasts
sharply with the behavior seen in bacterial colonies on agar commonly used to
study range expansion and evolution in biofilms. Our results show that biofilm
population dynamics can be controlled by patterning the surface, and
demonstrate how a better understanding of the physics of bacterial growth can
pave the way for new strategies in steering microbial evolution.


  In this article, we will establish a number of results concerning the
limiting behaviour of the longest edges in the genealogical tree generated by a
continuous-time Galton-Watson (GW) process. Separately, we consider the large
time behaviour of the longest pendant edges, the longest (strictly) interior
edges, and the longest of all the edges. These results extend the special case
of long pendant edges of birth-death processes established in Bocharov, Harris,
Kominek, Mooers, and Steel [1] .


  Retrosynthesis planning is a fundamental challenge in chemistry which aims at
designing reaction pathways from commercially available starting materials to a
target molecule. Each step in multi-step retrosynthesis planning requires
accurate prediction of possible precursor molecules given the target molecule
and confidence estimates to guide heuristic search algorithms. We model
single-step retrosynthesis planning as a distribution learning problem in a
discrete state space. First, we introduce the Markov Bridge Model, a generative
framework aimed to approximate the dependency between two intractable discrete
distributions accessible via a finite sample of coupled data points. Our
framework is based on the concept of a Markov bridge, a Markov process pinned
at its endpoints. Unlike diffusion-based methods, our Markov Bridge Model does
not need a tractable noise distribution as a sampling proxy and directly
operates on the input product molecules as samples from the intractable prior
distribution. We then address the retrosynthesis planning problem with our
novel framework and introduce RetroBridge, a template-free retrosynthesis
modeling approach that achieves state-of-the-art results on standard evaluation
benchmarks.


  Signalling pathways are conserved across different species, therefore making
yeast a model organism to study these via disruption of kinase activity. Yeast
has 159 genes that encode protein kinases and phosphatases, and 136 of these
have counterparts in humans. Therefore any insight in this model organism could
potentially offer indications of mechanisms of action in the human kinome. The
study utilises a Prolog-based approach, data from a yeast kinase deletions
strains study and publicly available kinase-protein associations. Prolog, a
programming language that is well-suited for symbolic reasoning is used to
reason over the data and infer compensatory kinase networks. This approach is
based on the idea that when a kinase is knocked out, other kinases may
compensate for this loss of activity. Background knowledge on kinases targeting
proteins is used to guide the analysis. This knowledge is used to infer the
potential compensatory interactions between kinases based on the changes in
phosphorylation observed in the phosphoproteomics data from the yeast study.
The results demonstrate the effectiveness of the Prolog-based approach in
analysing complex cell signalling mechanisms in yeast. The inferred
compensatory kinase networks provide new insights into the regulation of cell
signalling in yeast and may aid in the identification of potential therapeutic
targets for modulating signalling pathways in yeast and other organisms.


  We quantitatively assess the impacts of Downgrading Protected Areas (PAD) on
biodiversity in the U.S.. Results show that PAD events significantly reduce
biodiversity. The proximity to PAD events decreases the biodiversity by 26.0%
within 50 km compared with records of species further away from the PAD events.
We observe an overall 32.3% decrease in abundance after those nearest PAD
events are enacted. Abundance declines more in organisms living in contact with
water and non-mammals. Species abundance is more sensitive to the negative
impacts in areas where PAD events were later reversed, as well as in areas
close to protected areas belonging to the International Union for Conservation
of Nature (IUCN) category. The enacted PAD events between the period 1903 to
2018 in the U.S. lead to economic losses of approximately $689.95 million due
to decrease in abundance. Our results contribute to the understanding on the
impact of environmental interventions such as PAD events on biodiversity change
and provide important implications on biodiversity conservation policies.


  Molecular chaperones are vital proteins that maintain protein homeostasis by
assisting in protein folding, activation, degradation, and stress protection.
Among them, heat-shock protein 90 (Hsp90) stands out as an essential
proteostasis hub in eukaryotes, chaperoning hundreds of "clients" (substrates).
After decades of research, several "known unknowns" about the molecular
function of Hsp90 remain unanswered, hampering rational drug design for the
treatment of cancers, neurodegenerative and other diseases. We highlight three
fundamental open questions, reviewing the current state of the field for each,
and discuss new opportunities, including single-molecule technologies, to
answer the known unknowns of the Hsp90 chaperone.


  In recent decades, antibodies have emerged as indispensable therapeutics for
combating diseases, particularly viral infections. However, their development
has been hindered by limited structural information and labor-intensive
engineering processes. Fortunately, significant advancements in deep learning
methods have facilitated the precise prediction of protein structure and
function by leveraging co-evolution information from homologous proteins.
Despite these advances, predicting the conformation of antibodies remains
challenging due to their unique evolution and the high flexibility of their
antigen-binding regions. Here, to address this challenge, we present the
Bio-inspired Antibody Language Model (BALM). This model is trained on a vast
dataset comprising 336 million 40% non-redundant unlabeled antibody sequences,
capturing both unique and conserved properties specific to antibodies. Notably,
BALM showcases exceptional performance across four antigen-binding prediction
tasks. Moreover, we introduce BALMFold, an end-to-end method derived from BALM,
capable of swiftly predicting full atomic antibody structures from individual
sequences. Remarkably, BALMFold outperforms those well-established methods like
AlphaFold2, IgFold, ESMFold, and OmegaFold in the antibody benchmark,
demonstrating significant potential to advance innovative engineering and
streamline therapeutic antibody development by reducing the need for
unnecessary trials.


  Phylogenetic trees describe the relationships between species in the
evolutionary process, and provide information about the rates of
diversification. To understand the mechanisms behind macroevolution, we
consider a class of multitype branching processes called Markovian binary trees
(MBTs). MBTs allow for trait-based variation in diversification rates, and
provide a flexible and realistic probabilistic model for phylogenetic trees. We
use an approximate Bayesian computation (ABC) scheme to infer the rates of MBTs
by exploiting information in the shapes of phylogenetic trees. We evaluate the
accuracy of this inference method using simulation studies, and find that our
method is accurate and able to detect variation in the diversification rates.
In an application to a real-life phylogeny of squamata, we reinforce
conclusions drawn from earlier studies, in particular supporting the existence
of ovi-/viviparity transitions in both directions. This demonstrates the
potential for more complex models of evolution to be employed in phylogenetic
inference.


  As new COVID-19 variants continue to emerge globally, we develop an SISI-type
mathematical model to determine whether the transmission risks arising when
such a new variant enters Indonesia can be mitigated by solely accelerating the
country's ongoing four-dose vaccination programme. We begin by determining the
model's basic reproduction number, as well as the model's equilibria and their
stability. Subsequently, employing parameter values representing the country's
situation as of March 20, 2023, we conduct a numerical sensitivity analysis in
two simulated cases corresponding to two different levels of the new variant's
transmission. The results show that, a satisfactory mitigation relying solely
on vaccinations necessitates a drastic acceleration in the low-transmission
case, and proves unachievable in the high-transmission case. Accordingly, we
recommend that the acceleration of the ongoing four-dose vaccinations be
carried out in conjunction with other intervention measures, such as
improvements of the vaccine's efficacy and the disease's recovery rate.


  Plasmodium vivax is one of the most geographically widespread malaria
parasites in the world due to its ability to remain dormant in the human liver
as hypnozoites and subsequently reactivate after the initial infection (i.e.
relapse infections). More than 80% of P. vivax infections are due to hypnozoite
reactivation. Mathematical modelling approaches have been widely applied to
understand P. vivax dynamics and predict the impact of intervention outcomes.
In this article, we provide a scoping review of mathematical models that
capture P. vivax transmission dynamics published between January 1988 and May
2023 to provide a comprehensive summary of the mathematical models and
techniques used to model P. vivax dynamics. We aim to assist researchers
working on P. vivax transmission and other aspects of P. vivax malaria by
highlighting best practices in currently published models and highlighting
where future model development is required. We provide an overview of the
different strategies used to incorporate the parasite's biology, use of
multiple scales (within-host and population-level), superinfection, immunity,
and treatment interventions. In most of the published literature, the rationale
for different modelling approaches was driven by the research question at hand.
Some models focus on the parasites' complicated biology, while others
incorporate simplified assumptions to avoid model complexity. Overall, the
existing literature on mathematical models for P. vivax encompasses various
aspects of the parasite's dynamics. We recommend that future research should
focus on refining how key aspects of P. vivax dynamics are modelled, including
spatial heterogeneity in exposure risk, the accumulation of hypnozoite
variation, the interaction between P. falciparum and P. vivax, acquisition of
immunity, and recovery under superinfection.


  BRCA genes, comprising BRCA1 and BRCA2 play indispensable roles in preserving
genomic stability and facilitating DNA repair mechanisms. The presence of
germline mutations in these genes has been associated with increased
susceptibility to various cancers, notably breast and ovarian cancers. Recent
advancements in cost-effective sequencing technologies have revolutionized the
landscape of cancer genomics, leading to a notable rise in the number of
sequenced cancer patient genomes, enabling large-scale computational studies.
In this study, we delve into the BRCA mutations in the dbSNP, housing an
extensive repository of 41,177 and 44,205 genetic mutations for BRCA1 and
BRCA2, respectively. Employing meticulous computational analysis from an
umbrella perspective, our research unveils intriguing findings pertaining to a
number of critical aspects. Namely, we discover that the majority of BRCA
mutations in dbSNP have unknown clinical significance. We find that, although
exon 11 for both genes contains the majority of the mutations and may seem as
if it is a mutation hot spot, upon analyzing mutations per base pair, we find
that all exons exhibit similar levels of mutations. Investigating mutations
within introns, while we observe that the recorded mutations are generally
uniformly distributed, almost all of the pathogenic mutations in introns are
located close to splicing regions (at the beginning or the end). In addition to
the findings mentioned earlier, we have also made other discoveries concerning
mutation types and the level of confidence in observations within the dbSNP
database.


  Scientific computing has experienced a surge empowered by advancements in
technologies such as neural networks. However, certain important tasks are less
amenable to these technologies, benefiting from innovations to traditional
inference schemes. One such task is protein re-design. Recently a new re-design
algorithm, AOBB-K*, was introduced and was competitive with state-of-the-art
BBK* on small protein re-design problems. However, AOBB-K* did not scale well.
In this work we focus on scaling up AOBB-K* and introduce three new versions:
AOBB-K*-b (boosted), AOBB-K*-DH (with dynamic heuristics), and AOBB-K*-UFO
(with underflow optimization) that significantly enhance scalability.


  Molecular property prediction with deep learning has gained much attention
over the past years. Owing to the scarcity of labeled molecules, there has been
growing interest in self-supervised learning methods that learn generalizable
molecular representations from unlabeled data. Molecules are typically treated
as 2D topological graphs in modeling, but it has been discovered that their 3D
geometry is of great importance in determining molecular functionalities. In
this paper, we propose the Geometry-aware line graph transformer (Galformer)
pre-training, a novel self-supervised learning framework that aims to enhance
molecular representation learning with 2D and 3D modalities. Specifically, we
first design a dual-modality line graph transformer backbone to encode the
topological and geometric information of a molecule. The designed backbone
incorporates effective structural encodings to capture graph structures from
both modalities. Then we devise two complementary pre-training tasks at the
inter and intra-modality levels. These tasks provide properly supervised
information and extract discriminative 2D and 3D knowledge from unlabeled
molecules. Finally, we evaluate Galformer against six state-of-the-art
baselines on twelve property prediction benchmarks via downstream fine-tuning.
Experimental results show that Galformer consistently outperforms all baselines
on both classification and regression tasks, demonstrating its effectiveness.


  Single and collective cell migration are fundamental processes critical for
physiological phenomena ranging from embryonic development and immune response
to wound healing and cancer metastasis. To understand cell migration from a
physical perspective, a broad variety of models for the underlying physical
mechanisms that govern cell motility have been developed. A key challenge in
the development of such models is how to connect them to experimental
observations, which often exhibit complex stochastic behaviours. In this
review, we discuss recent advances in data-driven theoretical approaches that
directly connect with experimental data to infer dynamical models of stochastic
cell migration. Leveraging advances in nanofabrication, image analysis, and
tracking technology, experimental studies now provide unprecedented large
datasets on cellular dynamics. In parallel, theoretical efforts have been
directed towards integrating such datasets into physical models from the single
cell to the tissue scale with the aim of conceptualizing the emergent behaviour
of cells. We first review how this inference problem has been addressed in both
freely migrating and confined cells. Next, we discuss why these dynamics
typically take the form of underdamped stochastic equations of motion, and how
such equations can be inferred from data. We then review applications of
data-driven inference and machine learning approaches to heterogeneity in cell
behaviour, subcellular degrees of freedom, and to the collective dynamics of
multicellular systems. Across these applications, we emphasize how data-driven
methods can be integrated with physical active matter models of migrating
cells, and help reveal how underlying molecular mechanisms control cell
behaviour.


  The cornerstone of structural biology is the unique relationship between
protein sequence and the 3D structure at equilibrium. Although intrinsically
disordered proteins (IDPs) do not fold into a specific 3D structure, breaking
this paradigm, some IDPs exhibit large-scale organization, such as
liquid-liquid phase separation. In such cases, the structural plasticity has
the potential to form numerous self-assembled structures out of thermal
equilibrium. Here, we report that high-temperature incubation time is a
defining parameter for micro and nanoscale self-assembly of resilin-like IDPs.
Interestingly, high-resolution scanning electron microscopy micrographs reveal
that an extended incubation time leads to the formation of micron-size rods and
ellipsoids that depend on the amino acid sequence. More surprisingly, a
prolonged incubation time also induces amino acid composition-dependent
formation of short-range nanoscale order, such as periodic lamellar
nanostructures. We, therefore, suggest that regulating the period of
high-temperature incubation, in the one-phase regime, can serve as a unique
method of controlling the hierarchical self-assembly mechanism of structurally
disordered proteins.


  Forest fire spreading is a complex phenomenon characterized by a stochastic
behavior. Nowadays, the enormous quantity of georeferenced data and the
availability of powerful techniques for their analysis can provide a very
careful picture of forest fires opening the way to more realistic models. We
propose a stochastic spreading model continuous in space and time that is able
to use such data in their full power. The state of the forest fire is described
by the subprobability densities of the green trees and of the trees on fire
that can be estimated thanks to data coming from satellites and earth
detectors. The fire dynamics is encoded into a density probability kernel which
can take into account wind conditions, land slope, spotting phenomena and so
on, bringing to a system of integro-differential equations for the probability
densities. Existence and uniqueness of the solutions is proved by using
Banach's fixed point theorem. The asymptotic behavior of the model is analyzed
as well. Stochastic models based on cellular automata can be considered as
particular cases of the present model from which they can be derived by space
and/or time discretization. Suggesting a particular structure for the kernel,
we obtain numerical simulations of the fire spreading under different
conditions. For example, in the case of a forest fire evolving towards a river,
the simulations show that the probability density of the trees on fire is
different from zero beyond the river due to the spotting phenomenon.
Firefighters interventions and weather changes can be easily introduced into
the model.


  In an increasingly changing world, predicting the fate of species across the
globe has become a major concern. Understanding how the population dynamics of
various species and communities will unfold requires predictive tools that
experimental data alone can not capture. Here, we introduce our combinatorial
framework, Widespread Ecological Networks and their Dynamical Signatures
(WENDyS) which, using data on the relative strengths of interactions and growth
rates within a community of species predicts all possible long-term outcomes of
the community. To this end, WENDyS partitions the multidimensional parameter
space (formed by the strengths of interactions and growth rates) into a finite
number of regions, each corresponding to a unique set of coarse population
dynamics. Thus, WENDyS ultimately creates a library of all possible outcomes
for the community. On the one hand, our framework avoids the typical
``parameter sweeps'' that have become ubiquitous across other forms of
mathematical modeling, which can be computationally expensive for ecologically
realistic models and examples. On the other hand, WENDyS opens the opportunity
for interdisciplinary teams to use standard experimental data (i.e., strengths
of interactions and growth rates) to filter down the possible end states of a
community. To demonstrate the latter, here we present a case study from the
Indonesian Coral Reef. We analyze how different interactions between anemone
and anemonefish species lead to alternative stable states for the coral reef
community, and how competition can increase the chance of exclusion for one or
more species. WENDyS, thus, can be used to anticipate ecological outcomes and
test the effectiveness of management (e.g., conservation) strategies.


  We introduce a mathematical model for the mechanical behaviour of the
eukaryotic cell cytoskeleton. This discrete model involves a regular array of
pre-stressed protein filaments that exhibit resistance to enthalpic stretching,
joined at crosslinks to form a network. Assuming that the inter-crosslink
distance is much shorter than the lengthscale of the cell, we upscale the
discrete force balance to form a continuum system of governing equations and
deduce the corresponding macroscopic stress tensor. We use these discrete and
continuum models to analyse the imposed displacement of a bead placed in the
domain, characterising the cell rheology through the force-displacement curve.
We further derive an analytical approximation to the stress and strain fields
in the limit of small bead radius, predicting the net force required to
generate a given deformation and elucidating the dependency on the microscale
properties of the filaments. We apply these models to networks of the
intermediate filament vimentin and demonstrate good agreement between
predictions of the discrete, continuum and analytical approaches. In
particular, our model predicts that the network stiffness increases sublinearly
with the filament pre-stress and scales logarithmically with the bead size.


  There are multiple factors which can cause the phylogenetic inference process
to produce two or more conflicting hypotheses of the evolutionary history of a
set X of biological entities. That is: phylogenetic trees with the same set of
leaf labels X but with distinct topologies. This leads naturally to the goal of
quantifying the difference between two such trees T_1 and T_2. Here we
introduce the problem of computing a 'maximum relaxed agreement forest' (MRAF)
and use this as a proxy for the dissimilarity of T_1 and T_2, which in this
article we assume to be unrooted binary phylogenetic trees. MRAF asks for a
partition of the leaf labels X into a minimum number of blocks S_1, S_2, ...
S_k such that for each i, the subtrees induced in T_1 and T_2 by S_i are
isomorphic up to suppression of degree-2 nodes and taking the labels X into
account. Unlike the earlier introduced maximum agreement forest (MAF) model,
the subtrees induced by the S_i are allowed to overlap. We prove that it is
NP-hard to compute MRAF, by reducing from the problem of partitioning a
permutation into a minimum number of monotonic subsequences (PIMS).
Furthermore, we show that MRAF has a polynomial time O(log n)-approximation
algorithm where n=|X| and permits exact algorithms with single-exponential
running time. When at least one of the two input trees has a caterpillar
topology, we prove that testing whether a MRAF has size at most k can be
answered in polynomial time when k is fixed. We also note that on two
caterpillars the approximability of MRAF is related to that of PIMS. Finally,
we establish a number of bounds on MRAF, compare its behaviour to MAF both in
theory and in an experimental setting and discuss a number of open problems.


  DCAlign is a new alignment method able to cope with the conservation and the
co-evolution signals that characterize the columns of multiple sequence
alignments of homologous sequences. However, the pre-processing steps required
to align a candidate sequence are computationally demanding. We show in v1.0
how to dramatically reduce the overall computing time by including an empirical
prior over an informative set of variables mirroring the presence of insertions
and deletions.


  Recently, there has been interest in extending long-known results about the
multispecies coalescent tree to other models of gene trees. Results about the
gene duplication and loss model have mathematical proofs, including species
tree identifiability, estimability, and sample complexity of popular algorithms
like ASTRAL. Here, this work is continued by characterizing the anomaly zones
of uniformly sampled gene trees. The anomaly zone for species trees is the set
of parameters where some discordant gene tree occurs with the maximal
probability. The detection of anomalous gene trees is an important problem in
phylogenomics, as their presence renders effective estimation methods to being
positively misleading. Under the multispecies coalescent, anomaly zones are
known to exist for rooted species trees with as few as four species.
  The gene duplication and loss process is a generalization of the generalized
linear-birth death process to the rooted species tree, where each edge is
treated as a single timeline with exponential-rate duplication and loss. The
methods and results come from a detailed probabilistic analysis of trajectories
observed from this stochastic process. It is shown that anomaly zones do not
exist for rooted GDL balanced trees on four species, but they may exist for
rooted caterpillar trees, as with MSC.


  Biological sequence analysis relies on the ability to denoise the imprecise
output of sequencing platforms. We consider a common setting where a short
sequence is read out repeatedly using a high-throughput long-read platform to
generate multiple subreads, or noisy observations of the same sequence.
Denoising these subreads with alignment-based approaches often fails when too
few subreads are available or error rates are too high. In this paper, we
propose a novel method for blindly denoising sets of sequences without directly
observing clean source sequence labels. Our method, Self-Supervised Set
Learning (SSSL), gathers subreads together in an embedding space and estimates
a single set embedding as the midpoint of the subreads in both the latent and
sequence spaces. This set embedding represents the "average" of the subreads
and can be decoded into a prediction of the clean sequence. In experiments on
simulated long-read DNA data, SSSL methods denoise small reads of $\leq 6$
subreads with 17% fewer errors and large reads of $>6$ subreads with 8% fewer
errors compared to the best baseline. On a real dataset of antibody sequences,
SSSL improves over baselines on two self-supervised metrics, with a significant
improvement on difficult small reads that comprise over 60% of the test set. By
accurately denoising these reads, SSSL promises to better realize the potential
of high-throughput DNA sequencing data for downstream scientific applications.


  Plant-frugivore interactions play a central role for plant persistence and
spatial distribution by promoting the long-range dispersal of seeds by
frugivores. However, plant-frugivore interactions are increasingly being
threatened by anthropogenic activities. An important anthropogenic threat that
could expose plant-frugivore systems to extinction risk is fruit harvesting.
Here, we develop an individual-based and a pair-approximation model of
plant-frugivore-human interactions to elucidate the effects of human harvesting
of fruits on plant establishment, persistence, and spatial distribution. Our
results show that frugivores strongly affect global density of plants and
gradually shift their spatial distribution from aggregated to random, depending
on the attack rate and dispersal efficiency of frugivores. We find that, in the
absence of frugivores, plants experiencing intense fruit harvesting cannot
persist even if their fecundity is high. In the presence of frugivores, fruit
harvesting profoundly affects the global dispersal of seeds and thus changes
the spatial distributions of plants from random to aggregated, potentially
causing plant extinction. Our results demonstrate that sufficiently efficient
frugivores mitigate the negative impact of fruit harvesting on plant
populations and enable plant establishment precluded by harvesting. Taken
together, these results draw attention to previously underappreciated impacts
of fruit harvesting in plant-frugivore-human interactions.


  A collection of identical and independent rare event first passage times is
considered. The problem of finding the fastest out of $N$ such events to occur
is called an extreme first passage time. The rare event times are singular and
limit to infinity as a positive parameter scaling the noise magnitude is
reduced to zero. In contrast, previous work has shown that the mean of the
fastest event time goes to zero in the limit of an infinite number of walkers.
The combined limit is studied. In particular, the mean time and the most likely
path taken by the fastest random walker are investigated. Using techniques from
large deviation theory, it is shown that there is a distinguished limit where
the mean time for the fastest walker can take any positive value, depending on
a single proportionality constant. Furthermore, it is shown that the mean time
and most likely path can be approximated using the solution to a variational
problem related to the single-walker rare event.


  The ubiquity of information transmission via molecular communication between
cells is comprehensively documented on Earth; this phenomenon might even have
played a vital role in the origin(s) and early evolution of life. Motivated by
these considerations, a simple model for molecular communication entailing the
diffusion of signaling molecules from transmitter to receiver is elucidated.
The channel capacity $C$ (maximal rate of information transmission) and an
optimistic heuristic estimate of the actual information transmission rate
$\mathcal{I}$ are derived for this communication system; the two quantities,
especially the latter, are demonstrated to be broadly consistent with
laboratory experiments and more sophisticated theoretical models. The channel
capacity exhibits a potentially weak dependence on environmental parameters,
whereas the actual information transmission rate may scale with the
intercellular distance $d$ as $\mathcal{I} \propto d^{-4}$ and could vary
substantially across settings. These two variables are roughly calculated for
diverse astrobiological environments, ranging from Earth's upper oceans ($C
\sim 3.1 \times 10^3$ bits/s; $\mathcal{I} \sim 4.7 \times 10^{-2}$ bits/s) and
deep sea hydrothermal vents ($C \sim 4.2 \times 10^3$ bits/s; $\mathcal{I} \sim
1.2 \times 10^{-1}$ bits/s) to the hydrocarbon lakes and seas of Titan ($C \sim
3.8 \times 10^3$ bits/s; $\mathcal{I} \sim 2.6 \times 10^{-1}$ bits/s).


  In this paper, we consider an individual-based model with power law mutation
probability. In this setting, we use the large population limit with a
subsequent ``small mutations'' limit to derive the canonical equation of
adaptive dynamics. For a one-dimensional trait space this corresponds to well
established results and we can formulate a criterion for evolutionary branching
in the spirit of Champagnat and M\'el\'eard (2011). In higher dimensional trait
spaces, we find that the speed at which the solution of the canonical equation
moves through space is reduced due to mutations being restricted to the
underlying grid on the trait space. However, as opposed to the canonical
equation with rare mutations, we can explicitly calculate the path which the
dominant trait will take without having to solve the equation itself.


  In the last years algebraic tools have been proven to be useful in
phylogenetic reconstruction and model selection by means of the study of
phylogenetic invariants. However, up to now, the models studied from an
algebraic viewpoint are either too general or too restrictive (as group-based
models with a uniform stationary distribution) to be used in practice.
  In this paper we provide a new framework to work with time-reversible models,
which are the most widely used by biologists. In our approach we consider
algebraic time-reversible models on phylogenetic trees (as defined by Allman
and Rhodes) and introduce a new inner product to make all transition matrices
of the process diagonalizable through the same orthogonal eigenbasis. This
framework generalizes the Fourier transform widely used to work with
group-based models and recovers some of the well known results. As
illustration, we exploit the combination of our technique with algebraic
geometry tools to provide relevant phylogenetic invariants for trees evolving
under the Tamura-Nei model of nucleotide substitution.


  Climate change substantially impacts ecological systems. Marine species are
shifting their distribution because of climate change towards colder waters,
potentially compromising the benefits of currently established Marine Protected
Areas (MPAs). Therefore, we demonstrate how three case study regions will be
impacted by warming ocean waters to prepare stakeholders to understand how the
fisheries around the MPAs is predicted to change. We chose the case studies to
focus on large scale MPAs in i) a cold, polar region, ii) a tropical region
near the equator, and iii) a tropical region farther from the equator. We
quantify the biological impacts of shifts in species distribution due to
climate change for fishing communities that depend on the Palau National Marine
Sanctuary, the Great Barrier Reef Marine National Park Zone, and the North
Bering Sea Research Area MPAs. We find that fisheries sectors will be impacted
differently in different regions and show that all three regions can be
supported by this methodology for decision making that joins sector income and
species diversity.


  Mathematical and computational approaches in chemistry and biochemistry fill
a gap in respect to the analysis of the physicochemical features of compounds
and their functionality and provide an overview of known as well as yet
unknown, but hypothetically possible structures. Nitrogen containing
heterocycles such as aziridine, azetidine and pyrrolidine bear a high potential
in pharmacology, biotechnology and synthetic biology. Here, we present a
mathematical enumeration procedure for all possible azaheterocycles with at
least one substituent depending on the number of atoms in the ring, in the
sense of saturated and unsaturated congeners. One subgroup belonging to that
substance class is constituted by ring-shaped amino acids with a secondary
amino group such as proline. A recursion formula is derived, which results in a
modified Lucas number series. Moreover, an explicit formula for determining the
number of such substances based on the Golden Ratio is given and a second one,
based on binomial coefficients, is newly derived. This enumeration is a helpful
tool for construction or complementation of virtual compound databases and for
computer-assisted chemical synthesis route planning.


  We construct and analyze monomeric and multimeric models of the stochastic
disassembly of a single nucleosome. Our monomeric model predicts the time
needed for a number of histone-DNA contacts to spontaneously break, leading to
dissociation of a non-fragmented histone from DNA. The dissociation process can
be facilitated by DNA binding proteins or processing molecular motors that
compete with histones for histone-DNA contact sites. Eigenvalue analysis of the
corresponding master equation allows us to evaluate histone detachment times
under both spontaneous detachment and protein-facilitated processes. We find
that competitive DNA binding of remodeling proteins can significantly reduce
the typical detachment time but only if these remodelers have DNA-binding
affinities comparable to those of histone-DNA contact sites. In the presence of
processive motors, the histone detachment rate is shown to be proportional to
the product of the histone single-bond dissociation constant and the speed of
motor protein procession. Our simple intact-histone model is then extended to
allow for multimeric nucleosome kinetics that reveal additional pathways of
disassembly. In addition to a dependence of complete disassembly times on
subunit-DNA contact energies, we show how histone subunit concentrations in
bulk solution can mediate the disassembly process by rescuing partially
disassembled nucleosomes. Moreover, our kinetic model predicts that remodeler
binding can also bias certain pathways of nucleosome disassembly, with higher
remodeler binding rates favoring intact-histone detachment.


  There have been many studies to examine whether one trait is correlated with
another trait across a group of present-day species (for example, do species
with larger brains tend to have longer gestation times. Since the introduction
of the phylogenetic comparative method some authors have argued that it is
necessary to have a biologically realistic model to generate evolutionary trees
that incorporates information about the ecological niche occupied by species.
Price presented a simple model along these lines in 1997. He defined a
two-dimensional niche space formed by two continuous-valued traits, in which
new niches arise with trait values drawn from a bivariate normal distribution.
When a new niche arises, it is occupied by a descendant species of whichever
current species is closest in ecological niche space. In sequence, more species
are then evolved from already-existing species to which they are ecologically
closest.
  Here we explore ways of extending Price's adaptive radiation model. One
extension is to increase the dimensionality of the niche space by considering
more than two continuous traits. A second extension is to allow both extinction
of species (which may leave unoccupied niches) and removal of niches (which
causes species occupying them to go extinct). To model this problem, we
consider a continuous-time stochastic process which implicitly defines a
phylogeny. To explore if trees generated under such a model (or under different
parametrizations of the model) are realistic we can compute a variety of
summary statistics that can be compared to those of empirically observed
phylogenies. For example, there are existing statistics that aim to measure:
tree balance, the relative rate of diversification, and phylogenetic signal of
traits.


  In this work, we develop a stochastic model of gene gain and loss with the
aim of inferring when (if at all) in evolutionary history and association
between two genes arises. The data we consider is a species tree along with
information on the presence or absence of two genes in each of the species. The
biological motivation for our model is that if two genes are involved in the
same biochemical pathway, i.e. they are both required for some function, then
the rate of gain or loss of one gene in the pathway should depend upon the
presence or absence of the other gene in the pathway. However, if the two genes
are not functionally linked, then the rate of gain or loss of one gene should
be independent of the state of another gene.
  We simulate data under this model to determine under what conditions a shift
from the independent rates class to the dependent rates class can be detected.
For example, how large a tree is required and how large a shift in the rates is
needed before Akaike information criterion (AIC) supports a model with two rate
classes over a simpler model with just one rate class? If a model with two rate
classes is preferred, can it correctly detect where on the evolutionary tree
the shift occurred?


  Natural selection imply that any organisms including human being will evolve
to improve its fitness advantage and the selected genotype or phenotype in
equilibrium state will not vary over the time. However, evolutionary process of
biological organisms in reality is greatly affected by the environmental change
and historical accidents. In this research, we construct a co-evolutionary
system to investigate the impact of species-environment feedback. When we talk
about an invasion species or mutation, positive feedback is detrimental to the
success of the invasion because positive feedback benefits a large number of
individuals, whereas negative feedback benefits the invasion because negative
feedback disadvantages a large number of individuals. In the case of a
competition between two species with initially equal numbers of individuals,
both positive and negative feedback will favor the species with low fitness,
increasing its chances of taking over the whole population. The reason for this
is that feedback allows initially inferior species to have greater fitness than
initially dominating species in the early stages, emphasizing the importance of
early random accident. Our findings emphasize the significance of the
evolutionary path driven by species-environment feedback.


  Recent advances in Language Models have enabled the protein modeling
community with a powerful tool since protein sequences can be represented as
text. Specifically, by taking advantage of Transformers, sequence-to-property
prediction will be amenable without the need for explicit structural data. In
this work, inspired by recent progress in Large Language Models (LLMs), we
introduce PeptideBERT, a protein language model for predicting three key
properties of peptides (hemolysis, solubility, and non-fouling). The
PeptideBert utilizes the ProtBERT pretrained transformer model with 12
attention heads and 12 hidden layers. We then finetuned the pretrained model
for the three downstream tasks. Our model has achieved state of the art (SOTA)
for predicting Hemolysis, which is a task for determining peptide's potential
to induce red blood cell lysis. Our PeptideBert non-fouling model also achieved
remarkable accuracy in predicting peptide's capacity to resist non-specific
interactions. This model, trained predominantly on shorter sequences, benefits
from the dataset where negative examples are largely associated with insoluble
peptides. Codes, models, and data used in this study are freely available at:
https://github.com/ChakradharG/PeptideBERT


  With the fast-growing and evolving omics data, the demand for streamlined and
adaptable tools to handle the analysis continues to grow. In response to this
need, we introduce Auto Bioinformatics Analysis (AutoBA), an autonomous AI
agent based on a large language model designed explicitly for conventional
omics data analysis. AutoBA simplifies the analytical process by requiring
minimal user input while delivering detailed step-by-step plans for various
bioinformatics tasks. Through rigorous validation by expert bioinformaticians,
AutoBA's robustness and adaptability are affirmed across a diverse range of
omics analysis cases, including whole genome sequencing (WGS), RNA sequencing
(RNA-seq), single-cell RNA-seq, ChIP-seq, and spatial transcriptomics. AutoBA's
unique capacity to self-design analysis processes based on input data
variations further underscores its versatility. Compared with online
bioinformatic services, AutoBA deploys the analysis locally, preserving data
privacy. Moreover, different from the predefined pipeline, AutoBA has
adaptability in sync with emerging bioinformatics tools. Overall, AutoBA
represents a convenient tool, offering robustness and adaptability for complex
omics data analysis.


  We describe a web-based tool, MakeSBML (https://sys-bio.github.io/makesbml/),
that provides an installation-free application for creating, editing, and
searching the Biomodels repository for SBML-based models. MakeSBML is a
client-based web application that translates models expressed in human-readable
Antimony to the System Biology Markup Language (SBML) and vice-versa. Since
MakeSBML is a web-based application it requires no installation on the user's
part. Currently, MakeSBML is hosted on a GitHub page where the client-based
design makes it trivial to move to other hosts. This model for software
deployment also reduces maintenance costs since an active server is not
required. The SBML modeling language is often used in systems biology research
to describe complex biochemical networks and makes reproducing models much
easier. However, SBML is designed to be computer-readable, not human-readable.
We therefore employ the human-readable Antimony language to make it easy to
create and edit SBML models.


  The protein NLRP3 and its complexes are associated with an array of
inflammatory pathologies, among which neurodegenerative, autoimmune, and
metabolic diseases. Targeting the NLRP3 inflammasome represents a promising
strategy for easing the symptoms of pathologic neuroinflammation. When the
inflammasome is activated, NLRP3 undergoes a conformational change triggering
the production of pro-inflammatory cytokines IL-1\beta and IL-18, as well as
cell death by pyroptosis. NLRP3 nucleotide-binding and oligomerization (NACHT)
domain plays a crucial role in this function by binding and hydrolysing ATP and
is primarily responsible, together with conformational transitions involving
the PYD domain, for the complex-assembly process. Allosteric ligands proved
able to induce NLRP3 inhibition. Herein, we examine the origins of allosteric
inhibition of NLRP3. Through the use of molecular dynamics (MD) simulations and
advanced analysis methods, we provide molecular-level insights into how
allosteric binding affects protein structure and dynamics, remodelling of the
conformational ensembles populated by the protein, with key reverberations on
how NLRP3 is preorganized for assembly and ultimately function. The data are
used to develop a Machine Learning model to define the protein as Active or
Inactive, only based on the analysis of its internal dynamics. We propose this
model as a novel tool to select allosteric ligands.


  Hexapeptides are increasingly applied as model systems for studying the
amyloidogenecity properties of oligo- and polypeptides. It is possible to
construct 64 million different hexapeptides from the twenty proteinogenic amino
acid residues. Today's experimental amyloid databases contain only a fraction
of these annotated hexapeptides. For labeling all the possible hexapeptides as
"amyloidogenic" or "non-amyloidogenic" there exist several computational
predictors with good accuracies. It may be of interest to define and study a
simple graph structure on the 64 million hexapeptides as nodes when two
hexapeptides are connected by an edge if they differ by only a single residue.
For example, in this graph, HIKKLM is connected to AIKKLM, or HIKKNM, or
HIKKLC, but it is not connected with an edge to VVKKLM or HIKNPM. In the
present contribution, we consider our previously published artificial
intelligence-based tool, the Budapest Amyloid Predictor (BAP for short), and
demonstrate a spectacular property of this predictor in the graph defined
above. We show that for any two hexapeptides predicted to be "amyloidogenic" by
the BAP predictor, there exists an easily constructible path of length at most
6 that passes through neighboring hexapeptides all predicted to be
"amyloidogenic" by BAP. For example, the predicted amyloidogenic ILVWIW and
FWLCYL hexapeptides can be connected through the length-6 path
ILVWIW-IWVWIW-IWVCIW-IWVCIL-FWVCIL-FWLCIL-FWLCYL in such a way that the
neighbors differ in exactly one residue, and all hexapeptides on the path are
predicted to be amyloidogenic by BAP. The symmetric statement also holds for
non-amyloidogenic hexapeptides. It is noted that the mentioned property of the
Budapest Amyloid Predictor \url{https://pitgroup.org/bap} is not proprietary;
it is also true for any linear Support Vector Machine (SVM)-based predictors.


  Motivation: We explored how explainable AI (XAI) can help to shed light into
the inner workings of neural networks for protein function prediction, by
extending the widely used XAI method of integrated gradients such that latent
representations inside of transformer models, which were finetuned to Gene
Ontology term and Enzyme Commission number prediction, can be inspected too.
Results: The approach enabled us to identify amino acids in the sequences that
the transformers pay particular attention to, and to show that these relevant
sequence parts reflect expectations from biology and chemistry, both in the
embedding layer and inside of the model, where we identified transformer heads
with a statistically significant correspondence of attribution maps with ground
truth sequence annotations (e.g., transmembrane regions, active sites) across
many proteins. Availability and Implementation: Source code can be accessed at
https://github.com/markuswenzel/xai-proteins .


  Animal communication is frequently studied with conventional network
representations that link pairs of individuals who interact, for example,
through vocalisation. However, acoustic signals often have multiple
simultaneous receivers, or receivers integrate information from multiple
signallers, meaning these interactions are not dyadic. Additionally, non-dyadic
social structures often shape an individual's behavioural response to vocal
communication. Recently, major advances have been made in the study of these
non-dyadic, higher-order networks (e.g., hypergraphs and simplicial complexes).
Here, we show how these approaches can provide new insights into vocal
communication through three case studies that illustrate how higher-order
network models can: a) alter predictions made about the outcome of
vocally-coordinated group departures; b) generate different patterns of song
synchronisation than models that only include dyadic interactions; and c)
inform models of cultural evolution of vocal communication. Together, our
examples highlight the potential power of higher-order networks to study animal
vocal communication. We then build on our case studies to identify key
challenges in applying higher-order network approaches in this context and
outline important research questions these techniques could help answer.


  A ChatGPT-like system for drug compounds could be a game-changer in
pharmaceutical research, accelerating drug discovery, enhancing our
understanding of structure-activity relationships, guiding lead optimization,
aiding drug repurposing, reducing the failure rate, and streamlining clinical
trials. In this work, we make an initial attempt towards enabling ChatGPT-like
capabilities on drug molecule graphs, by developing a prototype system
DrugChat. DrugChat works in a similar way as ChatGPT. Users upload a compound
molecule graph and ask various questions about this compound. DrugChat will
answer these questions in a multi-turn, interactive manner. The DrugChat system
consists of a graph neural network (GNN), a large language model (LLM), and an
adaptor. The GNN takes a compound molecule graph as input and learns a
representation for this graph. The adaptor transforms the graph representation
produced by the GNN into another representation that is acceptable to the LLM.
The LLM takes the compound representation transformed by the adaptor and users'
questions about this compound as inputs and generates answers. All these
components are trained end-to-end. To train DrugChat, we collected instruction
tuning datasets which contain 10,834 drug compounds and 143,517 question-answer
pairs. The code and data is available at
\url{https://github.com/UCSD-AI4H/drugchat}