  Motivation: A pan-genome graph represents a collection of genomes and encodes
sequence variations between them. It is a powerful data structure for studying
multiple similar genomes. Sequence-to-graph alignment is an essential step for
the construction and the analysis of pan-genome graphs. However, existing
algorithms incur runtime proportional to the product of sequence length and
graph size, making them inefficient for aligning long sequences against large
graphs. Results: We propose the graph wavefront alignment algorithm (Gwfa), a
new method for aligning a sequence to a sequence graph. Although the worst-case
time complexity of Gwfa is the same as the existing algorithms, it is designed
to run faster for closely matching sequences, and its runtime in practice often
increases only moderately with the edit distance of the optimal alignment. On
four real datasets, Gwfa is up to four orders of magnitude faster than other
exact sequence-to-graph alignment algorithms. We also propose a graph pruning
heuristic on top of Gwfa, which can achieve an additional $\sim$10-fold speedup
on large graphs. Availability: Gwfa code is accessible at
https://github.com/lh3/gwfa.


  Compartmental models are an important quantitative tool in epidemiology,
enabling us to forecast the course of a communicable disease. However, the
model parameters, such as the infectivity rate of the disease, are riddled with
uncertainties, which has motivated the development and use of stochastic
compartmental models. Here, we first show that a common stochastic model, which
treats the uncertainties as white noise, is fundamentally flawed since it
erroneously implies that greater parameter uncertainties will lead to the
eradication of the disease. Then, we present a principled modeling of the
uncertainties based on reasonable assumptions on the contacts of each
individual. Using the central limit theorem and Doob's theorem on Gaussian
Markov processes, we prove that the correlated Ornstein-Uhlenbeck process is
the appropriate tool for modeling uncertainties in the infectivity rate. We
demonstrate our results using a compartmental model of the COVID-19 pandemic
and the available US data from the Johns Hopkins University COVID-19 database.
In particular, we show that the white noise stochastic model systematically
underestimates the severity of the Omicron variant of COVID-19, whereas the
Ornstein-Uhlenbeck model correctly forecasts the course of this variant.
Moreover, using an SIS model of sexually transmitted disease, we derive an
exact closed-form solution for the asymptotic distribution of infected
individuals. This analytic result shows that the white noise model
underestimates the severity of the pandemic because of unrealistic
noise-induced transitions. Our results strongly support the need for temporal
correlations in modeling of uncertainties in compartmental models of infectious
disease.


  Many countries have experienced at least two waves of the COVID-19 pandemic.
The second wave is far more dangerous as distinct strains appear more harmful
to human health, but it stems from the complacency about the first wave. This
paper introduces a parsimonious yet representative stochastic epidemic model
that simulates the uncertain spread of the disease regardless of the latency
and recovery time distributions. We also propose a Markov decision process to
seek an optimal trade-off between the usage of the healthcare system and the
economic costs of an epidemic. We apply the model to COVID-19 data from New
Delhi, India and simulate the epidemic spread with different policy review
times. The results show that the optimal policy acts swiftly to curb the
epidemic in the first wave, thus avoiding the collapse of the healthcare system
and the future costs of posterior outbreaks. An analysis of the recent collapse
of the healthcare system of India during the second COVID-19 wave suggests that
many lives could have been preserved if swift mitigation was promoted after the
first wave.


  In this communication we draw attention on serious flaws that plague recently
reported antioxidant properties of atorvastatin (ATV) in methanol. First and
foremost, we emphasize that the O-H bond dissociation energies (BDE) of about
400\,kcal/mol previously reported are completely wrong. Further, we present
results refuting the previous claim that the proton affinity (PA) of ATV is
smaller than that of the ascorbic acid. That unfounded claim relies on
  incorrect data for PA's ascorbic acid (which we correct here) circulated in
the literature. Further, we correct the values of the chemical reactivity
indices (e.g., chemical hardness, electrophilicity index, electroaccepting and
electrodonating powers), which were inadequately estimated previously via
Kohn-Sham HOMO and LUMO energies. Finally, our updated values for \ce{O-H} bond
dissociation enthalpy (BDE = 91.4\,kcal/mol) and electron transfer enthalpy
(ETE = 105.7\,kcal/mol) tentatively suggest that direct H-atom transfer (HAT)
and sequential proton loss electron transfer (SPLET) may coexist.


  It is difficult to observe a spontaneous translocation of cell-penetrating
peptides(CPPs) within a short time scale (e.g., a few hundred ns) in all-atom
molecular dynamics(MD) simulations because the time required for the
translocation of usual CPPs is on the order of minutes or so. In this work, we
report a spontaneous translocation of a single Arg$_9$(R9) across a
DOPC/DOPG(4:1) model membrane within an order of a few tens ns scale by using
the weighted ensemble(WE) method. We identify how water molecules and the
orientation of Arg$_9$ play a role in translocation. We also show how lipid
molecules are transported along with Arg$_9$. In addition, we present free
energy profiles of the translocation across the membrane using umbrella
sampling and show that a single Arg$_9$ translocation is energetically
unfavorable. We expect that the WE method can help study interactions of CPPs
with various model membranes within MD simulation approaches.


  A central challenge in every field of biology is to use existing measurements
to predict the outcomes of future experiments. In this work, we consider the
wealth of antibody inhibition data against variants of the influenza virus. Due
to this viru's genetic diversity and evolvability, the variants examined in one
study will often have little-to-no overlap with other studies, making it
difficult to discern common patterns or unify datasets for further analysis. To
that end, we develop a computational framework that predicts how an antibody or
serum would inhibit any variant from any other study. We use this framework to
greatly expand seven influenza datasets utilizing hemagglutination inhibition,
validating our method upon 200,000 existing measurements and predicting
2,000,000 new values along with their uncertainties. With these new values, we
quantify the transferability between seven vaccination and infection studies in
humans and ferrets, show that the serum potency is negatively correlated with
breadth, and present a tool for pandemic preparedness. This data-driven
approach does not require any information beyond each virus's name and
measurements, and even datasets with as few as 5 viruses can be expanded,
making this approach widely applicable. Future influenza studies using
hemagglutination inhibition can directly utilize our curated datasets to
predict newly measured antibody responses against ~80 H3N2 influenza viruses
from 1968-2011, whereas immunological studies utilizing other viruses or a
different assay only need a single partially-overlapping dataset to extend
their work. In essence, this approach enables a shift in perspective when
analyzing data from "what you see is what you get" into "what anyone sees is
what everyone gets."


  We propose an application of molecular information theory to analyze the
folding of single domain proteins. We analyze results from various areas of
protein science, such as sequence-based potentials, reduced amino acid
alphabets, backbone configurational entropy, secondary structure content,
residue burial layers, and mutational studies of protein stability changes. We
found that the average information contained in the sequences of evolved
proteins is very close to the average information needed to specify a fold ~2.2
$\pm$ 0.3 bits/(site operation). The effective alphabet size in evolved
proteins equals the effective number of conformations of a residue in the
compact unfolded state at around 5. We calculated an energy-to-information
conversion efficiency upon folding of around 50%, lower than the theoretical
limit of 70%, but much higher than human built macroscopic machines. We propose
a simple mapping between molecular information theory and energy landscape
theory and explore the connections between sequence evolution, configurational
entropy and the energetics of protein folding.


  Predicting the consensus structure of a set of aligned RNA homologs is a
convenient method to find conserved structures in an RNA genome, which has
applications in SARS-CoV-2 diagnostics and therapeutics. However, the
state-of-the-art algorithm for this task, RNAalifold, is prohibitively slow for
long sequences, due to a cubic scaling with the sequence length, and even
slower when analyzing many such sequences, due to a superlinear scaling with
the number of homologs, taking 4 days on 200 SARS-CoV variants. We present
LinearAlifold, an efficient algorithm for folding aligned RNA homologs that
scales linearly with both the sequence length and the number of sequences,
based on our recent work LinearFold that folds a single RNA in linear time. Our
work is orders of magnitude faster than RNAalifold (e.g., 0.5 hours on the
above 200 sequences or 316 times speedup) and achieves comparable accuracies
compared to a database of known structures. More interestingly, LinearAlifold's
prediction on SARS-CoV-2 correlates well with experimentally determined
structures, outperforming RNAalifold. Finally, LinearAlifold supports three
modes: minimum free energy (MFE), partition function, and stochastic sampling,
each of which takes under an hour for hundreds of SARS-CoV variants, while only
the MFE mode of RNAalifold works for them, taking days or weeks.


  Mixed microbial communities, usually composed of various bacterial and fungal
species, are fundamental in a plethora of environments, from soil to human gut
and skin. Their evolution is a paradigmatic example of intertwined dynamics,
where not just the relations among species plays a role, but also the
opportunities -- and possible harms -- that each species presents to the
others. These opportunities are in fact \textit{affordances}, which can be
seized by heritable variation and selection. In this paper, starting from a
systemic viewpoint of mixed microbial communities, we focus on the pivotal role
of affordances in evolution and we contrast it to the artificial evolution of
programs and robots. We maintain that the two realms are neatly separated, in
that natural evolution proceeds by extending the space of its possibilities in
a completely open way, while the latter is inherently limited by the
algorithmic framework it is defined. This discrepancy characterises also an
envisioned setting in which robots evolve in the physical world. We present
arguments supporting our claim and we propose an experimental setting for
assessing our statements. Rather than just discussing the limitations of the
artificial evolution of machines, the aim of this contribution is to emphasize
the tremendous potential of the evolution of the biosphere, beautifully
represented by the evolution of communities of microbes.


  Tracing potentially infected contacts of confirmed cases is important when
