  Wildlife camera traps and crowd-sourced image material provide novel
possibilities to monitor endangered animal species. However, massive image
volumes that these methods produce are overwhelming for researchers to go
through manually which calls for automatic systems to perform the analysis. The
analysis task that has gained the most attention is the re-identification of
individuals, as it allows, for example, to study animal migration or to
estimate the population size. The Saimaa ringed seal (Pusa hispida saimensis)
is an endangered subspecies only found in the Lake Saimaa, Finland, and is one
of the few existing freshwater seal species. Ringed seals have permanent pelage
patterns that are unique to each individual which can be used for the
identification of individuals. Large variation in poses further exacerbated by
the deformable nature of seals together with varying appearance and low
contrast between the ring pattern and the rest of the pelage makes the Saimaa
ringed seal re-identification task very challenging, providing a good benchmark
to evaluate state-of-the-art re-identification methods. Therefore, we make our
Saimaa ringed seal image (SealID) dataset (N=57) publicly available for
research purposes. In this paper, the dataset is described, the evaluation
protocol for re-identification methods is proposed, and the results for two
baseline methods HotSpotter and NORPPA are provided. The SealID dataset has
been made publicly available.


  The human proteome contains a vast network of interacting kinases and
substrates. Even though some kinases have proven to be immensely useful as
therapeutic targets, a majority are still understudied. In this work, we
present a novel knowledge graph representation learning approach to predict
novel interaction partners for understudied kinases. Our approach uses a
phosphoproteomic knowledge graph constructed by integrating data from iPTMnet,
Protein Ontology, Gene Ontology and BioKG. The representation of kinases and
substrates in this knowledge graph are learned by performing directed random
walks on triples coupled with a modified SkipGram or CBOW model. These
representations are then used as an input to a supervised classification model
to predict novel interactions for understudied kinases. We also present a
post-predictive analysis of the predicted interactions and an ablation study of
the phosphoproteomic knowledge graph to gain an insight into the biology of the
understudied kinases.


  Amidst the COVID-19 pandemic, travel restrictions have emerged as crucial
interventions for mitigating the spread of the virus. In this study, we enhance
the predictive capabilities of our model, Sequence-to-Sequence Epidemic
Attention Network (S2SEA-Net), by incorporating an attention module, allowing
us to assess the impact of distinct classes of travel distances on epidemic
dynamics. Furthermore, our model provides forecasts for new confirmed cases and
deaths. To achieve this, we leverage daily data on population movement across
various travel distance categories, coupled with county-level epidemic data in
the United States. Our findings illuminate a compelling relationship between
the volume of travelers at different distance ranges and the trajectories of
COVID-19. Notably, a discernible spatial pattern emerges with respect to these
travel distance categories on a national scale. We unveil the geographical
variations in the influence of population movement at different travel
distances on the dynamics of epidemic spread. This will contribute to the
formulation of strategies for future epidemic prevention and public health
policies.


  Analyzing the SARS-CoV-2 pandemic outbreak based on actual data while
reflecting the characteristics of the real city provides beneficial information
for taking reasonable infection control measures in the future. We demonstrate
agent-based modeling for Tokyo based on GPS information and official national
statistics and perform a spatiotemporal analysis of the infection situation in
Tokyo. As a result of the simulation during the first wave of SARS-CoV-2 in
Tokyo using real GPS data, the infection occurred in the service industry, such
as restaurants, in the city center, and then the infected people brought back
the virus to the residential area; the infection spread in each area in Tokyo.
This phenomenon clarifies that the spread of infection can be curbed by
suppressing going out or strengthening infection prevention measures in service
facilities. It was shown that pandemic measures in Tokyo could be achieved not
only by strong control, such as the lockdown of cities, but also by thorough
infection prevention measures in service facilities, which explains the curb
phenomena in real Tokyo.


  Recent advances in applying Graph Neural Networks (GNNs) to molecular science
have showcased the power of learning three-dimensional (3D) structure
representations with GNNs. However, most existing GNNs suffer from the
limitations of insufficient modeling of diverse interactions, computational
expensive operations, and ignorance of vectorial values. Here, we tackle these
limitations by proposing a novel GNN model, Physics-aware Multiplex Graph
Neural Network (PaxNet), to efficiently and accurately learn the
representations of 3D molecules for both small organic compounds and
macromolecule complexes. PaxNet separates the modeling of local and non-local
interactions inspired by molecular mechanics, and reduces the expensive
angle-related computations. Besides scalar properties, PaxNet can also predict
vectorial properties by learning an associated vector for each atom. To
evaluate the performance of PaxNet, we compare it with state-of-the-art
baselines in two tasks. On small molecule dataset for predicting quantum
chemical properties, PaxNet reduces the prediction error by 15% and uses 73%
less memory than the best baseline. On macromolecule dataset for predicting
protein-ligand binding affinities, PaxNet outperforms the best baseline while
reducing the memory consumption by 33% and the inference time by 85%. Thus,
PaxNet provides a universal, robust and accurate method for large-scale machine
learning of molecules. Our code is available at
https://github.com/zetayue/Physics-aware-Multiplex-GNN.


  Breeding for new crop characteristics and adjusting management practices are
critical avenues to mitigate yield loss and maintain yield stability under a
changing climate. However, identifying high-performing plant traits and
management options for different growing regions through traditional breeding
practices and agronomic field trials is often time and resource-intensive.
Mechanistic crop simulation models can serve as powerful tools to help
synthesize cropping information, set breeding targets, and develop adaptation
strategies to sustain food production. In this study, we develop a modeling
framework for a mechanistic crop model (MAIZSIM) to run many simulations within
a trait x environment x management landscape and demonstrate how such a
modeling framework could be used to identify ideal trait-management
combinations that maximize yield and yield stability for different agro-climate
regions in the US.


  SARS-COV-19 is the most prominent issue which many countries face today. The
frequent changes in infections, recovered and deaths represents the dynamic
nature of this pandemic. It is very crucial to predict the spreading rate of
this virus for accurate decision making against fighting with the situation of
getting infected through the virus, tracking and controlling the virus
transmission in the community. We develop a prediction model using statistical
time series models such as SARIMA and FBProphet to monitor the daily active,
recovered and death cases of COVID-19 accurately. Then with the help of various
details across each individual patient (like height, weight, gender etc.), we
designed a set of rules using Semantic Web Rule Language and some mathematical
models for dealing with COVID19 infected cases on an individual basis. After
combining all the models, a COVID-19 Ontology is developed and performs various
queries using SPARQL query on designed Ontology which accumulate the risk
factors, provide appropriate diagnosis, precautions and preventive suggestions
for COVID Patients. After comparing the performance of SARIMA and FBProphet, it
is observed that the SARIMA model performs better in forecasting of COVID
cases. On individual basis COVID case prediction, approx. 497 individual
samples have been tested and classified into five different levels of COVID
classes such as Having COVID, No COVID, High Risk COVID case, Medium to High
Risk case, and Control needed case.


  Network medicine has been pursued for Covid-19 drug repurposing. One such
approach adopts structural controllability, a theory for controlling a network
(the cell). Motivated to protect the cell from viral infections, we extended
this theory to total controllability and introduced a new concept of control
hubs. Perturbation to any control hub renders the cell uncontrollable by
exogenous stimuli, e.g., viral infections, so control hubs are ideal drug
targets. We developed an efficient algorithm for finding all control hubs and
applied it to the largest homogenous human protein-protein interaction network.
Our new method outperforms several popular gene-selection methods, including
that based on structural controllability. The final 65 druggable control hubs
are enriched with functions of cell proliferation, regulation of apoptosis, and
responses to cellular stress and nutrient levels, revealing critical pathways
induced by SARS-CoV-2. These druggable control hubs led to drugs in 4 major
categories: antiviral and anti-inflammatory agents, drugs on central nerve
systems, and dietary supplements and hormones that boost immunity. Their
functions also provided deep insights into the therapeutic mechanisms of the
drugs for Covid-19 therapy, making the new approach an explainable drug
repurposing method. A remarkable example is Fostamatinib that has been shown to
lower mortality, shorten the length of ICU stay, and reduce disease severity of
hospitalized Covid-19 patients. The drug targets 10 control hubs, 9 of which
are kinases that play key roles in cell differentiation and programmed death.
One such kinase is RIPK1 that directly interacts with viral protein nsp12, the
RdRp of the virus. The study produced many control hubs that were not targets
of existing drugs but were enriched with proteins on membranes and the
NF-$\kappa$B pathway, so are excellent candidate targets for new drugs.


  Mismatch repair is a critical step in DNA replication that occurs after base
selection and proofreading, significantly increasing fidelity. However, the
mechanism of mismatch recognition has not been established for any repair
enzyme. Speculations in this area mainly focus on exploiting thermodynamic
equilibrium and free energy. Nevertheless, non-equilibrium processes may play a
more significant role in enhancing mismatch recognition accuracy by utilizing
adenosine triphosphate (ATP). This study aimed to investigate this possibility.
Considering our limited knowledge of actual mismatch repair enzymes, we
proposed a hypothetical enzyme that operates as a quantum system with three
discrete energy levels. When the enzyme is raised to its highest energy level,
a quantum transition occurs, leading to one of two low-energy levels
representing potential recognition outcomes: a correct match or a mismatch. The
probabilities of the two outcomes are exponentially different, determined by
the energy gap between the two low energy levels. By flipping the energy gap,
discrimination between mismatches and correct matches can be achieved. Within a
framework that combines quantum mechanics with thermodynamics, we established a
relationship between energy cost and the recognition error.


  Undulation of infection levels, usually called waves, are not well
