  We present the MSA-to-protein transformer, a generative model of protein
sequences conditioned on protein families represented by multiple sequence
alignments (MSAs). Unlike existing approaches to learning generative models of
protein families, the MSA-to-protein transformer conditions sequence generation
directly on a learned encoding of the multiple sequence alignment,
circumventing the need for fitting dedicated family models. By training on a
large set of well-curated multiple sequence alignments in Pfam, our
MSA-to-protein transformer generalizes well to protein families not observed
during training and outperforms conventional family modeling approaches,
especially when MSAs are small. Our generative approach accurately models
epistasis and indels and allows for exact inference and efficient sampling
unlike other approaches. We demonstrate the protein sequence modeling
capabilities of our MSA-to-protein transformer and compare it with alternative
sequence modeling approaches in comprehensive benchmark experiments.


  We present a systematic study of the influence of cell geometry on the
orientational distribution of microtubules (MTs) nucleated from a single
microtubule organizing center (MTOC). For simplicity we consider an elliptical
cell geometry, a setting appropriate to a generic non-spherical animal cell.
Within this context we introduce four models of increasing complexity, in each
case introducing additional mechanisms that govern the interaction of the MTs
with the cell boundary. In order, we consider the cases: MTs that can bind to
the boundary with a fixed mean residence time (M0), force-producing MTs that
can slide on the boundary towards the cell poles (MS), MTs that interact with a
generic polarity factor that is transported and deposited at the boundary, and
which in turn stabilizes the MTs at the boundary (MP), and a final model in
which both sliding and stabilization by polarity factors is taken into account
(MSP). In the baseline model (M0), the exponential length distribution of MTs
causes most of the interactions at the cell boundary to occur along the shorter
transverse direction in the cell, leading to transverse biaxial order. MT
sliding (MS) is able to reorient the main axis of this biaxial order along the
longitudinal axis. The polarization mechanism introduced in MP and MSP
overrules the geometric bias towards bipolar order observed in M0 and MS, and
allows the establishment of unipolar order either along the short- (MP) or the
long cell axis (MSP). The behavior of the latter two models can be
qualitatively reproduced by a very simple toy model with discrete MT
orientations.


  The field of movement ecology has seen a rapid increase in high-resolution
data in recent years, leading to the development of numerous statistical and
numerical methods to analyse relocation trajectories. Data are often collected
at the level of the individual and for long periods that may encompass a range
of behaviours. Here, we use the power spectral density (PSD) to characterise
the random movement patterns of a black-winged kite (Elanus caeruleus) and a
white stork (Ciconia ciconia). The tracks are first segmented and clustered
into different behaviours (movement modes), and for each mode we measure the
PSD and the ageing properties of the process. For the foraging kite we find
$1/f$ noise, previously reported in ecological systems mainly in the context of
population dynamics, but not for movement data. We further suggest plausible
models for each of the behavioural modes by comparing both the measured PSD
exponents and the distribution of the single-trajectory PSD to known
theoretical results and simulations.


  In the last few years, de novo molecular design using machine learning has
made great technical progress but its practical deployment has not been as
successful. This is mostly owing to the cost and technical difficulty of
synthesizing such computationally designed molecules. To overcome such
barriers, various methods for synthetic route design using deep neural networks
have been studied intensively in recent years. However, little progress has
been made in designing molecules and their synthetic routes simultaneously.
Here, we formulate the problem of simultaneously designing molecules with the
desired set of properties and their synthetic routes within the framework of
Bayesian inference. The design variables consist of a set of reactants in a
reaction network and its network topology. The design space is extremely large
because it consists of all combinations of purchasable reactants, often in the
order of millions or more. In addition, the designed reaction networks can
adopt any topology beyond simple multistep linear reaction routes. To solve
this hard combinatorial problem, we present a powerful sequential Monte Carlo
algorithm that recursively designs a synthetic reaction network by sequentially
building up single-step reactions. In a case study of designing drug-like
molecules based on commercially available compounds, compared with heuristic
combinatorial search methods, the proposed method shows overwhelming
performance in terms of computational efficiency and coverage and novelty with
respect to existing compounds.


  A large number of inorganic and organic compounds are able to bind DNA and
form complexes, among which drug-related molecules are important. Chromatin
accessibility changes not only directly affects drug-DNA interactions, but also
promote or inhibit the expression of critical genes associated with drug
resistance by affecting the DNA binding capacity of TFs and transcriptional
regulators. However, Biological experimental techniques for measuring it are
expensive and time consuming. In recent years, several kinds of computational
methods have been proposed to identify accessible regions of the genome.
Existing computational models mostly ignore the contextual information of bases
in gene sequences. To address these issues, we proposed a new solution named
SemanticCAP. It introduces a gene language model which models the context of
gene sequences, thus being able to provide an effective representation of a
certain site in gene sequences. Basically, we merge the features provided by
the gene language model into our chromatin accessibility model. During the
process, we designed some methods to make feature fusion smoother. Compared
with other systems under public benchmarks, our model proved to have better
performance.


  Cancer is one of the deadliest diseases worldwide. Accurate diagnosis and
classification of cancer subtypes are indispensable for effective clinical
treatment. Promising results on automatic cancer subtyping systems have been
published recently with the emergence of various deep learning methods.
However, such automatic systems often overfit the data due to the high
dimensionality and scarcity. In this paper, we propose to investigate automatic
subtyping from an unsupervised learning perspective by directly constructing
the underlying data distribution itself, hence sufficient data can be generated
to alleviate the issue of overfitting. Specifically, we bypass the strong
Gaussianity assumption that typically exists but fails in the unsupervised
learning subtyping literature due to small-sized samples by vector
quantization. Our proposed method better captures the latent space features and
models the cancer subtype manifestation on a molecular basis, as demonstrated
by the extensive experimental results.


  Proteins are fundamental biological entities mediating key roles in cellular
function and disease. This paper introduces a multi-scale graph construction of
a protein -- HoloProt -- connecting surface to structure and sequence. The
surface captures coarser details of the protein, while sequence as primary
component and structure -- comprising secondary and tertiary components --
capture finer details. Our graph encoder then learns a multi-scale
representation by allowing each level to integrate the encoding from level(s)
below with the graph at that level. We test the learned representation on
different tasks, (i.) ligand binding affinity (regression), and (ii.) protein
function prediction (classification). On the regression task, contrary to
previous methods, our model performs consistently and reliably across different
dataset splits, outperforming all baselines on most splits. On the
classification task, it achieves a performance close to the top-performing
model while using 10x fewer parameters. To improve the memory efficiency of our
construction, we segment the multiplex protein surface manifold into molecular
superpixels and substitute the surface with these superpixels at little to no
performance loss.


  Deploying generative machine learning techniques to generate novel chemical
structures based on molecular fingerprint representation has been well
established in molecular design. Typically, sequential learning (SL) schemes
such as hidden Markov models (HMM) and, more recently, in the sequential deep
learning context, recurrent neural network (RNN) and long short-term memory
(LSTM) were used extensively as generative models to discover unprecedented
molecules. To this end, emission probability between two states of atoms plays
a central role without considering specific chemical or physical properties.
Lack of supervised domain knowledge can mislead the learning procedure to be
relatively biased to the prevalent molecules observed in the training data that
are not necessarily of interest. We alleviated this drawback by augmenting the
training data with domain knowledge, e.g. quantitative estimates of the
drug-likeness score (QEDs). As such, our experiments demonstrated that with
this subtle trick called enriched sequential learning (ESL), specific patterns
of particular interest can be learnt better, which led to generating de novo
molecules with ameliorated QEDs.


  Nucleation is the rate-determining step in the kinetics of many self-assembly
processes. However, the importance of nucleation in the kinetics of DNA-origami
self-assembly, which involves both the binding of staple strands and the
folding of the scaffold strand, is unclear. Here, using Monte Carlo simulations
of a lattice model of DNA origami, we find that some, but not all, designs can
have a nucleation barrier and that this barrier disappears at lower
temperatures, rationalizing the success of isothermal assembly. We show that
the height of the nucleation barrier depends primarily on the coaxial stacking
of staples that are adjacent on the same helix, a parameter that can be
modified with staple design. Creating a nucleation barrier to DNA-origami
assembly could be useful in optimizing assembly times and yields, while
eliminating the barrier may allow for fast molecular sensors that can
assemble/disassemble without hysteresis in response to changes in the
environment.


  Proteins in complex with small molecule ligands represent the core of
structure-based drug discovery. However, three-dimensional representations are
absent from most deep-learning-based generative models. We here present a
graph-based generative modeling technology that encodes explicit 3D
protein-ligand contacts within a relational graph architecture. The models
combine a conditional variational autoencoder that allows for activity-specific
molecule generation with putative contact generation that provides predictions
of molecular interactions within the target binding pocket. We show that
molecules generated with our 3D procedure are more compatible with the binding
pocket of the dopamine D2 receptor than those produced by a comparable
ligand-based 2D generative method, as measured by docking scores, expected
stereochemistry, and recoverability in commercial chemical databases. Predicted
protein-ligand contacts were found among highest-ranked docking poses with a
high recovery rate. This work shows how the structural context of a protein
target can be used to enhance molecule generation.


  To interpret molecular dynamics simulations of biomolecular systems,