fighting outbreaks of many infectious diseases. The COVID-19 pandemic has
motivated researchers to examine how different contact tracing strategies
compare in terms of effectiveness (ability to mitigate infections) and cost
efficiency (number of prevented infections per isolation). Two important
strategies are so-called forward contact tracing (tracing to whom disease
spreads) and backward contact tracing (tracing from whom disease spreads).
Recently, Kojaku and colleagues reported that backward contact tracing was
``profoundly more effective'' than forward contact tracing, that contact
tracing effectiveness ``hinges on reaching the `source' of infection'', and
that contact tracing outperformed case isolation in terms of cost efficiency.
Here we show that these conclusions are not true in general. They were based in
part on simulations that vastly overestimated the effectiveness and efficiency
of contact tracing. Our results show that the efficiency of contact tracing
strategies is highly contextual; faced with a disease outbreak, the disease
dynamics determine whether tracing infection sources or new cases is more
impactful. Our results also demonstrate the importance of simulating disease
spread and mitigation measures in parallel rather than sequentially.


  Co-evolutionary algorithms have a wide range of applications, such as in
hardware design, evolution of strategies for board games, and patching software
bugs. However, these algorithms are poorly understood and applications are
often limited by pathological behaviour, such as loss of gradient, relative
over-generalisation, and mediocre objective stasis. It is an open challenge to
develop a theory that can predict when co-evolutionary algorithms find
solutions efficiently and reliable.
  This paper provides a first step in developing runtime analysis for
population-based competitive co-evolutionary algorithms. We provide a
mathematical framework for describing and reasoning about the performance of
co-evolutionary processes. An example application of the framework shows a
scenario where a simple co-evolutionary algorithm obtains a solution in
polynomial expected time. Finally, we describe settings where the
co-evolutionary algorithm needs exponential time with overwhelmingly high
probability to obtain a solution.


  In this work, we present a quantum circuit model for inferring gene
regulatory networks (GRNs). The model is based on the idea of using qubit-qubit
entanglement to simulate interactions between genes. We provide preliminary
results that suggest our quantum GRN modeling method is competitive and
warrants further investigation. Specifically, we present the results derived
from the single-cell transcriptomic data of human cell lines, focusing on genes
in involving innate immunity regulation. We demonstrate that our quantum
circuit model can be used to predict the presence or absence of regulatory
interactions between genes and estimate the strength and direction of the
interactions, setting the stage for further investigations on how quantum
computing finds applications in data-driven life sciences and, more
importantly, to invite exploration of quantum algorithm design that takes
advantage of the single-cell data. The application of quantum computing on
single-cell transcriptomic data likewise contributes to a novel understanding
of GRNs, given that the relationship between fully interconnected genes can be
approached more effectively by quantum modeling than by statistical
correlations.


  The existence of an exponential growth phase during early stages of a
pandemic is often taken for granted. However, for the 2019 novel coronavirus
epidemic, the early exponential phase lasted only for about six days, while the
quadratic growth prevailed for forty days until it spread to other countries
and continued, again quadratically, but with a larger coefficient. Here we show
that this rapid phase is followed by a subsequent slow-down where the
coefficient is reduced to almost the original value at the outbreak. This can
be explained by the merging of previously disconnected sites that occurred
after the disease jumped (nonlocally) to a relatively small number of separated
sites. Subsequent variations in the slope with continued growth can
qualitatively be explained as a result of reinfections and changes in their
rate. We demonstrate that the observed behavior can be described by a standard
epidemiological model with spatial extent and reinfections included.
Time-dependent changes in the spatial diffusion coefficient can also model
corresponding variations in the slope.


  Given a rooted, binary phylogenetic network and a rooted, binary phylogenetic
tree, can the tree be embedded into the network? This problem, called
\textsc{Tree Containment}, arises when validating networks constructed by
phylogenetic inference methods.We present the first algorithm for (rooted)
\textsc{Tree Containment} using the treewidth $t$ of the input network $N$ as
parameter, showing that the problem can be solved in $2^{O(t^2)}\cdot|N|$ time
and space.


  The origin of life and the detection of alien life have historically been
treated as separate scientific research problems. However, they are not
strictly independent. Here, we discuss the need for a better integration of the
sciences of life detection and origins of life. Framing these dual problems
within the formalism of Bayesian hypothesis testing, we show via simple
examples how high confidence in life detection claims require either (1) a
strong prior hypothesis about the existence of life in a particular alien
environment, or conversely, (2) signatures of life that are not susceptible to
false positives. As a case study, we discuss the role of priors and hypothesis
testing in recent results reporting potential detection of life in the Venusian
atmosphere and in the icy plumes of Enceladus. While many current leading
biosignature candidates are subject to false positives because they are not
definitive of life, our analyses demonstrate why it is necessary to shift focus
to candidate signatures that are definitive. This indicates a necessity to
develop methods that lack false positives, by using observables for life that
rely on prior hypotheses with strong theoretical and empirical support in
identifying defining features of life. Abstract theories developed in pursuit
of understanding universal features of life are more likely to be definitive
and to apply to life-as-we-don't-know-it. In the absence of alien examples
these are best validated in origin of life experiments, substantiating the need
for better integration between origins of life and biosignature science
research communities.


  The rational design of novel molecules with desired bioactivity is a critical
but challenging task in drug discovery, especially when treating a novel target
family or understudied targets. Here, we propose PGMG, a pharmacophore-guided
deep learning approach for bioactivate molecule generation. Through the
guidance of pharmacophore, PGMG provides a flexible strategy to generate
bioactive molecules with structural diversity in various scenarios using a
trained variational autoencoder. We show that PGMG can generate molecules
matching given pharmacophore models while maintaining a high level of validity,
uniqueness, and novelty. In the case studies, we demonstrate the application of
PGMG to generate bioactive molecules in ligand-based and structure-based drug
de novo design, as well as in lead optimization scenarios. Overall, the
flexibility and effectiveness of PGMG make it a useful tool for accelerating
the drug discovery process.


  Phylogenetic networks which are, as opposed to trees, suitable to describe
processes like hybridization and horizontal gene transfer, play a substantial
role in evolutionary research. However, while non-treelike events need to be
taken into account, they are relatively rare, which implies that biologically
relevant networks are often assumed to be similar to trees in the sense that
they can be obtained by taking a tree and adding some additional edges. This
observation led to the concept of so-called tree-based networks, which recently
gained substantial interest in the literature. Unfortunately, though,
identifying such networks in the unrooted case is an NP-complete problem.
Therefore, classes of networks for which tree-basedness can be guaranteed are
of the utmost interest.
  The most prominent such class is formed by so-called edge-based networks,
which have a close relationship to generalized series-parallel graphs known
from graph theory. They can be identified in linear time and are in some
regards biologically more plausible than general tree-based networks. While
concerning the latter proximity measures for general networks have already been
introduced, such measures are not yet available for edge-basedness. This means
that for an arbitrary unrooted network, the "distance" to the nearest
edge-based network could so far not be determined. The present manuscript fills
this gap by introducing two classes of proximity measures for edge-basedness,
one based on the given network itself and one based on its so-called leaf
shrink graph (LS graph). Both classes contain four different proximity
measures, whose similarities and differences we study subsequently.


  Recent developments in artificial intelligence and automation support a new
drug design paradigm: autonomous drug design. Under this paradigm, generative
models can provide suggestions on thousands of molecules with specific
properties, and automated laboratories can potentially make, test and analyze
molecules with minimal human supervision. However, since still only a limited
number of molecules can be synthesized and tested, an obvious challenge is how
to efficiently select among provided suggestions in a closed-loop system. We
formulate this task as a stochastic multi-armed bandit problem with multiple
plays, volatile arms and similarity information. To solve this task, we adapt
previous work on multi-armed bandits to this setting, and compare our solution
with random sampling, greedy selection and decaying-epsilon-greedy selection
strategies. According to our simulation results, our approach has the potential
to perform better exploration and exploitation of the chemical space for
autonomous drug design.


  It has been established that Terahertz (THz) band signals can interact with
biomolecules through resonant modes. Specifically, of interest here, protein
activation. Our research goal is to show how directing the mechanical signaling
inside protein molecules using THz signals can control changes in their
structure and activate associated biochemical and biomechanical events. To
establish that, we formulate a selectivity metric that quantifies the system
performance and captures the capability of the nanoantenna to induce a
conformational change in the desired protein molecule/population. The metric
provides a score between -1 and 1 that indicates the degree of control we have
over the system to achieve targeted protein interactions. To develop the
selectivity measure, we first use the Langevin stochastic equation driven by an
external force to model the protein behavior. We then determine the probability
of protein folding by computing the steady-state energy of the driven protein
and then generalize our model to account for protein populations. Our numerical
analysis results indicate that a maximum selectivity score is attained when
only the targeted population experiences a folding behavior due to the
impinging THz signal. From the achieved selectivity values, we conclude that
the system response not only depends on the resonant frequency but also on the
system controlling parameters namely, the nanoantenna force, the damping
constant, and the abundance of each protein population. The presented work
sheds light on the potential associated with the electromagnetic-based control
of protein networks, which could lead to a plethora of applications in the
medical field ranging from bio-sensing to targeted therapy.


  Intercity travel is one of the most important parameters for combating a
pandemic. The ongoing COVID-19 pandemic has resulted in different computational
studies involving intercity connections. In this study, the effects of
intercity connections during an epidemic such as COVID-19 are evaluated using a
new network model. This model considers the actual geographic neighborhood and
population density data. This new model is applied to actual Turkish data by
the means of provincial connections and populations. A Monte Carlo algorithm
with a hybrid lattice model is applied to a lattice with 8,802 data points.
Results show that this model is quantitatively very efficient in modeling real
world COVID-19 epidemic data based on populations and geographical intercity
connections, by the means of estimating the number of deaths, disease spread,
and epidemic termination.


  RNA structure determination and prediction can promote RNA-targeted drug
development and engineerable synthetic elements design. But due to the
intrinsic structural flexibility of RNAs, all the three mainstream structure
determination methods (X-ray crystallography, NMR, and Cryo-EM) encounter
challenges when resolving the RNA structures, which leads to the scarcity of
the resolved RNA structures. Computational prediction approaches emerge as
complementary to the experimental techniques. However, none of the \textit{de
novo} approaches is based on deep learning since too few structures are
available. Instead, most of them apply the time-consuming sampling-based
strategies, and their performance seems to hit the plateau. In this work, we
develop the first end-to-end deep learning approach, E2Efold-3D, to accurately
perform the \textit{de novo} RNA structure prediction. Several novel components
are proposed to overcome the data scarcity, such as a fully-differentiable
end-to-end pipeline, secondary structure-assisted self-distillation, and
parameter-efficient backbone formulation. Such designs are validated on the
independent, non-overlapping RNA puzzle testing dataset and reach an average
sub-4 \AA{} root-mean-square deviation, demonstrating its superior performance
compared to state-of-the-art approaches. Interestingly, it also achieves
promising results when predicting RNA complex structures, a feat that none of
the previous systems could accomplish. When E2Efold-3D is coupled with the
experimental techniques, the RNA structure prediction field can be greatly
advanced.


  A closer look at nature has recently brought more interest in exploring and
utilizing intra-body communication networks composed of cells as intrinsic,
perfectly biocompatible infrastructures to deliver therapeutics. Naturally
occurring cell-to-cell communication systems are being manipulated to release,
navigate, and take up soluble cell-derived messengers that are either
therapeutic by nature or carry therapeutic molecular cargo in their structures.
One example of such structures is extracellular vesicles (EVs) which have been
recently proven to have favorable pharmacokinetic properties, opening new
avenues for developing the next generation biotherapeutics. In this paper, we
study theoretical aspects of the EV transfer within heart tissue as a case
study by utilizing an information and communication technology-like approach in
analyzing molecular communication systems. Our modeling implies the abstraction
of the EV releasing cells as transmitters, the extracellular matrix as the
channel, and the EV receiving cells as receivers. Our results, derived from the
developed analytical models, indicate that the release can be modulated using
external forces such as electrical signals, and the transfer and reception can
be affected by the extracellular matrix and plasma membrane properties,
respectively. The results can predict the EV biodistributions and contribute to
avoiding unplanned administration, often resulting in side- and adverse
effects.


  Metapopulation models have been a popular tool for the study of epidemic
spread over a network of highly populated nodes (cities, provinces, countries)
and have been extensively used in the context of the ongoing COVID-19 pandemic.
In the present work, we revisit such a model, bearing a particular case example
in mind, namely that of the region of Andalusia in Spain during the period of
the summer-fall of 2020 (i.e., between the first and second pandemic waves).
Our aim is to consider the possibility of incorporation of mobility across the
province nodes focusing on mobile-phone time dependent data, but also
discussing the comparison for our case example with a gravity model, as well as
with the dynamics in the absence of mobility. Our main finding is that mobility
is key towards a quantitative understanding of the emergence of the second wave
of the pandemic and that the most accurate way to capture it involves dynamic
(rather than static) inclusion of time-dependent mobility matrices based on
cell-phone data. Alternatives bearing no mobility are unable to capture the
trends revealed by the data in the context of the metapopulation model
considered herein.


  Prion and prion-like molecules are a type of self replicating aggregate
protein that have been implicated in a variety of neurodegenerative diseases.
Over recent decades the molecular dynamics of prions have been characterized
both empirically and through mathematical models, providing insights into the
epidemiology of prion diseases, and the impact of prions on the evolution of
cellular processes. At the same time, a variety of evidence indicates that
prions are themselves capable of a form of evolution, in which changes to their
structure that impact their rate of growth or fragmentation are replicated,
making such changes subject to natural selection. Here we study the role of
such selection in shaping the characteristics of prions under the nucleated
polymerization model (NPM). We show that fragmentation rates evolve to an
evolutionary stable value which balances rapid reproduction of $PrP^{Sc}$
aggregates with the need to produce stable polymers. We further show that this
evolved fragmentation rate differs in general from the rate that optimizes
transmission between cells. We find that under the NPM, prions that are both
evolutionary stable and optimized for transmission have a characteristic length
of $3n$, i.e three times the critical length below which they become unstable.
Finally we study the dynamics of inter-cellular competition between strains,
and show that the eco-evolutionary tradeoff between intra- and inter-cellular
competition favors coexistence.


  Monod's model describes the growth of microorganisms using a hyperbolic
function of extracellular resource concentration. Under fluctuating or limited
resource concentrations this model performs poorly against experimental data,
motivating the more complex Droop model with a time-varying internal storage
pool. We extend the Monod model to incorporate memory of past conditions,
adding a single parameter motivated by a fractional calculus analysis. We show
how to interpret the memory element in a biological context and describe its
connection to a resource storage pool. Under nitrogen starvation at
non-equilibrium conditions, we validate the model with simulations and
empirical data obtained from lab cultures of diatoms (T. pseudonana and T.
weissflogii) and prasinophytes (Micromonas sp. and O. tauri), globally
influential phytoplankton taxa. Using statistical analysis, we show that our
Monod-memory model estimates the growth rate, cell density, and resource
concentration as well as the Droop model while requiring one less state
variable. Our simple model may improve descriptions of phytoplankton dynamics
in complex earth system models at a lower computational cost than is presently
achievable.


  We propose in this article a brief description of the work, over almost a
decade, resulting from a collaboration between mathematicians and biologists
from four different research laboratories, identifiable as the co-authors of
the articles whose results are described here, and implicitely co-authors of
this article, under the signature of Camille No\^us. This modeling work is part
of population genetics, of which we give a definition below. It is therefore
essentially at the interface between mathematical tools, more particularly
probabilistic ones, and biological data, more specifically genetic ones.


  We consider the problem of sampling transition paths between two given
metastable states of a molecular system, e.g. a folded and unfolded protein or
products and reactants of a chemical reaction. Due to the existence of high
energy barriers separating the states, these transition paths are unlikely to
be sampled with standard Molecular Dynamics (MD) simulation. Traditional
methods to augment MD with a bias potential to increase the probability of the
transition rely on a dimensionality reduction step based on Collective
Variables (CVs). Unfortunately, selecting appropriate CVs requires chemical
intuition and traditional methods are therefore not always applicable to larger
systems. Additionally, when incorrect CVs are used, the bias potential might
not be minimal and bias the system along dimensions irrelevant to the
transition. Showing a formal relation between the problem of sampling molecular
transition paths, the Schr\"odinger bridge problem and stochastic optimal
control with neural network policies, we propose a machine learning method for
sampling said transitions. Unlike previous non-machine learning approaches our
method, named PIPS, does not depend on CVs. We show that our method successful
generates low energy transitions for Alanine Dipeptide as well as the larger
Polyproline and Chignolin proteins.


  In the context of Evolutionary Game Theory, one of the most noteworthy
mechanisms to support cooperation is spatial reciprocity, usually accomplished
by distributing players in a spatial structure allowing cooperators to cluster
together and avoid exploitation. This raises an important question: how is the
survival of cooperation affected by different topologies? Here, to address this
question, we explore the Focal Public Goods (FPGG) and classic Public Goods
Games (PGG), and the Prisoner's Dilemma (PD) on several regular lattices:
honeycomb, square (with von Neumann and Moore neighborhoods), kagome,
triangular, cubic, and 4D hypercubic lattices using both analytical methods and
agent-based Monte Carlo simulations. We found that for both Public Goods Games,
a consistent trend appears on all two-dimensional lattices: as the number of
first neighbors increases, cooperation is enhanced. However, this is only
visible by analysing the results in terms of the payoff's synergistic factor
normalized by the number of connections. Besides this, clustered topologies,
i.e., those that allow two connected players to share neighbors, are the most
beneficial to cooperation for the FPGG. The same is not always true for the
classic PGG, where having shared neighbors between connected players may or may
not benefit cooperation. We also provide a reinterpretation of the classic PGG
as a focal game by representing the lattice structure of this category of games
as a single interaction game with longer-ranged, weighted neighborhoods, an
approach valid for any regular lattice topology. Finally, we show that
depending on the payoff parametrization of the PD, there can be an equivalency
between the PD and the FPGG; when the mapping between the two games is
imperfect, the definition of an effective synergy parameter can still be useful
to show their similarities.


  In this paper, we show how to modify a compartmental epidemic model, without
changing the dimension, such that separable static heterogeneity is taken into
account. The derivation is based on the Kermack-McKendrick renewal equation.


  Modeling and analysis of the large scale Covid-19 pandemic data can yield
inferences about it's dynamics and characteristics of disease propagation.
These inferences can then be correlated with contextual factors like population
density, effects of strategic interventions, heterogeneous disease propagation
etc, and such set of validated inferences can serve as precedents for designing
of subsequent mitigation strategies. In this work, we present the analysis of
Covid-19 pandemic data in Indian context using growth functions fitting
procedure and harmonic analysis method. Our results of growth function fitting
to the data indicate that the growth function parameters are quite sensitive to
the growth of the infected population indicating positive impact of lockdown
strategy, identification of inflection point and nearly synchronous statistical
features of disease spreading. The harmonic analysis of the data shows the
countrywide synchronous incident features due to simultaneous implementation of
control strategies. However, if one analyzes the data from each state of the
India, one can see various forms of travelling waves in the countrywide wave
pattern. Hence, one needs to do these analysis from time to time to understand
the effectiveness of any control strategy and to closely look at the disease
propagation to devise the required type of mitigation strategies.


  A binary phylogenetic network on a taxon set $X$ is a rooted acyclic digraph
in which the degree of each nonleaf node is three and its leaves
(i.e.degree-one nodes) are uniquely labeled with the taxa of $X$. It is
tree-child if each nonleaf node has at least one child of indegree one. A set
of binary phylogenetic trees may or may not be simultaneously displayed in a
binary tree-child network. Necessary conditions for multiple phylogenetic trees
being simultaneously displayed in a tree-child network are given here. In
particular, it is proved that any two phylogenetic trees can always
simultaneously be displayed in some tree-child network on the same taxa set. It
is also proved that any set of multiple binary phylogenetic trees can always
simultaneously be displayed in some non-binary tree-child network on the same
taxa set, where each nonleaf node is of either indegree one and outdegree two
or indegree at least two and outdegree out.


  Molecular property prediction is essential in chemistry, especially for drug
discovery applications. However, available molecular property data is often
limited, encouraging the transfer of information from related data. Transfer
learning has had a tremendous impact in fields like Computer Vision and Natural
Language Processing signaling for its potential in molecular property
prediction. We present a pre-training procedure for molecular representation
learning using reaction data and use it to pre-train a SMILES Transformer. We
fine-tune and evaluate the pre-trained model on 12 molecular property
prediction tasks from MoleculeNet within physical chemistry, biophysics, and
physiology and show a statistically significant positive effect on 5 of the 12
tasks compared to a non-pre-trained baseline model.


  We discuss various models for epidemics on networks that rely on Markov
chains. Random walks on graphs are often used to predict epidemic spread and to
investigate possible control actions to mitigate them. In this study, we
demonstrate that they do not fully reflect the dynamics of epidemics, as they
overestimate infection times. Accordingly, we explain how Markov chains may
still be used to accurately model the virus spread, and to correctly predict
infection times. We also provide an algorithm that efficiently estimates
infection times via a sampling strategy. Finally, we present a novel indicator
based on infection times, and we compare its node ranking properties with other
centrality measures based on random walks.


  The coalescent is a stochastic process representing ancestral lineages in a
population undergoing neutral genetic drift. We introduce an abstract
formulation of the coalescent that applies to a broad class of neutral drift
models, allowing for arbitrary spatial structure and mating patterns. In this
formulation, the coalescent is a random sequence of mappings $\mathcal{C} =
\left(C_t\right)_{t=0}^\infty$ from a finite set $G$ to itself, representing
the ancestry of alleles at a single genetic locus. Here, $G$ represents the
``sites'' (in individuals, at particular spatial locations) at which these
alleles can live. The state $C_t:G \to G$ maps each site $g \in G$ to the site
containing $g$'s ancestor, $t$ time steps into the past. $\mathcal{C}$ is
obtained by iterated composition of ``parentage mappings", drawn independently
from a fixed probability distribution. Mutation is incorporated by augmenting
the state to include a set of mutated sites. Using this formalism, we define
and analyze biologically relevant quantities, including coalescence time,
coalescence branch length, reproductive value, number of mutations prior to
coalescence, and stationary probabilities of identity-by-descent and
identity-by-state. The limit of low mutation is analyzed in detail. We
illustrate this coalescent with applications to well-mixed and graph-structured
populations, as well as to a diploid population with arbitrary sex ratio.


  We proposed the SIkJalpha model at the beginning of the COVID-19 pandemic
(early 2020). Since then, as the pandemic evolved, more complexities were added
to capture crucial factors and variables that can assist with projecting
desired future scenarios. Throughout the pandemic, multi-model collaborative
efforts have been organized to predict short-term outcomes (cases, deaths, and
hospitalizations) of COVID-19 and long-term scenario projections. We have been
participating in five such efforts. This paper presents the evolution of the
SIkJalpha model and its many versions that have been used to submit to these
collaborative efforts since the beginning of the pandemic. Specifically, we
show that the SIkJalpha model is an approximation of a class of epidemiological
models. We demonstrate how the model can be used to incorporate various
complexities, including under-reporting, multiple variants, waning of immunity,
and contact rates, and to generate probabilistic outputs.


  Unknown-view tomography (UVT) reconstructs a 3D density map from its 2D
projections at unknown, random orientations. A line of work starting with Kam
(1980) employs the method of moments (MoM) with rotation-invariant Fourier
features to solve UVT in the frequency domain, assuming that the orientations
are uniformly distributed. This line of work includes the recent orthogonal
matrix retrieval (OMR) approaches based on matrix factorization, which, while
elegant, either require side information about the density that is not
available, or fail to be sufficiently robust. For OMR to break free from those
restrictions, we propose to jointly recover the density map and the orthogonal
matrices by requiring that they be mutually consistent. We regularize the
resulting non-convex optimization problem by a denoised reference projection
and a nonnegativity constraint. This is enabled by the new closed-form
expressions for spatial autocorrelation features. Further, we design an
easy-to-compute initial density map which effectively mitigates the
non-convexity of the reconstruction problem. Experimental results show that the
proposed OMR with spatial consensus is more robust and performs significantly
better than the previous state-of-the-art OMR approach in the typical low-SNR
scenario of 3D UVT.


  The interplay between the virus, infected cells and the immune responses to
SARS-CoV-2 is still under debate. Extending the basic model of viral dynamics
we propose here a formal approach to describe the neutralizing versus weakly
(or non-)neutralizing scenarios and compare with the possible effects of
antibody-dependent enhancement (ADE). The theoretical model is consistent with
data available from the literature; we show that weakly neutralizing antibodies
or ADE can both give rise to either final virus clearance or disease
progression, but the immuno-dynamic is different in each case. Given that a
significant part of the world population is already naturally immunized or
vaccinated, we also discuss the implications on secondary infections infections
following vaccination or in presence of immune system dysfunctions.


  Summary: The DynaSig-ML (Dynamical Signatures - Machine Learning) Python
package allows the efficient, user-friendly exploration of 3D dynamics-function
relationships in biomolecules, using datasets of experimental measures from
large numbers of sequence variants. The DynaSig-ML package is built around the
Elastic Network Contact Model (ENCoM), the first and only sequence-sensitive
coarse-grained NMA model, which is used to generate the input Dynamical
Signatures. Starting from in silico mutated structures, the whole pipeline can
be run with just a few lines of Python and modest computational resources. The
compute-intensive steps can also easily be parallelized in the case of either
large biomolecules or vast amounts of sequence variants. As an example
application, we use the DynaSig-ML package to predict the evolutionary fitness
of the bacterial enzyme VIM-2 lactamase from deep mutational scan data.
Availability and implementation: DynaSig-ML is open source software available
at https: //github.com/gregorpatof/dynasigml package Contact:
rafael.najmanovich@umontreal.ca


  The only animals other than bats reported to have been infected with
SARS-CoV-2-related coronaviruses (SARS2r-CoVs) prior to the COVID-19 pandemic
are pangolins. In early 2020 multiple papers reported the identification of two
clades of SARS2r-CoVs, GD and GX, infecting pangolins. However the RNA-Seq
datasets supporting pangolin genome assembly were widely contaminated,
contained synthetic vectors or were heavily enriched or filtered with little
but coronavirus sequences left in the datasets. Here we investigate two
pangolin fecal samples sequenced by Li et al. (2021) provided in support of GD
PCoV infection of pangolins in Guangdong and find the read distribution
consistent with PCR amplicon contamination and SARS-CoV-2 contamination, and
further identify the presence of synthetic plasmid sequences. We also build
upon our previous work to further analyze the dataset GX/P3B by Lam et al.
(2020), which is the only non enriched/heavily filtered pangolin tissue dataset
sequenced by Lam et al. (2020). We identify synthetic vectors and confirm human
genomic origin samples in the dataset. Finally, we find human mitochondrial
sequences in all pangolin organ datasets and mouse and tiger mitochondrial
sequences in selected pangolin organ datasets sequenced by Liu et al. (2019).
We infer that human and mouse genomic origin sequences were probably sourced
from contamination prior to sequencing, while tiger origin sequence
contamination may have occurred due to index hopping during sequencing. These
observations are problematic for attributing pangolins as SARS2r-CoV hosts in
the datasets examined. The forensic methods developed and used here can be
applied to examine any third party SRA data sets.


  Multi-parameter persistent homology naturally arises in applications of
persistent topology to data that come with extra information depending on
additional parameters, like for example time series data. We introduce the
concept of a Vietoris-Rips transformation, a method that reduces the
computation of the one-parameter persistent homology of pathwise subcomplexes
in multi-filtered flag complexes to the computation of the Vietoris-Rips
persistent homology of certain semimetric spaces. The corresponding pathwise
persistence barcodes track persistence features of the ambient multi-filtered
complex and can in particular be used to recover the rank invariant in
multi-parameter persistent homology. We present MuRiT, a scalable algorithm
that computes the pathwise persistence barcodes of multi-filtered flag
complexes by means of Vietoris-Rips transformations. Moreover, we provide an
efficient software implementation of the MuRiT algorithm which resorts to
Ripser for the actual computation of Vietoris-Rips persistence barcodes. To
demonstrate the applicability of MuRiT to real-world datasets, we establish
MuRiT as part of our CoVtRec pipeline for the surveillance of the convergent
evolution of the coronavirus SARS-CoV-2 in the current COVID-19 pandemic.


  In the recent COVID-19 pandemic we assisted at a sequence of epidemic waves
intertwined by anomalous fade-outs with periods of low but persistent epidemic
prevalence. These long-living epidemic states complicate epidemic control and
challenge current modeling approaches as classical epidemic models fail to
explain their emergence. Inspired by this phenomenon, we propose a simple
mechanism able to reproduce several features observed in real data.
Specifically, here we introduce a modification of the
Susceptible-Infected-Recovered (SIR) model in a meta-population framework where
a small inflow of infected individuals accounts for undetected internal or
imported cases. Focusing on a regime where this external seeding is so small
that cannot be detected from the analysis of epidemic curves, we find that
outbreaks of finite duration percolate in time resulting in overall low but
long-living epidemic states below and above the epidemic threshold. Using a
two-state description of the local dynamics, we can extract analytical
predictions for the phase space. The comparison with epidemic data demonstrates
that our model is able to reproduce some critical signatures observed in
COVID-19 spreading in England. Finally, our findings defy our understanding of
the concept of epidemic threshold and its relationship with outbreaks survival
for disease control.


  In this work, we extend our previously developed compartmental SEIQRD model
for SARS-CoV-2 in Belgium. We introduce SARS-CoV-2 variants of concern,
vaccines, and seasonality in our model, as their addition has proven necessary
for modelling SARS-CoV-2 transmission dynamics during the 2020-2021 COVID-19
pandemic in Belgium. The model is geographically stratified into eleven spatial
patches (provinces), and a telecommunication dataset provided by Belgium's
biggest operator is used to incorporate interprovincial mobility. We calibrate
the model using the daily number of hospitalisations in each province and
serological data. We find the model adequately describes these data, but the
addition of interprovincial mobility was not necessary to obtain an accurate
description of the 2020-2021 SARS-CoV-2 pandemic in Belgium. We further
demonstrate how our model can be used to help policymakers decide on the
optimal timing of the release of social restrictions. We find that adding
spatial heterogeneity by geographically stratifying the model results in more
uncertain model projections as compared to an equivalent nation-level model,
which has both communicative advantages and disadvantages. We finally discuss
the impact of imposing local mobility or social contact restrictions to contain
an epidemic in a given province and find that lowering social contact is a more
effective strategy than lowering mobility.


  We introduce an individual-based model of a complex ecological community with
random interactions. The model contains a large number of species, each with a
finite population of individuals, subject to discrete reproduction and death
events. The interaction coefficients determining the rates of these events is
chosen from an ensemble of random matrices, and is kept fixed in time. The
set-up is such that the model reduces to the known generalised Lotka-Volterra
equations with random interaction coefficients in the limit of an infinite
population for each species. Demographic noise in the individual-based model
means that species which would survive in the Lotka-Volterra model can become
extinct. These noise-driven extinctions are the focus of the paper. We find
that, for increasing complexity of interactions, ecological communities
generally become less prone to extinctions induced by demographic noise. An
exception are systems composed entirely of predator-prey pairs. These systems
are known to be stable in deterministic Lotka-Volterra models with random
interactions, but, as we show, they are nevertheless particularly vulnerable to
fluctuations.


  The prevention of biodiversity loss in agricultural landscapes to protect
ecosystem stability and functions is of major importance in itself and for the
maintenance of associated ecosystem services. Intense agriculture leads to a
loss in species richness and homogenization of species pools as well as the
fragmentation of natural habitats and groundwater pollution. Constructed
wetlands stand as nature-based solutions (NBS) to buffer the degradation of
water quality by intercepting the transfer of particles, nutrients and
pesticides between crops and surface waters. In karstic watersheds where
sinkholes short-cut surface water directly to groundwater increasing water
resource vulnerability, constructed wetlands are recommended to mitigate
agricultural pollutants. Constructed wetlands also have the potential to
improve landscape connectivity by providing refuge and breeding sites for
wildlife, especially for amphibians. We propose here a methodology to identify
optimal locations for water pollution mitigation using constructed wetlands
from the perspective of habitat connectivity. We use ecological niche modelling
at the regional scale to model the potential of habitat suitability for nine
amphibian species, and to infer how the landscape impedes species movements. We
combine those results to graph theory to identify connectivity priorities at
the operational scale of an agricultural catchment area. Our framework allowed
us to identify optimal areas from the point of view of the species, to analyze
the effect of multifunctional constructed wetlands aiming to both reduce water
pollution and to improve amphibian species habitat overall connectivity. More
generally, we show the potential of habitat connectivity assessment to improve
multifunctionality of NBS for pollution mitigation.


  A key problem in computational biology is discovering the gene expression
changes that regulate cell fate transitions, in which one cell type turns into
another. However, each individual cell cannot be tracked longitudinally, and
cells at the same point in real time may be at different stages of the
transition process. This can be viewed as a problem of learning the behavior of
a dynamical system from observations whose times are unknown. Additionally, a
single progenitor cell type often bifurcates into multiple child cell types,
further complicating the problem of modeling the dynamics. To address this
problem, we developed an approach called variational mixtures of ordinary
differential equations. By using a simple family of ODEs informed by the
biochemistry of gene expression to constrain the likelihood of a deep
generative model, we can simultaneously infer the latent time and latent state
of each cell and predict its future gene expression state. The model can be
interpreted as a mixture of ODEs whose parameters vary continuously across a
latent space of cell states. Our approach dramatically improves data fit,
latent time inference, and future cell state estimation of single-cell gene
expression data compared to previous approaches.


  De novo genome assembly, i.e., rebuilding the sequence of an unknown genome
from redundant and erroneous short sequences, is a key but computationally
intensive step in many genomics pipelines. The exponential growth of genomic
data is increasing the computational demand and requires scalable,
high-performance approaches. In this work, we present a novel
distributed-memory algorithm that, from a string graph representation of the
genome and using sparse matrices, generates the contig set, i.e., overlapping
sequences that form a map representing a region of a chromosome. Using matrix
abstraction, we mask branches in the string graph and compute the connected
component to group genomic sequences that belong to the same linear chain
(i.e., contig). Then, we perform multiway number partitioning to minimize the
load imbalance in local assembly, i.e., concatenation of sequences from a given
contig. Based on the assignment obtained by partitioning, we compute the induce
subgraph function to redistribute sequences between processes, resulting in a
set of local sparse matrices. Finally, we traverse each matrix using
depth-first search to concatenate sequences. Our algorithm shows good scaling
with parallel efficiency up to 80% on 128 nodes, resulting in uniform genome
coverage and showing promising results in terms of assembly quality. Our contig
generation algorithm localizes the assembly process to significantly reduce the
amount of computation spent on this step. Our work is a step forward for
efficient de novo long read assembly of large genomes in a distributed memory.


  BLAST is a standard tool in bioinformatics for creating local sequence
alignments using a "seed-and-extend" approach. Here we introduce an analogous
seed-and-extend algorithm that produces local network alignments: BLANT (Basic
Local Alignment of Network Topology). In Part 1, we introduced BLANT-seed,
which generates graphlet-based seeds using only topological information. Here,
in Part 2, we describe BLANT-extend, which "grows" seeds to larger local
alignments using only topological information. We allow the user to specify
bounds on several measures an alignment must satisfy, including the edge
density, edge commonality (i.e., aligned edges), and node-pair similarity if
such a measure is used; the latter allows the inclusion of sequence-based
similarity, if desired, as well as local topological constraints. BLANT-extend
is able to enumerate all possible alignments satisfying the bounds that can be
grown from each seed, within a specified CPU time or number of generated
alignments. While topology-driven local network alignment has a wide variety of
potential applications outside bioinformatics, here we focus on the alignment
of Protein-Protein Interaction (PPI) networks. We show that BLANT is capable of
finding large, high-quality local alignments when the networks are known to
have high topological similarity -- for example recovering hundreds of
orthologs between networks of the recent Integrated Interaction Database (IID).
Predictably, however, it performs less well when true topological similarity is
absent, as is the case in most current experimental PPI networks that are noisy
and have wide disparity in edge density which results in low common coverage.


  BLAST is a standard tool in bioinformatics for creating local sequence
alignments using a "seed-and-extend" approach. Here we introduce an analogous
seed-and-extend algorithm that produces local network alignments: BLANT, for
Basic Local Alignment of Network Topology. This paper introduces BLANT-seed:
given an input graph, BLANT-seed uses network topology alone to create a
limited, high-specificity index of k-node induced subgraphs called k-graphlets
(analogous to BLASTS's k-mers). The index is constructed so that, if
significant common network topology exists between two graphs, their indexes
are likely to overlap. BLANT-seed then queries the indexes of two networks to
generate a list of common k-graphlets which, when paired, form a seed pair. Our
companion paper (submitted elsewhere) describes BLANT-extend, which "grows"
these seeds to larger local alignments, again using only topological
information.


  Under certain cellular conditions, transcription and mRNA translation in
prokaryotes appear to be "coupled," in which the formation of mRNA transcript
and production of its associated protein are temporally correlated. Such
transcription-translation coupling (TTC) has been evoked as a mechanism that
speeds up the overall process, provides protection during the transcription,
and/or regulates the timing of transcript and protein formation. What molecular
mechanisms underlie ribosome-RNAP coupling and how they can perform these
functions have not been explicitly modeled. We develop and analyze a
continuous-time stochastic model that incorporates ribosome and RNAP elongation
rates, initiation and termination rates, RNAP pausing, and direct ribosome and
RNAP interactions (exclusion and binding). Our model predicts how distributions
of delay times depend on these molecular features of transcription and
translation. We also propose additional measures for TTC: a direct
ribosome-RNAP binding probability and the fraction of time the
translation-transcription process is "protected" from attack by
transcription-terminating proteins. These metrics quantify different aspects of
TTC and differentially depend on parameters of known molecular processes. We
use our metrics to reveal how and when our model can exhibit either
acceleration or deceleration of transcription, as well as protection from
termination. Our detailed mechanistic model provides a basis for designing new
experimental assays that can better elucidate the mechanisms of TTC.


  We study the generalised rock-paper-scissors game with five species whose
organisms face local epidemic outbreaks. As an evolutionary behavioural
survival strategy, organisms of one out of the species move in the direction
with more enemies of their enemies to benefit from protection against
selection. We consider that each organism scans the environment, performing
social distancing instead of agglomerating when perceiving that the density of
sick organisms is higher than a tolerable threshold. Running stochastic
simulations, we study the interference of the adaptive movement survival
strategy in spatial pattern formation, calculating the characteristic length
scale of the typical spatial domains inhabited by organisms of each species. We
compute how social distancing trigger impacts the chances of an individual
being killed in the cyclic game and contaminated by the disease. The outcomes
show that the species predominates in the cyclic game because of the organisms'
local adaptation. The territory occupied by the species grows with the
proportion of individuals learning to trigger the social distancing tactic. We
also show that organisms that perceive large distances more properly execute
the adaptive strategy, promptly triggering the social distancing tactic and
choosing the correct direction to move. Our findings may contribute to
understanding the role of adaptive behaviour when environmental changes
threaten biodiversity.


  We derived a new speed limit in population dynamics, which is a fundamental
limit on the evolutionary rate. By splitting the contributions of selection and
mutation to the evolutionary rate, we obtained the new bound on the speed of
arbitrary observables, named the selection bound, that can be tighter than the
conventional Cram\'{e}r--Rao bound. Remarkably, the selection bound can be much
tighter if the contribution of selection is more dominant than that of
mutation. This tightness can be geometrically characterized by the correlation
between the observable of interest and the growth rate. We also numerically
illustrate the effectiveness of the selection bound in the transient dynamics
of evolutionary processes and discuss how to test our speed limit
experimentally.


  Cancer survival prediction is important for developing personalized
treatments and inducing disease-causing mechanisms. Multi-omics data
integration is attracting widespread interest in cancer research for providing
information for understanding cancer progression at multiple genetic levels.
Many works, however, are limited because of the high dimensionality and
heterogeneity of multi-omics data. In this paper, we propose a novel method to
integrate multi-omics data for cancer survival prediction, called Stacked
AutoEncoder-based Survival Prediction Neural Network (SAEsurv-net). In the
cancer survival prediction for TCGA cases, SAEsurv-net addresses the curse of
dimensionality with a two-stage dimensionality reduction strategy and handles
multi-omics heterogeneity with a stacked autoencoder model. The two-stage
dimensionality reduction strategy achieves a balance between computation
complexity and information exploiting. The stacked autoencoder model removes
most heterogeneities such as data's type and size in the first group of
autoencoders, and integrates multiple omics data in the second autoencoder. The
experiments show that SAEsurv-net outperforms models based on a single type of
data as well as other state-of-the-art methods.


  We analyze algorithmic and computational aspects of biological phenomena,
such as replication and programmed death, in the context of machine learning.
We use two different measures of neuron efficiency to develop machine learning
algorithms for adding neurons to the system (i.e. replication algorithm) and
removing neurons from the system (i.e. programmed death algorithm). We argue
that the programmed death algorithm can be used for compression of neural
networks and the replication algorithm can be used for improving performance of
the already trained neural networks. We also show that a combined algorithm of
programmed death and replication can improve the learning efficiency of
arbitrary machine learning systems. The computational advantages of the
bio-inspired algorithms are demonstrated by training feedforward neural
networks on the MNIST dataset of handwritten images.


  The prevalence of COVID-19 has been the most serious health challenge of the
21th century to date, concerning national health systems on a daily basis,
since December 2019 when it appeared in Wuhan City. Nevertheless, most of the
proposed mathematical methodologies aiming to describe the dynamics of an
epidemic, rely on deterministic models that are not able to reflect the true
nature of its spread. In this paper, we propose a SEIHCRDV model - an
extension/improvement of the classic SIR compartmental model - which also takes
into consideration the populations of exposed, hospitalized, admitted in
intensive care units (ICU), deceased and vaccinated cases, in combination with
an unscented Kalman filter (UKF), providing a dynamic estimation of the time
dependent system's parameters. The stochastic approach is considered necessary,
as both observations and system equations are characterized by uncertainties.
Apparently, this new consideration is useful for examining various pandemics
more effectively. The reliability of the model is examined on the daily
recordings of COVID-19 in France, over a long period of 265 days. Two major
waves of infection are observed, starting in January 2021, which signified the
start of vaccinations in Europe providing quite encouraging predictive
performance, based on the produced NRMSE values. Special emphasis is placed on
proving the non-negativity of SEIHCRDV model, achieving a representative basic
reproductive number R0 and demonstrating the existence and stability of disease
equilibria according to the formula produced to estimate R0. The model
outperforms in predictive ability not only deterministic approaches but also
state-of-the-art stochastic models that employ Kalman filters.


  A hypercycle equation with infinitely many types of macromolecules is
formulated and studied both analytically and numerically. The resulting model
is given by an integro-differential equation of the mixed type. Sufficient
conditions for the existence, uniqueness, and non-negativity of solutions are
formulated and proved. Analytical evidence is provided for the existence of
non-uniform (with respect to the second variable) steady states. Finally,
numerical simulations strongly indicate the existence of a stable nonlinear
wave in the form of the wave train.


  We identify a mechanism for biological spatial pattern formation arising when
the signals that mediate interactions between individuals in a population have
pulsed character. Our general population-signal framework shows that while for
a slow signal-dynamics limit no pattern formation is observed for any values of
the model parameters, for a fast limit, on the contrary, pattern formation can
occur. Furthermore, at these limits, our framework reduces, respectively, to
reaction-diffusion and spatially nonlocal models, thus bridging these
approaches.


  The ability of a chemical reaction network to generate itself by catalyzed
reactions from constantly present environmental food sources is considered a
fundamental property in origin-of-life research. Based on Kaufmann's
autocatalytic sets, Hordijk and Steel have constructed the versatile formalism
of catalytic reaction systems (CRS) to model and to analyze such
self-generating networks, which they named reflexively autocatalytic and food
generated (RAF). Previously, it was established that the subsequent and
simultaenous catalytic functions of the chemicals of a CRS give rise to an
algebraic structure, termed a semigroup model. The semigroup model allows to
naturally consider the function of any subset of chemicals on the whole CRS.
This gives rise to a generative dynamics by iteratively applying the function
of a subset to the externally supplied food set. The fixed point of this
dynamics yields the maximal self-generating set of chemicals. Moreover, the
lattice of all functionally closed self-generating sets of chemicals is
discussed and a structure theorem for this lattice is proven. It is also shown
that a CRS which contains self-generating sets of chemicals cannot be nilpotent
and thus a useful link to the combinatorial theory of finite semigroups is
established. The main technical tool introduced and utilized in this work is
the representation of the semigroup elements as decorated rooted trees,
allowing to translate the generation of chemicals from a given set of resources
into the semigroup language.


  Accurate protein structure prediction can significantly accelerate the
development of life science. The accuracy of AlphaFold2, a frontier end-to-end
structure prediction system, is already close to that of the experimental
determination techniques. Due to the complex model architecture and large
memory consumption, it requires lots of computational resources and time to
implement the training and inference of AlphaFold2 from scratch. The cost of
running the original AlphaFold2 is expensive for most individuals and
institutions. Therefore, reducing this cost could accelerate the development of
life science. We implement AlphaFold2 using PaddlePaddle, namely HelixFold, to
improve training and inference speed and reduce memory consumption. The
performance is improved by operator fusion, tensor fusion, and hybrid
parallelism computation, while the memory is optimized through Recompute,
BFloat16, and memory read/write in-place. Compared with the original AlphaFold2
(implemented with Jax) and OpenFold (implemented with PyTorch), HelixFold needs
only 7.5 days to complete the full end-to-end training and only 5.3 days when
using hybrid parallelism, while both AlphaFold2 and OpenFold take about 11
days. HelixFold saves 1x training time. We verified that HelixFold's accuracy
could be on par with AlphaFold2 on the CASP14 and CAMEO datasets. HelixFold's
code is available on GitHub for free download:
https://github.com/PaddlePaddle/PaddleHelix/tree/dev/apps/protein_folding/helixfold,
and we also provide stable web services on
https://paddlehelix.baidu.com/app/drug/protein/forecast.


  G-Protein Coupled Receptors (GPCRs) are a big family of eukaryotic cell
transmembrane proteins, responsible for numerous biological processes. From a
practical viewpoint around 34\% of the drugs approved by the US Food and Drug
Administration target these receptors. They can be analyzed from their
simulated molecular dynamics, including the prediction of their behavior in the
presence of drugs. In this paper, the capability of Long Short-Term Memory
Networks (LSTMs) are evaluated to learn and predict the molecular dynamic
trajectories of a receptor. Several models were trained with the 3D position of
the amino acids of the receptor considering different transformations on the
position of the amino acid, such as their centers of mass, the geometric
centers and the position of the $\alpha$--carbon for each amino acid. The error
of the prediction of the position was evaluated by the mean average error (MAE)
and root-mean-square deviation (RMSD). The LSTM models show a robust
performance, with results comparable to the state-of-the-art in non-dynamic 3D
predictions. The best MAE and RMSD values were found for the mass center of the
amino acids with 0.078 {\AA} and 0.156 {\AA} respectively. This work shows the
potential of LSTM to predict the molecular dynamics of GPRCs.


  Voluntary vaccination is effective to prevent infectious diseases from
spreading. Both vaccination behavior and cognition of the vaccination risk play
important roles in individual vaccination decision making. However, it is not
clear how the co-evolution of the two shapes the population-wide vaccination
behavior. We establish a coupled dynamics of epidemic, vaccination behavior and
perceived vaccination risk with three different time scales. We assume that the
increase of vaccination level inhibits the rise of perceived vaccination risk,
and the increase of perceived vaccination risk inhibits the rise of vaccination
level. It is shown that the resulting vaccination behavior is similar to the
stag-hunt game, provided that the basic reproductive ratio is moderate and that
the epidemic dynamics evolves fast. This is in contrast with the previous view
that vaccination is a snowdrift like game. Furthermore, we find that epidemic
breaks out repeatedly and eventually leads to vaccine scares if these three
dynamics evolve on a similar time scale. And we propose some ways to promote
vaccination behavior, such as controlling side-effect bias and perceived
vaccination costs. Our work sheds light on epidemic control via vaccination by
taking into account the co-evolutionary dynamics of cognition and behavior.


  Extracting informative representations of molecules using Graph neural
networks (GNNs) is crucial in AI-driven drug discovery. Recently, the graph
research community has been trying to replicate the success of self-supervised
pretraining in natural language processing, with several successes claimed.
However, we find the benefit brought by self-supervised pretraining on small
molecular data can be negligible in many cases. We conduct thorough ablation
studies on the key components of GNN pretraining, including pretraining
objectives, data splitting methods, input features, pretraining dataset scales,
and GNN architectures, to see how they affect the accuracy of the downstream
tasks. Our first important finding is, self-supervised graph pretraining do not
always have statistically significant advantages over non-pretraining methods
in many settings. Secondly, although noticeable improvement can be observed
with additional supervised pretraining, the improvement may diminish with
richer features or more balanced data splits. Thirdly, hyper-parameters could
have larger impacts on accuracy of downstream tasks than the choice of
pretraining tasks, especially when the scales of downstream tasks are small.
Finally, we provide our conjectures where the complexity of some pretraining
methods on small molecules might be insufficient, followed by empirical
evidences on different pretraining datasets.


  Computational antibody design seeks to automatically create an antibody that
binds to an antigen. The binding affinity is governed by the 3D binding
interface where antibody residues (paratope) closely interact with antigen
residues (epitope). Thus, predicting 3D paratope-epitope complex (docking) is
the key to finding the best paratope. In this paper, we propose a new model
called Hierarchical Equivariant Refinement Network (HERN) for paratope docking
and design. During docking, HERN employs a hierarchical message passing network
to predict atomic forces and use them to refine a binding complex in an
iterative, equivariant manner. During generation, its autoregressive decoder
progressively docks generated paratopes and builds a geometric representation
of the binding interface to guide the next residue choice. Our results show
that HERN significantly outperforms prior state-of-the-art on paratope docking
and design benchmarks.


  Improved understanding of the relation between the behavior of RAS and RAF
proteins and the local lipid environment in the cell membrane is critical for
getting insights into the mechanisms underlying cancer formation. In this work,
we employ deep learning (DL) to learn this relationship by predicting protein
orientational states of RAS and RAS-RAF protein complexes with respect to the
lipid membrane based on the lipid densities around the protein domains from
coarse-grained (CG) molecular dynamics (MD) simulations. Our DL model can
predict six protein states with an overall accuracy of over 80%. The findings
of this work offer new insights into how the proteins modulate the lipid
environment, which in turn may assist designing novel therapies to regulate
such interactions in the mechanisms associated with cancer development.


  Composed of amino acid chains that influence how they fold and thus dictating
their function and features, proteins are a class of macromolecules that play a
central role in major biological processes and are required for the structure,
function, and regulation of the body's tissues. Understanding protein functions
is vital to the development of therapeutics and precision medicine, and hence
the ability to classify proteins and their functions based on measurable
features is crucial; indeed, the automatic inference of a protein's properties
from its sequence of amino acids, known as its primary structure, remains an
important open problem within the field of bioinformatics, especially given the
recent advancements in sequencing technologies and the extensive number of
known but uncategorized proteins with unknown properties. In this work, we
demonstrate and compare the performance of several deep learning frameworks,
including novel bi-directional LSTM and convolutional models, on widely
available sequencing data from the Protein Data Bank (PDB) of the Research
Collaboratory for Structural Bioinformatics (RCSB), as well as benchmark this
performance against classical machine learning approaches, including k-nearest
neighbors and multinomial regression classifiers, trained on experimental data.
Our results show that our deep learning models deliver superior performance to
classical machine learning methods, with the convolutional architecture
providing the most impressive inference performance.


  Allostery commonly refers to the mechanism that regulates protein activity
through the binding of a molecule at a different, usually distal, site from the
orthosteric site. The omnipresence of allosteric regulation in nature and its
potential for drug design and screening render the study of allostery
invaluable. Nevertheless, challenges remain as few computational methods are
available to effectively predict allosteric sites, identify signalling pathways
involved in allostery, or to aid with the design of suitable molecules
targeting such sites. Recently, bond-to-bond propensity analysis has been shown
successful at identifying allosteric sites for a large and diverse group of
proteins from knowledge of the orthosteric sites and its ligands alone by using
network analysis applied to energy-weighted atomistic protein graphs. To
address the identification of signalling pathways, we propose here a method to
compute and score paths of optimised propensity that link the orthosteric site
with the identified allosteric sites, and identifies crucial residues that
contribute to those paths. We showcase the approach with three well-studied
allosteric proteins: h-Ras, caspase-1, and 3-phosphoinositide-dependent
kinase-1 (PDK1). Key residues in both orthosteric and allosteric sites were
identified and showed agreement with experimental results, and pivotal
signalling residues along the pathway were also revealed, thus providing
alternative targets for drug design. By using the computed path scores, we were
also able to differentiate the activity of different allosteric modulators.


  The origin of life is shrouded in mystery, with few surviving clues, obscured
by evolutionary competition. Previous reviews have touched on the complementary
approaches of top-down and bottom-up synthetic biology to augment our
understanding of living systems. Here we point out the synergies between these
fields, especially between bottom-up synthetic biology and origin of life
research. We explore recent progress made in artificial cell compartmentation
in line with the crowded cell, its metabolism, as well as cycles of growth and
division, and how those efforts are starting to be combined. Though the
complexity of current life is among its most striking characteristics, none of
life's essential features require it, and they are unlikely to have emerged
thus complex from the beginning. Rather than recovering the one true origin
lost in time, current research converges towards reproducing the emergence of
minimal life, by teasing out how complexity and evolution may arise from a set
of essential components.


  The computer artificial intelligence system AlphaFold has recently predicted
previously unknown three-dimensional structures of thousands of proteins.
Focusing on the subset with high-confidence scores, we algorithmically analyze
these predictions for cases where the protein backbone exhibits rare
topological complexity, i.e. knotting. Amongst others, we discovered a
$7_1$-knot, the most topologically complex knot ever found in a protein, as
well several 6-crossing composite knots comprised of two methyltransferase or
carbonic anhydrase domains, each containing a simple trefoil knot. These deeply
embedded composite knots occur evidently by gene duplication and
interconnection of knotted dimers. Finally, we report two new five-crossing
knots including the first $5_1$-knot. Our list of analyzed structures forms the
basis for future experimental studies to confirm these novel knotted topologies
and to explore their complex folding mechanisms.


  A central issue in the study of large complex network systems, such as power
grids, financial networks, and ecological systems, is to understand their
response to dynamical perturbations. Recent studies recognize that many real
networks show nonnormality and that nonnormality can give rise to
reactivity--the capacity of a linearly stable system to amplify its response to
perturbations, oftentimes exciting nonlinear instabilities. Here, we identify
network structural properties underlying the pervasiveness of nonnormality and
reactivity in real directed networks, which we establish using the most
extensive data set of such networks studied in this context to date. The
identified properties are imbalances between incoming and outgoing network
links and paths at each node. Based on this characterization, we develop a
theory that quantitatively predicts nonnormality and reactivity and explains
the observed pervasiveness. We suggest that these results can be used to
design, upgrade, control, and manage networks to avoid or promote network
instabilities.


  The transport of particles in cells is influenced by the properties of
intracellular networks they traverse while searching for localized target
regions or reaction partners. Moreover, given the rapid turnover in many
intracellular structures, it is crucial to understand how temporal changes in
the network structure affect diffusive transport. In this work, we use network
theory to characterize complex intracellular biological environments across
scales. We develop an efficient computational method to compute the mean first
passage times for simulating a particle diffusing along two-dimensional planar
networks extracted from fluorescence microscopy imaging. We first benchmark
this methodology in the context of synthetic networks, and subsequently apply
it to live-cell data from endoplasmic reticulum tubular networks.


  We propose to use the mean time between infections (MTBI) metric as obtained
from a recently introduced non-homogeneous Markov stochastic model. Different
types of parameter calibration are performed. We estimate the MTBI using data
from different time windows and from the whole stage history and compare the
results. In order to detect waves and stages in the input data, a preprocessing
filtering technique is applied. The results of applying this indicator to the
COVID-19 reported data of infections from Argentina, Germany and the United
States are shown. We find that the MTBI behaves similarly with respect to the
different data inputs, whereas the model parameters completely change their
behaviour. Evolution over time of the parameters and the MTBI indicator is also
shown. We show evidence to support the claim that the MTBI is a rather good
indicator in order to measure the spreading speed of an epidemic, having
similar values whatever the input data size.


  Evolutionary game theory has proved to be a powerful tool to probe the
self-organisation of collective behaviour by considering frequency-dependent
fitness in evolutionary processes. It has shown that the stability of a
strategy depends not only on the payoffs received after each encounter but also
on the population's size. Here, we study $2\times2$ games in well-mixed finite
populations by analysing the fixation probabilities of single mutants as
functions of population size. We proved that 9 out of the 24 possible games
always lead to monotonically decreasing functions, similarly to fixed fitness
scenarios. However, fixation functions showed increasing regions under 12
distinct anti-coordination, coordination, and dominance games. Perhaps
counter-intuitively, this establishes that single-mutant strategies often
benefit from being in larger populations. Fixation functions that increase from
a global minimum to a positive asymptotic value are pervasive but may have been
easily concealed by the weak selection limit. We obtained sufficient conditions
to observe fixation increasing for small populations and three distinct ways
this can occur. Finally, we describe fixation functions with increasing regions
bounded by two extremes under intermediate population sizes. We associate their
occurrence with transitions from having one global extreme to other shapes.


  Metacells are disjoint and homogeneous groups of single-cell profiles,
representing discrete and highly granular cell states. Existing metacell
algorithms tend to use only one modality to infer metacells, even though
single-cell multi-omics datasets profile multiple molecular modalities within
the same cell. Here, we present \textbf{C}ross-M\textbf{O}dal
\textbf{E}mbedding for \textbf{M}etaCell Identification (COEM), which utilizes
an embedded space leveraging the information of both scATAC-seq and scRNA-seq
to perform aggregation, balancing the trade-off between fine resolution and
sufficient sequencing coverage. COEM outperforms the state-of-the-art method
SEACells by efficiently identifying accurate and well-separated metacells
across datasets with continuous and discrete cell types. Furthermore, COEM
significantly improves peak-to-gene association analyses, and facilitates
complex gene regulatory inference tasks.


  Homochirality is a generic and unique property of all biochemical life and is
considered a universal and agnostic biosignature. Upon interaction with
unpolarized light, homochirality induces fractional circular polarization in
the light that is scattered from it, which can be sensed remotely. As such, it
can be a prime candidate biosignature in the context of future life-detection
missions and observatories. The linear polarizance of vegetation is also
sometimes envisaged as a biosignature, although it does not share the molecular
origin as circular polarization. It is known that the linear polarization of
surfaces is strongly dependent on the phase angle. The relation between the
phase angle and circular polarization stemming from macromolecular assemblies,
such as in vegetation, however, remained unclear. We demonstrate in this study
using the average of 27 different species that the circular polarization phase
angle dependency of vegetation induces relatively small changes in spectral
shape and mostly affects the signal magnitude. With these results we underline
the use of circular spectropolarimetry as a promising agnostic biosignature
complementary to the use of linear spectropolarimetry and scalar reflectance.


  Testing for infection with SARS-CoV-2 is an important intervention in
reducing onwards transmission of COVID-19, particularly when combined with the
isolation and contact-tracing of positive cases. Many countries with the
capacity to do so have made use of lab-processed Polymerase Chain Reaction
(PCR) testing targeted at individuals with symptoms and the contacts of
confirmed cases. Alternatively, Lateral Flow Tests (LFTs) are able to deliver a
result quickly, without lab-processing and at a relatively low cost. Their
adoption can support regular mass asymptomatic testing, allowing earlier
detection of infection and isolation of infectious individuals. In this paper
we extend and apply the agent-based epidemic modelling framework Covasim to
explore the impact of regular asymptomatic testing on the peak and total number
of infections in an emerging COVID-19 wave. We explore testing with LFTs at
different frequency levels within a population with high levels of immunity and
with background symptomatic PCR testing, case isolation and contact tracing for
testing. The effectiveness of regular asymptomatic testing was compared with
`lockdown' interventions seeking to reduce the number of non-household contacts
across the whole population through measures such as mandating working from
home and restrictions on gatherings. Since regular asymptomatic testing
requires only those with a positive result to reduce contact, while lockdown
measures require the whole population to reduce contact, any policy decision
that seeks to trade off harms from infection against other harms will not
automatically favour one over the other. Our results demonstrate that, where
such a trade off is being made, at moderate rates of early exponential growth
regular asymptomatic testing has the potential to achieve significant infection
control without the wider harms associated with additional lockdown measures.


  Learning expressive molecular representations is crucial to facilitate the
