  The body plan of the fruit fly is determined by the expression of just a
handful of genes. We show that the spatial patterns of expression for several
of these genes scale precisely with the size of the embryo. Concretely,
discrete positional markers such as the peaks in striped patterns have absolute
positions along the anterior-posterior axis that are proportional to embryo
length, with better than 1% accuracy. Further, the information (in bits) that
graded patterns of expression provide about position can be decomposed into
information about fractional or scaled position and information about absolute
position or embryo length; all of the available information is about scaled
position, again with ~1% accuracy. These observations suggest that the
underlying genetic network exhibits scale invariance in a deeper mathematical
sense. Taking this mathematical statement seriously requires that the network
dynamics have a zero mode, which connects to many other observations on this
system.


  Learning algorithms based on backpropagation have enabled transformative
technological advances but alternatives based on local energy-based rules offer
benefits in terms of biological plausibility and decentralized training. A
broad class of such local learning rules involve \textit{contrasting} a clamped
configuration with the free, spontaneous behavior of the system. However,
comparisons of clamped and free configurations require explicit memory or
switching between Hebbian and anti-Hebbian modes. Here, we show how a simple
form of implicit non-equilibrium memory in the update dynamics of each
``synapse'' of a network naturally allows for contrastive learning. During
training, free and clamped behaviors are shown in sequence over time using a
sawtooth-like temporal protocol that breaks the symmetry between those two
behaviors when combined with non-equilibrium update dynamics at each synapse.
We show that the needed dynamics is implicit in integral feedback control,
broadening the range of physical and biological systems naturally capable of
contrastive learning. Finally, we show that non-equilibrium dissipation
improves learning quality and determine the Landauer energy cost of contrastive
learning through physical dynamics.


  The tasks of designing messenger RNAs and non-coding RNAs are discrete
optimization problems, and several versions of these problems are NP-hard. As
an alternative to commonly used local search methods, we formulate these
problems as continuous optimization and develop a general framework for this
optimization based on a new concept of "expected partition function". The basic
idea is to start with a distribution over all possible candidate sequences, and
extend the objective function from a sequence to a distribution. We then use
gradient descent-based optimization methods to improve the extended objective
function, and the distribution will gradually shrink towards a one-hot sequence
(i.e., a single sequence). We consider two important case studies within this
framework, the mRNA design problem optimizing for partition function (i.e.,
ensemble free energy) and the non-coding RNA design problem optimizing for
conditional (i.e., Boltzmann) probability. In both cases, our approach
demonstrate promising preliminary results. We make our code available at
https://github.com/KuNyaa/RNA_Design_codebase.


  Motivated by the question of how biological systems maintain homeostasis in
changing environments, Shinar and Feinberg introduced in 2010 the concept of
absolute concentration robustness (ACR). A biochemical system exhibits ACR in
some species if the steady-state value of that species does not depend on
initial conditions. Thus, a system with ACR can maintain a constant level of
one species even as the environment changes. Despite a great deal of interest
in ACR in recent years, the following basic question remains open: How can we
determine quickly whether a given biochemical system has ACR? Although various
approaches to this problem have been proposed, we show that they are
incomplete. Accordingly, we present new methods for deciding ACR, which harness
computational algebra. We illustrate our results on several biochemical
signaling networks.


  Understanding the evolution of cooperation is pivotal in biology and social
science. Public resources sharing is a common scenario in the real world. In
our study, we explore the evolutionary dynamics of cooperation on a regular
graph with degree $k$, introducing the presence of a third strategy, namely the
benevolence, who does not evolve over time, but provides a fixed benefit to all
its neighbors. We find that the presence of the benevolence can foster the
development of cooperative behavior and it follows a simple rule: $b/c > k -
p_S(k-1)$. Our results provide new insights into the evolution of cooperation
in structured populations.


  The time evolution of a contagious viral disease is modeled as the dynamic
progression of different classes of populations that interact pairwise, aiming
to improve their condition with respect to a given target. This evolutionary
mechanism is based on binary interactions between agents, designed to enhance
their utility mutually. To achieve this goal, we introduce kinetic equations of
Boltzmann type to describe the time evolution of the probability distributions
of the agent system undergoing binary interactions. The fundamental idea is to
describe these interactions using principles from price theory, particularly by
employing Cobb-Douglas utility functions for the binary exchange and the
Edgeworth box to depict the common exchange area where utility increases for
both agents.


  Walsh coefficients have been applied extensively to biallelic systems for
quantifying pairwise and higher order epistasis, in particular for
demonstrating the empirical importance of higher order interactions. Circuits,
or minimal dependence relations, and related approaches that use triangulations
of polytopes have also been applied to biallelic systems. Here we provide
biological interpretations of Walsh coefficients for several alleles, and
discuss circuits in the same general setting.


  Exploring methods and techniques of machine learning (ML) to address specific
challenges in various fields is essential. In this work, we tackle a problem in
the domain of Cheminformatics; that is, providing a suitable solution to aid in
predicting the activity of a chemical compound to the best extent possible. To
address the problem at hand, this study conducts experiments on 100 different
combinations of existing techniques. These solutions are then selected based on
a set of criteria that includes the G-means, F1-score, and AUC metrics. The
results have been tested on a dataset of about 10,000 chemical compounds from
PubChem that have been classified according to their activity


  To improve the environment of the ecosystem, China has implemented the
Green-forGrain Program for two decades, which has resulted in an imbalance
among ecology.economy and food. This study focuses on the "ecology-food"
imbalance problem.taking Sichuan-Chongqing Region as an example, to set up
future scenarios topredicate the distribution of ESs. We first forecast land
use/cover change in 2050under four different scenarios: Natural Development
Scenarios; Arable LandConservation Scenarios; Ecological Priority Scenarios;
Ecology-Arable LandHarmonization Scenarios. Then we assess changes in five ESs:
habitat quality ,cropproduction, soil conservation, water yield, and carbon
storage from 1990 to 2020 and2050. Finally, we reveal the spacial distribution
of ESs. The following conclusions areobtained: (1) From 1990-2020, CS, SC, and
HQ reveal an increasing trend with growthrates of 1.68%, 0.08%, and 0.46%: CP
reveals a reduce rate of 2.75% . (2) S4 has anincrease in arable land, and CP
has increased by 7.56% compared to S1, reversingthe trend of reduced CP under
S1. (3) The high-high anomalies area of CP under S4 isoasically the same as
that under S2, which proves that S4 is a scenario policy that canbe referred to
for future development.


  We apply a computational modelling approach to investigate the relative
effectiveness of general isolation practices for mitigation of COVID-19
outbreaks in residential care facilities. Our study focuses on policies
intended to reduce contact between residents, without regard to confirmed
infection status. Despite the ubiquity of such policies, and their
controversial association with adverse physical and mental health outcomes,
little evidence exists evaluating their effectiveness at mitigating outbreaks.
Through detailed simulations of COVID-19 outbreaks in residential care
facilities, our results demonstrate that general isolation of residents
provides little additional impact beyond what is achievable through isolation
of confirmed cases and deployment of personal protective equipment.


  Molecular property optimization (MPO) problems are inherently challenging
since they are formulated over discrete, unstructured spaces and the labeling
process involves expensive simulations or experiments, which fundamentally
limits the amount of available data. Bayesian optimization (BO) is a powerful
and popular framework for efficient optimization of noisy, black-box objective
functions (e.g., measured property values), thus is a potentially attractive
framework for MPO. To apply BO to MPO problems, one must select a structured
molecular representation that enables construction of a probabilistic surrogate
model. Many molecular representations have been developed, however, they are
all high-dimensional, which introduces important challenges in the BO process
-- mainly because the curse of dimensionality makes it difficult to define and
perform inference over a suitable class of surrogate models. This challenge has
been recently addressed by learning a lower-dimensional encoding of a SMILE or
graph representation of a molecule in an unsupervised manner and then
performing BO in the encoded space. In this work, we show that such methods
have a tendency to "get stuck," which we hypothesize occurs since the mapping
from the encoded space to property values is not necessarily well-modeled by a
Gaussian process. We argue for an alternative approach that combines numerical
molecular descriptors with a sparse axis-aligned Gaussian process model, which
is capable of rapidly identifying sparse subspaces that are most relevant to
modeling the unknown property function. We demonstrate that our proposed method
substantially outperforms existing MPO methods on a variety of benchmark and
real-world problems. Specifically, we show that our method can routinely find
near-optimal molecules out of a set of more than $>100$k alternatives within
100 or fewer expensive queries.


  Consumer-resource dynamics is central in determining biomass transport across
ecosystems. The assumptions of mass action, chemostatic conditions and
stationarity in stochastic feeding dynamics lead to Holling type II functional
responses, whose use is widespread in macroscopic models of population
dynamics. However, to be useful for parameter inference, stochastic population
models need to be identifiable, this meaning that model parameters can be
uniquely inferred from a large number of model observations. In this
contribution we study parameter identifiability in a multi-resource
consumer-resource model, for which we can obtain the steady-state and
out-of-equilibrium probability distributions of predator's abundances by
analytically solving the master equation. Based on these analytical results, we
can conduct in silico experiments by tracking the abundance of consumers that
are either searching for or handling prey, data then used for maximum
likelihood parameter estimation. We show that, when model observations are
recorded out of equilibrium, feeding parameters are truly identifiable, whereas
if sampling is done at stationarity, only ratios of rates can be inferred from
data. We discuss the implications of our results when inferring parameters of
general dynamical models.


  Minimizing data storage poses a significant challenge in large-scale
metagenomic projects. In this paper, we present a new method for improving the
encoding of FASTQ files generated by metagenomic sequencing. This method
incorporates metagenomic classification followed by a recursive filter for
clustering reads by DNA sequence similarity to improve the overall
reference-free compression. In the results, we show an overall improvement in
the compression of several datasets. As hypothesized, we show a progressive
compression gain for higher coverage depth and number of identified species.
Additionally, we provide an implementation that is freely available at
https://github.com/cobilab/mizar and can be customized to work with other FASTQ
compression tools.


  Design of new drugs is a challenging process: a candidate molecule should
