  During its first two years, the SARS-CoV-2 pandemic manifested as multiple
waves shaped by complex interactions between variants of concern,
non-pharmaceutical interventions, and the immunological landscape of the
population. Understanding how the age-specific epidemiology of SARS-CoV-2 has
evolved throughout the pandemic is crucial for informing policy decisions. We
developed an inference-based modelling approach to reconstruct the burden of
true infections and hospital admissions in children, adolescents and adults
over the seven waves of four variants (wild-type, Alpha, Delta, Omicron BA.1)
during the first two years of the pandemic, using the Netherlands as the
motivating example. We find that reported cases are a considerable
underestimate and a generally poor predictor of true infection burden,
especially because case reporting differs by age. The contribution of children
and adolescents to total infection and hospitalization burden increased with
successive variants and was largest during the Omicron BA.1 period. Before the
Delta period, almost all infections were primary infections occurring in naive
individuals. During the Delta and Omicron BA.1 periods, primary infections were
common in children but relatively rare in adults who experienced either
re-infections or breakthrough infections. Our approach can be used to
understand age-specific epidemiology through successive waves in other
countries where random community surveys uncovering true SARS-CoV-2 dynamics
are absent but basic surveillance and statistics data are available.


  Model reduction of fast-slow chemical reaction networks based on the
quasi-steady state approximation fails when the fast subsystem has first
integrals. We call these first integrals approximate conservation laws. In
order to define fast subsystems and identify approximate conservation laws, we
use ideas from tropical geometry. We prove that any approximate conservation
law evolves slower than all the species involved in it and therefore represents
a supplementary slow variable in an extended system. By elimination of some
variables of the extended system, we obtain networks without approximate
conservation laws, which can be reduced by standard singular perturbation
methods. The field of applications of approximate conservation laws covers the
quasi-equilibrium approximation, well known in biochemistry. We discuss
reductions of slow-fast as well as multiple timescale systems. Networks with
multiple timescales have hierarchical relaxation. At a given timescale, our
multiple timescale reduction method defines three subsystems composed of (i)
slaved fast variables satisfying algebraic equations, (ii) slow driving
variables satisfying reduced ordinary differential equations, and (iii)
quenched much slower variables that are constant. The algebraic equations
satisfied by fast variables define chains of nested normally hyberbolic
invariant manifolds. In such chains, faster manifolds are of higher dimension
and contain the slower manifolds. Our reduction methods are introduced
algorithmically for networks with monomial reaction rates and linear, monomial
or polynomial approximate conservation laws. We propose symbolic algorithms to
reshape and rescale the networks such that geometric singular perturbation
theory can be applied to them, test the applicability of the theory, and
finally reduce the networks. As a proof of concept, we apply this method to a
model of the TGF-b signaling pathway.


  Standardising the representation of biomedical knowledge among all