systematic dimensionality reduction methods are commonly employed. Among
others, this includes principal component analysis (PCA) and time-lagged
independent component analysis (TICA), which aim to maximize the variance and
the timescale of the first components, respectively. A crucial first step of
such an analysis is the identification of suitable and relevant input
coordinates (the so-called features), such as backbone dihedral angles and
interresidue distances. As typically only a small subset of those coordinates
is involved in a specific biomolecular process, it is important to discard the
remaining uncorrelated motions or weakly correlated noise coordinates. This is
because they may exhibit large amplitudes or long timescales and therefore will
be erroneously be considered important by PCA and TICA, respectively. To
discriminate collective motions underlying functional dynamics from
uncorrelated motions, the correlation matrix of the input coordinates is
block-diagonalized by a clustering method. This strategy avoids possible bias
due to presumed functional observables and conformational states or variation
principles that maximize variance or timescales. Considering several linear and
nonlinear correlation measures and various clustering algorithms, it is shown
that the combination of linear correlation and the Leiden community detection
algorithm yields excellent results for all considered model systems. These
include the functional motion of T4 lysozyme to demonstrate the successful
identification of collective motion, as well as the folding of villin headpiece
to highlight the physical interpretation of the correlated motions in terms of
a functional mechanism.


  DNA has been considered a promising medium for storing digital information.
As an essential step in the DNA-based data storage workflow, coding algorithms
are responsible to implement functions including bit-to-base transcoding, error
correction, etc. In previous studies, these functions are normally realized by
introducing multiple algorithms. Here, we report a graph-based architecture,
named SPIDER-WEB, providing an all-in-one coding solution by generating
customized algorithms automatically. SPIDERWEB is able to correct a maximum of
4% edit errors in the DNA sequences including substitution and
insertion/deletion (indel), with only 5.5% redundant symbols. Since no DNA
sequence pretreatment is required for the correcting and decoding processes,
SPIDER-WEB offers the function of real-time information retrieval, which is
305.08 times faster than the speed of single-molecule sequencing techniques.
Our retrieval process can improve 2 orders of magnitude faster compared to the
conventional one under megabyte-level data and can be scalable to fit
exabyte-level data. Therefore, SPIDER-WEB holds the potential to improve the
practicability in large-scale data storage applications.


  We introduce a variant of the voter model in which agents may have different
degrees of confidence on their opinions. Those with low confidence are normal
voters whose state can change upon a single contact with a different
neighboring opinion. However, confidence increases with opinion reinforcement
and, above a certain threshold, these agents become zealots that do not change
opinion. We show that both strategies, normal voters and zealots, may coexist,
leading to a competition between two different kinetic mechanisms:
curvature-driven growth and interfacial noise. The kinetically constrained
zealots are formed well inside the clusters, away from the different opinions
at the surfaces that help keep the confidence not so high. Normal voters
concentrate in a region around the interfaces and their number, that is related
with the distance between the surface and the zealotry bulk, depends on the
rate the confidence changes. Despite this interface being rough and fragmented,
typical of the voter model, the presence of zealots in the bulk of these
domains, induces a curvature-driven dynamics, similar to the low temperature
coarsening behavior of the non-conserved Ising model after a temperature
quench.


  Complex energy landscapes often arise in biological systems, e.g. for protein
folding, biochemical reactions or intracellular transport processes. Their
physical effects are often reflected in the first-passage times arising from
these energy landscapes. However, their calculation is notoriously challenging
and it is often difficult to identify the most relevant features of a given
energy landscape. Here we show how this can be achieved by coarse-graining the
Fokker-Planck equation to a master equation and decomposing its first-passage
times in an iterative process. We apply this method to the electrostatic
interaction between two rods of nonmuscle myosin II (NM2), which is the main
molecular motor for force generation in nonmuscle cells. Energy landscapes are
computed directly from the amino acid sequences of the three different
isoforms. Our approach allows us to identify the most relevant energy barriers
for their self-assembly into nonmuscle myosin II minifilaments and how they
change under force. In particular, we find that antiparallel configurations are
more stable than parallel ones, but also show more changes under mechanical
loading. Our work demonstrates the rich dynamics that can be expected for
NM2-assemblies under mechanical load and in general shows how one can identify
the most relevant energy barriers in complex energy landscapes.


  To understand the mechanisms underlying species coexistence, ecologists often
study invasion growth rates of theoretical and data-driven models. These growth
rates correspond to average per-capita growth rates of one species with respect
to an ergodic measure supporting other species. In the ecological literature,
coexistence often is equated with the invasion growth rates being positive.
Intuitively, positive invasion growth rates ensure that species recover from
being rare. To provide a mathematically rigorous framework for this approach,
we prove theorems that answer two questions: (i) When do the signs of the
invasion growth rates determine coexistence? (ii) When signs are sufficient,
which invasion growth rates need to be positive? We focus on deterministic
models and equate coexistence with permanence, i.e., a global attractor bounded
away from extinction. For models satisfying certain technical assumptions, we
introduce invasion graphs where vertices correspond to proper subsets of
species (communities) supporting an ergodic measure and directed edges
correspond to potential transitions between communities due to invasions by
missing species. These directed edges are determined by the signs of invasion
growth rates. When the invasion graph is acyclic (i.e. there is no sequence of
invasions starting and ending at the same community), we show that permanence
is determined by the signs of the invasion growth rates. In this case,
permanence is characterized by the invasibility of all -i communities, i.e.,
communities without species i where all other missing species having negative
invasion growth rates. We show that dissipative Lotka-Volterra models satisfy
our technical assumptions and computing their invasion graphs reduces to
solving systems of linear equations. We provide additional applications of the
results and discuss open problems.


  We propose a revisited version of the in vivo model of the cyanobacterial
circadian clock. Our aim is to address the lack of robustness predicted for the
mutant cyanobacteria without transcriptional regulation of the original model.
For this, we rely on an in vitro model of the clock describing explicitly the
hexameric structure of the core protein of the clock. Our model is able to
reproduce oscillatory behavior for the mutant, as observed experimentally,
without finely tuned parameters.


  The native conformation of structured proteins is stabilized by a complex
network of interactions. We analyzed the elementary patterns that constitute
such network and ranked them according to their importance in shaping protein
sequence design. To achieve this goal, we employed a cluster expansion of the
partition function in the space of sequences and evaluated numerically the
statistical importance of each cluster. An important feature of this procedure
is that it is applied to a dense, finite system. We found that patterns that
contribute most to the partition function are cycles with even numbers of
nodes, while cliques are typically detrimental. Each cluster also gives a
contribute to the sequence entropy, which is a measure of the evolutionary
designability of a fold. We compared the entropies associated with different
interaction patterns to their abundances in the native structures of real
proteins.


  In an extant population, how much information do extant individuals provide
on the pedigree of their ancestors? Recent work by Kim, Mossel, Ramnarayan and
Turner (2020) studied this question under a number of simplifying assumptions,
including random mating, fixed length inheritance blocks and sufficiently large
founding population. They showed that under these conditions if the average
number of offspring is a sufficiently large constant, then it is possible to
recover a large fraction of the pedigree structure and genetic content by an
algorithm they named REC-GEN.
  We are interested in studying the performance of REC-GEN on simulated data
generated according to the model. As a first step, we improve the running time
of the algorithm. However, we observe that even the faster version of the
algorithm does not do well in any simulations in recovering the pedigree beyond
2 generations. We claim that this is due to the inbreeding present in any
setting where the algorithm can be run, even on simulated data. To support the
claim we show that a main step of the algorithm, called ancestral
reconstruction, performs accurately in a idealized setting with no inbreeding
but performs poorly in random mating populations.
  To overcome the poor behavior of REC-GEN we introduce a Belief-Propagation
based heuristic that accounts for the inbreeding and performs much better in
our simulations.


  In the last two decades dengue cases increased significantly throughout the