satisfy multiple conditions to act properly and make the least side-effect --
perfect candidates selectively attach to and influence only targets, leaving
off-targets intact. The amount of experimental data about various properties of
molecules constantly grows, promoting data-driven approaches. However, the
applicability of typical predictive machine learning techniques can be
substantially limited by a lack of experimental data about a particular target.
For example, there are many known Thrombin inhibitors (acting as
anticoagulants), but a very limited number of known Protein C inhibitors
(coagulants). In this study, we present our approach to suggest new inhibitor
candidates by building an effective representation of chemical space. For this
aim, we developed a deep learning model -- autoencoder, trained on a large set
of molecules in the SMILES format to map the chemical space. Further, we
applied different sampling strategies to generate novel coagulant candidates.
Symmetrically, we tested our approach on anticoagulant candidates, where we
were able to predict their inhibition towards Thrombin. We also compare our
approach with MegaMolBART -- another deep learning generative model, but
exploiting similar principles of navigation in a chemical space.


  Embedding sequential computations in biochemical environments is challenging
because the computations are carried out by chemical reactions, which are
inherently disordered. In this paper we apply modular design to specific
calculations through chemical reactions and provide a design scheme of
biochemical oscillator models in order to generate periodical species for the
order regulation of these reaction modules. We take the case of arbitrary
multi-module regulation into consideration, analyze the main errors in the
regulation process under \textit{mass-action kinetics} and demonstrate our
design scheme under existing synthetic biochemical oscillator models.


  Despite its remoteness from other continents, the Antarctic region cannot
escape the aftermath of human activities as it is highly influenced by
anthropogenic impacts that occur both in the regional and global context.
Contamination by microplastics, mostly caused by the improper disposal of
plastic waste, is widely recognised as a serious environmental threat due to
its ubiquity. In recent years, most researchers have focused on microplastic
pollution in the marine ecosystem of Antarctica, while pollution in the
terrestrial environment continues to be neglected. This study was conducted to
investigate the ability of Antarctic soil bacteria to use polypropylene (PP)
microplastics as the sole carbon source. Bushnell Haas (BH) medium inoculated
with bacteria and supplemented PP-microplastics as the sole carbon source was
used in the utilisation test. In this study, the growth response of Dermacoccus
sp. strain AYDL3 was assessed after exposure to PP-microplastics in a basal
medium for 40 days. The weight reduction of the polymer was determined to
further support the growth response. The highest and lowest weight loss
percentages were observed on day 20 (23.0%) and day 10 (7.75%), respectively.
Fourier transforms infrared (FTIR) spectroscopy and scanning electron
microscopy (SEM) analyses were used to confirm the utilisation of
PP-microplastics by strain AYDL3. Results indicate that the soil bacteria
possess a mechanism for breaking down microplastics allowing them to utilise
plastics as energy sources without any pre-treatment. This emphasises the
significance of these soil bacteria to adapt and subsequently manage the
plastic fragments in the soil in the future.


  Anticancer peptides (ACPs) are a class of molecules that have gained
significant attention in the field of cancer research and therapy. ACPs are
short chains of amino acids, the building blocks of proteins, and they possess
the ability to selectively target and kill cancer cells. One of the key
advantages of ACPs is their ability to selectively target cancer cells while
sparing healthy cells to a greater extent. This selectivity is often attributed
to differences in the surface properties of cancer cells compared to normal
cells. That is why ACPs are being investigated as potential candidates for
cancer therapy. ACPs may be used alone or in combination with other treatment
modalities like chemotherapy and radiation therapy. While ACPs hold promise as
a novel approach to cancer treatment, there are challenges to overcome,
including optimizing their stability, improving selectivity, and enhancing
their delivery to cancer cells, continuous increasing in number of peptide
sequences, developing a reliable and precise prediction model. In this work, we
propose an efficient transformer-based framework to identify anticancer
peptides for by performing accurate a reliable and precise prediction model.
For this purpose, four different transformer models, namely ESM, ProtBert,
BioBERT, and SciBERT are employed to detect anticancer peptides from amino acid
sequences. To demonstrate the contribution of the proposed framework, extensive
experiments are carried on widely-used datasets in the literature, two versions
of AntiCp2, cACP-DeepGram, ACP-740. Experiment results show the usage of
proposed model enhances classification accuracy when compared to the
state-of-the-art studies. The proposed framework, ESM, exhibits 96.45 of
accuracy for AntiCp2 dataset, 97.66 of accuracy for cACP-DeepGram dataset, and
88.51 of accuracy for ACP-740 dataset, thence determining new state-of-the-art.


  Protein folding and evolution are intimately linked phenomena. Here, we
revisit the concept of exons as potential protein folding modules across 38
abundant and conserved protein families. Taking advantage of genomic
exon-intron organization and extensive protein sequence data, we explore exon
boundary conservation and assess their foldon-like behavior using energy
landscape theoretic measurements. We found deviations in exon size distribution
from exponential decay indicating selection in evolution. We describe that
there is a pronounced independent foldability of segments corresponding to
conserved exons, supporting the exon-foldon correspondence. We further develop
a systematic partitioning of protein domains using exon boundary hot spots,
unveiling minimal common exons consisting of uninterrupted alpha and/or beta
elements for the majority but not all of the studied families.


  Microbial ecosystems are commonly modeled by fixed interactions between
species in steady exponential growth states. However, microbes often modify
their environments so strongly that they are forced out of the exponential
state into stressed or non-growing states. Such dynamics are typical of
ecological succession in nature and serial-dilution cycles in the laboratory.
Here, we introduce a phenomenological model, the Community State model, to gain
insight into the dynamic coexistence of microbes due to changes in their
physiological states. Our model bypasses specific interactions (e.g., nutrient
starvation, stress, aggregation) that lead to different combinations of
physiological states, referred to collectively as "community states", and
modeled by specifying the growth preference of each species along a global
ecological coordinate, taken here to be the total community biomass density. We
identify three key features of such dynamical communities that contrast starkly
with steady-state communities: increased tolerance of community diversity to
fast growth rates of species dominating different community states, enhanced
community stability through staggered dominance of different species in
different community states, and increased requirement on growth dominance for
the inclusion of late-growing species. These features, derived explicitly for
simplified models, are proposed here to be principles aiding the understanding
of complex dynamical communities. Our model shifts the focus of ecosystem
dynamics from bottom-up studies based on idealized inter-species interaction to
top-down studies based on accessible macroscopic observables such as growth
rates and total biomass density, enabling quantitative examination of
community-wide characteristics.


  Cellular ingredient concentrations can be stabilized by adjusting generation
and consumption rates through multiple pathways. To explore the portion of
cellular metabolism equipped with multiple pathways, we categorize individual
metabolic reactions and compounds as viable or inviable: A compound is viable
if processed by two or more reactions, and a reaction is viable if all of its
substrates and products are viable. Using this classification, we identify the
maximal subnetwork of viable nodes, referred to as the {\it viable core}, in
bipartite metabolic networks across thousands of species. The obtained viable
cores are remarkably larger than those in degree-preserving randomized
networks, while their broad degree distributions commonly enable the viable
cores to shrink gradually as reaction nodes are deleted. We demonstrate that
the positive degree-degree correlations of the empirical networks may underlie
the enlarged viable cores compared to the randomized networks. By investigating
the relation between degree and cross-species frequency of metabolic compounds
and reactions, we elucidate the evolutionary origin of the correlations.


  Denoising diffusion models have shown great potential in multiple research
areas. Existing diffusion-based generative methods on de novo 3D molecule
generation face two major challenges. Since majority heavy atoms in molecules
allow connections to multiple atoms through single bonds, solely using
pair-wise distance to model molecule geometries is insufficient. Therefore, the
first one involves proposing an effective neural network as the denoising
kernel that is capable to capture complex multi-body interatomic relationships
and learn high-quality features. Due to the discrete nature of graphs,
mainstream diffusion-based methods for molecules heavily rely on predefined
rules and generate edges in an indirect manner. The second challenge involves
accommodating molecule generation to diffusion and accurately predicting the
existence of bonds. In our research, we view the iterative way of updating
molecule conformations in diffusion process is consistent with molecular
dynamics and introduce a novel molecule generation method named
Geometric-Facilitated Molecular Diffusion (GFMDiff). For the first challenge,
we introduce a Dual-Track Transformer Network (DTN) to fully excevate global
spatial relationships and learn high quality representations which contribute
to accurate predictions of features and geometries. As for the second
challenge, we design Geometric-Facilitated Loss (GFLoss) which intervenes the
formation of bonds during the training period, instead of directly embedding
edges into the latent space. Comprehensive experiments on current benchmarks
demonstrate the superiority of GFMDiff.


  This paper focuses on learning representation on the whole graph level in an
unsupervised manner. Learning graph-level representation plays an important
role in a variety of real-world issues such as molecule property prediction,
protein structure feature extraction, and social network analysis. The
mainstream method is utilizing contrastive learning to facilitate graph feature
extraction, known as Graph Contrastive Learning (GCL). GCL, although effective,
suffers from some complications in contrastive learning, such as the effect of
false negative pairs. Moreover, augmentation strategies in GCL are weakly
adaptive to diverse graph datasets. Motivated by these problems, we propose a
novel framework called Structure Knowledge Refinement (SKR) which uses data
structure to determine the probability of whether a pair is positive or
negative. Meanwhile, we propose an augmentation strategy that naturally
preserves the semantic meaning of the original data and is compatible with our
SKR framework. Furthermore, we illustrate the effectiveness of our SKR
framework through intuition and experiments. The experimental results on the
tasks of graph-level classification demonstrate that our SKR framework is
superior to most state-of-the-art baselines.


  De Bruijn graph is one of the most important data structures used in de-novo