researchers is an insurmountable task, hindering the effectiveness of many
computational methods. To facilitate harmonisation and interoperability despite
this fundamental challenge, we propose to standardise the framework of
knowledge graph creation instead. We implement this standardisation in
BioCypher, a FAIR (findable, accessible, interoperable, reusable) framework to
transparently build biomedical knowledge graphs while preserving provenances of
the source data. Mapping the knowledge onto biomedical ontologies helps to
balance the needs for harmonisation, human and machine readability, and ease of
use and accessibility to non-specialist researchers. We demonstrate the
usefulness of this framework on a variety of use cases, from maintenance of
task-specific knowledge stores, to interoperability between biomedical domains,
to on-demand building of task-specific knowledge graphs for federated learning.
BioCypher (https://biocypher.org) frees up valuable developer time; we
encourage further development and usage by the community.


  Virtual screening (VS) is an essential technique for understanding
biomolecular interactions, particularly, drug design and discovery. The
best-performing VS models depend vitally on three-dimensional (3D) structures,
which are not available in general but can be obtained from molecular docking.
However, current docking accuracy is relatively low, rendering unreliable VS
models. We introduce sequence-based virtual screening (SVS) as a new generation
of VS models for modeling biomolecular interactions. The SVS model utilizes
advanced natural language processing (NLP) algorithms and optimizes deep
$K$-embedding strategies to encode biomolecular interactions without invoking
3D structure-based docking. We demonstrate the state-of-art performance of SVS
for four regression datasets involving protein-ligand binding, protein-protein,
protein-nucleic acid binding, and ligand inhibition of protein-protein
interactions and five classification datasets for the protein-protein
interactions in five biological species. SVS has the potential to dramatically
change the current practice in drug discovery and protein engineering.


  Population heterogeneity, especially in individuals' contact networks, plays
an important role in transmission dynamics of infectious diseases. For
vaccine-preventable diseases, outstanding issues like vaccine hesitancy and
availability of vaccines further lead to nonuniform coverage among groups, not
to mention the efficacy of vaccines and the mixing pattern varying from one
group to another. As the ongoing COVID-19 pandemic transitions to endemicity,
it is of interest and significance to understand the impact of aforementioned
population heterogeneity on the emergence and persistence of epidemics. Here we
analyze epidemic thresholds and characterize bifurcation dynamics by accounting
for heterogeneity caused by group-dependent characteristics, including
vaccination rate and efficacy as well as disease transmissibility. Our analysis
shows that increases in the difference in vaccination coverage among groups can
render multiple equilibria of disease burden to exist even if the overall basic
reproductive ratio is below one (also known as backward bifurcation). The
presence of other heterogeneity factors such as differences in vaccine
efficacy, transmission, mixing pattern, and group size can each exhibit subtle
impacts on bifurcation. We find that heterogeneity in vaccine efficacy can
undermine the condition for backward bifurcations whereas homophily tends to
aggravate disease endemicity. Our results have practical implications for
improving public health efforts by addressing the role of population
heterogeneity in the spread and control of diseases.


  The secondary structure of ribonucleic acid (RNA) is more stable and
accessible in the cell than its tertiary structure, making it essential for
functional prediction. Although deep learning has shown promising results in
this field, current methods suffer from poor generalization and high
complexity. In this work, we present RFold, a simple yet effective RNA
secondary structure prediction in an end-to-end manner. RFold introduces a
decoupled optimization process that decomposes the vanilla constraint
satisfaction problem into row-wise and column-wise optimization, simplifying
the solving process while guaranteeing the validity of the output. Moreover,
RFold adopts attention maps as informative representations instead of designing
hand-crafted features. Extensive experiments demonstrate that RFold achieves
competitive performance and about eight times faster inference efficiency than
the state-of-the-art method. The code and Colab demo are available in
\href{http://github.com/A4Bio/RFold}{http://github.com/A4Bio/RFold}.


  Non-sharable sensitive data collection and analysis in large-scale consortia
for genomic research is complicated. Time consuming issues in installing
software arise due to different operating systems, software dependencies and
running the software. Therefore, easier, more standardized, automated protocols
and platforms can be a solution to overcome these issues. We have developed one
such solution for genomic data analysis using software container technologies.
The platform, COGEDAP, consists of different software tools placed into
Singularity containers with corresponding pipelines and instructions on how to
perform genome-wide association studies (GWAS) and other genomic data analysis
via corresponding tools. Using a provided helper script written in Python,
users can obtain auto-generated scripts to conduct the desired analysis both on
high-performance computing (HPC) systems and on personal computers. The
analyses can be done by running these auto-generated scripts with the software
containers. The helper script also performs minor re-formatting of the
input/output data, so that the end user can work with a unified file format
regardless of which genetic software is used for the analysis. COGEDAP is
actively being used by users from different countries/projects to conduct their
genomic data analyses. Thanks to this platform, users can easily run GWAS and
other genomic analyses without spending much effort on software installation,
data formats, and other technical requirements.


  Rapid advancements in collection and dissemination of multi-platform
molecular and genomics data has resulted in enormous opportunities to aggregate
such data in order to understand, prevent, and treat human diseases. While
significant improvements have been made in multi-omic data integration methods
to discover biological markers and mechanisms underlying both prognosis and
treatment, the precise cellular functions governing these complex mechanisms
still need detailed and data-driven de-novo evaluations. We propose a framework
called Functional Integrative Bayesian Analysis of High-dimensional
Multiplatform Genomic Data (fiBAG), that allows simultaneous identification of
upstream functional evidence of proteogenomic biomarkers and the incorporation
of such knowledge in Bayesian variable selection models to improve signal
detection. fiBAG employs a conflation of Gaussian process models to quantify
(possibly non-linear) functional evidence via Bayes factors, which are then
mapped to a novel calibrated spike-and-slab prior, thus guiding selection and
providing functional relevance to the associations with patient outcomes. Using
simulations, we illustrate how integrative methods with functional calibration
have higher power to detect disease related markers than non-integrative
approaches. We demonstrate the profitability of fiBAG via a pan-cancer analysis
of 14 cancer types to identify and assess the cellular mechanisms of
proteogenomic markers associated with cancer stemness and patient survival.


  Autocatalysis is thought to have played an important role in the earliest
stages of the origin of life. An autocatalytic cycle's (AC) constituent
chemicals can collectively catalyze their own recreation. When the reactions of
multiple, interacting ACs are active in a region of space, they form an
autocatalytic chemical ecosystem (ACE). Previous work demonstrated that, in
chemostats, interactions between ACs in ACEs can be framed as analogous to
those between species in biological ecosystems. Here, we extend this framework
to investigate the effects of surface adsorption, desorption, and diffusion on
ACE ecology. Simulating ACEs as particle-based stochastic reaction-diffusion
systems in spatial settings, including open, two-dimensional reaction-diffusion
systems and adsorptive mineral surfaces, we demonstrate that spatial structure
can support more coexisting ACs and expose new AC traits to selection.


  This paper presents a dual-aptamer scheme to cancel the signal drifts from
structure-switching aptamers during long-term monitoring. Electrochemical
aptamer-based (E-AB) biosensors recently demonstrated their great potential for
in vivo continuous monitoring. Nevertheless, the detection accuracy is often
limited by the signaling drifts. Conventionally, these drifts are removed by
the kinetic differential measurements (KDM) when coupled with square-wave
voltammetry. Yet we discover that KDM does not apply to every aptamer as the
responses at different SWV frequencies heavily depend on its
structure-switching characteristics and the redox reporters' electron transfer
(ET) kinetics. To this end, we present a "dual-aptamer" scheme that uses two
aptamers responding differentially to the same molecular target for drift
cancellation. We identify these paired aptamers through (1) screening from the
existing aptamers pool and (2) engineering the signaling behavior of the redox
reporters. We demonstrate their differential signaling to ampicillin and ATP
molecules and show that the aptamer pair bears common drifts in undilute goat
serum. Through cancellation, sensor drift is reduced by 370-fold. Benefiting
from the "differential" signaling, the recording throughput is also doubled
using differential readout electronics. The authors believe the proposed
technique is beneficial for long-term in vivo monitoring.


  We consider a model of antiviral immune response proposed in the works of
G.I. Marchuk. The model is described by a system of delay differential
equations. The asymptotic stability of a stationary solution to the system
corresponding to a completely healthy organism is studied. Estimates of the
attraction set of the given stationary solution are obtained and estimates of
solutions characterizing the stabilization rate at infinity are established.
The results are obtained using the Lyapunov-Krasovskii functional.


  We derive asymptotic formulae in the limit when population size N tends to
infinity for mean fixation times (conditional and unconditional) in a
population with two types of individuals, A and B, governed by the Moran
process. We consider only the case in which the fitness of the two types do not
depend on the population frequencies. Our results start with the important
cases in which the initial condition is a single individual of any type, but we
also consider the initial condition of a fraction x, 0<x<1, of A individuals,
where x is kept fixed and the total population size tends to infinity. In the
cases covered by Antal and Scheuring (Bull Math Biol 68(8):1923-1944, 2006),
i.e. conditional fixation times for a single individual of any type, it will
turn out that our formulae are much more accurate than the ones they found. As
quoted, our results include other situations not treated by them.


  For polynomial ODE models, we introduce and discuss the concepts of exact and
approximate conservation laws, which are the first integrals of the full and
truncated sets of ODEs. For fast-slow systems, truncated ODEs describe the fast
dynamics. We define compatibility classes as subsets of the state space,
obtained by equating the conservation laws to constants. A set of conservation
laws is complete when the corresponding compatibility classes contain a finite
number of steady states. Complete sets of conservation laws can be used for
model order reduction and for studying the multistationarity of the model. We
provide algorithmic methods for computing linear, monomial, and polynomial
conservation laws of polynomial ODE models and for testing their completeness.
The resulting conservation laws and their completeness are either independent
or dependent on the parameters. In the latter case, we provide parametric case
distinctions. In particular, we propose a new method to compute polynomial
conservation laws by comprehensive Gr\"obner systems and syzygies.
  Keywords: First integrals, chemical reaction networks, polynomial
conservation laws, syzygies, comprehensive Gr\"obner systems.


  As shown recently by the author, constant population SI(R)S models map to
Hethcote's classic endemic model originally proposed in 1973. This unifies a
whole class of models with up to 10 parameters, all being isomorphic to a
simple 2-parameter master model for endemic bifurcation. In this work this
procedure is extended to a 14-parameter SSISS Model, including social behavior
parameters, a (diminished) susceptibility of the R-compartment and unbalanced
constant per capita birth and death rates, thus covering many prominent models
in the literature. Under mild conditions, in the dynamics for fractional
variables in this model all vital parameters become redundant at the cost of
possibly negative incidence rates. There is a symmetry group G acting on
parameter space A, such that systems with G-equivalent parameters are
isomorphic and map to the same normalized system. Using (Xrep,I) as canonical
coordinates, Xrep the replacement number, normalization reduces to parameter
space A/G with 5 parameters only. This approach reveals unexpected relations
between various models in the literature. Part two of this work will analyze
equilibria, stability and backward bifurcation and part three will further
reduce the number of essential parameters from 5 to 3.


  Learning the underlying distribution of molecular graphs and generating
high-fidelity samples is a fundamental research problem in drug discovery and
material science. However, accurately modeling distribution and rapidly
generating novel molecular graphs remain crucial and challenging goals. To
accomplish these goals, we propose a novel Conditional Diffusion model based on
discrete Graph Structures (CDGS) for molecular graph generation. Specifically,
we construct a forward graph diffusion process on both graph structures and
inherent features through stochastic differential equations (SDE) and derive
discrete graph structures as the condition for reverse generative processes. We
present a specialized hybrid graph noise prediction model that extracts the
global context and the local node-edge dependency from intermediate graph
states. We further utilize ordinary differential equation (ODE) solvers for
efficient graph sampling, based on the semi-linear structure of the probability
flow ODE. Experiments on diverse datasets validate the effectiveness of our
framework. Particularly, the proposed method still generates high-quality
molecular graphs in a limited number of steps. Our code is provided in
https://github.com/GRAPH-0/CDGS.


  Phylogenetic trees (i.e. evolutionary trees, additive trees or X-trees) play
a key role in the processes of modeling and representing species evolution.
Genome evolution of a given group of species is usually modeled by a species
phylogenetic tree that represents the main patterns of vertical descent.
However, the evolution of each gene is unique. It can be represented by its own
gene tree which can differ substantially from a general species tree
representation. Consensus trees and supertrees have been widely used in
evolutionary studies to combine phylogenetic information contained in
individual gene trees. Nevertheless, if the available gene trees are quite
different from each other, then the resulting consensus tree or supertree can
either include many unresolved subtrees corresponding to internal nodes of high
degree or can simply be a star tree. This may happen if the available gene
trees have been affected by different reticulate evolutionary events, such as
horizontal gene transfer, hybridization or genetic recombination. Thus, the
problem of inferring multiple alternative consensus trees or supertrees, using
clustering, becomes relevant since it allows one to regroup in different
clusters gene trees having similar evolutionary patterns (e.g. gene trees
representing genes that have undergone the same horizontal gene transfer or
recombination events). We critically review recent advances and methods in the
field of phylogenetic tree clustering, discuss the methods' mathematical
properties, and describe the main advantages and limitations of multiple
consensus tree and supertree approaches. In the application section, we show
how the multiple supertree clustering approach can be used to cluster aaRS gene
trees according to their evolutionary patterns.


  Nitrogen (N) deposition is a primary driver of species loss in plant
communities globally. However, the mechanisms by which high N availability
causes species loss remain unclear. Many hypotheses for species loss with
increasing N availability highlight density-dependent mechanisms, i.e., changes
in species interactions. However, an alternative set of hypotheses highlights
density-independent detrimental effects of nitrogen (e.g., N toxicity). We
tested the role of density-dependent and density-independent mechanisms in
reducing species performance. For this aim, we used 120 experimental plant
communities comprised of annual species growing together in containers under
four fertilization treatments: (1) no nutrient addition(, (2) all nutrients
except N (P, K, and micronutrients), (3) Low N, and (4) high N. Each
fertilization treatment included two sowing densities to differentiate between
the effects of competition (N * density interactions) and other detrimental
effects of N. We focused on three performance attributes: the probability of
reaching the reproduction period, biomass growth, and population growth. We
found that individual biomass and population growth rates decreased with
increasing sowing density in all nutrient treatments, implying that species
interactions were predominantly negative. The common grass had a higher biomass
and population growth under N enrichment, regardless of sowing density. In
contrast, the legume showed a density-independent reduction in biomass growth
with increasing N. Lastly, the small forb showed a density-dependent reduction
in population growth, i.e., the decline occurred only under high density. Our
results demonstrate that density-dependent and density-independent mechanisms
operate simultaneously to reduce species performance under high N availability.
Yet, their relative importance varies among species and life stages.


  In a diverse population, where many species are present, competitors can
fight for surviving at individual and collective levels. In particular,
species, which would beat each other individually, may form a specific alliance
that ensures them stable coexistence against the invasion of an external
species. Our principal goal is to identify those general features of a
formation which determine its vitality. Therefore, we here study a traditional
Lotka-Volterra model of eight-species where two four-species cycles can fight
for space. Beside these formations, there are other solutions which may emerge
when invasion rates are varied. The complete range of parameters is explored
and we find that in most of the cases those alliances prevail which are formed
by equally strong members. Interestingly, there are regions where the symmetry
is broken and the system is dominated by a solution formed by seven species.
Our work also highlights that serious finite-size effects may emerge which
prevent observing the valid solution in a small system.


  Missing values are a common problem in data science and machine learning.
Removing instances with missing values can adversely affect the quality of
further data analysis. This is exacerbated when there are relatively many more
features than instances, and thus the proportion of affected instances is high.
Such a scenario is common in many important domains, for example, single
nucleotide polymorphism (SNP) datasets provide a large number of features over
a genome for a relatively small number of individuals. To preserve as much
information as possible prior to modeling, a rigorous imputation scheme is
acutely needed. While Denoising Autoencoders is a state-of-the-art method for
imputation in high-dimensional data, they still require enough complete cases
to be trained on which is often not available in real-world problems. In this
paper, we consider missing value imputation as a multi-label classification
problem and propose Chains of Autoreplicative Random Forests. Using multi-label
Random Forests instead of neural networks works well for low-sampled data as
there are fewer parameters to optimize. Experiments on several SNP datasets
show that our algorithm effectively imputes missing values based only on
information from the dataset and exhibits better performance than standard
algorithms that do not require any additional information. In this paper, the
algorithm is implemented specifically for SNP data, but it can easily be
adapted for other cases of missing value imputation.


  First passage time (FPT) is the time a particle, subject to some stochastic
process, hits or crosses a closed surface for the very first time.
$\tau$-leaping methods are a class of stochastic algorithms in which, instead
of simulating every single reaction, many reactions are ``leaped" over in order
to shorten the computing time. In this paper we developed a $\tau$-leaping
method for computing the FPT and position in arbitrary volumes for a Brownian
particle governed by the Langevin equation. The $\tau$-leaping method proposed
here works as follows. A sphere is inscribed within the volume of interest
(VOI) centered at the initial particle's location. On this sphere, the FPT is
sampled, as well as the position, which becomes the new initial position. Then,
another sphere, centered at this new location, is inscribed. This process
continues until the sphere becomes smaller than some minimal radius
$R_{\text{min}}$. When this occurs, the $\tau$-leaping switches to the
conventional Monte Carlo, which runs until the particle either crosses the
surface of the VOI or finds its way to a position where a sphere of radius
$>R_{\text{min}}$ can be inscribed. The switching between $\tau$-leaping and MC
continues until the particle crosses the surface of the VOI. The size of this
radius depends on the system parameters and on one's notion of accuracy: the
larger this radius the more accurate the $\tau$-leaping method, but also less
efficient. This trade off between accuracy and efficiency is discussed. For two
VOI, the $\tau$-leaping method is shown to be accurate and more efficient than
MC by at least a factor of 10 and up to a factor of about 110. However, while
MC becomes exponentially slower with increasing VOI, the efficiency of the
$\tau$-leaping method remains relatively unchanged. Thus, the $\tau$-leaping
method can potentially be many orders of magnitude more efficient than MC.


  The classical two-species non-linear Predator-Prey system, often used in
population dynamics modeling, is expressed in terms of a single positive
coupling parameter $\lambda$. Based on standard logarithmic transformations, we
derive a novel $\lambda$-\textit{invariant} Hamiltonian resulting in two
coupled first-order ODEs for ``hybrid-species'', \textit{albeit} with one being
\textit{linear}; we thus derive a new exact, closed-form, single quadrature
solution valid for any value of $\lambda$ and the system's energy. In the
particular case $\lambda = 1$ the ODE system completely uncouples and a new,
exact, energy-only dependent simple quadrature solution is derived. In the case
$\lambda \neq 1$ an accurate practical approximation uncoupling the non-linear
system is proposed and solutions are provided in terms of explicit quadratures
together with high energy asymptotic solutions. A novel, exact, closed-form
expression of the system's oscillation period valid for any value of $\lambda$
and orbital energy is also derived; two fundamental properties of the period
are established; for $\lambda = 1$ the period is expressed in terms of a
universal energy function and shown to be the shortest.


  Outbreaks are complex multi-scale processes that are impacted not only by
cellular dynamics and the ability of pathogens to effectively reproduce and
spread, but also by population-level dynamics and the effectiveness of
mitigation measures. A timely exchange of information related to the spread of
novel pathogens, stay-at-home orders, and other containment measures can be
effective at containing an infectious disease, particularly during in the early
stages when testing infrastructure, vaccines, and other medical interventions
may not be available at scale. Using a multiplex epidemic model that consists
of an information layer (modeling information exchange between individuals) and
a spatially embedded epidemic layer (representing a human contact network), we
study how random and targeted disruptions in the information layer (\eg, errors
and intentional attacks on communication infrastructure) impact outbreak
dynamics. We calibrate our model to the early outbreak stages of the SARS-CoV-2
pandemic in 2020. Mitigation campaign can still be effective under random
disruptions, such as failure of information channels between a few individuals.
However, targeted disruptions or sabotage of hub nodes that exchange
information with a large number of individuals can abruptly change outbreak
characteristics such as the time to reach the peak infection. Our results
emphasize the importance of using a robust communication infrastructure that
can withstand both random and targeted disruptions.


  Cataloging the complex behaviors of dynamical systems can be challenging,
even when they are well-described by a simple mechanistic model. If such a
system is of limited analytical tractability, brute force simulation is often
the only resort. We present an alternative, optimization-driven approach using
tools from machine learning. We apply this approach to a novel,
fully-optimizable, reaction-diffusion model which incorporates complex chemical
reaction networks (termed "Dense Reaction-Diffusion Network" or "Dense RDN").
This allows us to systematically identify new states and behaviors, including
pattern formation, dissipation-maximizing nonequilibrium states, and
replication-like dynamical structures.


  The reconstruction of phylogenetic networks is an important but challenging
problem in phylogenetics and genome evolution, as the space of phylogenetic
networks is vast and cannot be sampled well. One approach to the problem is to
solve the minimum phylogenetic network problem, in which phylogenetic trees are
first inferred, then the smallest phylogenetic network that displays all the
trees is computed. The approach takes advantage of the fact that the theory of
phylogenetic trees is mature and there are excellent tools available for
inferring phylogenetic trees from a large number of biomolecular sequences. A
tree-child network is a phylogenetic network satisfying the condition that
every non-leaf node has at least one child that is of indegree one. Here, we
develop a new method that infers the minimum tree-child network by aligning
lineage taxon strings in the phylogenetic trees. This algorithmic innovation
enables us to get around the limitations of the existing programs for
phylogenetic network inference. Our new program, named ALTS, is fast enough to
infer a tree-child network with a large number of reticulations for a set of up
to 50 phylogenetic trees with 50 taxa that have only trivial common clusters in
about a quarter of an hour on average.


  Biological systems and processes are networks of complex nonlinear regulatory
interactions between nucleic acids, proteins, and metabolites. A natural way in
which to represent these interaction networks is through the use of a graph. In
this formulation, each node represents a nucleic acid, protein, or metabolite
and edges represent intermolecular interactions (inhibition, regulation,
promotion, coexpression, etc.). In this work, a novel algorithm for the
discovery of latent graph structures given experimental data is presented.


  Here, we demonstrate how machine learning enables the prediction of
comonomers reactivity ratios based on the molecular structure of monomers. We
combined multi-task learning, multi-inputs, and Graph Attention Network to
build a model capable of predicting reactivity ratios based on the monomers
chemical structures.


  A common type of cytoskeletal morphology involves multiple converging
microbutubules with their minus ends collected and stabilized by a microtubule
organizing center (MTOC) in the interior of the cell. This arrangement enables
the ballistic transport of cargo bound to microtubules, both dynein mediated
transport towards the MTOC and kinesin mediated transport away from it,
interspersed with diffusion for unbound cargo-motor complexes. Spatial and
temporal positioning of the MTOC allows for bidirectional transport towards and
away from specific organelles and locations within the cell and also the
sequestering and subsequent dispersal of dynein transported cargo. The general
principles governing dynamics, efficiency and tunability of such transport in
the MTOC vicinity is not fully understood. To address this, we develop a
one-dimensional model that includes advective transport towards an attractor
(such as the MTOC), and diffusive transport that allows particles to reach
absorbing boundaries (such as cellular membranes). We calculated the mean first
passage time (MFPT) for cargo to reach the boundaries as a measure of the
effectiveness of sequestering (large MFPT) and diffusive dispersal (low MFPT).
The MFPT experiences a dramatic growth in magnitude, transitioning from a low
to high MFPT regime (dispersal to sequestering) over a window of cargo
attachment/detachment rates that is close to in vivo values. We find that
increasing either the attachment or detachment rate, while fixing the other,
can result in optimal dispersal when the attractor is placed asymmetrically.
Finally, we describe a rare event regime, where the escape location is
exponentially sensitive to the attractor positioning. Our results suggest that
structures such as the MTOC allow for the sensitive control of the spatial and
temporal features of transport and corresponding function under physiological
conditions.


  Regulatory networks depict promoting or inhibiting interactions between
molecules in a biochemical system. We introduce a category-theoretic formalism
for regulatory networks, using signed graphs to model the networks and signed
functors to describe occurrences of one network in another, especially
occurrences of network motifs. With this foundation, we establish functorial
mappings between regulatory networks and other mathematical models in
biochemistry. We construct a functor from reaction networks, modeled as Petri
nets with signed links, to regulatory networks, enabling us to precisely define
when a reaction network could be a physical mechanism underlying a regulatory
network. Turning to quantitative models, we associate a regulatory network with
a Lotka-Volterra system of differential equations, defining a functor from the
category of signed graphs to a category of parameterized dynamical systems. We
extend this result from closed to open systems, demonstrating that
Lotka-Volterra dynamics respects not only inclusions and collapsings of
regulatory networks, but also the process of building up complex regulatory
networks by gluing together simpler pieces. Formally, we use the theory of
structured cospans to produce a lax double functor from the double category of
open signed graphs to that of open parameterized dynamical systems. Throughout
the paper, we ground the categorical formalism in examples inspired by systems
biology.


  We run stochastic simulations of the spatial version of the
rock-paper-scissors game, considering that individuals use sensory abilities to
scan the environment to detect the presence of enemies. If the local
dangerousness level is above a tolerable threshold, individuals aggregate
instead of moving randomly on the lattice. We study the impact of the locally
adaptive aggregation on the organisms' spatial organisation by measuring the
characteristic length scale of the spatial domains occupied by organisms of a
single species. Our results reveal that aggregation is beneficial if triggered
when the local density of opponents does not exceed $30\%$; otherwise, the
behavioural strategy may harm individuals by increasing the average death risk.
We show that if organisms can perceive further distances, they can accurately
scan and interpret the signals from the neighbourhood, maximising the effects
of the locally adaptive aggregation on the death risk. Finally, we show that
the locally adaptive aggregation behaviour promotes biodiversity independently
of the organism's mobility. The coexistence probability rises if organisms join
conspecifics, even in the presence of a small number of enemies. We verify that
our conclusions hold for more complex systems by simulating the generalised
rock-paper-scissors models with five and seven species. Our discoveries may be
helpful to ecologists in understanding systems where organisms' self-defence
behaviour adapts to local environmental cues.


  There are various sources of ionizing radiation exposure, where medical
exposure for radiation therapy or diagnosis is the most common human-made
source. Understanding how gene expression is modulated after ionizing radiation
exposure and investigating the presence of any dose-dependent gene expression
patterns have broad implications for health risks from radiotherapy, medical
radiation diagnostic procedures, as well as other environmental exposure. In
this paper, we perform a comprehensive pathway-based analysis of gene
expression profiles in response to low-dose radiation exposure, in order to
examine the potential mechanism of gene regulation underlying such responses.
To accomplish this goal, we employ a statistical framework to determine whether
a specific group of genes belonging to a known pathway display coordinated
expression patterns that are modulated in a manner consistent with the
radiation level. Findings in our study suggest that there exist complex yet
consistent signatures that reflect the molecular response to radiation
exposure, which differ between low-dose and high-dose radiation.


  Effective representation of molecules is a crucial factor affecting the
performance of artificial intelligence models. This study introduces a
flexible, fragment-based, multiscale molecular representation framework called
t-SMILES (tree-based SMILES) with three code algorithms: TSSA (t-SMILES with
Shared Atom), TSDY (t-SMILES with Dummy Atom) and TSID (t-SMILES with ID). It
describes molecules using SMILES-type strings obtained by performing a
breadth-first search on a full binary tree formed from a fragmented molecular
graph. Systematic evaluations using JTVAE, BRICS, MMPA, and Scaffold show the
feasibility to construct a multi-code molecular description system, where
various descriptions complement each other, enhancing the overall performance.
Additionally, it exhibits impressive performance on low-resource datasets,
whether the model is original, data augmented, or pre-training fine-tuned. It
significantly outperforms classical SMILES, DeepSMILES, SELFIES and baseline
models in goal-directed tasks. Furthermore, it surpasses start-of-the-art
fragment, graph and SMILES based approaches on ChEMBL, Zinc, and QM9.


  Recent studies of genotype-phenotype (GP) maps have reported universally
enhanced phenotypic robustness to genotype mutations, a feature essential to
evolution. Virtually all of these studies make a simplifying assumption that
each genotype maps deterministically to a single phenotype. Here, we introduce
probabilistic genotype-phenotype (PrGP) maps, where each genotype maps to a
vector of phenotype probabilities, as a more realistic framework for
investigating robustness. We study three model systems to show that our
generalized framework can handle uncertainty emerging from various physical
sources: (1) thermal fluctuation in RNA folding, (2) external field disorder in
spin glass ground state finding, and (3) superposition and entanglement in
quantum circuits, which are realized experimentally on a 7-qubit IBM quantum
computer. In all three cases, we observe a novel biphasic robustness scaling
which is enhanced relative to random expectation for more frequent phenotypes
and approaches random expectation for less frequent phenotypes.


  The diversity of terrestrial vascular plants plays a key role in maintaining
the stability and productivity of ecosystems. Monitoring species compositional
diversity across large spatial scales is challenging and time consuming. The
advanced spectral and spatial specification of the recently launched DESIS (the
DLR Earth Sensing Imaging Spectrometer) instrument provides a unique
opportunity to test the potential for monitoring plant species diversity with
spaceborne hyperspectral data. This study provides a quantitative assessment on
the ability of DESIS hyperspectral data for predicting plant species richness
in two different habitat types in southeast Australia. Spectral features were
first extracted from the DESIS spectra, then regressed against on-ground
estimates of plant species richness, with a two-fold cross validation scheme to
assess the predictive performance. We tested and compared the effectiveness of
Principal Component Analysis (PCA), Canonical Correlation Analysis (CCA), and
Partial Least Squares analysis (PLS) for feature extraction, and Kernel Ridge
Regression (KRR), Gaussian Process Regression (GPR), Random Forest Regression
(RFR) for species richness prediction. The best prediction results were r=0.76
and RMSE=5.89 for the Southern Tablelands region, and r=0.68 and RMSE=5.95 for
the Snowy Mountains region. Relative importance analysis for the DESIS spectral
bands showed that the red-edge, red, and blue spectral regions were more
important for predicting plant species richness than the green bands and the
near-infrared bands beyond red-edge. We also found that the DESIS hyperspectral
data performed better than Sentinel-2 multispectral data in the prediction of
plant species richness. Our results provide a quantitative reference for future
studies exploring the potential of spaceborne hyperspectral data for plant
biodiversity mapping.


  An automata network with $n$ components over a finite alphabet $Q$ of size
$q$ is a discrete dynamical system described by the successive iterations of a
function $f:Q^n\to Q^n$. In most applications, the main parameter is the
interaction graph of $f$: the digraph with vertex set $[n]$ that contains an
arc from $j$ to $i$ if $f_i$ depends on input $j$. What can be said on the set
$\mathbb{G}(f)$ of the interaction graphs of the automata networks isomorphic
to $f$? It seems that this simple question has never been studied. Here, we
report some basic facts. First, we prove that if $n\geq 5$ or $q\geq 3$ and $f$
is neither the identity nor constant, then $\mathbb{G}(f)$ always contains the
complete digraph $K_n$, with $n^2$ arcs. Then, we prove that $\mathbb{G}(f)$
always contains a digraph whose minimum in-degree is bounded as a function of
$q$. Hence, if $n$ is large with respect to $q$, then $\mathbb{G}(f)$ cannot
only contain $K_n$. However, we prove that $\mathbb{G}(f)$ can contain only
dense digraphs, with at least $\lfloor n^2/4 \rfloor$ arcs.


  Machine Learning-guided solutions for protein learning tasks have made
significant headway in recent years. However, success in scientific discovery
tasks is limited by the accessibility of well-defined and labeled in-domain
data. To tackle the low-data constraint, recent adaptions of deep learning
models pretrained on millions of protein sequences have shown promise; however,
the construction of such domain-specific large-scale model is computationally
expensive. Here, we propose Representation Learning via Dictionary Learning
(R2DL), an end-to-end representation learning framework in which we reprogram
deep models for alternate-domain tasks that can perform well on protein
property prediction with significantly fewer training samples. R2DL reprograms
a pretrained English language model to learn the embeddings of protein
sequences, by learning a sparse linear mapping between English and protein
sequence vocabulary embeddings. Our model can attain better accuracy and
significantly improve the data efficiency by up to $10^5$ times over the
baselines set by pretrained and standard supervised methods. To this end, we
reprogram an off-the-shelf pre-trained English language transformer and
benchmark it on a set of protein physicochemical prediction tasks (secondary
structure, stability, homology, stability) as well as on a biomedically
relevant set of protein function prediction tasks (antimicrobial, toxicity,
antibody affinity).


  When searching for inhabited exoplanets, understanding the boundaries of the
habitable zone around the parent star is key. If life can strongly influence
its global environment, then we would expect the boundaries of the habitable
zone to be influenced by the presence of life. Here using a simple abstract
model of `tangled-ecology' where life can influence a global parameter,
labelled as temperature, we investigate the boundaries of the habitable zone of
our model system. As with other models of life-climate interactions, the
species act to regulate the temperature. However, the system can also
experience `punctuations', where the system's state jumps between different
equilibria. Despite this, an ensemble of systems still tends to sustain or even
improve conditions for life on average, a feature we call Entropic Gaia. The
mechanism behind this is sequential selection with memory which is discussed in
detail. With this modelling framework we investigate questions about how Gaia
can affect and ultimately extend the habitable zone to what we call the Gaian
habitable zone. This generates concrete predictions for the size of the
habitable zone around stars, suggests directions for future work on the
simulation of exoplanets and provides insight into the Gaian bottleneck
hypothesis and the habitability/inhabitance paradox.


  Advancements in technologies related to working with omics data require novel
computation methods to fully leverage information and help develop a better
understanding of human diseases. This paper studies the effects of introducing
graph contrastive learning to help leverage graph structure and information to
produce better representations for downstream classification tasks for
multi-omics datasets. We present a learnining framework named Multi-Omics Graph
Contrastive Learner(MOGCL) which outperforms several aproaches for integrating
multi-omics data for supervised learning tasks. We show that pre-training graph
models with a contrastive methodology along with fine-tuning it in a supervised
manner is an efficient strategy for multi-omics data classification.


  Throughout history, prevention and control of dengue transmission have
challenged public health authorities worldwide. In the last decades, the
interaction of multiple factors, such as environmental and climate variability,
has influenced increments in incidence and geographical spread of the virus. In
Costa Rica, a country characterized by multiple microclimates separated by
short distances, dengue has been endemic since its introduction in 1993.
Understanding the role of climatic and environmental factors in the seasonal
and inter-annual variability of disease spread is essential to develop
effective surveillance and control efforts. In this study, we conducted a
wavelet time series analysis of weekly climate, local environmental variables,
and dengue cases (2001-2019) from 32 cantons in Costa Rica to identify
significant periods (e.g., annual, biannual) in which climate and environmental
variables co-varied with dengue cases. Wavelet coherence analysis was used to
characterize seasonality, multi-year outbreaks, and relative delays between the
time series. Results show that dengue outbreaks occurring every 3 years in
cantons located in the country's Central, North, and South Pacific regions were
highly coherent with the Oceanic Ni\~no 3.4 and the Tropical North Caribbean
Index (TNA). Dengue cases were in phase with El Ni\~no 3.4 and TNA, with El
Ni\~no 3.4 ahead of dengue cases by roughly nine months and TNA ahead by less
than three months. Annual dengue outbreaks were coherent with local
environmental variables (NDWI, EVI, Evapotranspiration, and Precipitation) in
most cantons except those located in the Central, South Pacific, and South
Caribbean regions of the country. The local environmental variables were in
phase with dengue cases and were ahead by around three months.


  The integration of multi-modal data, such as pathological images and genomic
data, is essential for understanding cancer heterogeneity and complexity for
personalized treatments, as well as for enhancing survival predictions. Despite
the progress made in integrating pathology and genomic data, most existing
methods cannot mine the complex inter-modality relations thoroughly.
Additionally, identifying explainable features from these models that govern
preclinical discovery and clinical prediction is crucial for cancer diagnosis,
prognosis, and therapeutic response studies. We propose PONET- a novel
biological pathway-informed pathology-genomic deep model that integrates
pathological images and genomic data not only to improve survival prediction
but also to identify genes and pathways that cause different survival rates in
patients. Empirical results on six of The Cancer Genome Atlas (TCGA) datasets
show that our proposed method achieves superior predictive performance and
reveals meaningful biological interpretations. The proposed method establishes
insight into how to train biologically informed deep networks on multimodal
biomedical data which will have general applicability for understanding
diseases and predicting response and resistance to treatment.


  In this paper we introduce some simple models, based on rolling dice, to
explore mechanisms proposed to explain planetary habitability. The idea is to
study these selection mechanisms in an analytically tractable setting,
isolating their consequences from other details which can confound or obscure
their effect in more realistic models. We find that while the observable of
interest, the face value shown on the die, `improves' over time in all models,
for two of the more popular ideas: Selection by Survival and Sequential
Selection, this is down to sampling effects. A modified version of Sequential
Selection, Sequential Selection with Memory, implies a statistical tendency for
systems to improve over time. We discuss the implications of this and its
relationship to the ideas of the `inhabitance paradox' and the `Gaian
bottleneck'.


  An effective evasion strategy allows prey to survive encounters with
predators. Prey are generally thought to escape in a direction that is either
random or serves to maximize the minimum distance from the predator. Here we
introduce a comprehensive approach to determine the most likely evasion
strategy among multiple hypotheses and the role of biomechanical constraints on
the escape response of prey fish. Through a consideration of six strategies
with sensorimotor noise and previous kinematic measurements, our analysis shows
that zebrafish larvae generally escape in a direction orthogonal to the
predator's heading. By sensing only the predator's heading, this orthogonal
strategy maximizes the distance from fast-moving predators, and, when operating
within the biomechanical constraints of the escape response, it provides the
best predictions of prey behavior among all alternatives. This work
demonstrates a framework for resolving the strategic basis of evastion in
predator-prey interactions, which could be applied to a broad diversity of
animals.


  Discovery of the molecular candidates for applications in drug targets,
biomolecular systems, catalysts, photovoltaics, organic electronics, and
batteries, necessitates development of machine learning algorithms capable of
rapid exploration of the chemical spaces targeting the desired functionalities.
Here we introduce a novel approach for the active learning over the chemical
spaces based on hypothesis learning. We construct the hypotheses on the
possible relationships between structures and functionalities of interest based
on a small subset of data and introduce them as (probabilistic) mean functions
for the Gaussian process. This approach combines the elements from the symbolic
regression methods such as SISSO and active learning into a single framework.
The primary focus of constructing this framework is to approximate physical
laws in an active learning regime toward a more robust predictive performance,
as traditional evaluation on hold-out sets in machine learning doesn't account
for out-of-distribution effects and may lead to a complete failure on unseen
chemical space. Here, we demonstrate it for the QM9 dataset, but it can be
applied more broadly to datasets from both domains of molecular and solid-state
materials sciences.


  Current protein language models (pLMs) predominantly focus on single-chain
protein sequences and often have not accounted for constraints on generative
design imposed by protein-protein interactions. To address this gap, we present
paired Antibody T5 (pAbT5), an encoder-decoder model to generate complementary
heavy or light chain from its pairing partner. We show that our model respects
conservation in framework regions and variability in hypervariable domains,
demonstrated by agreement with sequence alignment and variable-length CDR
loops. We also show that our model captures chain pairing preferences through
the recovery of ground-truth chain type and gene families. Our results showcase
the potential of pAbT5 in generative antibody design, incorporating biological
constraints from chain pairing preferences.


  Recent genome-wide association studies (GWAS) have been successful in
identifying associations between genetic variants and simple cardiac parameters
derived from cardiac magnetic resonance (CMR) images. However, the emergence of
big databases including genetic data linked to CMR, facilitates investigation
of more nuanced patterns of shape variability. Here, we propose a new framework
for gene discovery entitled Unsupervised Phenotype Ensembles (UPE). UPE builds
a redundant yet highly expressive representation by pooling a set of phenotypes
learned in an unsupervised manner, using deep learning models trained with
different hyperparameters. These phenotypes are then analyzed via (GWAS),
retaining only highly confident and stable associations across the ensemble. We
apply our approach to the UK Biobank database to extract left-ventricular (LV)
geometric features from image-derived three-dimensional meshes. We demonstrate
that our approach greatly improves the discoverability of genes influencing LV
shape, identifying 11 loci with study-wide significance and 8 with suggestive
significance. We argue that our approach would enable more extensive discovery
of gene associations with image-derived phenotypes for other organs or image
modalities.


  Gene expression profiling technologies have been used in various applications
such as cancer biology. The development of gene expression profiling has
expanded the scope of target discovery in transcriptomic studies, and each
technology produces data with distinct characteristics. In order to guarantee
biologically meaningful findings using transcriptomic experiments, it is
important to consider various experimental factors in a systematic way through
statistical power analysis. In this paper, we review and discuss the power
analysis for three types of gene expression profiling technologies from a
practical standpoint, including bulk RNA-seq, single-cell RNA-seq, and
high-throughput spatial transcriptomics. Specifically, we describe the existing
power analysis tools for each research objective for each of the bulk RNA-seq
and scRNA-seq experiments, along with recommendations. On the other hand, since
there are no power analysis tools for high-throughput spatial transcriptomics
at this point, we instead investigate the factors that can influence power
analysis.


  The oldest forms of living organisms on Earth are about 3.5 billion years
old, and they are found in hydrothermal deposits, and it is often hypothesized
that life originated there. But the hydrothermal systems with a fairly strong
flow of chemical components are not the optimal place for the prebiological
self-assembly of biomolecules and for the emergence of homochirality. This
article examines the possibility for that the self-assembly of homochiral
molecules took place in an aqueous environment in the Earth's crust. Based on
the latest literature regarding the conditions in the lithosphere there are
several factors that point to that the crust could be the location for the
prebiological self-assembly of biomolecules, and there is nothing against it.
The crust and the mantle contain a substantial amount of water, and at the time
prior to the emergence of life the crust contained most likely the necessary
chemical substances for the synthesis of the biomolecules and an aqueous
environments where the homochirality could be established.


  To demonstrate their commitment, for instance during wartime, members of a
group will sometimes all engage in the same ruinous display. Such uniform,
high-cost signals are hard to reconcile with standard models of signaling. For
signals to be stable, they should honestly inform their audience; yet, uniform
signals are trivially uninformative. To explain this phenomenon, we design a
simple model, which we call the signal runaway game. In this game, senders can
express outrage at non-senders. Outrage functions as a second-order signal. By
expressing outrage at non-senders, senders draw attention to their own signal,
and benefit from its increased visibility. Using our model and a simulation, we
show that outrage can stabilize uniform signals, and can lead signal costs to
run away. Second-order signaling may explain why groups sometimes demand
displays of commitment from all their members, and why these displays can
entail extreme costs.


  Computational discovery of ideal lead compounds is a critical process for
modern drug discovery. It comprises multiple stages: hit screening, molecular
property prediction, and molecule optimization. Current efforts are disparate,
involving the establishment of models for each stage, followed by multi-stage
multi-model integration. However, this is non-ideal, as clumsy integration of
incompatible models increases research overheads, and may even reduce success
rates in drug discovery. Facilitating compatibilities requires establishing
inherent model consistencies across lead discovery stages. Towards that effect,
we propose an open deep graph learning (DGL) based pipeline: generative
adversarial feature subspace enhancement (GAFSE), which first unifies the
modeling of these stages into one learning framework. GAFSE also offers
standardized modular design and streamlined interfaces for future expansions
and community support. GAFSE combines adversarial/generative learning, graph
attention network, graph reconstruction network, and optimizes the
classification/regression loss, adversarial/generative loss, and reconstruction
loss simultaneously. Convergence analysis theoretically guarantees model
generalization performance. Exhaustive benchmarking demonstrates that the GAFSE
pipeline achieves excellent performance across almost all lead discovery
stages, while also providing valuable model interpretability. Hence, we believe
this tool will enhance the efficiency and productivity of drug discovery
researchers.


  In this paper we are interested in modeling the production of mucus by
diatoms under the constraint of a nutrient limitation. The initial questioning
comes from the observation of the so-called ''liga'' on the Aquitaine coast.
The biological origin of the phenomenon is presented and discussed based on the
existing litterature.A mathematical model is proposed and its theoretical
properties are analized: well-posedness and differentiability with respect to
the model parameters.Finally, numerical experiments are provided, investigating
the possibility of parameter identification for the model using chemostat-type
experiments.


  In varying environments it is beneficial for organisms to utilize available
cues to infer the conditions they may encounter and express potentially
favorable traits. However, external cues can be unreliable or too costly to
use. We consider an alternative strategy where organisms exploit internal
sources of information. Even without sensing environmental cues, their internal
states may become correlated with the environment as a result of selection,
which then form a memory that helps predict future conditions. To demonstrate
the adaptive value of such internal memory in varying environments, we revisit
the classic example of seed dormancy in annual plants. Previous studies have
considered the germination fraction of seeds and its dependence on
environmental cues. In contrast, we consider a model of germination fraction
that depends on the seed age, which is an internal state that can serve as a
memory. We show that, if the environmental variation has temporal structure,
then age-dependent germination fractions will allow the population to have an
increased long-term growth rate. The more organisms can remember through their
internal states, the higher growth rate a population can potentially achieve.
Our results suggest experimental ways to infer internal memory and its benefit
for adaptation in varying environments.


  Understanding the viral dynamics and immunizing antibodies of the severe
acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is crucial for devising
better therapeutic and prevention strategies for COVID-19. Here, we present a
Bayesian hierarchical model that jointly estimates the diagnostic RNA viral
load reflecting genomic materials of SARS-CoV-2, the subgenomic RNAs (sgRNA)
viral load reflecting active viral replication, and the rate and timing of
seroconversion reflecting presence of antibodies. Our proposed method accounts
for the dynamical relationship and correlation structure between the two types
of viral load, allows for borrowing of information between viral load and
antibody data, and identifies potential correlates of viral load
characteristics and propensity for seroconversion. We demonstrate the features
of the joint model through application to the COVID-19 PEP study and conduct a
cross-validation exercise to illustrate the model's ability to impute the sgRNA
viral trajectories for people who only had diagnostic viral load data.


  Recent studies at individual cell resolution have revealed phenotypic
heterogeneity in nominally clonal tumor cell populations. The heterogeneity
affects cell growth behaviors, which can result in departure from the idealized
uniform exponential growth of the cell population. Here we measured the
stochastic time courses of growth of an ensemble of populations of HL60
leukemia cells in cultures, starting with distinct initial cell numbers to
capture a departure from the {uniform exponential growth model for the initial
growth (``take-off'')}. Despite being derived from the same cell clone, we
observed significant variations in the early growth patterns of individual
cultures with statistically significant differences in growth dynamics, which
could be explained by the presence of inter-converting subpopulations with
different growth rates, and which could last for many generations. Based on the
hypothesis of existence of multiple subpopulations, we developed a branching
process model that was consistent with the experimental observations.


  Given several number sequences, determining the longest common subsequence is
a classical problem in computer science. This problem has applications in
bioinformatics, especially determining transposable genes. Nevertheless,
related works only consider how to find one longest common subsequence. In this
paper, we consider how to determine the uniqueness of the longest common
subsequence. If there are multiple longest common subsequences, we also
determine which number appears in all/some/none of the longest common
subsequences. We focus on four scenarios: (1) linear sequences without
duplicated numbers; (2) circular sequences without duplicated numbers; (3)
linear sequences with duplicated numbers; (4) circular sequences with
duplicated numbers. We develop corresponding algorithms and apply them to gene
sequencing data.


  Motivation: Detection of structural variants (SV) from the alignment of
sample DNA reads to the reference genome is an important problem in
understanding human diseases. Long reads that can span repeat regions, along
with an accurate alignment of these long reads play an important role in
identifying novel SVs. Long read sequencers such as nanopore sequencing can
address this problem by providing very long reads but with high error rates,
making accurate alignment challenging. Many errors induced by nanopore
sequencing have a bias because of the physics of the sequencing process and
proper utilization of these error characteristics can play an important role in
designing a robust aligner for SV detection problems. In this paper, we design
and evaluate HQAlign, an aligner for SV detection using nanopore sequenced
reads. The key ideas of HQAlign include (i) using basecalled nanopore reads
along with the nanopore physics to improve alignments for SVs (ii)
incorporating SV specific changes to the alignment pipeline (iii) adapting
these into existing state-of-the-art long read aligner pipeline, minimap2
(v2.24), for efficient alignments.
  Results: We show that HQAlign captures about 4%-6% complementary SVs across
different datasets which are missed by minimap2 alignments while having a
standalone performance at par with minimap2 for real nanopore reads data. For
the common SV calls between HQAlign and minimap2, HQAlign improves the start
and the end breakpoint accuracy for about 10%-50% of SVs across different
datasets. Moreover, HQAlign improves the alignment rate to 89.35% from minimap2
85.64% for nanopore reads alignment to recent telomere-to-telomere CHM13
assembly, and it improves to 86.65% from 83.48% for nanopore reads alignment to
GRCh37 human genome.


  Accurate in-silico prediction of conformational B-cell epitopes would lead to
major improvements in disease diagnostics, drug design and vaccine development.
A variety of computational methods, mainly based on machine learning
approaches, have been developed in the last decades to tackle this challenging
problem. Here, we rigorously benchmarked nine state-of-the-art conformational
B-cell epitope prediction webservers, including generic and antibody-specific
methods, on a dataset of over 250 antibody-antigen structures. The results of
our assessment and statistical analyses show that all the methods achieve very
low performances, and some do not perform better than randomly generated
patches of surface residues. In addition, we also found that commonly used
consensus strategies that combine the results from multiple webservers are at
best only marginally better than random. Finally, we applied all the predictors
to the SARS-CoV-2 spike protein as an independent case study, and showed that
they perform poorly in general, which largely recapitulates our benchmarking
conclusions. We hope that these results will lead to greater caution when using
these tools until the biases and issues that limit current methods have been
addressed, promote the use of state-of-the-art evaluation methodologies in
future publications, and suggest new strategies to improve the performance of
conformational B-cell epitope prediction methods.


  Biological cells encode information about their environment through
biochemical signaling networks that control their internal state and response.
This information is often encoded in the dynamical patterns of the signaling
molecules, rather than just their instantaneous concentrations. Here, we
analytically calculate the information contained in these dynamics for a number
of paradigmatic cases in the linear regime, for both static and time-dependent
input signals. When considering oscillatory output dynamics, we report the
emergence of synergy between successive measurements, meaning that the joint
information in two measurements exceeds the sum of the individual information.
We extend our analysis numerically beyond the scope of linear input encoding to
reveal synergetic effects in the cases of frequency or damping modulation, both
of which are relevant to classical biochemical signaling systems.


  Living cells can leverage correlations in environmental fluctuations to
predict the future environment and mount a response ahead of time. To this end,
cells need to encode the past signal into the output of the intracellular
network from which the future input is predicted. Yet, storing information is
costly while not all features of the past signal are equally informative on the
future input signal. Here, we show, for two classes of input signals, that
cellular networks can reach the fundamental bound on the predictive information
as set by the information extracted from the past signal: push-pull networks
can reach this information bound for Markovian signals, while networks that
take a temporal derivative can reach the bound for predicting the future
derivative of non-Markovian signals. However, the bits of past information that
are most informative about the future signal are also prohibitively costly. As
a result, the optimal system that maximizes the predictive information for a
given resource cost is, in general, not at the information bound. Applying our
theory to the chemotaxis network of Escherichia coli reveals that its adaptive
kernel is optimal for predicting future concentration changes over a broad
range of background concentrations, and that the system has been tailored to
predicting these changes in shallow gradients.


  Stationary distributions of continuous time Markov chains (CTMCs) are often a
main interest, but hard to find. We consider CTMCs modelling reaction networks,
and characterise complex balanced distributions (provided they exist) of
reaction networks with arbitrary transition functions through conditions on the
cycles of their corresponding digraph. The proof works by constructing a
dynamically equivalent reaction network of disjoint cycles. We further derive a
sufficient condition for the existence of a complex balanced distribution, and
give precise conditions on when it is necessary. The sufficient condition holds
for mass-action kinetics or if the digraph consists of only cyclic connected
components. Hence, we fully characterise the existence and form of complex
balanced distributions. Moreover, these results can be used to design complex
balanced reaction networks with a given stationary distribution.


  Protein folding neural networks (PFNNs) such as AlphaFold predict remarkably
accurate structures of proteins compared to other approaches. However, the
robustness of such networks has heretofore not been explored. This is
particularly relevant given the broad social implications of such technologies
and the fact that biologically small perturbations in the protein sequence do
not generally lead to drastic changes in the protein structure. In this paper,
we demonstrate that AlphaFold does not exhibit such robustness despite its high
accuracy. This raises the challenge of detecting and quantifying the extent to
which these predicted protein structures can be trusted. To measure the
robustness of the predicted structures, we utilize (i) the root-mean-square
deviation (RMSD) and (ii) the Global Distance Test (GDT) similarity measure
between the predicted structure of the original sequence and the structure of
its adversarially perturbed version. We prove that the problem of minimally
perturbing protein sequences to fool protein folding neural networks is
NP-complete. Based on the well-established BLOSUM62 sequence alignment scoring
matrix, we generate adversarial protein sequences and show that the RMSD
between the predicted protein structure and the structure of the original
sequence are very large when the adversarial changes are bounded by (i) 20
units in the BLOSUM62 distance, and (ii) five residues (out of hundreds or
thousands of residues) in the given protein sequence. In our experimental
evaluation, we consider 111 COVID-19 proteins in the Universal Protein resource
(UniProt), a central resource for protein data managed by the European
Bioinformatics Institute, Swiss Institute of Bioinformatics, and the US Protein
Information Resource. These result in an overall GDT similarity test score
average of around 34%, demonstrating a substantial drop in the performance of
AlphaFold.


  Bacterial cells are sensitive to a range of external signals used to learn
the environment. These incoming external signals are then processed using a
Gene Regulatory Network (GRN), exhibiting similarities to modern computing
algorithms. An in-depth analysis of gene expression dynamics suggests an
inherited Gene Regulatory Neural Network (GRNN) behavior within the GRN that
enables the cellular decision-making based on received signals from the
environment and neighbor cells. In this study, we extract a sub-network of
\textit{Pseudomonas aeruginosa} GRN that is associated with one virulence
factor: pyocyanin production as a use case to investigate the GRNN behaviors.
Further, using Graph Neural Network (GNN) architecture, we model a single
species biofilm to reveal the role of GRNN dynamics on ecosystem-wide
decision-making. Varying environmental conditions, we prove that the extracted
GRNN computes input signals similar to natural decision-making process of the
cell. Identifying of neural network behaviors in GRNs may lead to more accurate
bacterial cell activity predictive models for many applications, including
human health-related problems and agricultural applications. Further, this
model can produce data on causal relationships throughout the network, enabling
the possibility of designing tailor-made infection-controlling mechanisms. More
interestingly, these GRNNs can perform computational tasks for bio-hybrid
computing systems.


  Households play an important role in disease dynamics. Many infections
happening there due to the close contact, while mitigation measures mainly
target the transmission between households. Therefore, one can see households
as boosting the transmission depending on household size. To study the effect
of household size and size distribution, we differentiated the within and
between household reproduction rate. There are basically no preventive
measures, and thus the close contacts can boost the spread. We explicitly
incorporated that typically only a fraction of all household members are
infected. Thus, viewing the infection of a household of a given size as a
splitting process generating a new, small fully infected sub-household and a
remaining still susceptible sub-household we derive a compartmental ODE-model
for the dynamics of the sub-households. In this setting, the basic reproduction
number as well as prevalence and the peak of an infection wave in a population
with given households size distribution can be computed analytically. We
compare numerical simulation results of this novel household-ODE model with
results from an agent--based model using data for realistic household size
distributions of different countries. We find good agreement of both models
showing the catalytic effect of large households on the overall disease
dynamics.


  Classification algorithms using RNA-Sequencing (RNA-Seq) data as input are
used in a variety of biological applications. By nature, RNA-Seq data is
subject to uncontrolled fluctuations both within and especially across
datasets, which presents a major difficulty for a trained classifier to
generalize to an external dataset. Replacing raw gene counts with the rank of
gene counts inside an observation has proven effective to mitigate this
problem. However, the rank of a feature is by definition relative to all other
features, including highly variable features that introduce noise in the
ranking. To address this problem and obtain more robust ranks, we propose a
logistic regression model, optirank, which learns simultaneously the parameters
of the model and the genes to use as a reference set in the ranking. We show
the effectiveness of this method on simulated data. We also consider real
classification tasks, which present different kinds of distribution shifts
between train and test data. Those tasks concern a variety of applications,
such as cancer of unknown primary classification, identification of specific
gene signatures, and determination of cell type in single-cell RNA-Seq
datasets. On those real tasks, optirank performs at least as well as the
vanilla logistic regression on classical ranks, while producing sparser
solutions. In addition, to increase the robustness against dataset shifts, we
propose a multi-source learning scheme and demonstrate its effectiveness when
used in combination with rank-based classifiers.


  Biological systems are characterized by the ubiquitous roles of weak, that
is, non-covalent molecular interactions, small, often very small, numbers of
specific molecules per cell, and Brownian motion. These combine to produce
stochastic behaviors at all levels from the molecular and cellular to the
behavioral. That said, students are rarely introduced to the ubiquitous role of
stochastic processes in biological systems, and how they produce unpredictable
behaviors. Here I present the case that they need to be and provide some
suggestions as to how it might be approached.


  Opioid use disorder (OUD) continuously poses major public health challenges
and social implications worldwide with dramatic rise of opioid dependence
leading to potential abuse. Despite that a few pharmacological agents have been
approved for OUD treatment, the efficacy of said agents for OUD requires
further improvement in order to provide safer and more effective
pharmacological and psychosocial treatments. Preferable therapeutic treatments
of OUD rely on the advances in understanding the neurobiological mechanism of
opioid dependence. Proteins including mu, delta, kappa, nociceptin, and zeta
opioid receptors are the direct targets of opioids. Each receptor has a large
protein-protein interaction (PPI) network, that behaves differently when
subjected to various treatments, thus increasing the complexity in the drug
development process for an effective opioid addiction treatment. The report
below analyzes the work by presenting a PPI-network informed machine-learning
study of OUD. We have examined more than 500 proteins in the five opioid
receptor networks and subsequently collected 74 inhibitor datasets. Machine
learning models were constructed by pairing gradient boosting decision tree
(GBDT) algorithm with two advanced natural language processing (NLP)-based
molecular fingerprints. With these models, we systematically carried out
evaluations of screening and repurposing potential of drug candidates for four
opioid receptors. In addition, absorption, distribution, metabolism, excretion,
and toxicity (ADMET) properties were also considered in the screening of
potential drug candidates. Our study can be a valuable and promising tool of
pharmacological development for OUD treatments.


  Altruistic punishment, where individuals incur personal costs to punish
others who have harmed third parties, presents an evolutionary conundrum as it
undermines individual fitness. Resolving this puzzle is crucial for
understanding the emergence and maintenance of human cooperation. This study
investigates the role of an alternative strategy, the exit option, in
explaining altruistic punishment. We analyze a two-stage prisoner's dilemma
game in well-mixed and networked populations, considering both finite and
infinite scenarios. Our findings reveal that the exit option does not
significantly enhance altruistic punishment in well-mixed populations. However,
in networked populations, the exit option enables the existence of altruistic
punishment and gives rise to complex dynamics, including cyclic dominance and
bi-stable states. This research contributes to our understanding of costly
punishment and sheds light on the effectiveness of different voluntary
participation strategies in addressing the conundrum of punishment.


  The search for life elsewhere in the universe is one of the central aims of
science in the 21st century. While most of this work is aimed at planets
orbiting other stars, the search for life in our own Solar System is an
important part of this endeavour. Venus is often thought to have too harsh an
environment for life, but it may have been a more hospitable place in the
distant past. If life evolved there in the past then the cloud decks of Venus
are the only remaining niche where life as we know it might survive today. The
discovery of the molecule phosphine, PH$_3$, in these clouds has reinvigorated
research looking into the possibility of life in the clouds. In this review we
examine the background to studies of the possibility of life on Venus, discuss
the discovery of phosphine, review conflicting and confirming observations and
analyses, and then look forward to future observations and space missions that
will hopefully provide definitive answers as to the origin of phosphine on
Venus and to the question of whether life might exist there.


  Previous studies have shown that the topological properties of a complex
network, such as heterogeneity and average degree, affect the evolutionary game
dynamics on it. However, traditional numerical simulations are usually
time-consuming and demand a lot of computational resources. In this paper, we
propose the method of dynamical approximate master equations (DAMEs) to
accurately approximate the evolutionary outcomes on complex networks. We
demonstrate that the accuracy of DAMEs supersedes previous standard pairwise
approximation methods, and DAMEs require far fewer computational resources than
traditional numerical simulations. We use prisoner's dilemma and snowdrift game
on regular and scale-free networks to demonstrate the applicability of DAMEs.
Overall, our method facilitates the investigation of evolutionary dynamics on a
broad range of complex networks, and provides new insights into the puzzle of
cooperation.


  Globally, we are witnessing the rise of complex, non-communicable diseases
(NCDs) related to changes in our daily environments. Obesity, asthma,
cardiovascular disease, and type 2 diabetes are part of a long list of
"lifestyle" diseases that were rare throughout human history but are now
common. A key idea from anthropology and evolutionary biology--the evolutionary
mismatch hypothesis--seeks to explain this phenomenon. It posits that humans
evolved in environments that radically differ from the ones experienced by most
people today, and thus traits that were advantageous in past environments may
now be "mismatched" and disease-causing. This hypothesis is, at its core, a
genetic one: it predicts that loci with a history of selection will exhibit
"genotype by environment" (GxE) interactions and have differential health
effects in ancestral versus modern environments. Here, we discuss how this
concept could be leveraged to uncover the genetic architecture of NCDs in a
principled way. Specifically, we advocate for partnering with small-scale,
subsistence-level groups that are currently transitioning from environments
that are arguably more "matched" with their recent evolutionary history to
those that are more "mismatched". These populations provide diverse genetic
backgrounds as well as the needed levels and types of environmental variation
necessary for mapping GxE interactions in an explicit mismatch framework. Such
work would make important contributions to our understanding of environmental
and genetic risk factors for NCDs across diverse ancestries and sociocultural
contexts.


  The Wright--Fisher diffusion is important in population genetics in modelling
the evolution of allele frequencies over time subject to the influence of
biological phenomena such as selection, mutation, and genetic drift. Simulating
paths of the process is challenging due to the form of the transition density.
We present EWF, a robust and efficient sampler which returns exact draws for
the diffusion and diffusion bridge processes, accounting for general models of
selection including those with frequency-dependence. Given a configuration of
selection, mutation, and endpoints, EWF returns draws at the requested sampling
times from the law of the corresponding Wright--Fisher process. Output was
validated by comparison to approximations of the transition density via the
Kolmogorov--Smirnov test and QQ plots. All software is available at
https://github.com/JaroSant/EWF


  Porcine reproductive and respiratory syndrome virus (PRRSV) remains widely
distributed across the U.S. swine industry. Between-farm movement of animals
and transportation vehicles, along with local transmission are the primary
routes by which PRRSV is spread. Given the farm-to-farm proximity in high pig
production areas, local transmission is an important pathway in the spread of
PRRSV; however, there is limited understanding of the role local transmission
plays in the dissemination of PRRSV, specifically, the distance at which there
is increased risk for transmission from infected to susceptible farms. We used
a spatial and spatiotemporal kernel density approach to estimate PRRSV relative
risk and utilized a Bayesian spatiotemporal hierarchical model to assess the
effects of environmental variables, between-farm movement data, and on-farm
biosecurity features on PRRSV outbreaks. The maximum spatial distance
calculated through the kernel density approach was 15.3 km in 2018, 17.6 km in
2019, and 18 km in 2020. Spatiotemporal analysis revealed greater variability
throughout the study period, with significant differences between the different
farm types. Sow farms were consistently categorized as high risk farm types,
while downstream farms (i.e., finisher and nursery farms) had more susceptible
farms within areas of significant-high relative risk. Factors associated with
PRRSV outbreaks were farms with higher number of access points to barns, higher
numbers of outgoing movements of pigs, and higher number of days where
temperatures were between 4{\deg}C and 10{\deg}C. Results obtained from this
study may be used to guide the reinforcement of biosecurity and surveillance
strategies at farms and areas within the distance threshold of PRRSV positive
farms.


  The Spot Prawn trap fishery off the west coast of British Columbia (BC) is
managed using a fixed escapement strategy that aims to prevent recruitment
overfishing while maximizing expected long-term yield by closing the fishery
when the catch rate of spawners, projected to the following spring, drops below
1.7 spawners per trap (the de jure rule). We develop a management strategy
evaluation framework for BC's Spot Prawn fishery that examines the expected
performance of the management procedure implemented in practice (the de facto
rule), which was significantly more conservative than the de jure rule, usually
closing the fishery when spawner catch rates were at least twice as high as
specified by the de jure rule. Simulations indicate that the de facto spawner
index rule using average empirical March 31st targets from 2000 to 2019
maintains most stocks near or above 0.8 BMSY with or without accounting for
environmental effects and/or increasing future SST on recruitment. Abundance
indices were found to be strongly hyperstable, with fishing efficiency 1.5 to
3.0 times higher under low biomass than high biomass.


  Understanding the mutational history of tumor cells is a critical endeavor in
unraveling the mechanisms underlying cancer. Since the modeling of tumor cell
evolution employs labeled trees, researchers are motivated to develop different
methods to assess and compare mutation trees and other labeled trees. While the
Robinson-Foulds distance is a widely utilized metric for comparing phylogenetic
trees, its applicability to labeled trees reveals certain limitations. This
paper introduces the $k$-Robinson-Foulds dissimilarity measures, tailored to
address the challenges of labeled tree comparison. The Robinson-Foulds distance
is succinctly expressed as n-RF in the space of labeled trees with n nodes.
Like the Robinson-Foulds distance, the k-Robinson-Foulds is a pseudometric for
multiset-labeled trees and becomes a metric in the space of 1-labeled trees. By
setting k to a small value, the k-Robinson-Foulds dissimilarity can capture
analogous local regions in two labeled trees with different size or different
labels.


  Retrosynthesis is the cornerstone of organic chemistry, providing chemists in
material and drug manufacturing access to poorly available and brand-new
molecules. Conventional rule-based or expert-based computer-aided synthesis has
obvious limitations, such as high labor costs and limited search space. In
recent years, dramatic breakthroughs driven by artificial intelligence have
revolutionized retrosynthesis. Here we aim to present a comprehensive review of
recent advances in AI-based retrosynthesis. For single-step and multi-step
retrosynthesis both, we first list their goal and provide a thorough taxonomy
of existing methods. Afterwards, we analyze these methods in terms of their
mechanism and performance, and introduce popular evaluation metrics for them,
in which we also provide a detailed comparison among representative methods on
several public datasets. In the next part we introduce popular databases and
established platforms for retrosynthesis. Finally, this review concludes with a
discussion about promising research directions in this field.


  Tropical rainforests exhibit a rich repertoire of spatial patterns emerging
from the intricate relationship between the microscopic interaction between
species. In particular, the distribution of vegetation clusters can shed much
light on the underlying process regulating the ecosystem. Analyzing the
distribution of vegetation clusters at different resolution scales, we show the
first robust evidence of scale-invariant clusters of vegetation, suggesting the
coexistence of multiple intertwined scales in the collective dynamics of
tropical rainforests. We use field data and computational simulations to
confirm our hypothesis, proposing a predictor that could be particularly
interesting to monitor the ecological resilience of the world's 'green lungs'.


  Eukaryotic cells maintain their inner order by a hectic process of
distillation of molecular factors taking place on the surface of their lipid
membranes. To understand the properties of this molecular sorting process, a
physical model of the process has been recently proposed [arXiv:1811.06760],
based on (a) the phase separation of a single, initially dispersed molecular
species into spatially localized sorting domains on the lipid membrane, and (b)
domain-induced membrane bending leading to the nucleation of submicrometric
lipid vesicles, naturally enriched in the molecules of the engulfed sorting
domain. The analysis of the model has shown the existence of an optimal region
of the parameter space where sorting is most efficient. Here, the model is
extended to account for the simultaneous distillation of a pool of distinct
molecular species. We find that the mean time spent by sorted molecules on the
membrane increases with the heterogeneity of the pool (i.e., the number of
distinct molecular species sorted) according to a simple scaling law, and that
a large number of distinct molecular species can in principle be sorted in
parallel on a typical cell membrane region without significantly interfering
with each other. Moreover, sorting is found to be most efficient when the
distinct molecular species have comparable homotypic affinities. We also
consider how valence (i.e., the average number of interacting neighbors of a
molecule in a sorting domain) affects the sorting process, finding that
higher-valence molecules can be sorted with greater efficiency than
lower-valence molecules.


  Computational models are increasingly used in high-impact decision making in
science, engineering, and medicine. The National Aeronautics and Space
Administration (NASA) uses computational models to perform complex experiments
that are otherwise prohibitively expensive or require a microgravity
environment. Similarly, the Food and Drug Administration (FDA) and European
Medicines Agency (EMA) have began accepting models and simulations as form of
evidence for pharmaceutical and medical device approval. It is crucial that
computational models meet a standard of credibility when using them in
high-stakes decision making. For this reason, institutes including NASA, the
FDA, and the EMA have developed standards to promote and assess the credibility
of computational models and simulations. However, due to the breadth of models
these institutes assess, these credibility standards are mostly qualitative and
avoid making specific recommendations. On the other hand, modeling and
simulation in systems biology is a narrow domain and several standards are
already in place. As systems biology models increase in complexity and
influence, the development of a credibility assessment system is crucial. Here
we review existing standards in systems biology, credibility standards in other
science, engineering, and medical fields, and propose the development of a
credibility standard for systems biology models.


  Tumour heterogeneity is increasingly recognized as a major obstacle to
therapeutic success across neuro-oncology. Gliomas are characterised by
distinct combinations of genetic and epigenetic alterations, resulting in
complex interactions across multiple molecular pathways. Predicting disease
evolution and prescribing individually optimal treatment requires statistical
models complex enough to capture the intricate (epi)genetic structure
underpinning oncogenesis. Here, we formalize this task as the inference of
distinct patterns of connectivity within hierarchical latent representations of
genetic networks. Evaluating multi-institutional clinical, genetic, and outcome
data from 4023 glioma patients over 14 years, across 12 countries, we employ
Bayesian generative stochastic block modelling to reveal a hierarchical network
structure of tumour genetics spanning molecularly confirmed glioblastoma, IDH-
wildtype; oligodendroglioma, IDH-mutant and 1p/19q codeleted; and astrocytoma,
IDH- mutant. Our findings illuminate the complex dependence between features
across the genetic landscape of brain tumours, and show that generative network
models reveal distinct signatures of survival with better prognostic fidelity
than current gold standard diagnostic categories.


  Understanding and accurately predicting protein-ligand binding affinity are
essential in the drug design and discovery process. At present, machine
learning-based methodologies are gaining popularity as a means of predicting
binding affinity due to their efficiency and accuracy, as well as the
increasing availability of structural and binding affinity data for
protein-ligand complexes. In biomolecular studies, graph theory has been widely
applied since graphs can be used to model molecules or molecular complexes in a
natural manner. In the present work, we upgrade the graph-based learners for
the study of protein-ligand interactions by integrating extensive atom types
such as SYBYL and extended connectivity interactive features (ECIF) into
multiscale weighted colored graphs (MWCG). By pairing with the gradient
boosting decision tree (GBDT) machine learning algorithm, our approach results
in two different methods, namely $^\text{sybyl}\text{GGL}$-Score and
$^\text{ecif}\text{GGL}$-Score. Both of our models are extensively validated in
their scoring power using three commonly used benchmark datasets in the drug
design area, namely CASF-2007, CASF-2013, and CASF-2016. The performance of our
best model $^\text{sybyl}\text{GGL}$-Score is compared with other
state-of-the-art models in the binding affinity prediction for each benchmark.
While both of our models achieve state-of-the-art results, the SYBYL atom-type
model $^\text{sybyl}\text{GGL}$-Score outperforms other methods by a wide
margin in all benchmarks.


  One of the main challenges in drug discovery is to find molecules that bind
specifically and strongly to their target protein while having minimal binding
to other proteins. By predicting binding affinity, it is possible to identify
the most promising candidates from a large pool of potential compounds,
reducing the number of compounds that need to be tested experimentally.
Recently, deep learning methods have shown superior performance than
traditional computational methods for making accurate predictions on large
datasets. However, the complexity and time-consuming nature of these methods
have limited their usage and development. Quantum machine learning is an
emerging technology that has the potential to improve many classical machine
learning algorithms. In this work we present a hybrid quantum-classical
convolutional neural network, which is able to reduce by 20% the complexity of
the classical network while maintaining optimal performance in the predictions.
Additionally, it results in a significant time savings of up to 40% in the
training process, which means a meaningful speed up of the drug discovery
process.


  In this paper, we study a two stream species Lotka-Volterra competition patch
model with the patches aligned along a line. The two species are supposed to be
identical except for the diffusion rates. For each species, the diffusion rates
between patches are the same, while the drift rates vary. Our results show that
the convexity of the drift rates has a significant impact on the competition
outcomes: if the drift rates are convex, then the species with larger diffusion
rate wins the competition; if the drift rates are concave, then the species
with smaller diffusion rate wins the competition.


  In this paper, we study a three-patch two-species Lotka-Volterra competition
patch model over a stream network. The individuals are subject to both random
and directed movements, and the two species are assumed to be identical except
for the movement rates. The environment is heterogeneous, and the carrying
capacity is lager in upstream locations. We treat one species as a resident
species and investigate whether the other species can invade or not. Our
results show that the spatial heterogeneity of environment and the magnitude of
the drift rates have a large impact on the competition outcomes of the stream
species.