accurate prediction of molecular properties. Despite the significant
advancement of graph neural networks (GNNs) in molecular representation
learning, they generally face limitations such as neighbors-explosion,
under-reaching, over-smoothing, and over-squashing. Also, GNNs usually have
high computational costs because of the large-scale number of parameters.
Typically, such limitations emerge or increase when facing relatively
large-size graphs or using a deeper GNN model architecture. An idea to overcome
these problems is to simplify a molecular graph into a small, rich, and
informative one, which is more efficient and less challenging to train GNNs. To
this end, we propose a novel molecular graph coarsening framework named FunQG
utilizing Functional groups, as influential building blocks of a molecule to
determine its properties, based on a graph-theoretic concept called Quotient
Graph. By experiments, we show that the resulting informative graphs are much
smaller than the molecular graphs and thus are good candidates for training
GNNs. We apply the FunQG on popular molecular property prediction benchmarks
and then compare the performance of some popular baseline GNNs on the obtained
datasets with the performance of several state-of-the-art baselines on the
original datasets. By experiments, this method significantly outperforms
previous baselines on various datasets, besides its dramatic reduction in the
number of parameters and low computational costs. Therefore, the FunQG can be
used as a simple, cost-effective, and robust method for solving the molecular
representation learning problem.


  Conformational dynamics is crucial for ribonucleic acid (RNA) function.