genome assembly algorithms, especially for NGS data. There is a growing need
for parallel data structures and algorithms due to the increasing number of
cores in modern computers. The assembly task is an indispensable step in
sequencing genomes of new organisms and studying structural genomic changes. In
recent years, the dynamic development of next-generation sequencing (NGS)
methods raises hopes for making whole-genome sequencing a fast and reliable
tool used, for example, in medical diagnostics. However, this is hampered by
the slowness and computational requirements of the current processing
algorithms, which raises the need to develop more efficient algorithms. One
possible approach, still little explored, is the use of quantum computing.
  We created the lock-free version of the de Bruijn graph, as well as a
lock-free algorithm to build such graph from reads. Our algorithm and data
structures are developed to use parallel threads of execution and do not use
mutexes or other locking mechanisms, instead, we used only compare-and-swap
instruction and other atomic operations. It makes our algorithm very fast and
efficiently scaling.
  The presented article depicts the new lock-free de Bruijn graph data
structure with a graph build algorithm. We developed a C++ library and tested
its performance to depict its high speed and scalability compared to other
available tools.


  Large Language Models (LLMs), such as ChatGPT, have taken the world by storm
and have passed certain forms of the Turing test. However, LLMs are not limited
to human language and analyze sequential data, such as DNA, protein, and gene
expression. The resulting foundation models can be repurposed to identify the
complex patterns within the data, resulting in powerful, multi-purpose
prediction tools able to explain cellular systems. This review outlines the
different types of LLMs and showcases their recent uses in biology. Since LLMs
have not yet been embraced by the plant community, we also cover how these
models can be deployed for the plant kingdom.


  Species distribution modeling is a highly versatile tool for understanding
the intricate relationship between environmental conditions and species
occurrences. However, the available data often lacks information on confirmed
species absence and is limited to opportunistically sampled, presence-only
observations. To overcome this limitation, a common approach is to employ
pseudo-absences, which are specific geographic locations designated as negative
samples. While pseudo-absences are well-established for single-species
distribution models, their application in the context of multi-species neural
networks remains underexplored. Notably, the significant class imbalance
between species presences and pseudo-absences is often left unaddressed.
Moreover, the existence of different types of pseudo-absences (e.g., random and
target-group background points) adds complexity to the selection process.
Determining the optimal combination of pseudo-absences types is difficult and
depends on the characteristics of the data, particularly considering that
certain types of pseudo-absences can be used to mitigate geographic biases. In
this paper, we demonstrate that these challenges can be effectively tackled by
integrating pseudo-absences in the training of multi-species neural networks
through modifications to the loss function. This adjustment involves assigning
different weights to the distinct terms of the loss function, thereby
addressing both the class imbalance and the choice of pseudo-absence types.
Additionally, we propose a strategy to set these loss weights using spatial
block cross-validation with presence-only data. We evaluate our approach using
a benchmark dataset containing independent presence-absence data from six
different regions and report improved results when compared to competing
approaches.


  From a mass action law approach, we modelled and calibrated the distributions
of the seven-stripe patterns of Even-skipped (\textit{Eve}) and Fushi-tarazu
(\textit{Ftz}) proteins in the Drosophila embryo. The regulators of
\textit{Eve} are stripe-specific enhancers from the gap family of proteins. We
achieved remarkable data matching of the \textit{Eve} stripe pattern. We have
identified the putative repressive combinations for three \textit{Ftz}
enhancers, and we have explored the relationship between \textit{Eve} and
\textit{Ftz} for complementary patterns. Extended work to infer the Wingless
(\textit{Wg}) 14 stripe pattern from \textit{Eve} and \textit{Ftz} enhancers
have been proposed, given clues for the hierarchical structure of
\textit{Drosphila}'s genetic expression network during early development prior
to cellularisation.


  Artificial intelligence based on neural networks has made significant
progress. However, there are concerns about the reliability and security of
this approach due to its lack of transparency. This is the black box problem of
AI. Here we show how this problem can be solved using symbolic AI, which has a
transparent white box nature. The widespread use of symbolic AI is hindered by
the opacity of mathematical models and natural language terms, the lack of a
unified ontology, and the combinatorial explosion of search options. To solve
the AI black box problem and to implement general-purpose symbolic AI, we
propose to use deterministic logic cellular automata with rules based on first
principles of the general theory of the relevant domain. In this case, the
general theory of the relevant domain plays the role of a knowledge base for
the cellular automaton inference. A cellular automaton implements automatic
parallel logical inference at three levels of organization of a complex system.
Our verification of several ecological hypotheses provides a successful
precedent for the implementation of white-box AI. Finally, we discuss a program
for creating a general-purpose symbolic AI capable of processing knowledge and
ensuring the reliability and safety of automated decisions.


  Dimensional reduction techniques have long been used to visualize the
structure and geometry of high dimensional data. However, most widely used
techniques are difficult to interpret due to nonlinearities and opaque
optimization processes. Here we present a specific graph based construction for
dimensionally reducing continuous stochastic systems with multiplicative noise
moving under the influence of a potential. To achieve this, we present a
specific graph construction which generates the Fokker-Planck equation of the
stochastic system in the continuum limit. The eigenvectors and eigenvalues of
the normalized graph Laplacian are used as a basis for the dimensional
reduction and yield a low dimensional representation of the dynamics which can
be used for downstream analysis such as spectral clustering. We focus on the
use case of single cell RNA sequencing data and show how current diffusion map
implementations popular in the single cell literature fit into this framework.


  Molecular property prediction refers to the task of labeling molecules with
some biochemical properties, playing a pivotal role in the drug discovery and
design process. Recently, with the advancement of machine learning, deep
learning-based molecular property prediction has emerged as a solution to the
resource-intensive nature of traditional methods, garnering significant
attention. Among them, molecular representation learning is the key factor for
molecular property prediction performance. And there are lots of
sequence-based, graph-based, and geometry-based methods that have been
proposed. However, the majority of existing studies focus solely on one
modality for learning molecular representations, failing to comprehensively
capture molecular characteristics and information. In this paper, a novel
multi-modal representation learning model, which integrates the sequence,
graph, and geometry characteristics, is proposed for molecular property
prediction, called SGGRL. Specifically, we design a fusion layer to fusion the
representation of different modalities. Furthermore, to ensure consistency
across modalities, SGGRL is trained to maximize the similarity of
representations for the same molecule while minimizing similarity for different
molecules. To verify the effectiveness of SGGRL, seven molecular datasets, and
several baselines are used for evaluation and comparison. The experimental
results demonstrate that SGGRL consistently outperforms the baselines in most
cases. This further underscores the capability of SGGRL to comprehensively
capture molecular information. Overall, the proposed SGGRL model showcases its
potential to revolutionize molecular property prediction by leveraging
multi-modal representation learning to extract diverse and comprehensive
molecular insights. Our code is released at
https://github.com/Vencent-Won/SGGRL.


  During the COVID-19 pandemic, a major driver of new surges has been the
emergence of new variants. When a new variant emerges in one or more countries,
other nations monitor its spread in preparation for its potential arrival. The
impact of the variant and the timing of epidemic peaks in a country highly
depend on when the variant arrives. The current methods for predicting the
spread of new variants rely on statistical modeling, however, these methods
work only when the new variant has already arrived in the region of interest
and has a significant prevalence. The question arises: Can we predict when (and
if) a variant that exists elsewhere will arrive in a given country and reach a
certain prevalence? We propose a variant-dynamics-informed Graph Neural Network
(GNN) approach. First, We derive the dynamics of variant prevalence across
pairs of regions (countries) that applies to a large class of epidemic models.
The dynamics suggest that ratios of variant proportions lead to simpler
patterns. Therefore, we use ratios of variant proportions along with some
parameters estimated from the dynamics as features in a GNN. We develop a
benchmarking tool to evaluate variant emergence prediction over 87 countries
and 36 variants. We leverage this tool to compare our GNN-based approach
against our dynamics-only model and a number of machine learning models.
Results show that the proposed dynamics-informed GNN method retrospectively
outperforms all the baselines, including the currently pervasive framework of
Physics-Informed Neural Networks (PINNs) that incorporates the dynamics in the
loss function.


  RNA secondary structure is modeled with the novel arbitrary-order hidden
Markov model ({\alpha}-HMM). The {\alpha}-HMM extends over the traditional HMM
with capability to model stochastic events that may be in influenced by
historically distant ones, making it suitable to account for long-range
canonical base pairings between nucleotides, which constitute the RNA secondary
structure. Unlike previous heavy-weight extensions over HMM, the {\alpha}-HMM
has the flexibility to apply restrictions on how one event may influence
another in stochastic processes, enabling efficient prediction of RNA secondary
structure including pseudoknots.


  The years 2020 and 2021 were characterized by the covid-19 pandemic. The true
impact of the pandemic on populations health and life has still to be
completely discerned. Objective difficulties to account deaths due to covid-19,
determined the risk of underestimating such impact, hence official mortality
data started to be examined. The main objective of this work is to discern the
true impact of covid-19 in EU during the crucial years 2020 and 2021. The
mortality trends, already increasing in EU in pre-pandemic years, are modelled.
The excess mortality attributable to covid-19, on yearly basis, is estimated
via a novel strategy combining datasets from different official sources.
Considering demographic and geographic factors, new indexes are formulated,
which allow rankings of EU countries, and socio-political and economic
reflections. The new indexes here formulated provide a different perspective to
see the relative impact of the pandemic in EU. This work, also in line with the
conclusions of previously published authoritative papers on demographic studies
of excess mortality, represents an original methodology, which timely
implemented, can be at the service of public decision makers for future
emerging needs.


  Single-cell RNA sequencing (scRNA-seq) data are important for studying the