world. In several regions dengue re-emerged, particularly in Latin America,
where dengue cases not only increased but also occurred more frequently. It is
therefore necessary to understand the mechanisms that drive epidemic outbreaks
in non-endemic regions, to help in the design of control strategies. We develop
a stochastic model that includes climate variables, social structure, and
mobility between a non-endemic city and an endemic area. We choose as a case
study the non-endemic city of San Ram{\'o}n de la Nueva Or{\'a}n, located in
Northwest Argentina. Human mobility is intense through the border with Bolivia,
where dengue transmission is sustained during the whole year. City population
was modelled as a meta-population taking into account households and population
data for each patch. Climate variability was considered by including rainfall,
relative humidity and temperature time series into the models. Those climatic
variables were input of a mosquito population ecological model, which in turn
is coupled to an epidemiological model. Different hypotheses regarding people's
mobility between an endemic and non-endemic area are tested, taking into
account the local climatic variation, typical of the non-endemic city.
Simulations are qualitatively consistent with weekly clinical data reported
from 2009 to 2016. Our model results allow to explain the observed pattern of
outbreaks, that alternates large dengue epidemics and several years with
smaller outbreaks. We found that the number of vectors per host and an
effective reproductive number are proxies for large epidemics, both related
with climate variability such as rainfall and temperature, opening the
possibility to test these meteorological variables for forecast purposes.


  Understanding the coexistence of diverse species in a changing environment is
an important problem in community ecology. Bet-hedging is a strategy that helps
species survive in such changing environments. However, studies of bet-hedging
have often focused on the expected long-term growth rate of the species by
itself, neglecting competition with other coexisting species. Here we study the
extinction risk of a bet-hedging species in competition with others. We show
that there are three contributions to the extinction risk. The first is the
usual demographic fluctuation due to stochastic reproduction and selection
processes in finite populations. The second, due to the fluctuation of
population growth rate caused by environmental changes, may counterintuitively
reduce the extinction risk for small populations. Besides those two, we reveal
a third contribution, which is unique to bet-hedging species that diversify
into multiple phenotypes: The phenotype composition of the population will
fluctuate over time, resulting in increased extinction risk. We compare such
compositional fluctuation to the demographic and environmental contributions,
showing how they have different effects on the extinction risk depending on the
population size, generation overlap, and environmental correlation.


  Planktonic communities are extremely diverse and include a vast number of
rare species. The dynamics of these rare species is best described by
individual-based models. However, individual-based approaches to planktonic
diversity face substantial difficulties, due to the large number of individuals
required to make realistic predictions. In this paper, we study diversity of
planktonic communities by means of a spatial coalescence model, that
incorporates transport by oceanic currents. As a main advantage, our approach
requires simulating a number of individuals equal to the size of the sample one
is interested in, rather than the size of the entire community. By theoretical
analysis and simulations, we explore the conditions upon which our coalescence
model is equivalent to individual-based dynamics. As an application, we use our
model to predict the impact of chaotic advection by oceanic currents on
biodiversity. We conclude that the coalescent approach permits to simulate
marine microbial communities much more efficiently than with individual-based
models.


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


  Population dynamics is the study of temporal and spatial variation in the
size of populations of organisms and is a major part of population ecology. One
of the main difficulties in analyzing population dynamics is that we can only
obtain observation data with coarse time intervals from fixed-point
observations due to experimental costs or measurement constraints. Recently,
modeling population dynamics by using continuous normalizing flows (CNFs) and
dynamic optimal transport has been proposed to infer the sample trajectories
from a fixed-point observed population. While the sample behavior in CNFs is
deterministic, the actual sample in biological systems moves in an essentially
random yet directional manner. Moreover, when a sample moves from point A to
point B in dynamical systems, its trajectory typically follows the principle of
least action in which the corresponding action has the smallest possible value.
To satisfy these requirements of the sample trajectories, we formulate the
Lagrangian Schr\"odinger bridge (LSB) problem and propose to solve it
approximately by modeling the advection-diffusion process with regularized
neural SDE. We also develop a model architecture that enables faster
computation of the loss function. Experimental results show that the proposed
method can efficiently approximate the population-level dynamics even for
high-dimensional data and that using the prior knowledge introduced by the
Lagrangian enables us to estimate the sample-level dynamics with stochastic
behavior.


  We investigate the problem of using a learning-based strategy to stabilize a
synthetic toggle switch via an external control approach. To overcome the data
efficiency problem that would render the algorithm unfeasible for practical use
in synthetic biology, we adopt a sim-to-real paradigm where the policy is
learnt via training on a simplified model of the toggle switch and it is then
subsequently exploited to control a more realistic model of the switch
parameterized from in-vivo experiments. Our in-silico experiments confirm the
viability of the approach suggesting its potential use for in-vivo control
implementations.


  Infectious disease outbreaks can have a disruptive impact on public health
and societal processes. As decision making in the context of epidemic
mitigation is hard, reinforcement learning provides a methodology to
automatically learn prevention strategies in combination with complex epidemic
models. Current research focuses on optimizing policies w.r.t. a single
objective, such as the pathogen's attack rate. However, as the mitigation of
epidemics involves distinct, and possibly conflicting criteria (i.a.,
prevalence, mortality, morbidity, cost), a multi-objective approach is
warranted to learn balanced policies. To lift this decision-making process to
real-world epidemic models, we apply deep multi-objective reinforcement
learning and build upon a state-of-the-art algorithm, Pareto Conditioned
Networks (PCN), to learn a set of solutions that approximates the Pareto front
of the decision problem. We consider the first wave of the Belgian COVID-19
epidemic, which was mitigated by a lockdown, and study different deconfinement
strategies, aiming to minimize both COVID-19 cases (i.e., infections and
hospitalizations) and the societal burden that is induced by the applied
mitigation measures. We contribute a multi-objective Markov decision process
that encapsulates the stochastic compartment model that was used to inform
policy makers during the COVID-19 epidemic. As these social mitigation measures
are implemented in a continuous action space that modulates the contact matrix
of the age-structured epidemic model, we extend PCN to this setting. We
evaluate the solution returned by PCN, and observe that it correctly learns to
reduce the social burden whenever the hospitalization rates are sufficiently
low. In this work, we thus show that multi-objective reinforcement learning is
attainable in complex epidemiological models and provides essential insights to
balance complex mitigation policies.


  A fractional-order compartmental model was recently used to describe real
data of the first wave of the COVID-19 pandemic in Portugal [Chaos Solitons
Fractals 144 (2021), Art. 110652]. Here, we modify that model in order to
correct time dimensions and use it to investigate the third wave of COVID-19
that occurred in Portugal from December 2020 to February 2021, and that has
surpassed all previous waves, both in number and consequences. A new fractional
optimal control problem is then formulated and solved, with vaccination and
preventive measures as controls. A cost-effectiveness analysis is carried out,
and the obtained results are discussed.


  This work systematically conducts a data analysis based on the numbers of
both cumulative and daily confirmed COVID-19 cases and deaths in a time span
through April 2020 to June 2022 for over 200 countries around the world. Such
research feature aims to reveal the temporal and spatial evolution of the
country-level distribution observed in COVID-19 pandemic, and obtains some
interesting results as follows. (1) The distributions of the numbers for
cumulative confirmed cases and deaths obey power-law in early stages of
COVID-19 and stretched exponential function in subsequent course. (2) The
distributions of the numbers for daily confirmed cases and deaths obey
power-law in early and late stages of COVID-19 and stretched exponential
function in middle stages. The crossover region between power-law and stretched
exponential behaviour seems to depend on the evolution of "infection" event and
"death" event. Such observation implies a kind of important symmetry related to
the dynamics process of COVID-19 spreading. (3) The distributions of the
normalized numbers for each metric show a temporal scaling behaviour in 2-year
period, and are well described by stretched exponential function. The
observation of power-law and stretched exponential behaviour in such
country-level distributions suggests underlying intrinsic dynamics of a virus
spreading process in human interconnected society. And thus it is important for
understanding and mathematically modeling the COVID-19 pandemic.


  Dynamical systems modeling, particularly via systems of ordinary differential
equations, has been used to effectively capture the temporal behavior of
different biochemical components in signal transduction networks. Despite the
recent advances in experimental measurements, including sensor development and
'-omics' studies that have helped populate protein-protein interaction networks
in great detail, modeling in systems biology lacks systematic methods to
estimate kinetic parameters and quantify associated uncertainties. This is
because of multiple reasons, including sparse and noisy experimental
measurements, lack of detailed molecular mechanisms underlying the reactions,
and missing biochemical interactions. Additionally, the inherent nonlinearities
with respect to the states and parameters associated with the system of
differential equations further compound the challenges of parameter estimation.
In this study, we propose a comprehensive framework for Bayesian parameter
estimation and complete quantification of the effects of uncertainties in the
data and models. We apply these methods to a series of signaling models of
increasing mathematical complexity. Systematic analysis of these dynamical
systems showed that parameter estimation depends on data sparsity, noise level,
and model structure, including the existence of multiple steady states. These
results highlight how focused uncertainty quantification can enrich systems
biology modeling and enable additional quantitative analyses for parameter
estimation.


  The Hierarchical Redundancy Eliminated Tree Augmented Naive Bayes (HRE-TAN)
classifier is a semi-naive Bayesian model that learns a type of hierarchical
redundancy-free tree-like feature representation to estimate the data
distribution. In this work, we propose two new types of positive feature values
prioritized hierarchical redundancy eliminated tree augmented naive Bayes
classifiers that focus on features bearing positive instance values. The two
newly proposed methods are applied to 28 real-world bioinformatics datasets
showing better predictive performance than the conventional HRE-TAN classifier.


  In this paper we discuss reflection groups and root systems, in particular
non-crystallographic ones, and a Clifford algebra framework for both these
concepts. A review of historical as well as more recent work on viral capsid
symmetries motivates the focus on the icosahedral root system $H_3$. We discuss
a notion of affine extension for non-crystallographic groups with applications
to fullerenes and viruses. The icosahedrally ordered component of the nucleic
acid within the virus capsid and the interaction between the two have shed
light on the viral assembly process with interesting applications to antiviral
therapies, drug delivery and nanotechnology. The Clifford algebra framework is
very natural, as it uses precisely the structure that is already implicit in
this root system and reflection group context, i.e. a vector space with an
inner product. In addition, it affords a uniquely simple reflection formula, a
double cover of group transformations, and more insight into the geometry, e.g.
the geometry of the Coxeter plane. This approach made possible a range of root
system induction proofs, such as the constructions of $E_8$ from $H_3$ and the
exceptional 4D root systems from 3D root systems. This makes explicit various
connections between what Arnold called Trinities (sets of three exceptional
cases). In fact this generalises further since the induction construction
contains additional cases, namely two infinite families of cases. It therefore
actually yields $ADE$ correspondences between three sets of different
mathematical concepts that are usually thought of separately as polytopes,
subgroups of $SU(2)$ and $ADE$ Lie algebras. Here we connect them explicitly
and in a unified way by thinking of them as three different sets of root
systems.


  A significant proportion of the infections driving the current {SARS-CoV-2}
pandemic are transmitted asymptomatically. Here we introduce and study a simple
epidemic model with separate compartments comprising asymptomatic and
symptomatic infected individuals. The linear dynamics determining the outbreak
condition of the model is equivalent to a renewal theory approach with
exponential waiting time distributions. Exploiting a nontrivial conservation
law of the full nonlinear dynamics, we derive analytic bounds on the peak
number of infections in the absence and presence of mitigation through
isolation and testing. The bounds are compared to numerical solutions of the
differential equations.


  A long-standing goal of evolutionary biology is to decode how gene regulatory
processes contribute to organismal diversity, both within and between species.
This question has remained challenging to answer, due both to the difficulties
of predicting function from non-coding sequence, and to the technological
constraints of laboratory research with non-model taxa. However, a recent
methodological development in functional genomics, the massively parallel
reporter assay (MPRA), makes it possible to test thousands to millions of
sequences for regulatory activity in a single in vitro experiment. It does so
by combining traditional, single-locus episomal reporter assays (e.g.,
luciferase reporter assays) with the scalability of high-throughput sequencing.
In this perspective, we discuss the execution, advantages, and limitations of
MPRAs for research in evolutionary biology. We review recent studies that have
made use of this approach to address explicitly evolutionary questions,
highlighting study designs that we believe are particularly well-positioned to
gain from MPRA approaches. Additionally, we propose solutions for extending
these powerful assays to rare taxa and those with limited genomic resources. In
doing so, we underscore the broad potential of MPRAs to drive genome-scale
functional evolutionary genetics studies in non-traditional model organisms.


  1. Abrupt environmental changes can lead to evolutionary shifts in trait
evolution. Identifying these shifts is an important step in understanding the
evolutionary history of phenotypes.
  2. We propose an ensemble variable selection method (R package ELPASO) for
the evolutionary shift detection task and compare it with existing methods (R
packages l1ou and PhylogeneticEM) under several scenarios.
  3. The performances of methods are highly dependent on the selection
criterion. When the signal sizes are small, the methods using the Bayesian
information criterion (BIC) have better performances. And when the signal sizes
are large enough, the methods using the phylogenetic Bayesian information
criterion (pBIC) (Khabbazian et al., 2016) have better performance. Moreover,
the performance is heavily impacted by measurement error and tree
reconstruction error.
  4. Ensemble method + pBIC tends to perform less conservatively than l1ou +
pBIC, and Ensemble method + BIC is more conservatively than l1ou + BIC.
PhylogeneticEM is even more conservative with small signal sizes and falls
between l1ou + pBIC and Ensemble method + BIC with large signal sizes. The
results can differ between the methods, but none clearly outperforms the
others. By applying multiple methods to a single dataset, we can access the
robustness of each detected shift, based on the agreement among methods.


  Non-native conformations drive protein misfolding diseases, complicate
bioengineering efforts, and fuel molecular evolution. No current experimental
technique is well-suited for elucidating them and their phenotypic effects.
Especially intractable are the transient conformations populated by
intrinsically disordered proteins. We describe an approach to systematically
discover, stabilize, and purify native and non-native conformations, generated
in vitro or in vivo, and directly link conformations to molecular, organismal,
or evolutionary phenotypes. This approach involves high-throughput disulfide
scanning (HTDS) of the entire protein. To reveal which disulfides trap which
chromatographically resolvable conformers, we devised a deep-sequencing method
for double-Cys variant libraries of proteins that precisely and simultaneously
locates both Cys residues within each polypeptide. HTDS of the abundant E. coli
periplasmic chaperone HdeA revealed distinct classes of disordered hydrophobic
conformers with variable cytotoxicity depending on where the backbone was
cross-linked. HTDS can bridge conformational and phenotypic landscapes for many
proteins that function in disulfide-permissive environments.


  Machine learning (ML) is widely used in drug discovery to train models that
predict protein-ligand binding. These models are of great value to medicinal
chemists, in particular if they provide case-specific insight into the physical
interactions that drive the binding process. In this study we derive ML models
from over 50 fragment-screening campaigns to introduce two important elements
that we believe are absent in most -- if not all -- ML studies of this type
reported to date: First, alongside the observed hits we use to train our
models, we incorporate true misses and show that these experimentally validated
negative data are of significant importance to the quality of the derived
models. Second, we provide a physically interpretable and verifiable
representation of what the ML model considers important for successful binding.
This representation is derived from a straightforward attribution procedure
that explains the prediction in terms of the (inter-)action of chemical
environments. Critically, we validate the attribution outcome on a large scale
against prior annotations made independently by expert molecular modellers. We
find good agreement between the key molecular substructures proposed by the ML
model and those assigned manually, even when the model's performance in
discriminating hits from misses is far from perfect. By projecting the
attribution onto predefined interaction prototypes (pharmacophores), we show
that ML allows us to formulate simple rules for what drives fragment binding
against a target automatically from screening data.


  Genome-wide association studies (GWAS) have successfully identified a large
number of genetic variants associated with traits and diseases. However, it
still remains challenging to fully understand functional mechanisms underlying
many associated variants. This is especially the case when we are interested in
variants shared across multiple phenotypes. To address this challenge, we
propose graph-GPA 2.0 (GGPA 2.0), a novel statistical framework to integrate
GWAS datasets for multiple phenotypes and incorporate functional annotations
within a unified framework. We conducted simulation studies to evaluate GGPA
2.0. The results indicate that incorporating functional annotation data using
GGPA 2.0 does not only improve detection of disease-associated variants, but
also allows to identify more accurate relationships among diseases. We analyzed
five autoimmune diseases and five psychiatric disorders with the functional
annotations derived from GenoSkyline and GenoSkyline-Plus and the prior disease
graph generated by biomedical literature mining. For autoimmune diseases, GGPA
2.0 identified enrichment for blood, especially B cells and regulatory T cells
across multiple diseases. Psychiatric disorders were enriched for brain,
especially prefrontal cortex and inferior temporal lobe for bipolar disorder
(BIP) and schizophrenia (SCZ), respectively. Finally, GGPA 2.0 successfully
identified the pleiotropy between BIP and SCZ. These results demonstrate that
GGPA 2.0 can be a powerful tool to identify associated variants associated with
each phenotype or those shared across multiple phenotypes, while also promoting
understanding of functional mechanisms underlying the associated variants.


  This study deals with an existing mathematical model of asymmetrically
interacting agents. We analyze the following two previously unfocused features
of the model: (i) synchronization of growth rates and (ii) initial value
dependence of damped oscillation. By applying the techniques of variable
transformation and time-scale separation, we perform the stability analysis of
a diverging solution. We find that (i) all growth rates synchronize to the same
value that is as small as the smallest growth rate and (ii) oscillatory
dynamics appear if the initial value of the slowest-growing agent is
sufficiently small. Furthermore, our analytical method proposes a way to apply
stability analysis to an exponentially diverging solution, which we believe is
also a contribution of this study. Although the employed model is originally
proposed as a model of infectious disease, we do not discuss its biological
relevance but merely focus on the technical aspects.


  Mounting evidence shows that oscillatory activity is widespread in cell
signaling. Here we review some of this recent evidence, focusing on both the
molecular mechanisms that potentially underlie such dynamical behavior, and the
potential advantages that signaling oscillations might have in cell function.
The biological processes considered include tissue maintenance, embryonic
development and wound healing. With the aid of mathematical modeling, we show
that a common principle, namely delayed negative feedback, underpins this wide
variety of phenomena.


  AlphaFold2 (AF) is a promising tool, but is it accurate enough to predict
single mutation effects? Here, we report that the localized structural
deformation between protein pairs differing by only 1-3 mutations -- as
measured by the effective strain -- is correlated across \num{3901}
experimental and AF-predicted structures. Furthermore, analysis of ${\sim}
11000$ proteins shows that the local structural change correlates with various
phenotypic changes. These findings suggest that AF can predict the range and
magnitude of single-mutation effects on average, and we propose a method to
improve precision of AF predictions and to indicate when predictions are
unreliable.


  Computational models starting from large ensembles of evolutionarily related
protein sequences capture a representation of protein families and learn
constraints associated to protein structure and function. They thus open the
possibility for generating novel sequences belonging to protein families.
Protein language models trained on multiple sequence alignments, such as MSA
Transformer, are highly attractive candidates to this end. We propose and test
an iterative method that directly employs the masked language modeling
objective to generate sequences using MSA Transformer. We demonstrate that the
resulting sequences score as well as natural sequences, for homology,
coevolution and structure-based measures. For large protein families, our
synthetic sequences have similar or better properties compared to sequences
generated by Potts models, including experimentally-validated ones. Moreover,
for small protein families, our generation method based on MSA Transformer
outperforms Potts models. Our method also more accurately reproduces the
higher-order statistics and the distribution of sequences in sequence space of
natural data than Potts models. MSA Transformer is thus a strong candidate for
protein sequence generation and protein design.


  By the peak of COVID-19 restrictions on April 8, 2020, up to 1.5 billion
students across 188 countries were by the suspension of physical attendance in
schools. Schools were among the first services to reopen as vaccination
campaigns advanced. With the emergence of new variants and infection waves, the
question now is to find safe protocols for the continuation of school
activities. We need to understand how reliable these protocols are under
different levels of vaccination coverage, as many countries have a meager
fraction of their population vaccinated, including Uganda where the coverage is
about 8\%. We investigate the impact of face-to-face classes under different
protocols and quantify the surplus number of infected individuals in a city.
Using the infection transmission when schools were closed as a baseline, we
assess the impact of physical school attendance in classrooms with poor air
circulation. We find that (i) resuming school activities with people only
wearing low-quality masks leads to a near fivefold city-wide increase in the
number of cases even if all staff is vaccinated, (ii) resuming activities with
students wearing good-quality masks and staff wearing N95s leads to about a
threefold increase, (iii) combining high-quality masks and active monitoring,
activities may be carried out safely even with low vaccination coverage. These
results highlight the effectiveness of good mask-wearing. Compared to ICU
costs, high-quality masks are inexpensive and can help curb the spreading.
Classes can be carried out safely, provided the correct set of measures are
implemented.


  The terrestrial water cycle links the soil and atmosphere moisture reservoirs
through four fluxes: precipitation, evaporation, runoff and atmospheric
moisture convergence. Each of these fluxes is essential for human and ecosystem
well-being. However, predicting how the water cycle responds to changes in
vegetation cover, remains a challenge (Lawrence and Vandecar, 2015; Ellison et
al., 2017; te Wierik et al., 2021). Recently, rainfall was shown to decrease
disproportionally with reduced forest transpiration following deforestation
(Baudena et al., 2021). Here, combining these findings with the law of matter
conservation, we show that in a sufficiently wet atmosphere forest
transpiration can control atmospheric moisture convergence such that increased
transpiration enhances atmospheric moisture import. Conversely, in a drier
atmosphere increased transpiration reduces atmospheric moisture convergence and
runoff. This previously unrecognized dichotomy can explain the seemingly random
observations of runoff and soil moisture sometimes increasing and sometimes
reducing in response to re-greening (e.g., Zheng et al., 2021). Evaluating the
transition between the two regimes is crucial both for characterizing the risk
posed by deforestation as well as for motivating and guiding global ecosystem
restoration.


  The absorption, distribution, metabolism, excretion, and toxicity (ADMET)
properties are important in drug discovery as they define efficacy and safety.
In this work, we applied an ensemble of features, including fingerprints and
descriptors, and a tree-based machine learning model, extreme gradient
boosting, for accurate ADMET prediction. Our model performs well in the
Therapeutics Data Commons ADMET benchmark group. For 22 tasks, our model is
ranked first in 18 tasks and top 3 in 21 tasks. The trained machine learning
models are integrated in ADMETboost, a web server that is publicly available at
https://ai-druglab.smu.edu/admet.


  Ecological management problems often involve navigating from an initial to a
desired community state. We ask whether navigation is possible without
brute-force additions and deletions of species, using actions of varying costs:
adding/deleting a small number of individuals of a species, changing the
environment, and waiting. Navigation can yield direct paths (single sequence of
actions) or shortcut paths (multiple sequences of actions with lower cost than
a direct path). We ask (1) when is non-brute-force navigation possible?; (2) do
shortcuts exist and what are their properties?; and (3) what heuristics predict
shortcut existence? Using several empirical datasets, we show that (1)
non-brute-force navigation is only possible between some state pairs, (2)
shortcuts exist between many state pairs; and (3) changes in abundance and
richness are the strongest predictors of shortcut existence, independent of
dataset and algorithm choices. State diagrams thus unveil hidden strategies for
efficiently shifting between states.


  We examine how spatial heterogeneity combines with mobility network structure
to influence vector-borne disease dynamics. Specifically, we consider a
Ross-Macdonald-type disease model on $n$ spatial locations that are coupled by
host movement on a strongly connected, weighted, directed graph. We derive a
closed form approximation to the domain reproduction number using a Laurent
series expansion, and use this approximation to compute sensitivities of the
basic reproduction number to model parameters. To illustrate how these results
can be used to help inform mitigation strategies, as a case study we apply
these results to malaria dynamics in Namibia, using published cell phone data
and estimates for local disease transmission. Our analytical results are
particularly useful for understanding drivers of transmission when mobility
sinks and transmission hot spots do not coincide.


  RNA virus (e.g., SARS-CoV-2) evolves in a complex manner. Studying RNA virus
evolution is vital for understanding molecular evolution and medicine
development. Scientists lack, however, general frameworks to characterize the
dynamics of RNA virus evolution directly from empirical data and identify
potential physical laws. To fill this gap, we present a theory to characterize
the RNA virus evolution as a physical system with absorbing states and
avalanche behaviors. This approach maps accessible biological data (e.g.,
phylogenetic tree and infection) to a general stochastic process of RNA virus
infection and evolution, enabling researchers to verify potential
self-organized criticality underlying RNA virus evolution. We apply our
framework to SARS-CoV-2, the virus accounting for the global epidemic of
COVID-19. We find that SARS-CoV-2 exhibits scale-invariant avalanches as
mean-field theory predictions. The observed scaling relation, universal
collapse, and slowly decaying auto-correlation suggest a self-organized
critical dynamics of SARS-CoV-2 evolution. Interestingly, the lineages that
emerge from critical evolution processes coincidentally match with threatening
lineages of SARS-CoV-2 (e.g., the Delta virus). We anticipate our approach to
be a general formalism to portray RNA virus evolution and help identify
potential virus lineages to be concerned.


  The hybridization and dehybridization of DNA subject to tension is relevant
to fundamental genetic processes and to the design of DNA-based mechanobiology
assays. While strong tension accelerates DNA melting and decelerates DNA
annealing, the effects of tension weaker than 5 pN are less clear. In this
study, we developed a DNA bow assay, which uses the bending rigidity of
double-stranded DNA (dsDNA) to exert weak tension on a single-stranded DNA
(ssDNA) target in the range of 2 pN to 6 pN. Combining this assay with
single-molecule FRET, we measured the hybridization and dehybridization
kinetics between a 15 nt ssDNA under tension and a 8-9 nt oligo, and found that
both the hybridization and dehybridization rates monotonically increase with
tension for various nucleotide sequences tested. These findings suggest that
the nucleated duplex in its transition state is more extended than the pure
dsDNA or ssDNA counterpart. Our simulations using the coarse-grained oxDNA2
model indicate that the increased extension of the transition state is due to
exclusion interactions between unpaired ssDNA regions in close proximity to one
another. This study highlights an example where the ideal worm-like chain
models fail to explain the kinetic behavior of DNA in the low force regime.


  Maintenance and regeneration of adult tissues rely on the self-renewal of
stem cells. Regeneration without over-proliferation requires precise regulation
of the stem cell proliferation and differentiation rates. The nature of such
regulatory mechanisms in different tissues, and how to incorporate them in
models of stem cell population dynamics, is incompletely understood. The
critical birth-death (CBD) process is widely used to model stem cell
populations, capturing key phenomena, such as scaling laws in clone size
distributions. However, the CBD process neglects regulatory mechanisms. Here,
we propose the birth-death process with volume exclusion (vBD), a variation of
the birth-death process that considers crowding effects, such as may arise due
to limited space in a stem cell niche. While the deterministic rate equations
predict a single non-trivial attracting steady state, the master equation
predicts extinction and transient distributions of stem cell numbers with three
possible behaviours: long-lived quasi-steady state, and short-lived bimodal or
unimodal distributions. In all cases, we approximate solutions to the vBD
master equation using a renormalised system-size expansion, quasi-steady state
approximation and the WKB method. Our study suggests that the size distribution
of a stem cell population bears signatures that are useful to detect negative
feedback mediated via volume exclusion.


  The discovery of novel inhibitor molecules for emerging drug-target proteins
is widely acknowledged as a challenging inverse design problem: Exhaustive
exploration of the vast chemical search space is impractical, especially when
the target structure or active molecules are unknown. Here we validate
experimentally the broad utility of a deep generative framework trained
at-scale on protein sequences, small molecules, and their mutual interactions
-- that is unbiased toward any specific target. As demonstrators, we consider
two dissimilar and relevant SARS-CoV-2 targets: the main protease and the spike
protein (receptor binding domain, RBD). To perform target-aware design of novel
inhibitor molecules, a protein sequence-conditioned sampling on the generative
foundation model is performed. Despite using only the target sequence
information, and without performing any target-specific adaptation of the
generative model, micromolar-level inhibition was observed in in vitro
experiments for two candidates out of only four synthesized for each target.
The most potent spike RBD inhibitor also exhibited activity against several
variants in live virus neutralization assays. These results therefore establish
that a single, broadly deployable generative foundation model for accelerated
hit discovery is effective and efficient, even in the most general case where
neither target structure nor binder information is available.


  A fundamental problem in drug discovery is to design molecules that bind to
specific proteins. To tackle this problem using machine learning methods, here
we propose a novel and effective framework, known as GraphBP, to generate 3D
molecules that bind to given proteins by placing atoms of specific types and
locations to the given binding site one by one. In particular, at each step, we
first employ a 3D graph neural network to obtain geometry-aware and chemically
informative representations from the intermediate contextual information. Such
context includes the given binding site and atoms placed in the previous steps.
Second, to preserve the desirable equivariance property, we select a local
reference atom according to the designed auxiliary classifiers and then
construct a local spherical coordinate system. Finally, to place a new atom, we
generate its atom type and relative location w.r.t. the constructed local
coordinate system via a flow model. We also consider generating the variables
of interest sequentially to capture the underlying dependencies among them.
Experiments demonstrate that our GraphBP is effective to generate 3D molecules
with binding ability to target protein binding sites. Our implementation is
available at https://github.com/divelab/GraphBP.


  Emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)