understood. In this paper we propose a mathematical model that exhibits
undulation and decay towards a stable state. The model is a re-interpretation
of the original SIR-model obtained by postulating different constitutive
relations whereby classical logistic growth with recovery is obtained. The
recovery relation is based on the premise that infectiousness only lasts for
some time. This leads to a differential-difference (delay) equation which
intrinsically exhibits periodicity in its solutions but not necessarily decay
to asymptotically stable equilibrium. Limit cycles can indeed occur. An
appropriate linearization of the governing equation provides a firm basis for
heuristic reasoning as well as confidence in numerical calculations.


  In recent years, aligning a sequence to a pangenome has become a central
problem in genomics and pangenomics. A fast and accurate solution to this
problem can serve as a toolkit to many crucial tasks such as read-correction,
Multiple Sequences Alignment (MSA), genome assemblies, variant calling, just to
name a few. In this paper we propose a new, fast and exact method to align a
string to a D-string, the latter possibly representing an MSA, a pan-genome or
a partial assembly. An implementation of our tool dsa is publicly available at
https://github.com/urbanslug/dsa


  In this work we study Petri nets associated to ordinary differential
equations coming from epidemiological models. We propose a geometric procedure
to obtain the basic reproduction number. This number is defined as the average
number of secondary individuals during their entire infectious lifetime.


  Designing accurate deep learning models for molecular property prediction
plays an increasingly essential role in drug and material discovery. Recently,
due to the scarcity of labeled molecules, self-supervised learning methods for
learning generalizable and transferable representations of molecular graphs
have attracted lots of attention. In this paper, we argue that there exist two
major issues hindering current self-supervised learning methods from obtaining
desired performance on molecular property prediction, that is, the ill-defined
pre-training tasks and the limited model capacity. To this end, we introduce
Knowledge-guided Pre-training of Graph Transformer (KPGT), a novel
self-supervised learning framework for molecular graph representation learning,
to alleviate the aforementioned issues and improve the performance on the
downstream molecular property prediction tasks. More specifically, we first
introduce a high-capacity model, named Line Graph Transformer (LiGhT), which
emphasizes the importance of chemical bonds and is mainly designed to model the
structural information of molecular graphs. Then, a knowledge-guided
pre-training strategy is proposed to exploit the additional knowledge of
molecules to guide the model to capture the abundant structural and semantic
information from large-scale unlabeled molecular graphs. Extensive
computational tests demonstrated that KPGT can offer superior performance over
current state-of-the-art methods on several molecular property prediction
tasks.


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


  Single-cell RNA sequencing (scRNAseq) is rapidly advancing our understanding
of cellular composition within complex tissues and organisms. A major
limitation in most scRNAseq analysis pipelines is the reliance on manual
annotations to determine cell identities, which are time consuming, subjective,
and require expertise. Given the surge in cell sequencing, supervised
methods-especially deep learning models-have been developed for automatic cell
type identification (ACTI), which achieve high accuracy and scalability.
However, all existing deep learning frameworks for ACTI lack interpretability
and are used as "black-box" models. We present N-ACT (Neural-Attention for Cell
Type identification): the first-of-its-kind interpretable deep neural network
for ACTI utilizing neural-attention to detect salient genes for use in
cell-type identification. We compare N-ACT to conventional annotation methods
on two previously manually annotated data sets, demonstrating that N-ACT
accurately identifies marker genes and cell types in an unsupervised manner,
while performing comparably on multiple data sets to current state-of-the-art
model in traditional supervised ACTI.


  Construction of a scaffold structure that supports a desired motif,
conferring protein function, shows promise for the design of vaccines and
enzymes. But a general solution to this motif-scaffolding problem remains open.
Current machine-learning techniques for scaffold design are either limited to
unrealistically small scaffolds (up to length 20) or struggle to produce
multiple diverse scaffolds. We propose to learn a distribution over diverse and
longer protein backbone structures via an E(3)-equivariant graph neural
network. We develop SMCDiff to efficiently sample scaffolds from this
distribution conditioned on a given motif; our algorithm is the first to
theoretically guarantee conditional samples from a diffusion model in the
large-compute limit. We evaluate our designed backbones by how well they align
with AlphaFold2-predicted structures. We show that our method can (1) sample
scaffolds up to 80 residues and (2) achieve structurally diverse scaffolds for
a fixed motif.


  Much evidence from biological theory and empirical data indicates that, gene
tree, phylogenetic trees reconstructed from different genes (loci), do not have
to have exactly the same tree topologies. Such incongruence between gene trees
might be caused by some ``unusual'' evolutionary events, such as meiotic sexual
recombination in eukaryotes or horizontal transfers of genetic material in
prokaryotes. However, most of gene trees are constrained by the tree topology
of its species tree, that is, the phylogenetic tree of a given species
following their evolutionary history. In order to discover ``outlying'' gene
trees which do not follow the ``main distribution(s)'' of trees, we propose to
apply the ``tropical metric'' with the max-plus algebra from tropical geometry
to a non-parametric estimation of gene trees over the space of phylogenetic
trees. In this research we apply the ``tropical metric,'' a well-defined metric
over the space of phylogenetic trees under the max-plus algebra, to
non-parametric estimation of gene trees distribution over the tree space.
Kernel density estimator (KDE) is one of the most popular non-parametric
estimation of a distribution from a given sample, and we propose an analogue of
the classical KDE in the setting of tropical geometry with the tropical metric
which measures the length of an intrinsic geodesic between trees over the tree
space. We estimate the probability of an observed tree by empirical frequencies
of nearby trees, with the level of influence determined by the tropical metric.
Then, with simulated data generated from the multispecies coalescent model, we
show that the non-parametric estimation of gene tree distribution using the
tropical metric performs better than one using the Billera-Holmes-Vogtmann
(BHV) metric developed by Weyenberg et al. in terms of computational times and
accuracy. We then apply it to Apicomplexa data.


  This paper presents the foundation for a decomposition theory for Boolean
networks, a type of discrete dynamical system that has found a wide range of
applications in the life sciences, engineering, and physics. Given a Boolean
network satisfying certain conditions, there is a unique collection of
subnetworks so that the network can be reconstructed from these subnetworks by
an extension operation. The main result of the paper is that this structural
decomposition induces a corresponding decomposition of the network dynamics.
The theory is motivated by the search for a mathematical framework to formalize
the hypothesis that biological systems are modular, widely accepted in the life
sciences, but not well-defined and well-characterized. As an example of how
dynamic modularity could be used for the efficient identification of phenotype
control, the control strategies for the network can be found by identifying
controls in its modules, one at a time.


  Imaging biomarkers offer a non-invasive way to predict the response of
immunotherapy prior to treatment. In this work, we propose a novel type of deep
radiomic features (DRFs) computed from a convolutional neural network (CNN),
which capture tumor characteristics related to immune cell markers and overall
survival. Our study uses four MRI sequences (T1-weighted, T1-weighted
post-contrast, T2-weighted and FLAIR) with corresponding immune cell markers of
151 patients with brain tumor. The proposed method extracts a total of 180 DRFs
by aggregating the activation maps of a pre-trained 3D-CNN within labeled tumor
regions of MRI scans. These features offer a compact, yet powerful
representation of regional texture encoding tissue heterogeneity. A
comprehensive set of experiments is performed to assess the relationship
between the proposed DRFs and immune cell markers, and measure their
association with overall survival. Results show a high correlation between DRFs
and various markers, as well as significant differences between patients
grouped based on these markers. Moreover, combining DRFs, clinical features and
immune cell markers as input to a random forest classifier helps discriminate
between short and long survival outcomes, with AUC of 72\% and
p=2.36$\times$10$^{-5}$. These results demonstrate the usefulness of proposed
DRFs as non-invasive biomarker for predicting treatment response in patients
with brain tumors.


  We describe a kernel of size 9k-8 for the NP-hard problem of computing the
Tree Bisection and Reconnect (TBR) distance k between two unrooted binary
phylogenetic trees. We achieve this by extending the existing portfolio of
reduction rules with three novel new reduction rules. Two of the rules are
based on the idea of topologically transforming the trees in a
distance-preserving way in order to guarantee execution of earlier reduction
rules. The third rule extends the local neighbourhood approach introduced in
(Kelk and Linz, Annals of Combinatorics 24(3), 2020) to more global structures,
allowing new situations to be identified when deletion of a leaf definitely
reduces the TBR distance by one. The bound on the kernel size is tight up to an
additive term. Our results also apply to the equivalent problem of computing a
Maximum Agreement Forest (MAF) between two unrooted binary phylogenetic trees.
We anticipate that our results will be more widely applicable for computing
agreement-forest based dissimilarity measures.


  A set of discrete individual points located in an embedding continuum space
can be seen as percolating or non-percolating, depending on the radius of the
discs/spheres associated with each of them. This problem is relevant in
theoretical ecology to analyze, e.g., the spatial percolation of a tree species
in a tropical forest or a savanna. Here, we revisit the problem of aggregating
random points in continuum systems (from $2$ to $6-$dimensional Euclidean
spaces) to analyze the nature of the corresponding percolation transition in
spatial point processes. This problem finds a natural description in terms of
the canonical ensemble but not in the usual grand-canonical one, customarily
employed to describe percolation transitions. This leads us to analyze the
question of ensemble equivalence and study whether the resulting canonical
continuum percolation transition shares its universal properties with standard
percolation transitions, analyzing diverse homogeneous and heterogeneous
spatial point processes. We, therefore, provide a powerful tool to characterize
and classify a vast class of natural point patterns, revealing their
fundamental properties based on percolation phase transitions.


  All-atom and coarse-grained molecular dynamics are two widely used