biology of development or diseases at single-cell level. To better understand
the properties of the data, to build controlled benchmark data for testing
downstream methods, and to augment data when collecting sufficient real data is
challenging, generative models have been proposed to computationally generate
synthetic scRNA-seq data. However, the data generated with current models are
not very realistic yet, especially when we need to generate data with
controlled conditions. In the meantime, the Diffusion models have shown their
power in generating data in computer vision at high fidelity, providing a new
opportunity for scRNA-seq generation.
  In this study, we developed scDiffusion, a diffusion-based model to generate
high-quality scRNA-seq data with controlled conditions. We designed multiple
classifiers to guide the diffusion process simultaneously, enabling scDiffusion
to generate data under multiple condition combinations. We also proposed a new
control strategy called Gradient Interpolation. This strategy allows the model
to generate continuous trajectories of cell development from a given cell
state.
  Experiments showed that scDiffusion can generate single-cell gene expression
data closely resembling real scRNA-seq data, surpassing state-of-the-art models
in multiple metrics. Also, scDiffusion can conditionally produce data on
specific cell types including rare cell types. Furthermore, we could use the
multiple-condition generation of scDiffusion to generate cell type that was out
of the training data. Leveraging the Gradient Interpolation strategy, we
generated a continuous developmental trajectory of mouse embryonic cells. These
experiments demonstrate that scDiffusion is a powerful tool for augmenting the
real scRNA-seq data and can provide insights into cell fate research.


  During the founding of the field of quantitative genetics, Fisher formulated
in 1918 his infinitesimal model that provided a novel mathematical framework to
describe the Mendelian transmission of quantitative traits. If the
infinitely-many genes in that model are assumed to segregate independently
during reproduction, corresponding to having no linkage, directional selection
asymptotically leads to a constant genetic gain at each generation. However, in
reality genes are subject to linkage because they lie on chromosomes and thus
segregate in a correlated way. Various approximations have been used in the
past to study that more realistic case of the infinitesimal model and they lead
to the prediction that the asymptotic gain per generation is modestly
decreased. Here we provide an exact solution to the problem where genes lie on
continuous chromosomes. Surprisingly, the consequences of genetic linkage are
in fact rather singular, changing the nature of the long term gain per
generation: the asymptotic gain vanishes rather than being simply decreased.
Nevertheless, the rate at which the per-generation gain tends to zero is slow
enough for the total gain, accumulated over generations, to be unbounded.


  The Boltzmann distribution of a protein provides a roadmap to all of its
functional states. Normalizing flows are a promising tool for modeling this
distribution, but current methods are intractable for typical pharmacological
targets; they become computationally intractable due to the size of the system,
heterogeneity of intra-molecular potential energy, and long-range interactions.
To remedy these issues, we present a novel flow architecture that utilizes
split channels and gated attention to efficiently learn the conformational
distribution of proteins defined by internal coordinates. We show that by
utilizing a 2-Wasserstein loss, one can smooth the transition from maximum
likelihood training to energy-based training, enabling the training of
Boltzmann Generators for macromolecules. We evaluate our model and training
strategy on villin headpiece HP35(nle-nle), a 35-residue subdomain, and protein
G, a 56-residue protein. We demonstrate that standard architectures and
training strategies, such as maximum likelihood alone, fail while our novel
architecture and multi-stage training strategy are able to model the
conformational distributions of protein G and HP35.


  Consider the evolutionary dynamics of an ensemble of genetic variations
within a given population of organisms as a function of time. Such a system can
be described microscopically in terms of individual biological sequences or
macroscopically in terms of macroscopic or thermodynamic quantities such as
entropy of sequences, population size, etc. While the microscopic approach is
more useful for studying fine-grained details of evolution, such as individual
mutations, the macroscopic approach is better suited for modeling
coarse-grained or thermodynamic properties of evolving systems. In this paper,
we develop a macroscopic description by following the temporal dynamics of the
total Shannon entropy of sequences, denoted by S, and the average Hamming
distance between them, denoted by H. We argue that a biological system can
persist in the so-called quasi-equilibrium state for an extended period,
characterized by strong correlations between S and H, before undergoing a phase
transition to another quasi-equilibrium state. To demonstrate the results, we
conducted a statistical analysis of SARS-CoV-2 data from the United Kingdom
during the period between March, 2020 and December, 2023. In this analysis, we
identified eight quasi-equilibrium states separated from each other by seven
phase transitions. To study the temporal dynamics in the evolution of the
virus, the alignment was done separately for each month, using Mafft alignment
program, and on only relatively clean sequences, i.e. for which the fraction of
ambiguous characters is smaller than one per cent. From a purely theoretical
perspective, this allows us to systematically study various types of phase
transitions described by a discontinuous change in the thermodynamic
parameters. From a more practical point of view, the analysis can be used, for
example, as an early warning system for pandemics.


  The characterization of interactions between autophagy modifiers (Atg8-family
proteins) and their natural ligands (peptides and proteins) or small molecules
is important for a detailed understanding of selective autophagy mechanisms and
for design of potential Atg8 inhibitors that affect the autophagy processes in
cells. The fluorescence polarization (FP) assay is a rapid, cost-effective and
robust method which provides affinity and selectivity information of small
molecules and peptide ligands targeting human Atg8 proteins. This chapter
introduces the basic principles of FP assays are introduced. In addition, a
case study for FP small molecule and peptide interactions with human Atg8
proteins (LC3/GABARAPs) is described. Finally, data analysis and quality
control of FP studies are discussed for the proper calculation of Ki values of
the measured compounds.


  The success of drug discovery and development relies on the precise
prediction of molecular activities and properties. While in silico molecular
property prediction has shown remarkable potential, its use has been limited so
far to assays for which large amounts of data are available. In this study, we
use a fine-tuned large language model to integrate biological assays based on
their textual information, coupled with Barlow Twins, a Siamese neural network
using a novel self-supervised learning approach. This architecture uses both
assay information and molecular fingerprints to extract the true molecular
information. TwinBooster enables the prediction of properties of unseen
bioassays and molecules by providing state-of-the-art zero-shot learning tasks.
Remarkably, our artificial intelligence pipeline shows excellent performance on
the FS-Mol benchmark. This breakthrough demonstrates the application of deep
learning to critical property prediction tasks where data is typically scarce.
By accelerating the early identification of active molecules in drug discovery
and development, this method has the potential to help streamline the
identification of novel therapeutics.


  Epidemiological dynamics are affected by the spatial and demographic
structure of the host population. Households and neighbourhoods are known to be
important demographic structures but little is known about the epidemiological
interplay between them. Here we present a multi-scale model consisting of
neighbourhoods of households. In our analysis we focus on key parameters which
control household size, the importance of transmission within households
relative to outside of them, and the degree to which the non-household
transmission is localised within neighbourhoods. We construct the household
reproduction number $R_*$ over all neighbourhoods and derive the analytic
probability of an outbreak occurring from a single infected individual in a
specific neighbourhood. We find that reduced localisation of transmission
within neighbourhoods reduces $R_*$ when household size differs between
neighbourhoods. This effect is amplified by larger differences between
household sizes and larger divergence between transmission rates within
households and outside of them. However, the impact of neighbourhoods with
larger household sizes is mainly limited to these neighbourhoods. We consider
various surveillance scenarios and show that household size information from
the initial infectious cases is often more important than neighbourhood
information while household size and neighbourhood localisation influences the
sequence of neighbourhoods in which an outbreak is observed.


  The method developed by Michaelis and Menten was foundational in the
development of our understanding of biochemical reaction kinetics. Extended
models of metabolism encapsulated by reaction rate theory, stochastic reaction
models, and dynamic flux estimation, amongst others, address aspects of this
fundamental idea. The limitations of these approaches are well understood, and
efforts to overcome those issues so far have been plentiful but with limited
success. The known issues can be summarised as the sole dependent relation with
substrate concentration, the encapsulation of rate in a single relevant scalar,
and the subsequent lack of functional control that results from this
assumption. The Rate Control of Chaos (RCC) is a nonlinear control method that
has been shown to be effective in controlling the dynamic state of biological
oscillators based on the concept of rate limitation of the exponential growth
in chaotic systems. Extending RCC with allosteric properties allows robust
control of the enzymatic process, and replicates the Michaelis-Menten kinetics.
The emergent dynamics is robust to perturbations and noise but susceptible to
regulatory adjustments. This control method adapts the control parameters
dynamically in the presence of a ligand, and permits introduction of energy
relations into the control function. The dynamic nature of the control
eliminates the steady-state requirements and allows the modelling of
large-scale dynamic behaviour, potentially addressing issues in metabolic
disorder and failure of metabolic control.


  Food web topology and energy flow rates across food web linkages can
influence ecosystem properties such as stability. Stability predictions from
current models of energy flow are often sensitive to details in their
formulation, and their complexity makes it difficult to elucidate underlying
mechanisms of general phenomena. Here, within the maximum information entropy
inference framework (MaxEnt), we derive a simple formula for the energy flow
carried by each linkage between two adjacent trophic layers. Inputs to the
model are the topological structure of the food web and aggregate energy fluxes
entering or exiting each species node. For ecosystems with interactions
dominated by consumer-resource interactions between two trophic layers, we
construct a model of species dynamics based on the energy flow predictions from
the MaxEnt model. Mathematical analyses and simulations of the model show that
a food web topology with a higher matrix dipole moment promotes stability
against small perturbations in population sizes, where the \textit{matrix
dipole moment} is a simple nestedness metric that we introduce. Since nested
bipartite subnetworks arise naturally in food webs, our result provides an
explanation for the stability of natural communities.


  During photoheterotrophic growth on organic substrates, purple nonsulfur