variants have compromised existing vaccines and posed a grand challenge to
coronavirus disease 2019 (COVID-19) prevention, control, and global economic
recovery. For COVID-19 patients, one of the most effective COVID-19 medications
is monoclonal antibody (mAb) therapies. The United States Food and Drug
Administration (U.S. FDA) has given the emergency use authorization (EUA) to a
few mAbs, including those from Regeneron, Eli Elly, etc. However, they are also
undermined by SARS-CoV-2 mutations. It is imperative to develop effective
mutation-proof mAbs for treating COVID-19 patients infected by all emerging
variants and/or the original SARS-CoV-2. We carry out a deep mutational
scanning to present the blueprint of such mAbs using algebraic topology and
artificial intelligence (AI). To reduce the risk of clinical trial-related
failure, we select five mAbs either with FDA EUA or in clinical trials as our
starting point. We demonstrate that topological AI-designed mAbs are effective
to variants of concerns and variants of interest designated by the World Health
Organization (WHO), as well as the original SARS-CoV-2. Our topological AI
methodologies have been validated by tens of thousands of deep mutational data
and their predictions have been confirmed by results from tens of experimental
laboratories and population-level statistics of genome isolates from hundreds
of thousands of patients.


  This paper focuses on three critical problems on protein classification.
Firstly, Carbohydrate-active enzyme (CAZyme) classification can help people to
understand the properties of enzymes. However, one CAZyme may belong to several
classes. This leads to Multi-label CAZyme classification. Secondly, to capture
information from the secondary structure of protein, protein classification is
modeled as graph classification problem. Thirdly, compound-protein interactions
prediction employs graph learning for compound with sequential embedding for
protein. This can be seen as classification task for compound-protein pairs.
This paper proposes three models for protein classification. Firstly, this
paper proposes a Multi-label CAZyme classification model using CNN-LSTM with
Attention mechanism. Secondly, this paper proposes a variational graph
autoencoder based subspace learning model for protein graph classification.
Thirdly, this paper proposes graph isomorphism networks (GIN) and
Attention-based CNN-LSTM for compound-protein interactions prediction, as well
as comparing GIN with graph convolution networks (GCN) and graph attention
networks (GAT) in this task. The proposed models are effective for protein
classification. Source code and data are available at
https://github.com/zshicode/GNN-AttCL-protein. Besides, this repository
collects and collates the benchmark datasets with respect to above problems,
including CAZyme classification, enzyme protein graph classification,
compound-protein interactions prediction, drug-target affinities prediction and
drug-drug interactions prediction. Hence, the usage for evaluation by benchmark
datasets can be more conveniently.


  Given its status as a classic problem and its importance to both