Techniques such as nuclear magnetic resonance, cryo-electron microscopy, small-
and wide-angle X-ray scattering, chemical probing, single-molecule F\"orster
resonance energy transfer or even thermal or mechanical denaturation
experiments probe RNA dynamics at different time and space resolutions. Their
combination with accurate atomistic molecular dynamics (MD) simulations paves
the way for quantitative and detailed studies of RNA dynamics. First,
experiments provide a quantitative validation tool for MD simulations. Second,
available data can be used to refine simulated structural ensembles to match
experiments. Finally, comparison with experiments allows for improving MD force
fields that are transferable to new systems for which data is not available.
Here we review the recent literature and provide our perspective on this field.


  Could nature be harnessing quantum mechanics in cilia to optimize the
sensitivity of the mechanism of left-right symmetry breaking during development
in vertebrates? I evaluate whether mechanosensing -- i.e., the detection of a
left-right asymmetric signal through mechanical stimulation of sensory cilia,
as opposed to biochemical signalling -- might be functioning in the embryonic
left-right organizer of the vertebrate bodyplan through quantum mechanics. I
conclude that there is a possible role for quantum biology in mechanosensing in
cilia. The system may not be limited by classical thermal noise, but instead by
quantum noise, with an amplification process providing active cooling.


  Traditional Chinese medicine (TCM) relies on natural medical products to