photosynthetic bacteria like Rhodospirillum rubrum can acquire electrons by
multiple means, including oxidation of organic substrates, oxidation of
inorganic electron donors (e.g. H$_2$), and by reverse electron flow from the
photosynthetic electron transport chain. These electrons are stored in the form
of reduced electron-carrying cofactors (e.g. NAD(P)H and ferredoxin). The ratio
of oxidized to reduced redox cofactors (e.g. ratio of NAD(P)+:NAD(P)H), or
'redox poise` is difficult to understand or predict, as are the the cellular
processes for dissipating these reducing equivalents. Using physics-based
models that capture mass action kinetics consistent with the thermodynamics of
reactions and pathways, a range of redox conditions for heterophototrophic
growth are evaluated, from conditions in which the NADP+/NADPH levels
approached thermodynamic equilibrium to conditions in which the NADP+/NADPH
ratio is far above the typical physiological values. Modeling results together
with experimental measurements of macro molecule levels (DNA, RNA, proteins and
fatty acids) indicate that the redox poise of the cell results in large-scale
changes in the activity of biosynthetic pathways. Phototrophic growth is less
coupled than expected to producing reductant, NAD(P)H, by reverse electron flow
from the quinone pool. Instead, it primarily functions for ATP production
(photophosphorylation), which drives reduction even when NADPH levels are
relatively low compared to NADP+. The model, in agreement with experimental
measurements of macromolecule ratios of cells growing on different carbon
substrates, indicate that the dynamics of nucleotide versus lipid and protein
production is likely a significant mechanism of balancing oxidation and
reduction in the cell.


  Liquid-liquid phase separation (LLPS) involving intrinsically disordered
protein regions (IDRs) is a major physical mechanism for biological
membraneless compartmentalization. The multifaceted electrostatic effects in
these biomolecular condensates are exemplified here by the different salt- and
ATP-dependent LLPSs of an IDR of messenger RNA-regulating protein Caprin1 and
its phosphorylated variant pY-Caprin1, exhibiting, e.g., reentrant behaviors in
some instances but not others. Physical modeling by analytical theory,
molecular dynamics, and polymer field-theoretic simulations indicates that, in
general, interchain salt bridges enhance LLPS of polyelectrolytes such as
Caprin1 and that the high valence of ATP-magnesium is a significant factor for
its colocalization with the condensed phases, as similar trends are observed
for several other IDRs. It is probable, therefore, that high concentrations of
ATP driven largely by electrostatics into some biomolecular condensates may
facilitate intracompartmental functional biochemical processes.


  Antimicrobial peptides (AMPs) play important roles in cancer, autoimmune
diseases, and aging. A critical aspect of AMP functionality is their targeted
interaction with pathogen membranes, which often possess altered lipid
compositions. Designing AMPs with enhanced therapeutic properties relies on a
nuanced understanding of these interactions, which are believed to trigger a
rearrangement of these peptides from random coil to alpha-helical
conformations, essential for their lytic action. Traditional supercomputing has
consistently encountered difficulties in accurately modeling these structural
changes, especially within membrane environments, thereby opening an
opportunity for more advanced approaches. This study extends an existing
quantum computing algorithm to address the complexities of antimicrobial
peptide interactions at interfaces. Our approach enables the prediction of the
optimal conformation of peptides located in the transition region between
hydrophilic and hydrophobic phases, akin to lipid membranes. The new method has
been applied to model the structure of three 10-amino-acid-long peptides, each
exhibiting hydrophobic, hydrophilic, or amphipathic properties in different
media and at interfaces between solvents of different polarity. Notably, our
approach does not demand a higher number of qubits compared to simulations in
homogeneous media, making it more feasible with current quantum computing
resources. Despite existing limitations in computational power and qubit
accessibility, our findings demonstrate the significant potential of quantum
computing in accurately characterizing complex biomolecular processes,
particularly the folding of AMPs at membrane models. This research paves the
way for future advances in quantum computing to enhance the accuracy and
applicability of biomolecular simulations.


  Coarse-grained models have emerged as valuable tools to simulate long DNA
molecules while maintaining computational efficiency. These models aim at
preserving interactions among coarse-grained variables in a manner that mirrors
the underlying atomistic description. We explore here a method for testing
coarse-grained vs. all-atom models using stiffness matrices in Fourier space
($q$-stiffnesses), which are particularly suited to probe DNA elasticity at
different length scales. We focus on a class of coarse-grained rigid base DNA
models known as cgDNA and its most recent version cgDNA+. Our analysis shows
that while cgDNA+ follows closely the $q$-stiffnesses of the all-atom model,
the original cgDNA shows some deviations for twist and bending variables which
are rather strong in the $q \to 0$ (long length scale) limit. The consequence
is that while both cgDNA and cgDNA+ give a suitable description of local
elastic behavior, the former misses some effects which manifest themselves at
longer length scales. In particular, cgDNA performs poorly on the twist
stiffness with a value much lower than expected for long DNA molecules.
Conversely, the all-atom and cgDNA+ twist is strongly length scale dependent:
DNA is torsionally soft at a few base pair distances, but becomes more rigid at
distances of a few dozens base pairs. Our analysis shows that the bending
persistence length in all-atom and cgDNA+ is somewhat overestimated.


  The field of antibody-based therapeutics has grown significantly in recent
years, with targeted antibodies emerging as a potentially effective approach to
personalized therapies. Such therapies could be particularly beneficial for
complex, highly individual diseases such as cancer. However, progress in this
field is often constrained by the extensive search space of amino acid
sequences that form the foundation of antibody design. In this study, we
introduce a novel reinforcement learning method specifically tailored to
address the unique challenges of this domain. We demonstrate that our method
can learn the design of high-affinity antibodies against multiple targets in
silico, utilizing either online interaction or offline datasets. To the best of
our knowledge, our approach is the first of its kind and outperforms existing
methods on all tested antigens in the Absolut! database.


  Accurate RNA secondary structure prediction is vital for understanding
cellular regulation and disease mechanisms. Deep learning (DL) methods have
surpassed traditional algorithms by predicting complex features like
pseudoknots and multi-interacting base pairs. However, traditional distance
measures can hardly deal with such tertiary interactions and the currently used
evaluation measures (F1 score, MCC) have limitations. We propose the
Weisfeiler-Lehman graph kernel (WL) as an alternative metric. Embracing
graph-based metrics like WL enables fair and accurate evaluation of RNA
structure prediction algorithms. Further, WL provides informative guidance, as
demonstrated in an RNA design experiment.


  In this work, we introduce AutoFragDiff, a fragment-based autoregressive
diffusion model for generating 3D molecular structures conditioned on target
protein structures. We employ geometric vector perceptrons to predict atom
types and spatial coordinates of new molecular fragments conditioned on
molecular scaffolds and protein pockets. Our approach improves the local
geometry of the resulting 3D molecules while maintaining high predicted binding
affinity to protein targets. The model can also perform scaffold extension from
user-provided starting molecular scaffold.


  The post-2020 global biodiversity framework needs ambitious, research-based
targets. Estimating the accelerated extinction risk due to climate change is
critical. The International Union for Conservation of Nature (IUCN) measures
the extinction risk of species. Automatic methods have been developed to
provide information on the IUCN status of under-assessed taxa. However, these
compensatory methods are based on current species characteristics, mainly
geographical, which precludes their use in future projections. Here, we
evaluate a novel method for classifying the IUCN status of species benefiting
from the generalisation power of species distribution models based on deep
learning. Our method matches state-of-the-art classification performance while
relying on flexible SDM-based features that capture species' environmental
preferences. Cross-validation yields average accuracies of 0.61 for status
classification and 0.78 for binary classification. Climate change will reshape
future species distributions. Under the species-environment equilibrium
hypothesis, SDM projections approximate plausible future outcomes. Two extremes
of species dispersal capacity are considered: unlimited or null. The projected
species distributions are translated into features feeding our IUCN
classification method. Finally, trends in threatened species are analysed over
time and i) by continent and as a function of average ii) latitude or iii)
altitude. The proportion of threatened species is increasing globally, with
critical rates in Africa, Asia and South America. Furthermore, the proportion
of threatened species is predicted to peak around the two Tropics, at the
Equator, in the lowlands and at altitudes of 800-1,500 m.


  Recently, there has been a growing interest in the relationships between
unrooted and rooted phylogenetic networks. In this context, a natural question
to ask is if an unrooted phylogenetic network U can be oriented as a rooted
phylogenetic network such that the latter satisfies certain structural
properties. In a recent preprint, Bulteau et al. claim that it is computational
hard to decide if U has a funneled (resp. funneled tree-child) orientation, for
when the internal vertices of U have degree at most 5. Unfortunately, the proof
of their funneled tree-child result appears to be incorrect. In this paper, we
present a corrected proof and show that hardness remains for other popular
classes of rooted phylogenetic networks such as funneled normal and funneled
reticulation-visible. Additionally, our results hold regardless of whether U is
rooted at an existing vertex or by subdividing an edge with the root.


  We study how Cas9, a central component of the CRISPR/Cas9 system, searches
for a target sequence on the DNA. We propose a model that includes as key
ingredients 3D diffusion, 1D sliding along the DNA, and the effect of short
binding sequences preceding the target (protospacer adjacent sequences --
PAMs). This latter aspect constitutes the main difference with traditional
facilitated diffusion of transcription factors. We solve our model, obtaining
an expression for the average search time of Cas9 for its target. We find that
experimentally measured kinetic parameters are close to the values yielding an
optimal search time. Our results rationalize the role of PAMs in guiding the
search process, and show that Cas9 searches for its targets in a nearly optimal
way.


  Generative models of macromolecules carry abundant and impactful implications