theoreticians and practitioners, edit distance provides an excellent lens
through which to understand how the theoretical analysis of algorithms impacts
practical implementations. From an applied perspective, the goals of
theoretical analysis are to predict the empirical performance of an algorithm
and to serve as a yardstick to design novel algorithms that perform well in
practice. In this paper, we systematically survey the types of theoretical
analysis techniques that have been applied to edit distance and evaluate the
extent to which each one has achieved these two goals. These techniques include
traditional worst-case analysis, worst-case analysis parametrized by edit
distance or entropy or compressibility, average-case analysis, semi-random
models, and advice-based models. We find that the track record is mixed. On one
hand, two algorithms widely used in practice have been born out of theoretical
analysis and their empirical performance is captured well by theoretical
predictions. On the other hand, all the algorithms developed using theoretical
analysis as a yardstick since then have not had any practical relevance. We
conclude by discussing the remaining open problems and how they can be tackled.


  We investigate a tritrophic system whose cyclic dominance is modelled by the
rock-paper-scissors game. We consider that organisms of one or two species are
affected by movement limitations, which unbalances the cyclic spatial game.
Performing stochastic simulations, we show that mobility unevenness controls
the population dynamics. In the case of one slow species, the predominant
species depends on the level of mobility restriction, with the slow species
being preponderant if the mobility limitations are substantial. If two species
face mobility limitations, our outcomes show that being higher dispersive does
not constitute an advantage in terms of population growth. On the contrary, if
organisms move with higher mobility, they expose themselves to enemies more
frequently, being more vulnerable to being eliminated. Finally, our findings
show that biodiversity benefits in regions where species are slowed.
Biodiversity loss for high mobility organisms, common to cyclic systems, may be
avoided with coexistence probability being higher for robust mobility
limitations. Our results may help biologists understand the dynamics of
unbalanced spatial systems where organisms' dispersal is fundamental to
biodiversity conservation.


  Free-energy landscapes for short peptides -- specifically for variants of the