computational tools to study the conformational states of proteins. Yet, these
two simulation methods suffer from the fact that without access to
supercomputing resources, the time and length scales at which these states
become detectable are difficult to achieve. One alternative to such methods is
based on encoding the atomistic trajectory of molecular dynamics as a shorthand
version devoid of physical particles, and then learning to propagate the
encoded trajectory through the use of artificial intelligence. Here we show
that a simple textual representation of the frames of molecular dynamics
trajectories as vectors of Ramachandran basin classes retains most of the
structural information of the full atomistic representation of a protein in
each frame, and can be used to generate equivalent atom-less trajectories
suitable to train different types of generative neural networks. In turn, the
trained generative models can be used to extend indefinitely the atom-less
dynamics or to sample the conformational space of proteins from their
representation in the models latent space. We define intuitively this
methodology as molecular dynamics without molecules, and show that it enables
to cover physically relevant states of proteins that are difficult to access
with traditional molecular dynamics.


  Retrosynthesis is a procedure where a target molecule is transformed into
potential reactants and thus the synthesis routes can be identified. Recently,
computational approaches have been developed to accelerate the design of
synthesis routes. In this paper, we develop a generative framework
$\mathsf{G^2Retro}$ for one-step retrosynthesis prediction. $\mathsf{G^2Retro}$
imitates the reversed logic of synthetic reactions. It first predicts the
reaction centers in the target molecules (products), identifies the synthons
needed to assemble the products, and transforms these synthons into reactants.
$\mathsf{G^2Retro}$ defines a comprehensive set of reaction center types, and
learns from the molecular graphs of the products to predict potential reaction
centers. To complete synthons into reactants, $\mathsf{G^2Retro}$ considers all
the involved synthon structures and the product structures to identify the
optimal completion paths, and accordingly attaches small substructures
sequentially to the synthons. Here we show that $\mathsf{G^2Retro}$ is able to
better predict the reactants for given products in the benchmark dataset than
the state-of-the-art methods.


  The topic of this paper is prevalence estimation from the perspective of
active information. Prevalence among tested individuals has an upward bias
under the assumption that individuals' willingness to be tested for the disease
increases with the strength of their symptoms. Active information due to
testing bias quantifies the degree at which the willingness to be tested
correlates with infection status. Interpreting incomplete testing as a missing
data problem, the missingness mechanism impacts the degree at which the bias of
the original prevalence estimate can be removed. The reduction in prevalence,
when testing bias is adjusted for, translates into an active information due to
bias correction, with opposite sign to active information due to testing bias.
Prevalence and active information estimates are asymptotically normal, a
behavior also illustrated through simulations.


  This paper carries out an analysis of the global properties of solutions of
an in-host model of hepatitis C for general values of its parameters. A
previously unknown stable steady state on the boundary of the positive orthant
is exhibited. It is proved that the model exhibits Hopf bifurcations and hence
periodic solutions. A general parametrization of positive steady states is
given and it is determined when the number of steady states is odd or even,
according to the value of a certain basic reproductive ratio. This implies, in
particular, that when this reproductive ratio is greater than one there always
exists at least one positive steady state. A positive steady state which
bifurcates from an infection-free state when the reproductive ratio passes
through one is always stable, i.e. no backward bifurcation occurs in this
model.


  Stemmatology is a subfield of philology where one approach to understand the
copy-history of textual variants of a text (witnesses of a tradition) is to
generate an evolutionary tree. Computational methods are partly shared between
the sister discipline of phylogenetics and stemmatology. In 2022, a surveypaper
in nature communications found that Deep Learning (DL), which otherwise has
brought about major improvements in many fields (Krohn et al 2020) has had only
minor successes in phylogenetics and that "it is difficult to conceive of an
end-to-end DL model to directly estimate phylogenetic trees from raw data in
the near future"(Sapoval et al. 2022, p.8). In stemmatology, there is to date
no known DL approach at all. In this paper, we present a new DL approach to
placement of manuscripts on a stemma and demonstrate its potential. This could
be extended to phylogenetics where the universal code of DNA might be an even
better prerequisite for the method using sequence to sequence based neural
networks in order to retrieve tree distances.


  Cell crawling on flat substrates is based on intracellular flows of the actin
cytoskeleton that are driven by both actin polymerization at the front and
myosin contractility at the back. The new experimental tool of optogenetics
makes it possible to spatially control contraction and thereby possibly also
cell migration. Here we analyze this situation theoretically using a
one-dimensional active gel model in which the excluded volume interactions of
myosin and their aggregation into minifilaments is modeled by a supercritical
van der Waals fluid. This physically simple and transparent, but nonlinear and
thermodynamically rigorous model predicts bistability between sessile and
motile solutions. We then show that one can switch between these two states at
realistic parameter ranges via optogenetic activation or inhibition of
contractility, in agreement with recent experiments. We also predict the
required activation strengths and initiation times.


  This paper presents the methods that have participated in the SHREC 2022
contest on protein-ligand binding site recognition. The prediction of
protein-ligand binding regions is an active research domain in computational
biophysics and structural biology and plays a relevant role for molecular
docking and drug design. The goal of the contest is to assess the effectiveness
of computational methods in recognizing ligand binding sites in a protein based
on its geometrical structure. Performances of the segmentation algorithms are
analyzed according to two evaluation scores describing the capacity of a
putative pocket to contact a ligand and to pinpoint the correct binding region.
Despite some methods perform remarkably, we show that simple
non-machine-learning approaches remain very competitive against data-driven
algorithms. In general, the task of pocket detection remains a challenging
learning problem which suffers of intrinsic difficulties due to the lack of
negative examples (data imbalance problem).


  Maps have played an important role in epidemiology and public health since
the beginnings of these disciplines. With the advent of geographical
information systems and advanced information visualization techniques,
interactive maps have become essential tools for the analysis of geographical
patterns of disease incidence and prevalence, as well as communication of
public health knowledge, as dramatically illustrated by the proliferation of
web-based maps and disease surveillance ``dashboards'' during the COVID-19
pandemic. While such interactive maps are usually effective in supporting
static spatial analysis, support for spatial epidemiological visualization and
modelling involving distributed and dynamic data sources, and support for
analysis of temporal aspects of disease spread have proved more challenging.
Combining these two aspects can be crucial in applications of interactive maps
in epidemiology and public health work. In this paper, we discuss these issues
in the context of support for disease surveillance in remote regions, including
tools for distributed data collection, simulation and analysis, and enabling
multidisciplinary collaboration.


  Finding cancer driver genes has been a focal theme of cancer research and
clinical studies. One of the recent approaches is based on network structural
controllability that focuses on finding a control scheme and driver genes that
can steer the cell from an arbitrary state to a designated state. While
theoretically sound, this approach is impractical for many reasons, e.g., the
control scheme is often not unique and half of the nodes may be driver genes
for the cell. We developed a novel approach that transcends structural
controllability. Instead of considering driver genes for one control scheme, we
considered control hub genes that reside in the middle of a control path of
every control scheme. Control hubs are the most vulnerable spots for
controlling the cell and exogenous stimuli on them may render the cell
uncontrollable. We adopted control hubs as cancer-keep genes (CKGs) and applied
them to a gene regulatory network of bladder cancer (BLCA). All the genes on
the cell cycle and p53 singling pathways in BLCA are CKGs, confirming the
importance of these genes and the two pathways in cancer. A smaller set of 35
sensitive CKGs (sCKGs) for BLCA was identified by removing network links. Six
sCKGs (RPS6KA3, FGFR3, N-cadherin (CDH2), EP300, caspase-1, and FN1) were
subjected to small-interferencing-RNA knockdown in four cell lines to validate
their effects on the proliferation or migration of cancer cells. Knocking down
RPS6KA3 in a mouse model of BLCA significantly inhibited the growth of tumor
xenografts in the mouse model. Combined, our results demonstrated the value of
CKGs as therapeutic targets for cancer therapy and the potential of CKGs as an
effective means for studying and characterizing cancer etiology.


  Solving population balance equations, we derive analytical steady-state cell
size distributions for single-lineage experiments, such as the mother machine.
These experiments are fundamentally different from batch cultures where
populations of cells grow freely, and the statistical bias between them is
obtained by comparing our results to cell size distributions measured in
population. For exponential single-cell growth, characterizing most bacteria,
the lineage-population bias is obtained explicitly. In addition, if volume is
evenly split between the daughter cells at division, we show that cells are on
average smaller in populations. For more general power-law growth rates and
deterministic volume partitioning, both symmetric and asymmetric, we derive the
exact lineage distribution. This solution is in good agreement with E. Coli
mother machine data, and can be used to infer cell cycle parameters, such as
the strength of the size control and the asymmetry of the division. When
introducing stochastic volume partitioning, we derive the large-size and
small-size tails of the lineage size distributions, and show that they
respectively do not depend on the partitioning of volume and on the size
control strength. Finally, we show that introducing noise, either on the volume
partitioning or on the single cell growth rate, can cancel the
lineage-population bias.


  Reproducibility and replicability of research findings are central to the