for industrial and biomedical efforts in protein engineering. However, existing
methods are currently limited to modeling protein structures or sequences,
independently or jointly, without regard to the interactions that commonly
occur between proteins and other macromolecules. In this work, we introduce
MMDiff, a generative model that jointly designs sequences and structures of
nucleic acid and protein complexes, independently or in complex, using joint
SE(3)-discrete diffusion noise. Such a model has important implications for
emerging areas of macromolecular design including structure-based transcription
factor design and design of noncoding RNA sequences. We demonstrate the utility
of MMDiff through a rigorous new design benchmark for macromolecular complex
generation that we introduce in this work. Our results demonstrate that MMDiff
is able to successfully generate micro-RNA and single-stranded DNA molecules
while being modestly capable of joint modeling DNA and RNA molecules in
interaction with multi-chain protein complexes. Source code:
https://github.com/Profluent-Internships/MMDiff.


  De novo drug design is a pivotal issue in pharmacology and a new area of
focus in AI for science research. A central challenge in this field is to
generate molecules with specific properties while also producing a wide range
of diverse candidates. Although advanced technologies such as transformer
models and reinforcement learning have been applied in drug design, their
potential has not been fully realized. Therefore, we propose MolRL-MGPT, a
reinforcement learning algorithm with multiple GPT agents for drug molecular
generation. To promote molecular diversity, we encourage the agents to
collaborate in searching for desirable molecules in diverse directions. Our
algorithm has shown promising results on the GuacaMol benchmark and exhibits
efficacy in designing inhibitors against SARS-CoV-2 protein targets. The codes
are available at: https://github.com/HXYfighter/MolRL-MGPT.


  We propose a new large-scale molecular model, named AdaMR, which stands for
Adjustable Molecular Representation for Unified Pre-training Strategy. Unlike
recent large-scale molecular models that use a single molecular encoding, AdaMR
employs a granularity-adjustable molecular encoder, learning molecular
representations at both the atomic and substructure levels. For the
pre-training process, we designed a task for molecular canonicalization, which
involves transforming ltiple generic molecular representations into canonical
representations. By adjusting the granularity of molecular encoding, the
trained model can improve the effects on multiple downstream tasks, such as
model attribute prediction and molecule generation. Substructure-level
molecular representation retains information of specific atom groups or
arrangements that determine chemical properties and have similar functions,
which is beneficial for tasks like property prediction. Meanwhile, atomic-level
representation, combined with generative molecular canonicalization
pre-training tasks, enhances the validity, novelty, and uniqueness in
generative tasks. These features of AdaMR demonstrate its strong performance in
numerous downstream tasks. We use different molecular properties prediction
tasks on six different datasets on MoleculeNet and two generative tasks on
ZINC250K dataset to evaluate our proposed molecular encoding and pre-training
methods, and obtain state-of-the-art (SOTA) results on five of these tasks.


  Numerous prognostic factors are currently assessed histopathologically in
biopsies of canine mast cell tumors to evaluate clinical behavior. In addition,
PCR analysis of the c-Kit exon 11 mutational status is often performed to
evaluate the potential success of a tyrosine kinase inhibitor therapy. This
project aimed at training deep learning models (DLMs) to identify the c-Kit-11
mutational status of MCTs solely based on morphology without additional
molecular analysis. HE slides of 195 mutated and 173 non-mutated tumors were
stained consecutively in two different laboratories and scanned with three
different slide scanners. This resulted in six different datasets
(stain-scanner variations) of whole slide images. DLMs were trained with single
and mixed datasets and their performances was assessed under scanner and
staining domain shifts. The DLMs correctly classified HE slides according to
their c-Kit 11 mutation status in, on average, 87% of cases for the best-suited
stain-scanner variant. A relevant performance drop could be observed when the
stain-scanner combination of the training and test dataset differed.
Multi-variant datasets improved the average accuracy but did not reach the
maximum accuracy of algorithms trained and tested on the same stain-scanner
variant. In summary, DLM-assisted morphological examination of MCTs can predict
c-Kit-exon 11 mutational status of MCTs with high accuracy. However, the
recognition performance is impeded by a change of scanner or staining protocol.
Larger data sets with higher numbers of scans originating from different
laboratories and scanners may lead to more robust DLMs to identify c-Kit
mutations in HE slides.


  While modern biotechnologies allow synthesizing new proteins and function
measurements at scale, efficiently exploring a protein sequence space and
engineering it remains a daunting task due to the vast sequence space of any
given protein. Protein engineering is typically conducted through an iterative
process of adding mutations to the wild-type or lead sequences, recombination
of mutations, and running new rounds of screening. To enhance the efficiency of
such a process, we propose a tree search-based bandit learning method, which
expands a tree starting from the initial sequence with the guidance of a bandit
machine learning model. Under simplified assumptions and a Gaussian Process
prior, we provide theoretical analysis and a Bayesian regret bound,
demonstrating that the combination of local search and bandit learning method
can efficiently discover a near-optimal design. The full algorithm is
compatible with a suite of randomized tree search heuristics, machine learning
models, pre-trained embeddings, and bandit techniques. We test various
instances of the algorithm across benchmark protein datasets using simulated
screens. Experiment results demonstrate that the algorithm is both
sample-efficient and able to find top designs using reasonably small mutation
counts.


  When hybridization or other forms of lateral gene transfer have occurred,
evolutionary relationships of species are better represented by phylogenetic
networks than by trees. While inference of such networks remains challenging,
several recently proposed methods are based on quartet concordance factors --
the probabilities that a tree relating a gene sampled from the species displays
the possible 4-taxon relationships. Building on earlier results, we investigate
what level-1 network features are identifiable from concordance factors under
the network multispecies coalescent model. We obtain results on both
topological features of the network, and numerical parameters, uncovering a
number of failures of identifiability related to 3-cycles in the network.


  Translation of proteins is a fundamental part of gene expression that is
mediated by ribosomes. As ribosomes significantly contribute to both cellular
mass and energy consumption, achieving efficient management of the ribosome
population is also crucial to metabolism and growth. Inspired by biological
evidence for nutrient-dependent mechanisms that control both ribosome active
degradation and genesis, we introduce a dynamical model of protein production,
that includes the dynamics of resources and control over the ribosome
population. Under the hypothesis that active degradation and biogenesis are
optimal for maximizing and maintaining protein production, we aim to
qualitatively reproduce empirical observations of the ribosome population
dynamics. Upon formulating the associated optimization problem, we first
analytically study the stability and global behaviour of solutions under
constant resource input, and characterize the extent of oscillations and
convergence rate to a global equilibrium. We further use these results to
simplify and solve the problem under a quasi-static approximation. Using
biophysical parameter values, we find that optimal control solutions lead to
both control mechanisms and the ribosome population switching between periods
of feeding and fasting, suggesting that the intense regulation of ribosome
population observed in experiments allows to maximize and maintain protein
production. Finally, we find some range for the control values over which such
a regime can be observed, depending on the intensity of fasting.


  This paper studies a two microbial species model in competition for a single
resource in the chemostat including general interspecific density-dependent
growth rates with distinct removal rates for each species. We give the
necessary and sufficient conditions of existence, uniqueness, and local
stability of all steady states. We show that a positive steady state, if it
exists, then it is unique and unstable. In this case, the system exhibits a
bi-stability where the behavior of the process depends on the initial
condition. Our mathematical analysis proves that at most one species can
survive which confirms the competitive exclusion principle. We conclude that
adding only interspecific competition in the classical chemostat model is not
sufficient to show the coexistence of two species even considering mortality in
the dynamics of two species. Otherwise, we focus on the study, theoretically
and numerically, of the operating diagram which depicts the existence and the
stability of each steady state according to the two operating parameters of the
process which are the dilution rate and the input concentration of the
substrate. Using our mathematical analysis, we construct analytically the
operating diagram by plotting the curves that separate their various regions.
Our numerical method using MATCONT software validates these theoretical results
but it reveals new bifurcations that occur by varying two parameters as
Bogdanov-Takens and Zero-Hopf bifurcations. The bifurcation analysis shows that
all steady states can appear or disappear only through transcritical
bifurcations.


  The emergence and spread of deadly pandemics has repeatedly occurred
throughout history, causing widespread infections and loss of life. The rapid
spread of pandemics have made governments across the world adopt a range of
actions, including non-pharmaceutical measures to contain its impact. However,
the dynamic nature of pandemics makes selecting intervention strategies
challenging. Hence, the development of suitable monitoring and forecasting
tools for tracking infected cases is crucial for designing and implementing
effective measures. Motivated by this, we present a hybrid pandemic infection
forecasting methodology that integrates compartmental model and learning-based
approaches. In particular, we develop a compartmental model that includes
time-varying infection rates, which are the key parameters that determine the
pandemic's evolution. To identify the time-dependent infection rates, we
establish a hybrid methodology that combines the developed compartmental model
and tools from optimization and neural networks. Specifically, the proposed
methodology estimates the infection rates by fitting the model to available
data, regarding the COVID-19 pandemic in Cyprus, and then predicting their
future values through either a) extrapolation, or b) feeding them to neural
networks. The developed approach exhibits strong accuracy in predicting
infections seven days in advance, achieving low average percentage errors both
using the extrapolation (9.90%) and neural network (5.04%) approaches.


  In recent years, deep learning has demonstrated promising results in de novo
drug design. However, the proposed techniques still lack an efficient
exploration of the large chemical space. Most of these methods explore a small
fragment of the chemical space of known drugs, if the desired molecules were
not found, the process ends. In this work, we introduce a curiosity-driven
method to force the model to navigate many parts of the chemical space,
therefore, achieving higher desirability and diversity as well. At first, we
train a recurrent neural network-based general molecular generator (G), then we
fine-tune G to maximize curiosity and desirability. We define curiosity as the
Tanimoto similarity between two generated molecules, a first molecule generated
by G, and a second one generated by a copy of G (Gcopy). We only backpropagate
the loss through G while keeping Gcopy unchanged. We benchmarked our approach
against two desirable chemical properties related to drug-likeness and showed
that the discovered chemical space can be significantly expanded, thus,
discovering a higher number of desirable molecules with more diversity and
potentially easier to synthesize. All Code and data used in this paper are
available at https://github.com/amine179/Curiosity-RL-for-Drug-Design.


  Recent developments in single-cell omics technologies have enabled the