pH Low Insertion Peptide (pHLIP) -- in the heterogeneous environment of a lipid
bilayer or cell membrane are constructed, taking into account a set of dominant
interactions and the conformational preferences of the peptide backbone. Our
methodology interprets broken internal H-bonds along the backbone of a
polypeptide as statistically interacting quasiparticles, activated from the
helix reference state. The favored conformation depends on the local
environment (ranging from polar to nonpolar), specifically on the availability
of external H-bonds (with $\mathrm{H_2O}$ molecules or lipid headgroups) to
replace internal H-bonds. The dominant side-chain contribution is accounted for
by residue-specific transfer free energies between polar and nonpolar
environments. The free-energy landscape is sensitive to the level of pH in the
aqueous environment surrounding the membrane. For high pH, we identify pathways
of descending free energy that suggest a coexistence of membrane-adsorbed
peptides with peptides in solution. A drop in pH raises the degree of
protonation of negatively charged residues and thus increases the
hydrophobicity of peptide segments near the C terminus. For low pH, we identify
insertion pathways between the membrane-adsorbed state and a stable
trans-membrane state with the C terminus having crossed the membrane.


  Predicting gene functions is a challenge for biologists in the post genomic
era. Interactions among genes and their products compose networks that can be
used to infer gene functions. Most previous studies adopt a linkage assumption,
i.e., they assume that gene interactions indicate functional similarities
between connected genes. In this study, we propose to use a gene's context
graph, i.e., the gene interaction network associated with the focal gene, to
infer its functions. In a kernel-based machine-learning framework, we design a
context graph kernel to capture the information in context graphs. Our
experimental study on a testbed of p53-related genes demonstrates the advantage
of using indirect gene interactions and shows the empirical superiority of the
proposed approach over linkage-assumption-based methods, such as the algorithm
to minimize inconsistent connected genes and diffusion kernels.


  Gene/protein interactions provide critical information for a thorough
understanding of cellular processes. Recently, considerable interest and effort
has been focused on the construction and analysis of genome-wide gene networks.
The large body of biomedical literature is an important source of gene/protein
interaction information. Recent advances in text mining tools have made it
possible to automatically extract such documented interactions from free-text
literature. In this paper, we propose a comprehensive framework for
constructing and analyzing large-scale gene functional networks based on the
gene/protein interactions extracted from biomedical literature repositories
using text mining tools. Our proposed framework consists of analyses of the
network topology, network topology-gene function relationship, and temporal
network evolution to distill valuable information embedded in the gene
functional interactions in literature. We demonstrate the application of the
proposed framework using a testbed of P53-related PubMed abstracts, which shows
that literature-based P53 networks exhibit small-world and scale-free
properties. We also found that high degree genes in the literature-based
networks have a high probability of appearing in the manually curated database
and genes in the same pathway tend to form local clusters in our
literature-based networks. Temporal analysis showed that genes interacting with
many other genes tend to be involved in a large number of newly discovered
interactions.


  We present a three-stage probabilistic model for the progression of Chronic
Myeloid Leukemia (CML), as manifested by the leukemic stem cells, progenitor
cells and mature leukemic cells. This progression is captured through the
process of cell division and cell mutation, with probabilities of occurrence
being assigned to both of them. The key contributions of this study include,
the determination of the expected number of the leukemic stem cells, progenitor
cells, mature leukemic cells, as well as total number of these cells (in terms
of probabilities, and contingent on the initial cell count), expected time to
reach a threshold level of total and injurious leukemic cells, as well as the
critical time when the disease changes its phases, the probability of
extinction of CML, and the dynamics of CML evolution consequent to primary
therapy. Finally, the results obtained are demonstrated with numerical
illustrations.


  Structured populations are ubiquitous across the biological sciences.
Mathematical models of these populations allow us to understand how individual
physiological traits drive the overall dynamics in aggregate. For example,
linear age- or age-and-size-structured models establish constraints on
individual growth under which the age- or age-and-size- distribution
stabilizes, even as the population continues to grow without bound. However,
individuals in real-world populations exhibit far more structural features than
simply age and size. Notably, cyanobacteria contain carboxysome organelles
which are central to carbon fixation and can be older (if inherited from parent
cells) or younger (if created after division) than the enveloping cell.
Motivated by a desire to understand how carboxysome age impacts growth at the
colony level, we develop a multi-structured model which allows for an arbitrary
(but finite) number of structure variables. We then derive and solve the
renewal equation for cell division to obtain an asymptotic solution, and show
that, under certain conditions, a stable structural distribution is reached.


  Identifying and characterizing mutational paths is an important issue in
evolutionary biology and in bioengineering. We here introduce a generic
description of mutational paths in terms of the goodness of sequences and of
the mutational dynamics (how sequences change) along the path. We first propose
an algorithm to sample mutational paths, which we benchmark on exactly solvable
models of proteins in silico, and apply to data-driven models of natural
proteins learned from sequence data with Restricted Boltzmann Machines. We then
use mean-field theory to characterize the properties of mutational paths for
different mutational dynamics of interest, and show how it can be used to
extend Kimura's estimate of evolutionary distances to sequence-based epistatic
models of selection.


  Generative models for structure-based molecular design hold significant