treat symptoms and diseases. While clinical data have demonstrated the
effectiveness of selected TCM-based treatments, the mechanistic root of how TCM
herbs treat diseases remains largely unknown. More importantly, current
approaches focus on single herbs or prescriptions, missing the high-level
general principles of TCM. To uncover the mechanistic nature of TCM on a system
level, in this work we establish a generic network medicine framework for TCM
from the human protein interactome. Applying our framework reveals a network
pattern between symptoms (diseases) and herbs in TCM. We first observe that
genes associated with a symptom are not distributed randomly in the
interactome, but cluster into localized modules; furthermore, a short network
distance between two symptom modules is indicative of the symptoms'
co-occurrence and similarity. Next, we show that the network proximity of a
herb's targets to a symptom module is predictive of the herb's effectiveness in
treating the symptom. We validate our framework with real-world hospital
patient data by showing that (1) shorter network distance between symptoms of
inpatients correlates with higher relative risk (co-occurrence), and (2)
herb-symptom network proximity is indicative of patients' symptom recovery rate
after herbal treatment. Finally, we identified novel herb-symptom pairs in
which the herb's effectiveness in treating the symptom is predicted by network
and confirmed in hospital data, but previously unknown to the TCM community.
These predictions highlight our framework's potential in creating herb
discovery or repurposing opportunities. In conclusion, network medicine offers
a powerful novel platform to understand the mechanism of traditional medicine
and to predict novel herbal treatment against diseases.


  The present paper is devoted to the study of the long term dynamics of