quantification of molecular profiles in individual cells at an unparalleled
resolution. Deep learning, a rapidly evolving sub-field of machine learning,
has instilled a significant interest in single-cell omics research due to its
remarkable success in analysing heterogeneous high-dimensional single-cell
omics data. Nevertheless, the inherent multi-layer nonlinear architecture of
deep learning models often makes them `black boxes' as the reasoning behind
predictions is often unknown and not transparent to the user. This has
stimulated an increasing body of research for addressing the lack of
interpretability in deep learning models, especially in single-cell omics data
analyses, where the identification and understanding of molecular regulators
are crucial for interpreting model predictions and directing downstream
experimental validations. In this work, we introduce the basics of single-cell
omics technologies and the concept of interpretable deep learning. This is
followed by a review of the recent interpretable deep learning models applied
to various single-cell omics research. Lastly, we highlight the current
limitations and discuss potential future directions. We anticipate this review
to bring together the single-cell and machine learning research communities to
foster future development and application of interpretable deep learning in
single-cell omics research.


  Edge-based percolation methods can be used to analyze disease transmission on
complex social networks. This allows us to include complex social heterogeneity
in our models while maintaining tractability. Here we review the seminal works
on this field by Newman et al (2001); Newman (2002, 2003), and Miller et al
(2012). We present a systematic discussion of the theoretical background behind
these models, including an extensive derivation of the major results. We also
connect these results relate back to the classical literature in random graph
theory Molloy and Reed (1995, 1998). Finally, we also present an accompanying R
package that takes epidemic and network parameters as input and generates
estimates of the epidemic trajectory and final size. This manuscript and the R
package was developed to help researchers easily understand and use network
models to investigate the interaction between different community structures
and disease transmission.


  This work discloses an epidemiological mathematical model to predict an
empirical treatment for dogs infected by Pseudomonas aeruginosa. This dangerous
pathogen is one of the leading causes of multi-resistant infections and can be
transmitted from dogs to humans. Numerical simulations and appropriated codes
were developed using Matlab software to gather information concerning long-time
dynamics of the susceptible, infected and recovered individuals. All data
compiled from the mathematical model was used to provide an appropriated
antibiotic sensitivity panel for this specific infection. In this study,
several variables have been included in this model to predict which treatment
should be prescribed in emergency cases, when there is no time to perform an
antibiogram or the cost of it could not be assumed. In particular, we highlight
the use of this model aiming to become part of the convenient toolbox of Public
Health research and decision-making in the design of the mitigation strategy of
bacterial pathogens.


  The energy landscape theory finds its both extensive and intensive
application in studying stochastic dynamics of physical and biological systems.
Although the weighted summation of the Gaussian approximation (WSGA) approach
has been proposed for quantifying the energy landscape in multistable systems
by solving the diffusion equation approximately from moment equations, we are
still lacking an accurate approach for quantifying the energy landscape of the
periodic oscillatory systems. To address this challenge, we propose an
approach, called the diffusion decomposition of the Gaussian approximation
(DDGA). Using typical oscillatory systems as examples, we demonstrate the
efficacy of the proposed DDGA in quantifying the energy landscape of
oscillatory systems and corresponding stochastic dynamics, in comparison with
existing approaches. By further applying the DDGA to a high-dimensional cell
cycle network, we are able to uncover more intricate biological mechanisms in
cell cycle, which cannot be discerned using previously developed approaches.


  Chemical reaction networks (CRN) comprise an important class of models to
understand biological functions such as cellular information processing, the
robustness and control of metabolic pathways, circadian rhythms, and many more.
However, any CRN describing a certain function does not act in isolation but is
a part of a much larger network and as such is constantly subject to external
changes. In [Shinar, Alon, and Feinberg. "Sensitivity and robustness in
chemical reaction networks." SIAM J App Math (2009): 977-998.], the responses
of CRN to changes in the linear conserved quantities, called sensitivities,
were studied in and the question of how to construct absolute, i.e.,
basis-independent, sensitivities was raised. In this article, by applying
information geometric methods, such a construction is provided. The idea is to
track how concentration changes in a particular chemical propagate to changes
of all the other chemicals within a steady state. This is encoded in the matrix
of absolute sensitivites. A linear algebraic characterization of the matrix of
absolute sensitivities for quasi-thermostatic CRN is derived via a Cramer-Rao
bound for CRN, which is based on the the analogy between quasi-thermostatic
steady states and the exponential family of probability distributions.


  Cell polarization is a critical process that separates molecules into two
distinct regions in prokaryotic and eukaryotic cells, guiding biological
processes such as cell division and cell differentiation. Although several
underlying antagonistic reaction-diffusion networks capable of setting up cell
polarization have been identified experimentally and theoretically, our
understanding of how to manipulate pattern stability and asymmetry remains
incomplete, especially when only a subset of network components are known. Here
we present numerical results to show that the polarized pattern of an
antagonistic 2-node network collapses into a homogeneous state when subjected
to single-sided self-regulation, single-sided additional regulation, or unequal
system parameters. However, polarity can be restored through a combination of
two modifications that have opposing effects. Additionally, spatially
inhomogeneous parameters favoring respective domains stabilize their interface
at designated locations. To connect our findings to cell polarity studies of
the nematode Caenorhabditis elegans zygote, we reconstituted a 5-node network
where a 4-node circuit with full mutual inhibitions between anterior and
posterior is modified by a mutual activation in the anterior and an additional
mutual inhibition between the anterior and the posterior. Once again, a generic
set of kinetic parameters moves the interface towards either the anterior or
posterior end, yet a polarized pattern can be stabilized through spatial tuning
of one or more parameters coupled to intracellular or extracellular cues. A
user-friendly software, PolarSim, is introduced to facilitate the exploration
of networks with alternative node numbers, parameter values, and regulatory
pathways.


  This manuscript delves into the intersection of genomics and phenotypic
prediction, focusing on the statistical innovation required to navigate the
complexities introduced by noisy covariates and confounders. The primary
emphasis is on the development of advanced robust statistical models tailored
for genomic prediction from single nucleotide polymorphism (SNP) data collected
from genome-wide association studies (GWAS) in plant and animal breeding and
multi-field trials. The manuscript explores the limitations of traditional
marker-assisted recurrent selection, highlighting the significance of
incorporating all estimated effects of marker loci into the statistical
framework and aiming to reduce the high dimensionality of GWAS data while
preserving critical information. This paper introduces a new robust statistical
framework for genomic prediction, employing one-stage and two-stage linear
mixed model analyses along with utilizing the popular robust minimum density
power divergence estimator (MDPDE) to estimate genetic effects on phenotypic
traits. The study illustrates the superior performance of the proposed
MDPDE-based genomic prediction and associated heritability estimation
procedures over existing competitors through extensive empirical experiments on
artificial datasets and application to a real-life maize breeding dataset. The
results showcase the robustness and accuracy of the proposed MDPDE-based
approaches, especially in the presence of data contamination, emphasizing their
potential applications in improving breeding programs and advancing genomic
prediction of phenotyping traits.


  One of the main goals of developmental biology is to reveal the gene
regulatory networks (GRNs) underlying the robust differentiation of multipotent
progenitors into precisely specified cell types. Most existing methods to infer
GRNs from experimental data have limited predictive power as the inferred GRNs
merely reflect gene expression similarity or correlation. Here, we demonstrate,
how physics-informed neural networks (PINNs) can be used to infer the
parameters of predictive, dynamical GRNs that provide mechanistic understanding
of biological processes. Specifically we study GRNs that exhibit bifurcation
behavior and can therefore model cell differentiation. We show that PINNs
outperform regular feed-forward neural networks on the parameter inference task
and analyze two relevant experimental scenarios: 1. a system with cell
communication for which gene expression trajectories are available and 2.
snapshot measurements of a cell population in which cell communication is
absent. Our analysis will inform the design of future experiments to be
analyzed with PINNs and provides a starting point to explore this powerful
class of neural network models further.


  AI for drug discovery has been a research hotspot in recent years, and
SMILES-based language models has been increasingly applied in drug molecular
design. However, no work has explored whether and how language models
understand the chemical spatial structure from 1D sequences. In this work, we
pre-train a transformer model on chemical language and fine-tune it toward drug
design objectives, and investigate the correspondence between high-frequency
SMILES substrings and molecular fragments. The results indicate that language
models can understand chemical structures from the perspective of molecular
fragments, and the structural knowledge learned through fine-tuning is
reflected in the high-frequency SMILES substrings generated by the model.


  Much progress has been made in elucidating the inner workings of
voltage-gated ion channels, but less understood is the influence of lipid rafts
on gating kinetics. Here we propose that state-dependent channel affinity for
different lipid species provides a unified explanation for the experimentally
observed behaviors of clustering, cooperativity, and hysteresis. We develop
models of diffusing lipids and channels engaged in Ising-like interactions to
investigate the collective behaviors driven by raft formation in critical
membranes close to the demixing transition. The model channels demonstrate
lipid-mediated long-range interactions, activation curve steepening, and
long-term memory in ionic currents. These behaviors likely play a role in
channel-mediated cellular signaling and suggest a universal mechanism for
self-organization of biomolecular assemblies.


  Sponges have long been recognized as a rich source of bioactive natural
products. Various studies suggest that many of these compounds are produced by
symbiotic bacteria. However, substance supplies and functional insights about
the producers remain limited because cultivation remains unsuccessful. To
identify alternative, sustainable sources of sponge-derived polyketides, we
computationally analyzed 5289 characterized and orphan trans-acyltransferase
polyketide synthases, enzymes with widespread roles in polyketide biosynthesis
by bacterial symbionts. The analytical workflow predicted marine animal-derived
Acidobacteria of the family Acanthopleuribacteraceae with large sets of
biosynthetic gene clusters to be enriched in sponge-type chemistry. Targeted
compound isolation from a chiton-associated strain yielded new congeners of the
phorboxazoles and calyculins, potent and scarce cytotoxins exclusively known
from sponges. These first natural products of Acidobacteria and new coral
metagenomic data on a third family member suggest animal-associated
Acanthopleuribacteraceae as a rich source of sponge-type as well as novel
metabolites


  Cells contain multiple condensates which spontaneously form due to the
heterotypic interactions between their components. Although the proteins and
disordered region sequences that are responsible for condensate formation have
been extensively studied, the rule of interactions between the components that
allow demixing, i.e., the coexistence of multiple condensates, is yet to be
elucidated. Here we construct an effective theory of the interaction between
heteropolymers by fitting it to the molecular dynamics simulation results
obtained for more than 200 sequences sampled from the disordered regions of
human proteins. We find that the sum of amino acid pair interactions across two
heteropolymers predicts the Boyle temperature qualitatively well, which can be
quantitatively improved by the dimer pair approximation, where we incorporate
the effect of neighboring amino acids in the sequences. The improved theory,
combined with the finding of a metric that captures the effective interaction
strength between distinct sequences, allowed the selection of up to three
disordered region sequences that demix with each other in multicomponent
simulations, as well as the generation of artificial sequences that demix with
a given sequence. The theory points to a generic sequence design strategy to
demix or hypermix thanks to the low dimensional nature of the space of the
interactions that we identify. As a consequence of the geometric arguments in
the space of interactions, we find that the number of distinct sequences that
can demix with each other is strongly constrained, irrespective of the choice
of the coarse-grained model. Altogether, we construct a theoretical basis for
methods to estimate the effective interaction between heteropolymers, which can
be utilized in predicting phase separation properties as well as rules of
assignment in the localization and functions of disordered proteins.


  High-throughput omics profiling advancements have greatly enhanced cancer
patient stratification. However, incomplete data in multi-omics integration
presents a significant challenge, as traditional methods like sample exclusion
or imputation often compromise biological diversity and dependencies.
Furthermore, the critical task of accurately classifying new patients with
partial omics data into existing subtypes is commonly overlooked. To address
these issues, we introduce IntegrAO (Integrate Any Omics), an unsupervised
framework for integrating incomplete multi-omics data and classifying new
samples. IntegrAO first combines partially overlapping patient graphs from
diverse omics sources and utilizes graph neural networks to produce unified
patient embeddings. Our systematic evaluation across five cancer cohorts
involving six omics modalities demonstrates IntegrAO's robustness to missing
data and its accuracy in classifying new samples with partial profiles. An
acute myeloid leukemia case study further validates its capability to uncover
biological and clinical heterogeneity in incomplete datasets. IntegrAO's
ability to handle heterogeneous and incomplete data makes it an essential tool
for precision oncology, offering a holistic approach to patient
characterization.


  N6-methyladenosine (m6A) is the most abundant mRNA modification within
mammalian cells, holding pivotal significance in the regulation of mRNA
stability, translation, and splicing. Furthermore, it plays a critical role in
the regulation of RNA degradation by primarily recruiting the YTHDF2 reader
protein. However, the selective regulation of mRNA decay of the m6A-methylated
mRNA through YTHDF2 binding is poorly understood. To improve our understanding,
we developed m6A-BERT-Deg, a BERT model adapted for predicting YTHDF2-mediated
degradation of m6A-methylated mRNAs. We meticulously assembled a high-quality
training dataset by integrating multiple data sources for the HeLa cell line.
To overcome the limitation of small training samples, we employed a
pre-training-fine-tuning strategy by first performing a self-supervised
pre-training of the model on 427,760 unlabeled m6A site sequences. The test
results demonstrated the importance of this pre-training strategy in enabling
m6A-BERT-Deg to outperform other benchmark models. We further conducted a
comprehensive model interpretation and revealed a surprising finding that the
presence of co-factors in proximity to m6A sites may disrupt YTHDF2-mediated
mRNA degradation, subsequently enhancing mRNA stability. We also extended our
analyses to the HEK293 cell line, shedding light on the context-dependent
YTHDF2-mediated mRNA degradation.


  Background: Men and women with a migration background comprise an increasing
proportion of incident HIV cases across Western Europe. Several studies
indicate a substantial proportion acquire HIV post-migration.
  Methods: We used partial HIV consensus sequences with linked demographic and
clinical data from the opt-out ATHENA cohort of people with HIV in the
Netherlands to quantify population-level sources of transmission to Dutch-born
and foreign-born Amsterdam men who have sex with men (MSM) between 2010-2021.
We identified phylogenetically and epidemiologically possible transmission
pairs in local transmission chains and interpreted these in the context of
estimated infection dates, quantifying transmission dynamics between
sub-populations by world region of birth.
  Results: We estimate the majority of Amsterdam MSM who acquired their
infection locally had a Dutch-born Amsterdam MSM source (56% [53-58%]).
Dutch-born MSM were the predominant source population of infections among
almost all foreign-born Amsterdam MSM sub-populations. Stratifying by two-year
intervals indicated shifts in transmission dynamics, with a majority of
infections originating from foreign-born MSM since 2018, although uncertainty
ranges remained wide.
  Conclusions: In the context of declining HIV incidence among Amsterdam MSM,
our data suggest whilst native-born MSM have predominantly driven transmissions
in 2010-2021, the contribution from foreign-born MSM living in Amsterdam is
increasing.


  Background: The secure management of sensitive clinical data, particularly
human genomics data, has become a critical requirement in modern biomedical
research. Although the necessary software and algorithms are readily available,
their use by non-IT experts poses significant challenges.
  Methods: We developed dabih, an open-source web application specifically
designed to facilitate user-friendly encrypted data management. dabih enables
web-based uploading, storing, sharing, and downloading of sensitive data in any
format. Its approach to data security involves a two-stage envelope encryption
process. We combine symmetric-key encryption for data and public-key encryption
as key encapsulation mechanism. The private key necessary for decrypting the
data remains exclusively on the owner's device. Thus, accessing data is
impossible without explicit permission from the keyholder.
  Results: dabih is available open-source on GitHub
https://github.com/spang-lab/dabih, as ready to use containers on docker hub
and includes a command line interface and a graphical bulk upload tool as
pre-built binaries. Documentation is available as part of the web application.
  Conclusions: dabih enables everyone to use strong cryptography for their
data, while being just as simple to use as other, non-encrypted, data storage
solutions. All the cryptography occurs seamlessly in the background as users
interact with a secure web portal, simply by dragging and dropping files.


  This study introduces the Supervised Magnitude-Altitude Scoring (SMAS)
methodology, a machine learning-based approach, for analyzing gene expression
data obtained from nonhuman primates (NHPs) infected with Ebola virus (EBOV).
We utilize a comprehensive dataset of NanoString gene expression profiles from
Ebola-infected NHPs, deploying the SMAS system for nuanced host-pathogen
interaction analysis. SMAS effectively combines gene selection based on
statistical significance and expression changes, employing linear classifiers
such as logistic regression to accurately differentiate between RT-qPCR
positive and negative NHP samples. A key finding of our research is the
identification of IFI6 and IFI27 as critical biomarkers, demonstrating
exceptional predictive performance with 100% accuracy and Area Under the Curve
(AUC) metrics in classifying various stages of Ebola infection. Alongside IFI6
and IFI27, genes, including MX1, OAS1, and ISG15, were significantly
upregulated, highlighting their essential roles in the immune response to EBOV.
Our results underscore the efficacy of the SMAS method in revealing complex
genetic interactions and response mechanisms during EBOV infection. This
research provides valuable insights into EBOV pathogenesis and aids in
developing more precise diagnostic tools and therapeutic strategies to address
EBOV infection in particular and viral infection in general.


  Measures of tree balance play an important role in different research areas
such as mathematical phylogenetics or theoretical computer science. The balance
of a tree is usually quantified in a single number, called a balance or
imbalance index, and several such indices exist in the literature. Here, we
focus on the stairs2 balance index for rooted binary trees, which was first
introduced in the context of viral phylogenetics but has not been fully
analyzed from a mathematical viewpoint yet. While it is known that the
caterpillar tree uniquely minimizes the stairs2 index for all leaf numbers and
the fully balanced tree uniquely maximizes the stairs2 index for leaf numbers
that are powers of two, understanding the maximum value and maximal trees for
arbitrary leaf numbers is an open problem in the literature. In this note, we
fill this gap by showing that for all leaf numbers, there is a unique rooted
binary tree maximizing the stairs2 index. Additionally, we obtain recursive and
closed expressions for the maximum value of the stairs2 index of a rooted
binary tree with $n$ leaves.


  The study of rigid protein-protein docking plays an essential role in a
variety of tasks such as drug design and protein engineering. Recently, several
learning-based methods have been proposed for the task, exhibiting much faster
docking speed than those computational methods. In this paper, we propose a
novel learning-based method called ElliDock, which predicts an elliptic
paraboloid to represent the protein-protein docking interface. To be specific,
our model estimates elliptic paraboloid interfaces for the two input proteins
respectively, and obtains the roto-translation transformation for docking by
making two interfaces coincide. By its design, ElliDock is independently
equivariant with respect to arbitrary rotations/translations of the proteins,
which is an indispensable property to ensure the generalization of the docking
process. Experimental evaluations show that ElliDock achieves the fastest
inference time among all compared methods and is strongly competitive with
current state-of-the-art learning-based models such as DiffDock-PP and Multimer
particularly for antibody-antigen docking.