promise for drug discovery, with the potential to speed up the hit-to-lead
development cycle, while improving the quality of drug candidates and reducing
costs. Data sparsity and bias are, however, two main roadblocks to the
development of 3D-aware models. Here we propose a first-in-kind training
protocol based on multi-level contrastive learning for improved bias control
and data efficiency. The framework leverages the large data resources available
for 2D generative modelling with datasets of ligand-protein complexes. The
result are hierarchical generative models that are topologically unbiased,
explainable and customizable. We show how, by deconvolving the generative
posterior into chemical, topological and structural context factors, we not
only avoid common pitfalls in the design and evaluation of generative models,
but furthermore gain detailed insight into the generative process itself. This
improved transparency significantly aids method development, besides allowing
fine-grained control over novelty vs familiarity.


  The linear sequence of amino acids determines protein structure and function.
Protein design, known as the inverse of protein structure prediction, aims to
obtain a novel protein sequence that will fold into the defined structure.
Recent works on computational protein design have studied designing sequences
for the desired backbone structure with local positional information and
achieved competitive performance. However, similar local environments in
different backbone structures may result in different amino acids, indicating
that protein structure's global context matters. Thus, we propose the
Global-Context Aware generative de novo protein design method (GCA), consisting
of local and global modules. While local modules focus on relationships between
neighbor amino acids, global modules explicitly capture non-local contexts.
Experimental results demonstrate that the proposed GCA method outperforms
state-of-the-arts on de novo protein design. Our code and pretrained model will
be released.


  We consider a general class of birth-and-death processes with state space
$\{0,1,2,3,\ldots\}$ which describes the size of a population going eventually
to extinction with probability one. We obtain the complete spectrum of the
generator of the process killed at $0$ in the large population limit, that is,
we scale the process by a parameter $K$, and take the limit $K\to+\infty$. We
assume that the differential equation $\mathrm{d} x/\mathrm{d} t=b(x)-d(x)$
describing the infinite population limit (in any finite-time interval) has a
repulsive fixed point at $0$, and an attractive fixed point $x_*>0$. We prove
that, asymptotically, the spectrum is the superposition of two spectra. One is
the spectrum of the generator of an Ornstein-Uhlenbeck process, which is
$n(b'(x_*)-d'(x_*))$, $n\ge 0$. The other one is the spectrum of a
continuous-time binary branching process conditioned on non-extinction, and is
given by $n(d'(0)-b'(0))$, $n\ge 1$. A major difficulty is that different
scales and function spaces are involved. We work at the level of the
eigenfunctions that we split over different regions, and study their asymptotic
dependence on $K$ in each region. In particular, we prove that the spectral gap
goes to $\min\big\{b'(0)-d'(0),\,d'(x_*)-b'(x_*)\big\}$. This work complements
a previous work of ours in which we studied the approximation of the
quasi-stationary distribution and of the mean time to extinction.


  Reputations provide a powerful mechanism to sustain cooperation, as
individuals cooperate with those of good social standing. But how should moral
reputations be updated as we observe social behavior, and when will a
population converge on a common norm of moral assessment? Here we develop a
mathematical model of cooperation conditioned on reputations, for a population
that is stratified into groups. Each group may subscribe to a different social
norm for assessing reputations, and so norms compete as individuals choose to
move from one group to another. We show that a group initially comprising a
minority of the population may nonetheless overtake the entire
population--especially if it adopts the Stern Judging norm, which assigns a bad
reputation to individuals who cooperate with those of bad standing. When
individuals do not change group membership, stratifying reputation information
into groups tends to destabilize cooperation, unless individuals are strongly
insular and favor in-group social interactions. We discuss the implications of
our results for the structure of information flow in a population and the
evolution of social norms of moral judgment.


  The Covid-19 outbreak of 2020 has required many governments to develop and
adopt mathematical-statistical models of the pandemic for policy and planning
purposes. To this end, this work provides a tutorial on building a
compartmental model using Susceptible, Exposed, Infected, Recovered, Deaths and
Vaccinated (SEIRDV) status through time. The proposed model uses interventions
to quantify the impact of various government attempts made to slow the spread
of the virus. Furthermore, a vaccination parameter is also incorporated in the
model, which is inactive until the time the vaccine is deployed. A Bayesian
framework is utilized to perform both parameter estimation and prediction.
Predictions are made to determine when the peak Active Infections occur. We
provide inferential frameworks for assessing the effects of government
interventions on the dynamic progression of the pandemic, including the impact
of vaccination. The proposed model also allows for quantification of number of
excess deaths averted over the study period due to vaccination.


  Objective: To predict structure and function of Glutamine synthetase (GS)
from Pseudoalteromonas sp. by bioinformatics technology, and to provide a
theoretical basis for further study. Methods: Open reading frame (ORF) of GS
sequence from Pseudoalteromonas sp. was obtained by ORF finder and was
translated into amino acid residue. The structure domain was analyzed by Blast.
By the method of analysis tools: Protparam, ProtScale, SignalP-4.0, TMHMM,
SOPMA, SWISS-MODEL, NCBI SMART-BLAST and MAGA 7.0, the structure and function
of the protein were predicted and analyzed. Results: The results showed that
the sequence was GS with 468 amino acid residues, theoretical molecular weight
was 51986.64 Da. The protein has the closest evolutionary status with
Shewanella oneidensis. Then it had no signal peptide site and transmembrane
domain. Secondary structure of GS contained 35.04% alpha-helix, 16.67% Extended
chain, 5.34% beta-turn, 42.95% RandomCoil. Conclusions: This GU was a variety
of biological functions of protein that may be used as a molecular samples of
microbial nitrogen metabolism in extreme environments.


  We derive rigorous estimates on the speed of invasion of an advantageous
trait in a spatially advancing population in the context of a system of
one-dimensional F-KPP equations. The model was introduced and studied
heuristically and numerically in a paper by Venegas-Ortiz et al. In that paper,
it was noted that the speed of invasion by the mutant trait is faster faster
when the resident population ist expanding in space compared to the speed when
the resident population is already present everywhere. We use the Feynman-Kac
representation to provide rigorous estimates that confirm these predictions.


  We study the following model for a diploid population of constant size $N$:
Every individual carries a random number of (genetic) elements. Upon a
reproduction event each of the two parents passes each element independently
with probability $\tfrac 12$ on to the offspring. We study the process $X^N =
(X^N(1), X^N(2),...)$, where $X_t^N(k)$ is the frequency of individuals at time
$t$ that carry $k$ elements, and prove convergence (in some weak sense) of
$X^N$ jointly with its empirical first moment $Z^N$ to the ``slow-fast'' system
$(Z,X)$, where $X_t = \text{Poi}(Z_t)$ and $Z$ evolves according to a critical
Feller branching process. We discuss heuristics explaining this finding and
some extensions and limitations.


  We propose a Markovian stochastic approach to model the spread of a
SARS-CoV-2-like infection within a closed group of humans. The model takes the
form of a Partially Observable Markov Decision Process (POMDP), whose states
are given by the number of subjects in different health conditions. The model
also exposes the different parameters that have an impact on the spread of the
disease and the various decision variables that can be used to control it (e.g,
social distancing, number of tests administered to single out infected
subjects). The model describes the stochastic phenomena that underlie the
spread of the epidemic and captures, in the form of deterministic parameters,
some fundamental limitations in the availability of resources (hospital beds
and test swabs). The model lends itself to different uses. For a given control
policy, it is possible to verify if it satisfies an analytical property on the
stochastic evolution of the state (e.g., to compute probability that the
hospital beds will reach a fill level, or that a specified percentage of the
population will die). If the control policy is not given, it is possible to
apply POMDP techniques to identify an optimal control policy that fulfils some
specified probabilistic goals. Whilst the paper primarily aims at the model
description, we show with numeric examples some of its potential applications.


  Seagrasses are key structural elements in coastal ecosystems, and studying
how temperature affects these species is crucial to anticipate the implications
of global warming. In this work, we use an empirically-based numerical model to
study the combined dynamics of {\it Posidonia oceanica} and {\it Cymodocea
nodosa} and their resilience to sea warming. The model is parametrised using
seagrass growth rates measured at the Western Mediterranean Sea. Under
favorable growth conditions, our simulations predict the emergence of a
coexistence region at the front between mono-specific meadows. This region can
be characterised by its width and local shoot densities, which are found to
depend on the coupling parameter between {\it Posidonia oceanica} and {\it
Cymodocea nodosa} species. Such regions have been empirically observed in Ses
Olles de Son Saura (Balearic Islands, Western Mediterranean Sea). A comparison
between the field measurements at the study site with the model predictions has
been used to fit the value of the coupling parameter. Field data also relates
the width of the coexistence region to the average length of {\it Posidonia
oceanica} leaves at the front. Remarkably, a linear relationship is found
between the coupling parameter and the leaf length. In the presence of sea
warming, the model predicts an exponential decay in the population of {\it
Posidonia oceanica}, which is highly sensitive to temperature. This behaviour
is a direct consequence of the clonal nature of the plant and can be
characterised by the model parameters. Considering a scenario of high
greenhouse emissions, our model forecasts that {\it Posidonia oceanica} meadows
will lose 70$\%$ of their population by the year $2050$. {\it Cymodocea
nodosa}, with higher thermal resilience, acts as an opportunistic species
conquering the space left by the degraded {\it Posidonia oceanica}.


  With the present highly infectious dominant SARS-CoV-2 strain of B1.1.529 or