scientific integrity of epidemiology. In addition, many research questions
require combiningdata from multiple sources to achieve adequate statistical
power. However, barriers related to confidentiality, costs, and incentives
often limit the extent and speed of sharing resources, both data and code.
Epidemiological practices that follow FAIR principles can address these
barriers by making resources (F)indable with the necessary metadata ,
(A)ccessible to authorized users and (I)nteroperable with other data, to
optimize the (R)e-use of resources with appropriate credit to its creators. We
provide an overview of these principles and describe approaches for
implementation in epidemiology. Increasing degrees of FAIRness can be achieved
by moving data and code from on-site locations to the Cloud, using
machine-readable and non-proprietary files, and developing open-source code.
Adoption of these practices will improve daily work and collaborative analyses,
and facilitate compliance with data sharing policies from funders and
scientific journals. Achieving a high degree of FAIRness will require funding,
training, organizational support, recognition, and incentives for sharing
resources. But these costs are amply outweighed by the benefits of making
research more reproducible, impactful, and equitable by facilitating the re-use
of precious research resources by the scientific community.


  Cumulative cultural evolution occurs when adaptive innovations are passed
down to consecutive generations through social learning. This process has
shaped human technological innovation, but also occurs in non-human species.
While it is traditionally argued that cumulative culture relies on
high-fidelity social transmission and advanced cognitive skills, here I show
that a much simpler system suffices. Cumulative culture spontaneously emerged
in artificial agents who navigate with a minimal cognitive architecture of
goal-direction, social proximity, and route memory. Within each generation,
naive individuals benefitted from being paired with experienced navigators
because they could follow previously established routes. Crucially, experienced
navigators also benefitted from the presence of naive individuals through
regression to the goal. As experienced agents followed their memorised path,
their naive counterparts (unhindered by route memory) were more likely to err
towards than away from the goal, and thus biased the pair in that direction.
This improved route efficiency within each generation. In control experiments,
cumulative culture was attenuated when agents' social proximity or route memory
were lesioned, whereas eliminating goal-direction only reduced efficiency.
These results demonstrate that cumulative cultural evolution occurs even in the
absence of sophisticated communication or thought. One interpretation of this
finding is that current definitions are too loose, and should be narrowed. An
alternative conclusion is that rudimentary cumulative culture is an emergent
property of systems that seek social proximity and have an imprecise memory
capacity, providing a flexible complement to traditional evolutionary
mechanisms.


  Biological systems are very robust to morphological damage, but artificial
systems (robots) are currently not. In this paper we present a system based on
neural cellular automata, in which locomoting robots are evolved and then given
the ability to regenerate their morphology from damage through gradient-based
training. Our approach thus combines the benefits of evolution to discover a
wide range of different robot morphologies, with the efficiency of supervised
training for robustness through differentiable update rules. The resulting
neural cellular automata are able to grow virtual robots capable of regaining
more than 80\% of their functionality, even after severe types of morphological
damage.


  Dry-land ecosystem has turned into a matter of grave concern, due to growing
threat of land degradation and bioproductivity-loss. Self-organized vegetation
patterns are a remarkable characteristic of these ecosystems; apart from being
visually captivating, patterns modulate the system-response to increasing
environmental stress. Empirical studies hinted that herbivory is one the key
regulatory mechanism behind pattern formation and overall ecosystem
functioning. However most of the mathematical models have taken a mean-field
strategy to grazing; foraging has been considered to be independent of spatial
distribution of vegetation. To this end, an extended version of the celebrated
plant-water model due to Klausmeier, has been taken as the base here. To
encompass the effect of heterogeneous vegetation distribution on foraging
intensity and subsequent impact on entire ecosystem, grazing is considered here
to depend on spatially weighted average vegetation density, instead of density
at a particular point. Moreover, varying influence of vegetation at any
location over gazing elsewhere, is incorporated by choosing suitable averaging
function. A comprehensive analysis demonstrates that inclusion of spatial
non-locality, alters the understanding of system dynamics significantly. The
grazing ecosystem is found to be more resilient to increasing aridity than it
was anticipated to be in earlier studies on non-local grazing. The
system-response to rising environmental pressure is also observed to vary
depending on the grazer. Obtained results also suggest possibility of
multi-stability due the history-dependence of system-response. Overall, this
work indicates that the spatial heterogeneity in grazing intensity has a
decisive role to play in the functioning of water-limited ecosystems.


  A system of 4 nonlinearly-coupled Ordinary Differential Equations has been
recently introduced to investigate the evolution of human respiratory virus
epidemics. In this paper we point out that some explicit solutions of that
system can be obtained by algebraic operations, provided the parameters of the
model satisfy certain constraints.


  Efficient and effective drug-target binding affinity (DTBA) prediction is a
challenging task due to the limited computational resources in practical
applications and is a crucial basis for drug screening. Inspired by the good
representation ability of graph neural networks (GNNs), we propose a
simple-structured GNN model named SS-GNN to accurately predict DTBA. By
constructing a single undirected graph based on a distance threshold to
represent protein-ligand interactions, the scale of the graph data is greatly
reduced. Moreover, ignoring covalent bonds in the protein further reduces the
computational cost of the model. The GNN-MLP module takes the latent feature
extraction of atoms and edges in the graph as two mutually independent
processes. We also develop an edge-based atom-pair feature aggregation method
to represent complex interactions and a graph pooling-based method to predict
the binding affinity of the complex. We achieve state-of-the-art prediction
performance using a simple model (with only 0.6M parameters) without
introducing complicated geometric feature descriptions. SS-GNN achieves
Pearson's Rp=0.853 on the PDBbind v2016 core set, outperforming
state-of-the-art GNN-based methods by 5.2%. Moreover, the simplified model
structure and concise data processing procedure improve the prediction
efficiency of the model. For a typical protein-ligand complex, affinity
prediction takes only 0.2 ms. All codes are freely accessible at
https://github.com/xianyuco/SS-GNN.


  Social insects are ecologically and evolutionarily most successful organisms
on earth, which can achieve robust collective behaviors through local
interactions among group members. Colony migration has been considered as a
leading example of collective decision-making in social insects. In this paper,
a piecewise colony migration system with recruitment switching is proposed to
explore underlying mechanisms and synergistic effects of colony size and quorum
on the outcomes of collective decision. The completed dynamical analysis for
the non-smooth system (including the dynamics on subsystems, switching surface,
and full system) is performed, and the sufficient conditions significantly
related to colony size for the stability of equilibria are provided. The
theoretical results suggest that large colonies are more likely to emigrate to
a new site. More interesting findings include but not limit to: (a) the system
may exhibit oscillation when the colony size is below a critical level; (b) the
system may also exhibit a bistable state, and colonies migrate to a new site or
the old nest depending on their initial sizes of recruiters. Bifurcation
analysis shows that the variations of colony size and quorum threshold greatly
impact the dynamics. The results suggest that it is important to distinguish
two populations of recruiters in modeling. This work may provide important
insights on how simple and local interactions achieve the collective migrating
activity in social insects.


  Sea lice are a threat to the health of both wild and farmed salmon and an
economic burden for salmon farms. With a free living larval stage, sea lice can
disperse tens of kilometers in the ocean between salmon farms, leading to
connected sea lice populations that are difficult to control in isolation. In
this paper we develop a simple analytical model for the dispersal of sea lice
between two salmon farms. From the model we calculate the arrival time
distribution of sea lice dispersing between farms, as well as the level of
cross-infection of sea lice. We also use numerical flows from a hydrodynamic
model, coupled with a particle tracking model, to directly calculate the
arrival time of sea lice dispersing between two farms in the Broughton
Archipelago, BC, in order to fit our analytical model and find realistic
parameter estimates. Using the parametrized analytical model we show that there
is often an intermediate inter-farm spacing that maximizes the level of cross
infection between farms, and that increased temperatures will lead to increased
levels of cross infection.


  Genomic selection has increased genetic gain in several livestock species,
but due to the complicated genetics and reproduction biology not yet in honey
bees. Recently, 2 970 queens were genotyped to gather a reference population.
For the application of genomic selection in honey bees, this study analyses the
predictive ability and bias of pedigree-based and genomic breeding values for
honey yield, three workability traits and two traits for resistance against the
parasite Varroa destructor. For breeding value estimation, we use a honey
bee-specific model with maternal and direct effects, to account for the
contributions of the workers and the queen of a colony to the phenotypes. We
conducted a validation for the last generation and a five-fold
cross-validation. In the validation for the last generation, the predictive
ability of pedigree-based estimated breeding values was 0.06 for honey yield,
and ranged from 0.2 to 0.41 for the workability traits. The inclusion of
genomic marker data improved these predictive abilities to 0.11 for honey
yield, and a range from 0.22 to 0.44 for the workability traits. The inclusion
of genomic data did not improve the predictive ability for the disease related
traits. Traits with high heritability for maternal effects compared to the
heritability for direct effects showed the most promising results. Across all
traits, the bias with genomic methods was close to the bias with pedigree-based
BLUP. The results show that genomic selection can successfully be applied to
honey bees.


  Deriving emergent patterns from models of biological processes is a core
concern of mathematical biology. In the context of partial differential
equations (PDEs), these emergent patterns sometimes appear as local minimisers
of a corresponding energy functional. Here we give methods for determining the
qualitative structure of local minimum energy states of a broad class of
multi-species nonlocal advection-diffusion models, recently proposed for
modelling the spatial structure of ecosystems. We show that when each pair of
species respond to one another in a symmetric fashion (i.e. via mutual
avoidance or mutual attraction, with equal strength), the system admits an
energy functional that decreases in time and is bounded below. This suggests
that the system will eventually reach a local minimum energy steady state,
rather than fluctuating in perpetuity. We leverage this energy functional to
develop tools, including a novel application of computational algebraic
geometry, for making conjectures about the number and qualitative structure of
local minimum energy solutions. These conjectures give a guide as to where to
look for numerical steady state solutions, which we verify through numerical
analysis. Our technique shows that even with two species, multi-stability with
up to four classes of local minimum energy state can emerge. The associated
dynamics include spatial sorting via aggregation and repulsion both within and
between species. The emerging spatial patterns include a mixture of
territory-like segregation as well as narrow spike-type solutions. Overall, our
study reveals a general picture of rich multi-stability in systems of moving
and interacting species.


  A well-known limitation of existing molecular generative models is that the