diffusion processes modelling a single species that experiences both
demographic and environmental stochasticity. In our setting, the long term
dynamics of the diffusion process in the absence of demographic stochasticity
is determined by the sign of $\Lambda_0$, the external Lyapunov exponent, as
follows: $\Lambda_0<0$ implies (asymptotic) extinction and $\Lambda_0>0$
implies convergence to a unique positive stationary distribution $\mu_0$. If
the system is of size $\frac{1}{\epsilon^{2}}$ for small $\epsilon>0$ (the
intensity of demographic stochasticity), demographic effects will make the
extinction time finite almost surely. This suggests that to understand the
dynamics one should analyze the quasi-stationary distribution (QSD)
$\mu_\epsilon$ of the system. The existence and uniqueness of the QSD is
well-known under mild assumptions.
  We look at what happens when the population size is sent to infinity, i.e.,
when $\epsilon\to 0$. We show that the external Lyapunov exponent still plays a
key role: 1) If $\Lambda_0<0$, then $\mu_\epsilon\to \delta_0$, the mean
extinction time is of order $|\ln \epsilon|$ and the extinction rate associated
with the QSD $\mu_{\epsilon}$ has a lower bound of order
$\frac{1}{|\ln\epsilon|}$; 2) If $\Lambda_0>0$, then $\mu_\epsilon\to \mu_0$,
the mean extinction time is polynomial in $\frac{1}{\epsilon^{2}}$ and the
extinction rate is polynomial in $\epsilon^{2}$. Furthermore, when
$\Lambda_0>0$ we are able to show that the system exhibits multiscale dynamics:
at first the process quickly approaches the QSD $\mu_\epsilon$ and then, after
spending a polynomially long time there, it relaxes to the extinction state. We
give sharp asymptotics in $\epsilon$ for the time spent close to
$\mu_\epsilon$.


  The rapid spread of the COVID-19 pandemic has resulted in an unprecedented