Omicron spreading around the globe, there is concern that the COVID-19 pandemic
will not end soon and that it will be a race against time until a more
contagious and virulent variant emerges. One of the most promising approaches
for preventing virus propagation is to maintain continuous high vaccination
efficacy among the population, thereby strengthening the population protective
effect and preventing the majority of infection in the vaccinated population,
as is known to occur with the Omicron variant frequently. Countries must
structure vaccination programs in accordance with their populations'
susceptibility to infection, optimizing vaccination efforts by delivering
vaccines progressively enough to protect the majority of the population. We
present a feasibility study proposal for maintaining optimal continuous
vaccination by assessing the susceptible population, the decline of vaccine
efficacy in the population, and advising booster dosage deployment to maintain
the population's protective efficacy through the use of a predictive model.
Numerous studies have been conducted in the direction of analyzing vaccine
utilization; however, very little study has been conducted to substantiate the
optimal deployment of booster dosage vaccination with the help of a predictive
model based on machine learning algorithms.


  We investigate how initial and boundary conditions influence the competition
dynamics and outcome in dispersal-structured populations. The study is carried
out through numerical modeling of the heterogeneous Brownian bugs model, in
which individuals are characterized by diversified diffusion coefficients and
compete for resources within a finite interaction range. We observed that if
instead of being distributed randomly across the domain the population is
initially localized within a small region, the dynamics and the stationary
state of the system differ significantly. The outcome of the natural selection
is determined by different competition processes, fluctuations, and spreading
of the organisms.


  The function of a protein is defined by its interaction partners. Thus,
topology-driven network alignment of the protein-protein interaction (PPI)
networks of two species should uncover similar interaction patterns and allow
identification of functionally similar proteins. Howver, few of the fifty or
more algorithms for PPI network alignment have demonstrated a significant link
between network topology and functional similarity, and none have recovered
orthologs using network topology alone. We find that the major contributing
factors to this failure are: (i) edge densities in current PPI networks are too
low to expect topological network alignment to succeed; (ii) when edge
densities are high enough, some measures of topological similarity easily
uncover functionally similar proteins while others do not; and (iii) most
network alignment algorithms fail to optimize their own topological objective
functions, hampering their ability to use topology effectively. We demonstrate
that SANA-the Simulated Annealing Network Aligner-significantly outperforms
existing aligners at optimizing their own objective functions, even achieving
near-optimal solutions when optimal solution is known. We offer the first
demonstration of global network alignments based on topology alone that align
functionally similar proteins with p-values in some cases below 1e-300. We
predict that topological network alignment has a bright future as edge
densities increase towards the value where good alignments become possible. We
demonstrate that when enough common topology is present at high enough edge
densities-for example in the recent, partly synthetic networks of the
Integrated Interaction Database-topological network alignment easily recovers
most orthologs, paving the way towards high-throughput functional prediction
based on topology-driven network alignment.


  Decent social fairness is highly desired both for socio-economic activities
and individuals, as it is one of the cornerstones of our social welfare and
sustainability. How to effectively promote the level of fairness thus becomes a
significant issue to be addressed. Here, by adopting a pinning control
procedure, we find that when a very small fraction of individuals are pinned to
be fair players in the Ultimatum Game, the whole population unexpectedly
evolves into the full fairness level. The basic observations are quite robust
in homogeneous networks, but the converging time as a function of the pinning
number shows different laws for different underlying topologies. For
heterogeneous networks, this leverage effect is even more pronounced that one
hub node is sufficient for the aim, and a periodic on-off control procedure can
be applied to further save the control cost. Intermittent failures are seen
when the pinning control is marginally strong, our statistical analysis
indicates some sort of criticality. Our work suggests that the pinning control
procedure could potentially be a good strategy to promote the social fairness
for some real scenarios when necessary.


  In a highly simplified view, a disease can be seen as the phenotype emerging
from the interplay of genetic predisposition and fluctuating environmental
stimuli. We formalize this situation in a minimal model, where a network
(representing cellular regulation) serves as an interface between an input
layer (representing environment) and an output layer (representing functional
phenotype). Genetic predisposition for a disease is represented as a loss of
function of some network nodes. Reduced, but non-zero, output indicates
disease. The simplicity of this genetic disease model and its deep relationship
to percolation theory allows us to understand the interplay between disease,
network topology and the location and clusters of affected network nodes. We
find that our model generates two different characteristics of diseases, which
can be interpreted as chronic and acute diseases. In its stylized form, our
model provides a new view on the relationship between genetic mutations and the
type and severity of a disease.


  Microbes have a profound impact on our health and environment, but our
understanding of the diversity and function of microbial communities is
severely limited. Through DNA sequencing of microbial communities
(metagenomics), DNA fragments (reads) of the individual microbes can be
obtained, which through assembly graphs can be combined into long contiguous
DNA sequences (contigs). Given the complexity of microbial communities, single
contig microbial genomes are rarely obtained. Instead, contigs are eventually
clustered into bins, with each bin ideally making up a full genome. This
process is referred to as metagenomic binning.
  Current state-of-the-art techniques for metagenomic binning rely only on the
local features for the individual contigs. These techniques therefore fail to
exploit the similarities between contigs as encoded by the assembly graph, in
which the contigs are organized. In this paper, we propose to use Graph Neural
Networks (GNNs) to leverage the assembly graph when learning contig
representations for metagenomic binning. Our method, VaeG-Bin, combines
variational autoencoders for learning latent representations of the individual
contigs, with GNNs for refining these representations by taking into account
the neighborhood structure of the contigs in the assembly graph. We explore
several types of GNNs and demonstrate that VaeG-Bin recovers more high-quality
genomes than other state-of-the-art binners on both simulated and real-world
datasets.


  An epidemiological study of deaths is carried out in a dozen countries by
analyzing the first wave of the COVID-19 pandemic. These countries are among
those most affected by the first wave, i.e. where daily-death data series may
closely resemble a solution of the basic SIR equations. The SIR equations are
solved parametrically using the proper time as parameter. Some general
properties of the SIR solutions are studied such as time-scaling and asymmetry.
Additionally, we use approximations to the SIR solutions through Gumbel
functions, which present a very similar behavior. The parameters of the SIR
model and the Gumbel function are extracted from the data and compared for the
different countries. It is found that ten of the selected countries are very
well described by the solutions of the SIR model, with a basic reproduction
number between 3 and 8.


  Due to the lack of a method to efficiently represent the multimodal
information of a protein, including its structure and sequence information,
predicting compound-protein binding affinity (CPA) still suffers from low
accuracy when applying machine learning methods. To overcome this limitation,
in a novel end-to-end architecture (named FeatNN), we develop a coevolutionary
strategy to jointly represent the structure and sequence features of proteins
and ultimately optimize the mathematical models for predicting CPA.
Furthermore, from the perspective of data-driven approach, we proposed a
rational method that can utilize both high- and low-quality databases to
optimize the accuracy and generalization ability of FeatNN in CPA prediction
tasks. Notably, we visually interpret the feature interaction process between
sequence and structure in the rationally designed architecture. As a result,
FeatNN considerably outperforms the state-of-the-art (SOTA) baseline in virtual
drug screening tasks, indicating the feasibility of this approach for practical
use. FeatNN provides an outstanding method for higher CPA prediction accuracy
and better generalization ability by efficiently representing multimodal
information of proteins via a coevolutionary strategy.


  The COVID-19 pandemic has presented many challenges that have spurred
biotechnological research to address specific problems. Diagnostics is one area
where biotechnology has been critical. Diagnostic tests play a vital role in
managing a viral threat by facilitating the detection of infected and/or
recovered individuals. From the perspective of what information is provided,
these tests fall into two major categories, molecular and serological.
Molecular diagnostic techniques assay whether a virus is present in a
biological sample, thus making it possible to identify individuals who are
currently infected. Additionally, when the immune system is exposed to a virus,
it responds by producing antibodies specific to the virus. Serological tests
make it possible to identify individuals who have mounted an immune response to
a virus of interest and therefore facilitate the identification of individuals
who have previously encountered the virus. These two categories of tests
provide different perspectives valuable to understanding the spread of
SARS-CoV-2. Within these categories, different biotechnological approaches
offer specific advantages and disadvantages. Here we review the categories of
tests developed for the detection of the SARS-CoV-2 virus or antibodies against
SARS-CoV-2 and discuss the role of diagnostics in the COVID-19 pandemic.


  Topological network alignment aims to align two networks node-wise in order
to maximize the observed common connection (edge) topology between them. The
topological alignment of two Protein-Protein Interaction (PPI) networks should
thus expose protein pairs with similar interaction partners allowing, for
example, the prediction of common Gene Ontology (GO) terms. Unfortunately, no
network alignment algorithm based on topology alone has been able to achieve
this aim, though those that include sequence similarity have seen some success.
We argue that this failure of topology alone is due to the sparsity and
incompleteness of the PPI network data of almost all species, which provides
the network topology with a small signal-to-noise ratio that is effectively
swamped when sequence information is added to the mix. Here we show that the
weak signal can be detected using multiple stochastic samples of "good"
topological network alignments, which allows us to observe regions of the two
networks that are robustly aligned across multiple samples. The resulting
Network Alignment Frequency (NAF) strongly correlates with GO-based Resnik
semantic similarity and enables the first successful cross-species predictions
of GO terms based on topology-only network alignments. Our best predictions
have an AUPR of about 0.4, which is competitive with state-of-the-art
algorithms, even when there is no observable sequence similarity and no known
homology relationship. While our results provide only a "proof of concept" on
existing network data, we hypothesize that predicting GO terms from
topology-only network alignments will become increasingly practical as the
volume and quality of PPI network data increase.