generated molecules highly resemble those in the training set. To generate
truly novel molecules that may have even better properties for de novo drug
discovery, more powerful exploration in the chemical space is necessary. To
this end, we propose Molecular Out-Of-distribution Diffusion(MOOD), a
score-based diffusion scheme that incorporates out-of-distribution (OOD)
control in the generative stochastic differential equation (SDE) with simple
control of a hyperparameter, thus requires no additional costs. Since some
novel molecules may not meet the basic requirements of real-world drugs, MOOD
performs conditional generation by utilizing the gradients from a property
predictor that guides the reverse-time diffusion process to high-scoring
regions according to target properties such as protein-ligand interactions,
drug-likeness, and synthesizability. This allows MOOD to search for novel and
meaningful molecules rather than generating unseen yet trivial ones. We
experimentally validate that MOOD is able to explore the chemical space beyond
the training distribution, generating molecules that outscore ones found with
existing methods, and even the top 0.01% of the original training pool. Our
code is available at https://github.com/SeulLee05/MOOD.


  Topological protection has emerged as an organizing principle for
understanding and engineering robust collective behavior in electronic and
material systems. Recent work suggests that topology may also play a role in
organizing stochastic processes relevant to biology and self-assembly. Here, we
show that topological protection in chemical networks can be understood
entirely in terms of nonequilibrium thermodynamics. We illustrate these ideas
using simple examples inspired by the literature.


  A common theme in causal inference is learning causal relationships between
observed variables, also known as causal discovery. This is usually a daunting
task, given the large number of candidate causal graphs and the combinatorial
nature of the search space. Perhaps for this reason, most research has so far
focused on relatively small causal graphs, with up to hundreds of nodes.
However, recent advances in fields like biology enable generating experimental
data sets with thousands of interventions followed by rich profiling of
thousands of variables, raising the opportunity and urgent need for large
causal graph models. Here, we introduce the notion of factor directed acyclic
graphs (f-DAGs) as a way to restrict the search space to non-linear low-rank
causal interaction models. Combining this novel structural assumption with
recent advances that bridge the gap between causal discovery and continuous
optimization, we achieve causal discovery on thousands of variables.
Additionally, as a model for the impact of statistical noise on this estimation
procedure, we study a model of edge perturbations of the f-DAG skeleton based
on random graphs and quantify the effect of such perturbations on the f-DAG
rank. This theoretical analysis suggests that the set of candidate f-DAGs is
much smaller than the whole DAG space and thus may be more suitable as a search
space in the high-dimensional regime where the underlying skeleton is hard to
assess. We propose Differentiable Causal Discovery of Factor Graphs (DCD-FG), a
scalable implementation of -DAG constrained causal discovery for
high-dimensional interventional data. DCD-FG uses a Gaussian non-linear
low-rank structural equation model and shows significant improvements compared
to state-of-the-art methods in both simulations as well as a recent large-scale
single-cell RNA sequencing data set with hundreds of genetic interventions.


  The ability of biomolecules to exert forces on their surroundings or resist
compression from the environment is essential in a variety of biologically
relevant contexts. As has been understood for centuries, slender rods can only
be compressed so far until they buckle, adopting an intrinsically bent state
that may be unable to bear a compressive load. In the low-temperature limit and
for a constant compressive force, Euler buckling theory predicts a sudden
transition from a compressed to a bent state in these slender rods. In this
paper, we use a mean-field theory to show that if a semiflexible chain is
compressed at a finite temperature with a fixed end-to-end distance (permitting
fluctuations in the compressive forces), it exhibits a continuous phase
transition to a buckled state at a critical level of compression, and we
determine a quantitatively accurate prediction of the transverse position
distribution function of the midpoint of the chain that indicates the
transition. We find the mean compressive forces are non-monotonic as the
extension of the filament varies, consistent with the observation that strongly
buckled filaments are less able to bear an external load. We also find that for
the fixed extension (isometric) ensemble that the buckling transition does not
coincide with the local minimum of the mean force (in contrast to Euler
buckling). We also show the theory is highly sensitive to fluctuations in
length, and that the buckling transition can still be accurately recovered by
accounting for those fluctuations. These predictions may be useful in
understanding the behavior of filamentous biomolecules compressed by
fluctuating forces, relevant in a variety of biological contexts.


  Systems driven by Brownian motion are ubiquitous. A prevailing challenge is
inferring, from data, the diffusion and kinetic parameters that describe these
stochastic processes. In this work, we investigate a multi-state diffusion
process that arises in the context of single particle tracking (SPT), wherein
the motion of a particle is governed by a discrete set of diffusive states, and
the tendency of the particle to switch between these states is modeled as a
random process. We consider two models for this behavior: a mixture model and a
hidden Markov model (HMM). For both, we adopt a Bayesian approach to sample the
distributions of the underlying parameters and implement a Markov Chain Monte
Carlo (MCMC) scheme to compute the posterior distributions, as in Das, Cairo,
Coombs (2009). The primary contribution of this work is a study of the
robustness of this method to infer parameters of a three-state HMM, and a
discussion of the challenges and degeneracies that arise from considering three
states. Finally, we investigate the problem of determining the number of
diffusive states using model selection criteria. We present results from
simulated data that demonstrate proof of concept, as well as apply our method
to experimentally measured single molecule diffusion trajectories of monomeric
phosphoinositide-dependent kinase-1 (PDK1) on a synthetic target membrane where
it can associate with its binding partner protein kinase C alpha isoform
(PKC$\alpha$) to form a heterodimer detected by its significantly lower
diffusivity.
  All matlab software is available here:
\url{https://github.com/MathBioCU/SingleMolecule}


  Bayesian inference for phylogenetics is a gold standard for computing
distributions of phylogenies. It faces the challenging problem of. moving
throughout the high-dimensional space of trees. However, hyperbolic space
offers a low dimensional representation of tree-like data. In this paper, we
embed genomic sequences into hyperbolic space and perform hyperbolic Markov
Chain Monte Carlo for Bayesian inference. The posterior probability is computed
by decoding a neighbour joining tree from proposed embedding locations. We
empirically demonstrate the fidelity of this method on eight data sets. The
sampled posterior distribution recovers the splits and branch lengths to a high
degree. We investigated the effects of curvature and embedding dimension on the
Markov Chain's performance. Finally, we discuss the prospects for adapting this
method to navigate tree space with gradients.


  Accurate and efficient prediction of the molecular properties of drugs is one
of the fundamental problems in drug research and development. Recent
advancements in representation learning have been shown to greatly improve the
performance of molecular property prediction. However, due to limited labeled
data, supervised learning-based molecular representation algorithms can only
search limited chemical space, which results in poor generalizability. In this
work, we proposed a self-supervised representation learning framework for
large-scale unlabeled molecules. We developed a novel molecular graph
augmentation strategy, referred to as attention-wise graph mask, to generate
challenging positive sample for contrastive learning. We adopted the graph
attention network (GAT) as the molecular graph encoder, and leveraged the
learned attention scores as masking guidance to generate molecular augmentation
graphs. By minimization of the contrastive loss between original graph and
masked graph, our model can capture important molecular structure and
higher-order semantic information. Extensive experiments showed that our
attention-wise graph mask contrastive learning exhibit state-of-the-art
performance in a couple of downstream molecular property prediction tasks.


  This paper studies the features of a homopolymer translocating through a
flexible pore. The channel is modeled as a monolayer tube composed by monomers
with two elastic parameters: spring-like two body interaction and bending three
body recall interaction. In order to guarantee the stability of the system, the
membrane is compounded by a lipid bilayer structure having hydrophobic body
(internal), while the pore is hydrophilic in both edges. The polymer is
end-pulled from the cis-side to the trans-side by a cantilever, to which is
connected through a spring able to measure the force acting on the polymer
during the translocation. All the structure reacts to the impacts of the
monomers of the polymer with vibrations generated by the movement of its
constituent bodies. In these conditions, the work done by the cantilever shows
a nonmonotonic behavior with the elastic constant, revealing a resonant-like
behavior in a parameter region. Moreover, the force spectroscopy registered as
a function of time, is able to record the main kinetics of the polymer
progression inside the pore.


  This paper uses new and existing methods to study collective trends across
countries throughout the pandemic, with a focus on the multivariate time series
of reproduction numbers and vaccine proliferation. We begin with a time-varying
analysis of the collective nature of infectivity, where we evaluate the
eigenspectrum and collective magnitude of reproduction number time series on a
country-by-country basis. Next, we study the topology of this eigenspectrum,
measuring the deviation between all points in time, and introduce a
graph-theoretic methodology to reveal a clear partition in global infectivity
dynamics. Then, we compare countries' vaccine rollouts with economic indicators
such as their GDP and HDI in a collective fashion. We investigate time-varying
consistency and determine points in time where there is the greatest
discrepancy between these indicators as a whole. Our two primary findings are a
considerable increase in collective infectivity in the latter half of the
period, and a concave-up (``down then up'') pattern in the collective
consistency between vaccine coverage and economic/development indicators across
countries.


  Multisite phosphorylation is a signaling mechanism well known to give rise to
multiple steady states, a property termed multistationarity. When
phosphorylation occurs in a sequential and distributive manner, we obtain a
family of networks indexed by the number of phosphorylation sites $n$. This
work addresses the problem of understanding the parameter region where this
family of networks displays multistationarity, by focusing on the projection of
this region onto the set of kinetic parameters. The problem is phrased in the
context of real algebraic geometry and reduced to studying whether a
polynomial, defined as the determinant of a parametric matrix of size three,
attains negative values over the positive orthant. The coefficients of the
polynomial are functions of the kinetic parameters. For any $n$, we provide
sufficient conditions for the polynomial to be positive and hence, preclude
multistationarity, and also sufficient conditions for it to attain negative
values and hence, enable multistationarity. These conditions are derived by
exploiting the structure of the polynomial, its Newton polytope, and employing
circuit polynomials. A relevant consequence of our results is that the set of
kinetic parameters that enable or preclude multistationarity are both connected
for all $n$.


  Diversification models describe the random growth of evolutionary trees,
modeling the historical relationships of species through speciation and
extinction events. One class of such models allows for independently changing
traits, or types, of the species within the tree, upon which speciation and
extinction rates depend. Although identifiability of parameters is necessary to
justify parameter estimation with a model, it has not been formally established
for these models, despite their adoption for inference. This work establishes
generic identifiability up to label swapping for the parameters of one of the
simpler forms of such a model, a multitype pure birth model of speciation, from
an asymptotic distribution derived from a single tree observation as its depth
goes to infinity. Crucially for applications to available data, no observation
of types is needed at any internal points in the tree, nor even at the leaves.


  We show that a number of models in virus dynamics, epidemiology and plant
biology can be presented as ``damped" versions of the Lotka-Volterra
predator-prey model, by analogy to the damped harmonic oscillator. The analogy
deepens with the use of Lyapunov functions, which allow us to characterize
their dynamics and even make some estimates.


  The Turing completeness of continuous chemical reaction networks (CRNs)
states that any computable real function can be computed by a continuous CRN on
a finite set of molecular species, possibly restricted to elementary reactions,
i.e. with at most two reactants and mass action law kinetics. In this paper, we
introduce a notion of online analog computation for the CRNs that stabilize the
concentration of their output species to the result of some function of the
concentration values of their input species, whatever changes are operated on
the inputs during the computation. We prove that the set of real functions
stabilized by a CRN with mass action law kinetics is precisely the set of real
algebraic functions.


  Bacteriophages (or phages), which infect bacteria, have two distinct
lifestyles: virulent and temperate. Predicting the lifestyle of phages helps
decipher their interactions with their bacterial hosts, aiding phages'
applications in fields such as phage therapy. Because experimental methods for
annotating the lifestyle of phages cannot keep pace with the fast accumulation
of sequenced phages, computational method for predicting phages' lifestyles has
become an attractive alternative. Despite some promising results, computational
lifestyle prediction remains difficult because of the limited known annotations
and the sheer amount of sequenced phage contigs assembled from metagenomic
data. In particular, most of the existing tools cannot precisely predict
phages' lifestyles for short contigs. In this work, we develop PhaTYP (Phage
TYPe prediction tool) to improve the accuracy of lifestyle prediction on short
contigs. We design two different training tasks, self-supervised and
fine-tuning tasks, to overcome lifestyle prediction difficulties. We rigorously
tested and compared PhaTYP with four state-of-the-art methods: DeePhage,
PHACTS, PhagePred, and BACPHLIP. The experimental results show that PhaTYP
outperforms all these methods and achieves more stable performance on short
contigs. In addition, we demonstrated the utility of PhaTYP for analyzing the
phage lifestyle on human neonates' gut data. This application shows that PhaTYP
is a useful means for studying phages in metagenomic data and helps extend our
understanding of microbial communities.


  Accurate prediction of Drug-Target Affinity (DTA) is of vital importance in