amount of sequence data of the SARS-CoV-2 genome -- millions of sequences and
counting. This amount of data, while being orders of magnitude beyond the
capacity of traditional approaches to understanding the diversity, dynamics,
and evolution of viruses is nonetheless a rich resource for machine learning
(ML) approaches as alternatives for extracting such important information from
these data. It is of hence utmost importance to design a framework for testing
and benchmarking the robustness of these ML models.
  This paper makes the first effort (to our knowledge) to benchmark the
robustness of ML models by simulating biological sequences with errors. In this
paper, we introduce several ways to perturb SARS-CoV-2 genome sequences to
mimic the error profiles of common sequencing platforms such as Illumina and
PacBio. We show from experiments on a wide array of ML models that some
simulation-based approaches are more robust (and accurate) than others for
specific embedding methods to certain adversarial attacks to the input
sequences. Our benchmarking framework may assist researchers in properly
assessing different ML models and help them understand the behavior of the
SARS-CoV-2 virus or avoid possible future pandemics.


  We propose two SIR models which incorporate sociological behavior of groups
of individuals. It is these differences in behaviors which impose different
infection rates on the individual susceptible populations, rather than
biological differences. We compute the basic reproduction number for each
model, as well as analyze the sensitivity of $R_0$ to changes in sociological
parameter values.