early-stage drug discovery, facilitating the identification of drugs that can
effectively interact with specific targets and regulate their activities. While
wet experiments remain the most reliable method, they are time-consuming and
resource-intensive, resulting in limited data availability that poses
challenges for deep learning approaches. Existing methods have primarily
focused on developing techniques based on the available DTA data, without
adequately addressing the data scarcity issue. To overcome this challenge, we
present the SSM-DTA framework, which incorporates three simple yet highly
effective strategies: (1) A multi-task training approach that combines DTA
prediction with masked language modeling (MLM) using paired drug-target data.
(2) A semi-supervised training method that leverages large-scale unpaired
molecules and proteins to enhance drug and target representations. This
approach differs from previous methods that only employed molecules or proteins
in pre-training. (3) The integration of a lightweight cross-attention module to
improve the interaction between drugs and targets, further enhancing prediction
accuracy. Through extensive experiments on benchmark datasets such as
BindingDB, DAVIS, and KIBA, we demonstrate the superior performance of our
framework. Additionally, we conduct case studies on specific drug-target
binding activities, virtual screening experiments, drug feature visualizations,
and real-world applications, all of which showcase the significant potential of
our work. In conclusion, our proposed SSM-DTA framework addresses the data
limitation challenge in DTA prediction and yields promising results, paving the
way for more efficient and accurate drug discovery processes. Our code is
available at $\href{https://github.com/QizhiPei/SSM-DTA}{Github}$.


  Transfer learning has emerged as a powerful technique in many application
problems, such as computer vision and natural language processing. However,
this technique is largely ignored in application to genetic data analysis. In
this paper, we combine transfer learning technique with a neural network based
method(expectile neural networks). With transfer learning, instead of starting
the learning process from scratch, we start from one task that have been
learned when solving a different task. We leverage previous learnings and avoid
starting from scratch to improve the model performance by passing information
gained in different but related task. To demonstrate the performance, we run
two real data sets. By using transfer learning algorithm, the performance of
expectile neural networks is improved compared to expectile neural network
without using transfer learning technique.


  Objective: The reservoir of human immunodeficiency virus (HIV) latently
infected cells is the major obstacle for eradication of acquired
immunodeficiency syndrome (AIDS). Due to the noisy environment and multiple
influencing factors in the organism, current dynamical models cannot reach a
common understanding of the molecular mechanism of HIV latency. In this work,
through a new dynamical structure decomposition, the deterministic part of the
equation can be separated from the stochastic noise. Thus, the fixed-point
analysis of ordinary differential equation is enough to obtain the different
steady states of the system. Methods: We established a dynamical model of HIV
transcription process by using continuous stochastic differential equations,
which simplifies the dimensions of equations needed to describe the system and
increases the explorable space of the model. Different states between latency
and activation of virus and their relations can be intuitively represented by
potential functions and probability distribution functions. Results: Based on
our model, the influence of different dynamical parameters on stability is
quantitatively analyzed, the parameter ranges of the system in bistable and
monostable states are obtained respectively. The theoretical basis of this work
is verified by comparing the effects of different factors on dynamic
bifurcation with the results of biological experiments. Conclusion: This paper
goes beyond previous discrete stochastic methods, and can quantitatively
analyze the dynamic mechanism of HIV transcriptional regulation through
ordinary differential equations, which is beneficial to the promotion to deal
with the high-dimensional situation, and further study the occurrence and
development of AIDS in vivo, so as to guide the design of experiments and
search for clinical treatment.


  Infectious diseases remain among the top contributors to human illness and
death worldwide, among which many diseases produce epidemic waves of infection.
The unavailability of specific drugs and ready-to-use vaccines to prevent most
of these epidemics makes the situation worse. These force public health
officials and policymakers to rely on early warning systems generated by
reliable and accurate forecasts of epidemics. Accurate forecasts of epidemics
can assist stakeholders in tailoring countermeasures, such as vaccination
campaigns, staff scheduling, and resource allocation, to the situation at hand,
which could translate to reductions in the impact of a disease. Unfortunately,
most of these past epidemics exhibit nonlinear and non-stationary
characteristics due to their spreading fluctuations based on seasonal-dependent
variability and the nature of these epidemics. We analyse a wide variety of
epidemic time series datasets using a maximal overlap discrete wavelet
transform (MODWT) based autoregressive neural network and call it EWNet model.
MODWT techniques effectively characterize non-stationary behavior and seasonal
dependencies in the epidemic time series and improve the nonlinear forecasting
scheme of the autoregressive neural network in the proposed ensemble wavelet
network framework. From a nonlinear time series viewpoint, we explore the
asymptotic stationarity of the proposed EWNet model to show the asymptotic
behavior of the associated Markov Chain. We also theoretically investigate the
effect of learning stability and the choice of hidden neurons in the proposal.
From a practical perspective, we compare our proposed EWNet framework with
several statistical, machine learning, and deep learning models. Experimental
results show that the proposed EWNet is highly competitive compared to the
state-of-the-art epidemic forecasting methods.


  Identifiability and sloppiness are investigated in this paper for the
parameters of a descriptor system based on its frequency response samples. Two
metrics are suggested respectively for measuring absolute and relative
sloppiness of the parameter vector at a prescribed value. In this descriptor
system, system matrices are assumed to depend on its parameters through a
linear fractional transformation (LFT). When an associated transfer function
matrix (TFM) is of full normal row rank, a matrix rank based necessary and
sufficient condition is derived for parameter identifiability with a set of
finitely many frequency responses. This condition can be verified recursively
which is computationally quite appealing, especially when the system is of a
large scale. From this condition, an algorithm is suggested to find a set of
frequencies with which the frequency responses of the system are capable to
uniquely determine its parameters. An ellipsoid approximation is given for the
set consisting of all the parameter values with which the associated descriptor
system has a frequency response that deviates within a prescribed distance,
from that corresponding to a globally identifiable parameter vector value.
Explicit formulas are also derived for the suggested absolute and relative
sloppiness metrics.


  We develop a spatially realistic model of mutualistic metacommunities that
exploits the joint structure of spatial and interaction networks. This model
exhibits a sharp transition between a stable non-null equilibrium state and a
global extinction state. This behaviour allows defining a threshold on
colonisation/extinction parameters for the long-term metacommunity persistence.
This threshold, the 'metacommunity capacity', extends the metapopulation
capacity concept and can be calculated from the spatial and interaction
networks without simulating the whole dynamics. In several applications we
illustrate how the joint structure of the spatial and the interaction networks
affects metacommunity capacity. It results that a weakly modular spatial
network and a power-law degree distribution of the interaction network provide
the most favourable configuration for the long-term persistence of a
mutualistic metacommunity. Our model that encodes several explicit ecological
assumptions should pave the way for a larger exploration of spatially realistic
metacommunity models involving multiple interaction types.


  Aggregation is a common behavior by which groups of organisms arrange into
cohesive groups. Whether suspended in the air (like honey bee clusters), built
on the ground (such as army ant bridges), or immersed in water (such as sludge
worm blobs), these collectives serve a multitude of biological functions, from
protection against predation to the ability to maintain a relatively desirable
local environment despite a variable ambient environment. In this review, we
survey dense aggregations of a variety of insects, other arthropods, and worms
from a soft matter standpoint. An aggregation can be orders of magnitude larger
than its individual organisms, consisting of tens to hundreds of thousands of
individuals, and yet functions as a coherent entity. Understanding how
aggregating organisms coordinate with one another to form a superorganism
requires an interdisciplinary approach. We discuss how the physics of the
aggregation can yield additional insights to those gained from ecological and
physiological considerations, given that the aggregating individuals exchange
information, energy, and matter continually with the environment and one
another. While the connection between animal aggregations and the physics of
non-living materials has been proposed since the early 1900s, the recent advent
of physics of behavior studies provides new insights into social interactions
governed by physical principles. Current efforts focus on eusocial insects;
however, we show that these may just be the tip of an iceberg of superorganisms
that take advantage of physical interactions and simple behavioral rules to
adapt to changing environments. By bringing attention to a wide range of
invertebrate aggregations, we wish to inspire a new generation of scientists to
explore collective dynamics and bring a deeper understanding of the physics of
dense living aggregations.


  We study optimal two-sector (vegetative and reproductive) allocation models
of annual plants in temporally variable environments, that incorporate effects
of density dependent lifetime variability and juvenile mortality in a fitness
function whose expected value is maximized. Only special cases of arithmetic
and geometric mean maximizers have previously been considered in the
literature, and we also allow a wider range of production functions with
diminishing returns. The model predicts that the time of maturity is pushed to
an earlier date as the correlation between individual lifetimes increases, and
while optimal schedules are bang-bang at the extremes, the transition is
mediated by schedules where vegetative growth is mixed with reproduction for a
wide intermediate range. The mixed growth lasts longer when the production
function is less concave allowing for better leveraging of plant size when
generating seeds. Analytic estimates are obtained for the power means that
interpolate between arithmetic and geometric mean and correspond to partially
correlated lifetime distributions.


  Machine learning (ML) offers powerful methods for detecting and modeling
associations often in data with large feature spaces and complex associations.
Many useful tools/packages (e.g. scikit-learn) have been developed to make the
various elements of data handling, processing, modeling, and interpretation
accessible. However, it is not trivial for most investigators to assemble these
elements into a rigorous, replicatable, unbiased, and effective data analysis
pipeline. Automated machine learning (AutoML) seeks to address these issues by
simplifying the process of ML analysis for all. Here, we introduce STREAMLINE,
a simple, transparent, end-to-end AutoML pipeline designed as a framework to
easily conduct rigorous ML modeling and analysis (limited initially to binary
classification). STREAMLINE is specifically designed to compare performance
between datasets, ML algorithms, and other AutoML tools. It is unique among
other autoML tools by offering a fully transparent and consistent baseline of
comparison using a carefully designed series of pipeline elements including:
(1) exploratory analysis, (2) basic data cleaning, (3) cross validation
partitioning, (4) data scaling and imputation, (5) filter-based feature
importance estimation, (6) collective feature selection, (7) ML modeling with
`Optuna' hyperparameter optimization across 15 established algorithms
(including less well-known Genetic Programming and rule-based ML), (8)
evaluation across 16 classification metrics, (9) model feature importance
estimation, (10) statistical significance comparisons, and (11) automatically
exporting all results, plots, a PDF summary report, and models that can be
easily applied to replication data.


  Proteins are essential component of human life and their structures are
important for function and mechanism analysis. Recent work has shown the
potential of AI-driven methods for protein structure prediction. However, the
development of new models is restricted by the lack of dataset and benchmark
training procedure. To the best of our knowledge, the existing open source
datasets are far less to satisfy the needs of modern protein sequence-structure
related research. To solve this problem, we present the first million-level
protein structure prediction dataset with high coverage and diversity, named as
PSP. This dataset consists of 570k true structure sequences (10TB) and 745k
complementary distillation sequences (15TB). We provide in addition the
benchmark training procedure for SOTA protein structure prediction model on
this dataset. We validate the utility of this dataset for training by
participating CAMEO contest in which our model won the first place. We hope our
PSP dataset together with the training benchmark can enable a broader community
of AI/biology researchers for AI-driven protein related research.


  Microorganisms self-organize in very large communities exhibiting complex
fluctuations. Despite recent advances, still the mechanism by which these
systems are able to exhibit large variability at the one hand and dynamical
robustness on the other, is not fully explained. With that motivation, here we
analyze three aspects of the dynamics of the microbiota and plankton: the
density fluctuations, the correlation structure and the avalanching dynamics.
In all communities under study we find that the results exhibits scale-free
density fluctuations, anomalous variance' scaling, scale-free abundance
correlations and stationary scale-free avalanching dynamics. These behaviors,
typical in systems exhibiting critical dynamics, suggest criticality as a
potential mechanism to explain both the robustness and (paradoxical) high
irregularity of processes observed in very large microbial communities.


  Scientific and technological advances in medicine and systems biology have
unequivocally shown that health and disease must be viewed in the context of
the interplay among multiple molecular and environmental factors. Understanding
the effects of cellular interconnection on disease progression may lead to the
identification of novel disease genes and pathways, and hence influence
precision diagnostics and therapeutics. To accomplish this goal, the emerging
field of network medicine applies network science approaches to investigate
disease pathogenesis, integrating information from relevant Omics databases,
including protein-protein interaction, correlation-based, gene regulatory, and
Bayesian networks. However, this requires analysing and computing large amounts
of data. Moreover, if we are to efficiently search for new drugs and new drug
combinations, there is a pressing need for computational methods that could
allow us to access the immense chemical compound space until now largely
unexplored. Finally, at the microscopic level, drug-target chemistry simulation
is ultimately a quantum problem, and hence it requires a quantum solution. As
we will discuss, quantum computing may be a key ingredient in enabling the full
potential of network medicine. We propose to combine network medicine and
quantum algorithms in a novel research field, quantum network medicine, to lay
the foundations of a new era of disease prevention and drug design.


  Molecular optimization is a fundamental goal in the chemical sciences and is
of central interest to drug and material design. In recent years, significant
progress has been made in solving challenging problems across various aspects
of computational molecular optimizations, emphasizing high validity, diversity,
and, most recently, synthesizability. Despite this progress, many papers report
results on trivial or self-designed tasks, bringing additional challenges to
directly assessing the performance of new methods. Moreover, the sample
efficiency of the optimization--the number of molecules evaluated by the
oracle--is rarely discussed, despite being an essential consideration for
realistic discovery applications.
  To fill this gap, we have created an open-source benchmark for practical
molecular optimization, PMO, to facilitate the transparent and reproducible
evaluation of algorithmic advances in molecular optimization. This paper
thoroughly investigates the performance of 25 molecular design algorithms on 23
tasks with a particular focus on sample efficiency. Our results show that most
"state-of-the-art" methods fail to outperform their predecessors under a
limited oracle budget allowing 10K queries and that no existing algorithm can
efficiently solve certain molecular optimization problems in this setting. We
analyze the influence of the optimization algorithm choices, molecular assembly
strategies, and oracle landscapes on the optimization performance to inform
future algorithm development and benchmarking. PMO provides a standardized
experimental setup to comprehensively evaluate and compare new molecule
optimization methods with existing ones. All code can be found at
https://github.com/wenhao-gao/mol_opt.


  We propose a model of the evolution of a matrix along a phylogenetic tree, in
which transformations affect either entire rows or columns of the matrix. This
represents the change of both lexical and phonological aspects of linguistic
data, by allowing for new words to appear and for systematic phonological
changes to affect the entire vocabulary. We implement a Sequential Monte Carlo
method to sample from the posterior distribution, and infer jointly the
phylogeny, model parameters, and latent variables representing cognate births
and phonological transformations. We successfully apply this method to
synthetic and real data of moderate size.


  Posttranslational modifications (PTMs) are an integral component to how cells
respond to perturbation. While experimental advances have enabled improved PTM
identification capabilities, the same throughput for characterizing how
structural changes caused by PTMs equate to altered physiological function has
not been maintained. In this Perspective, we cover the history of computational
modeling and molecular dynamics simulations which have characterized the
structural implications of PTMs. We distinguish results from different
molecular dynamics studies based upon the timescales simulated and analysis
approaches used for PTM characterization. Lastly, we offer insights into how
opportunities for modern research efforts on in silico PTM characterization may
proceed given current state-of-the-art computing capabilities and
methodological advancements.


  When investigating epidemic dynamics through differential models, the
parameters needed to understand the phenomenon and to simulate forecast
scenarios require a delicate calibration phase, often made even more
challenging by the scarcity and uncertainty of the observed data reported by
official sources. In this context, Physics-Informed Neural Networks (PINNs), by
embedding the knowledge of the differential model that governs the physical
phenomenon in the learning process, can effectively address the inverse and
forward problem of data-driven learning and solving the corresponding epidemic
problem. In many circumstances, however, the spatial propagation of an
infectious disease is characterized by movements of individuals at different
scales governed by multiscale PDEs. This reflects the heterogeneity of a region
or territory in relation to the dynamics within cities and in neighboring
zones. In presence of multiple scales, a direct application of PINNs generally
leads to poor results due to the multiscale nature of the differential model in
the loss function of the neural network. To allow the neural network to operate
uniformly with respect to the small scales, it is desirable that the neural
network satisfies an Asymptotic-Preservation (AP) property in the learning
process. To this end, we consider a new class of AP Neural Networks (APNNs) for
multiscale hyperbolic transport models of epidemic spread that, thanks to an
appropriate AP formulation of the loss function, is capable to work uniformly
at the different scales of the system. A series of numerical tests for
different epidemic scenarios confirms the validity of the proposed approach,
highlighting the importance of the AP property in the neural network when
dealing with multiscale problems especially in presence of sparse and partially
observed systems.


  We improve a recent mathematical model for cholera by adding a time delay
that represents the time between the instant at which an individual becomes
infected and the instant at which he begins to have symptoms of cholera
disease. We prove that the delayed cholera model is biologically meaningful and
analyze the local asymptotic stability of the equilibrium points for positive
time delays. An optimal control problem is proposed and analyzed, where the
goal is to obtain optimal treatment strategies, through quarantine, that
minimize the number of infective individuals and the bacterial concentration,
as well as treatment costs. Necessary optimality conditions are applied to the
delayed optimal control problem, with a $L^1$ type cost functional. We show
that the delayed cholera model fits better the cholera outbreak that occurred
in the Department of Artibonite -- Haiti, from 1 November 2010 to 1 May 2011,
than the non-delayed model. Considering the data of the cholera outbreak in
Haiti, we solve numerically the delayed optimal control problem and propose
solutions for the outbreak control and eradication.


  This work provides an overview on deterministic and stochastic models that
have previously been proposed by us to study the transmission dynamics of the
Coronavirus Disease 2019 (COVID-19) in Europe and USA. Briefly, we describe
realistic deterministic and stochastic models for the evolution of the COVID-19
pandemic, subject to the lockdown and quarantine measures, which take into
account the time-delay for recovery or death processes. Realistic dynamic
equations for the entire process have been derived by adopting the so-called
"kinetic-type reactions approach". The lockdown and the quarantine measures are
modelled by some kind of inhibitor reactions where susceptible and infected
individuals can be "trapped" into inactive states. The dynamics for the
recovered people is obtained by accounting people who are only traced back to
hospitalised infected people. To model the role of the Hospitals we take
inspiration from the Michaelis-Menten's enzyme-substrate reaction model (the
so-called "MM reaction") where the "enzyme" is associated to the "available
hospital beds", the "substrate" to the "infected people", and the "product" to
the "recovered people", respectively. The statistical properties of the models,
in particular the relevant correlation functions and the probability density
functions, have duly been evaluated. We validate our theoretical predictions
with a large series of experimental data for Italy, Germany, France, Belgium
and United States, and we also compare data for Italy and Belgium with the
theoretical predictions of the logistic model. We found that our predictions
are in good agreement with the real world since the onset of COVID 19, contrary
to the the logistics model that only applies in the first days of the pandemic.
In the final part of the work, we can find the (theoretical) relationships that
should be satisfied to obtain the disappearance of the virus.


  Measuring gene expression simultaneously in both hosts and symbionts offers a
powerful approach to explore the biology underlying species interactions. Such
dual or simultaneous RNAseq approaches have primarily been used to gain insight
into gene function in model systems, but there is opportunity to expand and
apply these tools in new ways to understand ecological and evolutionary
questions. By incorporating genetic diversity in both hosts and symbionts and
studying how gene expression is correlated between partner species, we can gain
new insight into host-symbiont coevolution and the ecology of species
interactions. In this perspective, we explore how these relatively new tools
could be applied to study such questions. We review the mechanisms that could
be generating patterns of cross-species gene coexpression, including indirect
genetic effects and selective filters, how these tools could be applied across
different biological and temporal scales, and outline other methodological
considerations and experiment possibilities.


  Predicting the adaptation of populations to a changing environment is crucial
to assess the impact of human activities on biodiversity. Many theoretical
studies have tackled this issue by modeling the evolution of quantitative
traits subject to stabilizing selection around an optimum phenotype, whose
value is shifted continuously through time. In this context, the population
fate results from the equilibrium distribution of the trait, relative to the
moving optimum. Such a distribution may vary with the shape of selection, the
system of reproduction, the number of loci, the mutation kernel or their
interactions. Here, we develop a methodology that provides quantitative
measures of population maladaptation and potential of survival directly from
the entire profile of the phenotypic distribution, without any a priori on its
shape.We investigate two different models of reproduction (asexual and
infinitesimal sexual models of inheritance), with general forms of selection.In
particular, we recover that fitness functions such that selection weakens away
from the optimum lead to evolutionary tipping points, with an abrupt collapse
of the population when the speed of environmental change is too high. Our
unified framework furthermore allows highlighting the underlying mechanisms
that lead to this phenomenon. More generally, it allows discussing similarities
and discrepancies between the two reproduction models, the latter being
ultimately explained by different constraints on the evolution of the
phenotypic variance. We demonstrate that the mean fitness in the population
crucially depends on the shape of the selection function in the sexual
infinitesimal model, in contrast with the asexual model.In the asexual model,
we also investigate the effect of the mutation kernel and we show that kernels
with higher kurtosis tend to reduce maladaptation and improve fitness,
especially in fast changing environments.


  In a family of random variables, Taylor's law or Taylor's power law
offluctuation scaling is a variance function that gives the variance
$\sigma^{2}>0$ of a random variable (rv) $X$ with expectation $\mu >0$ as a
powerof $\mu$: $\sigma ^{2}=A\mu ^{b}$ for finite real $A>0,\ b$ that are
thesame for all rvs in the family. Equivalently, TL holds when $\log
\sigma^{2}=a+b\log \mu ,\ a=\log A$, for all rvs in some set. Here we analyze
thepossible values of the TL exponent $b$ in five families of
infinitelydivisible two-parameter distributions and show how the values of $b$
dependon the parameters of these distributions. The five families
areTweedie-Bar-Lev-Enis, negative binomial, compound Poisson-geometric,compound
geometric-Poisson (or P\'{o}lya-Aeppli), and gamma distributions.These families
arise frequently in empirical data and population models, and they are limit
laws of Markov processes that we exhibit in each case.


  The machine learning (ML) and deep learning (DL) techniques are widely
recognized to be powerful tools for virtual drug screening. The recently
reported ML- or DL-based scoring functions have shown exciting performance in
predicting protein-ligand binding affinities with fruitful application
prospects. However, the differentiation between highly similar ligand
conformations, including the native binding pose (the global energy minimum
state), remains challenging which could greatly enhance the docking. In this
work, we propose a fully differentiable framework for ligand pose optimization
based on a hybrid scoring function (SF) combined with a multi-layer perceptron
(DeepRMSD) and the traditional AutoDock Vina SF. The DeepRMSD+Vina, which
combines (1) the root mean square deviation (RMSD) of the docking pose with
respect to the native pose and (2) the AutoDock Vina score, is fully
differentiable thus is capable of optimizing the ligand binding pose to the
energy-lowest conformation. Evaluated by the CASF-2016 docking power dataset,
the DeepRMSD+Vina reaches a success rate of 95.4%, which is by far the best
reported SF to date. Based on this SF, an end-to-end ligand pose optimization
framework was implemented to improve the docking pose quality. We demonstrated
that this method significantly improves the docking success rate (by 15%) in
redocking and crossdocking tasks, revealing the high potentialities of this
framework in drug design and discovery.


  The study of a receptor-ligand system generally relies on the analysis of its
dose-response (or concentration-effect) curve, which quantifies the relation
between ligand concentration and the biological effect (or cellular response)
induced when binding its specific cell surface receptor. Mathematical models of
receptor-ligand systems have been developed to compute a dose-response curve
under the assumption that the biological effect is proportional to the number
of ligand-bound receptors. Given a dose-response curve, two quantities (or
metrics) have been defined to characterise the properties of the
ligand-receptor system under consideration: amplitude and potency (or
half-maximal effective concentration, and denoted by EC$_{50}$). Both the
amplitude and the EC$_{50}$ are key quantities commonly used in
pharmaco-dynamic modelling, yet a comprehensive mathematical investigation of
the behaviour of these two metrics is still outstanding; for a large (and
important) family of receptors, called cytokine receptors, we still do not know
how amplitude and EC$_{50}$ depend on receptor copy numbers. Here we make use
of algebraic approaches (Gr\"obner basis) to study these metrics for a large
class of receptor-ligand models, with a focus on cytokine receptors. In
particular, we introduce a method, making use of two motivating examples based
on the interleukin-7 (IL-7) receptor, to compute analytic expressions for the
amplitude and the EC$_{50}$. We then extend the method to a wider class of
receptor-ligand systems, sequential receptor-ligand systems with extrinsic
kinase, and provide some examples.