  The actin cytoskeleton of cells is in continuous motion due to both
polymerization of new filaments and their contraction by myosin II molecular
motors. Through adhesion to the substrate, such intracellular flow can be
converted into cell migration. Recently, optogenetics has emerged as a new
powerful experimental method to control both actin polymerization and myosin II
contraction. While optogenetic control of polymerization can initiate cell
migration by generating protrusion, it is less clear if and how optogenetic
control of contraction can also affect cell migration. Here we analyze the
latter situation using a minimal variant of active gel theory into which we
include optogenetic activation as a spatiotemporally constrained perturbation.
The model can describe the symmetrical flow of the actomyosin system observed
in optogenetic experiments, but not the long-lasting polarization required for
cell migration. Motile solutions become possible if cytoskeletal polymerization
is included through the boundary conditions. Optogenetic activation of
contraction can then initiate locomotion in a symmetrically spreading cell and
strengthen motility in an asymmetrically polymerizing one. If designed
appropriately, it can also arrest motility even for protrusive boundaries.


  Large self-supervised models pretrained on millions of protein sequences have
recently gained popularity in generating embeddings of protein sequences for
protein function prediction. However, the absence of random baselines makes it
difficult to conclude whether pretraining has learned useful information for
protein function prediction. Here we show that one-hot encoding and random
embeddings, both of which do not require any pretraining, are strong baselines
for protein function prediction across 14 diverse sequence-to-function tasks.


  Traditional Indigenous marriage rules have been studied extensively since the
mid 1800s. Despite this, they have historically been cast aside as having very
little utility. This is, in large part, due to a focus on trying to understand
broad-stroke marriage restrictions or how they may evolve. Here, taking the
Gamilaraay system as a case study, we instead ask how relatedness may be
distributed under such a system. We show, remarkably, that this system
dynamically trades off kin avoidance to minimise incidence of recessive
diseases against expected levels of cooperation, as understood formally through
Hamilton's rule.


  This paper deals with an important subject in classification problems
addressed by machine learning techniques: the evaluation of the influence of
each of the features on the classification of individuals. Specifically, a
measure of that influence is introduced using the Shapley value of cooperative
games. In addition, an axiomatic characterisation of the proposed measure is
provided based on properties of efficiency and balanced contributions.
Furthermore, some experiments have been designed in order to validate the
appropriate performance of such measure. Finally, the methodology introduced is
applied to a sample of COVID-19 patients to study the influence of certain
demographic or risk factors on various events of interest related to the
evolution of the disease.


  Cell polarization underlies many cellular processes, such as differentiation,
migration, and budding. Many living cells, such as budding yeast and fission
yeast, use cytoskeletal structures to actively transport proteins to one
location on the membrane and create a high density spot of membrane-bound
proteins. Yet, the thermodynamic constraints on filament-based cell
polarization remain unknown. We show by mathematical modeling that cell
polarization requires detailed balance to be broken, and we quantify the
free-energy cost of maintaining a polarized state of the cell. Our study
reveals that detailed balance cannot only be broken via the active transport of
proteins along filaments, but also via a chemical modification cycle, allowing
detailed balance to be broken by the shuttling of proteins between the
filament, membrane, and cytosol. Our model thus shows that cell polarization
can be established via two distinct driving mechanisms, one based on active
transport and one based on non-equilibrium binding. Furthermore, the model
predicts that the driven binding process dissipates orders of magnitude less
free-energy than the transport-based process to create the same membrane spot.
Active transport along filaments may be sufficient to create a polarized
distribution of membrane-bound proteins, but an additional chemical
modification cycle of the proteins themselves is more efficient and less
sensitive to the physical exclusion of proteins on the transporting filaments,
providing insight in the design principles of the Pom1/Tea1/Tea4 system in
fission yeast and the Cdc42 system in budding yeast.


  The outbreak of a novel coronavirus causing severe acute respiratory syndrome
in December 2019 has escalated into a worldwide pandemic. In this work, we
propose a compartmental model to describe the dynamics of transmission of
infection and use it to obtain the optimal vaccination control. The model
accounts for the various stages of the vaccination and the optimisation is
focused on minimising the infections to protect the population and relieve the
healthcare system. As a case study we selected the Republic of Ireland. We use
data provided by Ireland's COVID-19 Data-Hub and simulate the evolution of the
pandemic with and without the vaccination in place for two different scenarios,
one representative of a national lockdown situation and the other indicating
looser restrictions in place. One of the main findings of our work is that the
optimal approach would involve a vaccination programme where the older
population is vaccinated in larger numbers earlier while simultaneously part of
the younger population also gets vaccinated to lower the risk of transmission
between groups. We compare our simulated results with that of the vaccination
policy taken by the Irish government to explore the advantages of our
optimisation method. Our comparison suggests that a similar reduction in cases
may have been possible even with a reduced set of vaccinations being available
for use.


  Homeostasis is a mechanism by which a feature can remain invariant with
change in external parameters. We adopt the definition of homeostasis in the
context of singularity theory. We make a connection between homeostasis and the
theory of injective reaction networks. In particular, we show that a reaction
network cannot exhibit homeostasis if a modified reaction network (which we
call the homeostasis-associated reaction network) is injective. We provide
examples of reaction networks which can or cannot exhibit homeostasis by
analyzing the injectivity of the homeostasis-associated reaction network.


  A reaction system exhibits "absolute concentration robustness" (ACR) in some
species if the positive steady-state value of that species does not depend on
initial conditions. Mathematically, this means that the positive part of the
variety of the steady-state ideal lies entirely in a hyperplane of the form
$x_i=c$, for some $c>0$. Deciding whether a given reaction system -- or those
arising from some reaction network -- exhibits ACR is difficult in general, but
here we show that for many simple networks, assessing ACR is straightforward.
Indeed, our criteria for ACR can be performed by simply inspecting a network or
its standard embedding into Euclidean space. Our main results pertain to
networks with many conservation laws, so that all reactions are parallel to one
other. Such "one-dimensional" networks include those networks having only one
species. We also consider networks with only two reactions, and show that ACR
is characterized by a well-known criterion of Shinar and Feinberg. Finally, up
to some natural ACR-preserving operations -- relabeling species, lengthening a
reaction, and so on -- only three families of networks with two reactions and
two species have ACR. Our results are proven using algebraic and combinatorial
techniques.


  The Susceptible-Exposed-Infectious-Recovered (SEIR) model is applied in
several countries to ascertain the spread of the coronavirus disease 2019
(COVID-19). We consider discrete-time SEIR epidemic model in a closed system
which does not account for births or deaths, total population size under
consideration is constant. This dynamical system generated by a non-linear
evolution operator depending on four parameters. Under some conditions on
parameters we reduce the evolution operator to a quadratic stochastic operator
(QSO) which maps 3-dimensional simplex to itself. We show that the QSO has
uncountable set of fixed points (all laying on the boundary of the simplex). It
is shown that all trajectories of the dynamical system (generated by the QSO)
of the SEIR model are convergent (i.e. the QSO is regular). Moreover, we
discuss the efficiency of the model for Uzbekistan.


  Persistent homology has undergone significant development in recent years.
However, one outstanding challenge is to build a coherent statistical inference
procedure on persistent diagrams. In this paper, we first present a new lattice
path representation for persistent diagrams. We then develop a new exact
statistical inference procedure for lattice paths via combinatorial
enumerations. The lattice path method is applied to the topological
characterization of the protein structures of the COVID-19 virus. We
demonstrate that there are topological changes during the conformational change
of spike proteins.


  Mitochondrial diseases are largely caused by dysfunction in mitochondrial
proteins. However, annotations of human mitochondrial proteins are scattered
across various public databases and individual studies. To facilitate research
aimed at elucidating mitochondrial functions, we constructed the
MEGADOCK-Web-Mito database as a protein-protein interaction (PPI) prediction
data archive, including prediction results for exhaustive protein pairs of 654
mitochondria-related human proteins. MEGADOCK-Web-Mito enables users to search
for all PPI prediction results efficiently and comprehensively. In particular,
we linked functional annotations to each human mitochondrial protein. The
comprehensive and specialized human mitochondrial PPI prediction results and
searching function of MEGADOCK-Web-Mito will support further research on
mitochondria and mitochondrial diseases.


  Identifying the mutations that drive cancer growth is key in clinical
decision making and precision oncology. As driver mutations confer selective
advantage and thus have an increased likelihood of occurrence, frequency-based
statistical models are currently favoured. These methods are not suited to
rare, low frequency, driver mutations. The alternative approach to address this
is through functional-impact scores, however methods using this approach are
highly prone to false positives. In this paper, we propose a novel combination
method for driver mutation identification, which uses the power of both
statistical modelling and functional-impact based methods. Initial results show
this approach outperforms the state-of-the-art methods in terms of precision,
and provides comparable performance in terms of area under receiver operating
characteristic curves (AU-ROC). We believe that data-driven systems based on
machine learning, such as these, will become an integral part of precision
oncology in the near future.


  Adversarial defenses protect machine learning models from adversarial
attacks, but are often tailored to one type of model or attack. The lack of
information on unknown potential attacks makes detecting adversarial examples
challenging. Additionally, attackers do not need to follow the rules made by
the defender. To address this problem, we take inspiration from the concept of
Applicability Domain in cheminformatics. Cheminformatics models struggle to
make accurate predictions because only a limited number of compounds are known
and available for training. Applicability Domain defines a domain based on the
known compounds and rejects any unknown compound that falls outside the domain.
Similarly, adversarial examples start as harmless inputs, but can be
manipulated to evade reliable classification by moving outside the domain of
the classifier. We are the first to identify the similarity between
Applicability Domain and adversarial detection. Instead of focusing on unknown
attacks, we focus on what is known, the training data. We propose a simple yet
robust triple-stage data-driven framework that checks the input globally and
locally, and confirms that they are coherent with the model's output. This
framework can be applied to any classification model and is not limited to
specific attacks. We demonstrate these three stages work as one unit,
effectively detecting various attacks, even for a white-box scenario.


  In the new paradigm of health-centric governance, policy makers are in a
constant need for appropriate metrics and estimates in order to determine the
best policies in a non-arbitrary fashion. Thus, in this paper, a
compartmentalized model for the transmission of COVID-19 is developed to
facilitate policy making. A socially distanced compartment is added to the
model and its utility in quantifying the magnitude of voluntary social
distancing is illustrated. Modifications are made to incorporate inter-region
migration, and suitable metrics are proposed to quantify the impact of
migration on the rise of cases. The healthcare capacity is modeled and a method
is developed to study the consequences of the saturation of the healthcare
system. The model and related measures are used to study the nature of the
transmission and spread of COVID-19 in India, and appropriate insights are
drawn.


  Effects like selection in evolution as well as fertility inheritance in the
development of populations can lead to a higher degree of asymmetry in
evolutionary trees than expected under a null hypothesis. To identify and
quantify such influences, various balance indices were proposed in the
phylogenetic literature and have been in use for decades.
  However, so far no balance index was based on the number of \emph{symmetry
nodes}, even though symmetry nodes play an important role in other areas of
mathematical phylogenetics and despite the fact that symmetry nodes are a quite
natural way to measure balance or symmetry of a given tree.
  The aim of this manuscript is thus twofold: First, we will introduce the
\emph{symmetry nodes index} as an index for measuring balance of phylogenetic
trees and analyze its extremal properties. We also show that this index can be
calculated in linear time. This new index turns out to be a generalization of a
simple and well-known balance index, namely the \emph{cherry index}, as well as
a specialization of another, less established, balance index, namely
\emph{Rogers' $J$ index}. Thus, it is the second objective of the present
manuscript to compare the new symmetry nodes index to these two indices and to
underline its advantages. In order to do so, we will derive some extremal
properties of the cherry index and Rogers' $J$ index along the way and thus
complement existing studies on these indices. Moreover, we used the programming
language \textsf{R} to implement all three indices in the software package
\textsf{symmeTree}, which has been made publicly available.


  Ants are social insects. When the existing nest of an ant colony becomes
uninhabitable, the hunt for a new suitable location for migration of the colony
begins. Normally, multiple sites may be available as the potential new nest
site. Distinct sites may be chosen by different scout ants based on their own
assessments. Since the individual assessment is error prone, many ants may
choose inferior site(s). But, the collective decision that emerges from the
sequential and decentralized decision making process is often far better. We
develop a model for this multi-stage decision making process. A stochastic
drift-diffusion model (DDM) captures the sequential information accumulation by
individual scout ants for arriving at their respective individual choices. The
subsequent tandem runs of the scouts, whereby they recruit their active
nestmates, is modelled in terms of suitable adaptations of the totally
asymmetric simple exclusion processes (TASEP). By a systematic analysis of the
model we explore the conditions that determine the speed of the emergence of
the collective decision and the quality of that decision. More specifically, we
demonstrate that collective decision of the colony is much less error-prone
that the individual decisions of the scout ants. We also compare our
theoretical predictions with experimental data.


  How cooperation emerges in human societies is both an evolutionary enigma,
and a practical problem with tangible implications for societal health.
Population structure has long been recognized as a catalyst for cooperation
because local interactions enable reciprocity. Analysis of this phenomenon
typically assumes bi-directional social interactions, even though real-world
interactions are often uni-directional. Uni-directional interactions -- where
one individual has the opportunity to contribute altruistically to another, but
not conversely -- arise in real-world populations as the result of
organizational hierarchies, social stratification, popularity effects, and
endogenous mechanisms of network growth. Here we expand the theory of
cooperation in structured populations to account for both uni- and
bi-directional social interactions. Even though directed interactions remove
the opportunity for reciprocity, we find that cooperation can nonetheless be
favored in directed social networks and that cooperation is provably maximized
for networks with an intermediate proportion of directed interactions, as
observed in many empirical settings. We also identify two simple structural
motifs that allow efficient modification of interaction directionality to
promote cooperation by orders of magnitude. We discuss how our results relate
to the concepts of generalized and indirect reciprocity.


  Thanks to the increasing availability of genomics and other biomedical data,
many machine learning approaches have been proposed for a wide range of
therapeutic discovery and development tasks. In this survey, we review the
literature on machine learning applications for genomics through the lens of
therapeutic development. We investigate the interplay among genomics,
compounds, proteins, electronic health records (EHR), cellular images, and
clinical texts. We identify twenty-two machine learning in genomics
applications across the entire therapeutics pipeline, from discovering novel
targets, personalized medicine, developing gene-editing tools all the way to
clinical trials and post-market studies. We also pinpoint seven important
challenges in this field with opportunities for expansion and impact. This
survey overviews recent research at the intersection of machine learning,
genomics, and therapeutic development.


  Biclustering is a data mining technique which searches for local patterns in
numeric tabular data with main application in bioinformatics. This technique
has shown promise in multiple areas, including development of biomarkers for
cancer, disease subtype identification, or gene-drug interactions among others.
In this paper we introduce EBIC.JL - an implementation of one of the most
accurate biclustering algorithms in Julia, a modern highly parallelizable
programming language for data science. We show that the new version maintains
comparable accuracy to its predecessor EBIC while converging faster for the
majority of the problems. We hope that this open source software in a
high-level programming language will foster research in this promising field of
bioinformatics and expedite development of new biclustering methods for big
data.


  Cardiac muscle contraction is initiated by an elementary Ca signal (called Ca
spark) which is achieved by collective action of Ca release channels in a
cluster. The mechanism of this synchronization remains uncertain. This paper
approaches Ca spark activation as an emergent phenomenon of an interactive
system of release channels. We construct a Markov chain that applies an Ising
model formalism to such release channel clusters and realistic open channel
configurations to demonstrate that spark activation is described as a system
transition from a metastable to an absorbing state, analogous to the pressure
required to overcome surface tension in bubble formation. This yields
quantitative estimates of the spark generation probability as a function of
various system parameters. Our model of the release channel system yields
similar results for the sarcoplasmic reticulum Ca concentration threshold for
spark activation as previous experimental results, providing a mechanistic
explanation of the spark initiation. Additionally, we perform numerical
simulations to find spark probabilities as a function of sarcoplasmic reticulum
Ca concentration obtaining similar values for spark activation threshold as our
analytic model, as well as those reported in experimental studies.


  This review is about statistical genetics, an interdisciplinary topic between
statistical physics and population biology. The focus is on the phase of
quasi-linkage equilibrium (QLE). Our goals here are to clarify under which
conditions the QLE phase can be expected to hold in population biology and how
the stability of the QLE phase is lost. The QLE state, which has many
similarities to a thermal equilibrium state in statistical mechanics, was
discovered by M Kimura for a two-locus two-allele model, and was extended and
generalized to the global genome scale by (Neher and Shraiman, 2011). What we
will refer to as the Kimura-Neher-Shraiman (KNS) theory describes a population
evolving due to the mutations, recombination, natural selection and possibly
genetic drift. A QLE phase exists at sufficiently high recombination rate
and/or mutation rates with respect to selection strength. We show how in QLE it
is possible to infer the epistatic parameters of the fitness function from the
knowledge of the (dynamical) distribution of genotypes in a population. We
further consider the breakdown of the QLE regime for high enough selection
strength. We review recent results for the selection-mutation and
selection-recombination dynamics. Finally, we identify and characterize a new
phase which we call the non-random coexistence (NRC) where variability persists
in the population without either fixating or disappearing.


  Supertree methods are tree reconstruction techniques that combine several
smaller gene trees (possibly on different sets of species) to build a larger
species tree. The question of interest is whether the reconstructed supertree
converges to the true species tree as the number of gene trees increases (that
is, the consistency of supertree methods). In this paper, we are particularly
interested in the convergence rate of the maximum likelihood supertree.
Previous studies on the maximum likelihood supertree approach often formulate
the question of interest as a discrete problem and focus on reconstructing the
correct topology of the species tree. Aiming to reconstruct both the topology
and the branch lengths of the species tree, we propose an analytic approach for
analyzing the convergence of the maximum likelihood supertree method.
Specifically, we consider each tree as one point of a metric space and prove
that the distance between the maximum likelihood supertree and the species tree
converges to zero at a polynomial rate under some mild conditions. We further
verify these conditions for the popular exponential error model of gene trees.


  In this note we estimate herd immunity levels for the Covid-19 epidemic using
a standard SEIR system that models the disease dynamics in homogeneous
populations. The results obtained are indicative of values between 80% and 90%
for unprotected, fully susceptible populations. Basic protective measures such
as hand hygiene and mask wearing may be effective to lower herd immunity levels
down to values between 50% and 60% of the total population.


  Motivation: Quality control of genomic data is an essential but complicated
multi-step procedure, often requiring separate installation and expert
familiarity with a combination of disparate bioinformatics tools. Results: To
provide an automated solution that retains comprehensive quality checks and
flexible workflow architecture, we have developed snpQT, a scalable,
stand-alone software pipeline, offering some 36 discrete quality filters or
correction steps, with plots before-and-after user-modifiable thresholding.
This includes build conversion, population stratification against 1,000 Genomes
data, population outlier removal, and built-in imputation with its own pre- and
post- quality controls. Common input formats are used and users need not be
superusers nor have any prior coding experience. A comprehensive online
tutorial and installation guide is provided through to GWAS
(https://snpqt.readthedocs.io/en/latest/), introducing snpQT using a synthetic
demonstration dataset and a real-world Amyotrophic Lateral Sclerosis SNP-array
dataset. Availability: snpQT is open source and freely available at
https://github.com/nebfield/snpQT Contact: Vasilopoulou-C@ulster.ac.uk,
w.duddy@ulster.ac.uk


  Deep learning based methods have been widely applied to predict various kinds
of molecular properties in the pharmaceutical industry with increasingly more
success. Solvation free energy is an important index in the field of organic
synthesis, medicinal chemistry, drug delivery, and biological processes.
However, accurate solvation free energy determination is a time-consuming
experimental process. Furthermore, it could be useful to assess solvation free
energy in the absence of a physical sample. In this study, we propose two novel
models for the problem of free solvation energy predictions, based on the Graph
Neural Network (GNN) architectures: Message Passing Neural Network (MPNN) and
Graph Attention Network (GAT). GNNs are capable of summarizing the predictive
information of a molecule as low-dimensional features directly from its graph
structure without relying on an extensive amount of intra-molecular
descriptors. As a result, these models are capable of making accurate
predictions of the molecular properties without the time consuming process of
running an experiment on each molecule. We show that our proposed models
outperform all quantum mechanical and molecular dynamics methods in addition to
existing alternative machine learning based approaches in the task of solvation
free energy prediction. We believe such promising predictive models will be
applicable to enhancing the efficiency of the screening of drug molecules and
be a useful tool to promote the development of molecular pharmaceutics.


  Biomolecular condensates in cells are often rich in catalytically-active
enzymes. This is particularly true in the case of the large enzymatic complexes
known as metabolons, which contain different enzymes that participate in the
same catalytic pathway. One possible explanation for this self-organization is
the combination of the catalytic activity of the enzymes and a chemotactic
response to gradients of their substrate, which leads to a substrate-mediated
effective interaction between enzymes. These interactions constitute a purely
non-equilibrium effect and show exotic features such as non-reciprocity. Here,
we analytically study a model describing the phase separation of a mixture of
such catalytically-active particles. We show that a Michaelis-Menten-like
dependence of the particles' activities manifests itself as a screening of the
interactions, and that a mixture of two differently-sized active species can
exhibit phase separation with transient oscillations. We also derive a rich
stability phase diagram for a mixture of two species with both
concentration-dependent activity and size dispersity. This work highlights the
variety of possible phase separation behaviours in mixtures of
chemically-active particles, which provides an alternative pathway to the
passive interactions more commonly associated with phase separation in cells.
Our results highlight non-equilibrium organizing principles that can be
important for biologically relevant liquid-liquid phase separation.


  India has been hit by a huge second wave of Covid-19 that started in
mid-February 2021. Mumbai was amongst the first cities to see the increase. In
this report, we use our agent based simulator to computationally study the
second wave in Mumbai. We build upon our earlier analysis, where projections
were made from November 2020 onwards. We use our simulator to conduct an
extensive scenario analysis - we play out many plausible scenarios through
varying economic activity, reinfection levels, population compliance,
infectiveness, prevalence and lethality of the possible variant strains, and
infection spread via local trains to arrive at those that may better explain
the second wave fatality numbers. We observe and highlight that timings of peak
and valley of the fatalities in the second wave are robust to many plausible
scenarios, suggesting that they are likely to be accurate projections for
Mumbai. During the second wave, the observed fatalities were low in February
and mid-March and saw a phase change or a steep increase in the growth rate
after around late March. We conduct extensive experiments to replicate this
observed sharp convexity. This is not an easy phenomena to replicate, and we
find that explanations such as increased laxity in the population, increased
reinfections, increased intensity of infections in Mumbai transportation,
increased lethality in the virus, or a combination amongst them, generally do a
poor job of matching this pattern. We find that the most likely explanation is
presence of small amount of extremely infective variant on February 1 that
grows rapidly thereafter and becomes a dominant strain by Mid-March. From a
prescriptive view, this points to an urgent need for extensive and continuous
genome sequencing to establish existence and prevalence of different virus
strains in Mumbai and in India, as they evolve over time.


  Perturbations are ubiquitous in metabolism. A central tool to understand and
control their influence on metabolic networks is sensitivity analysis, which
investigates how the network responds to external perturbations. We follow here
a structural approach: the analysis is based on the network stoichiometry only
and it does not require any quantitative knowledge of the reaction rates. We
consider perturbations of reaction rates and metabolite concentrations, at
equilibrium, and we investigate the responses in the network. For general
metabolic systems, this paper focuses on the sign of the responses, i.e.
whether a response is positive, negative or whether its sign depends on the
parameters of the system. In particular, we identify and describe the
subnetworks that are the main players in the sign description. These
subnetworks are associated to certain kernel vectors of the stoichiometric
matrix and are thus independent from the chosen kinetics.


  The eukaryotic cell's cytoskeleton is a prototypical example of an active
material: objects embedded within it are driven by molecular motors acting on
the cytoskeleton, leading to anomalous diffusive behavior. Experiments tracking
the behavior of cell-attached objects have observed anomalous diffusion with a
distribution of displacements that is non-Gaussian, with heavy tails. This has
been attributed to "cytoquakes" or other spatially extended collective effects.
We show, using simulations and analytical theory, that a simple continuum
active gel model driven by fluctuating force dipoles naturally creates heavy
power-law tails in cytoskeletal displacements. We predict that this power law
exponent should depend on the geometry and dimensionality of where force
dipoles are distributed through the cell; we find qualitatively different
results for force dipoles in a 3D cytoskeleton and a quasi-two-dimensional
cortex. We then discuss potential applications of this model both in cells and
in synthetic active gels.


  Modeling and simulation (M&S) has revolutionized the way modern engineered
products are designed, tested and evaluated. Yet modeling and simulation is
much less frequently applied in the study of natural biological and active
matter systems. Two of the greatest challenges in applying M\&S approaches to
natural biological systems are (1) difficulty in specifying a model, and
developing a simulation from it, and (2) tuning and optimizing parameters. Here
we address the first challenge directly, and develop a software library that
can lead to progress in the second.


  The adaptation of biological species to their environment depends on their
traits. When various biological processes occur (survival, reproduction,
migration, etc.), the trait distribution may change with respect to time and
space. In the context of invasions, when considering the evolution of a
heritable trait that encodes the dispersive ability of individuals, the trait
distribution develops a particular spatial structure that leads to the
acceleration of the front propagation. That phenomenon is known as spatial
sorting. Many biological examples can be cited like the bush cricket in
Britain, the cane toad invasion in Australia or the common myna one in South
Africa. Adopting this framework, recent mathematical studies have led to
highlight the influence of the reproductive mode on the front propagation.
Asexual populations have been shown to spread with an asymptotic rate of t 3/2
in a minimal reactiondiffusion model, whereas the analogous rate for sexual
populations is of t 5/4 (where t denotes the time). However, the precise
description of the behaviour of the front propagation in the sexual case is
still an open question. The aim of this paper is to give precise approximations
for large times of its position, as well as some features of the local trait
distribution at the front. To do so, we solve explicitly the asymptotic problem
derived formally. Numerical simulations are shown to confirm these
calculations.


  Genetic association data from national biobanks and large-scale association
studies have provided new prospects for understanding the genetic evolution of
complex traits and diseases in humans. In turn, genomes from ancient human
archaeological remains are now easier than ever to obtain, and provide a direct
window into changes in frequencies of trait-associated alleles in the past.
This has generated a new wave of studies aiming to analyse the genetic
component of traits in historic and prehistoric times using ancient DNA, and to
determine whether any such traits were subject to natural selection. In humans,
however, issues about the portability and robustness of complex trait inference
across different populations are particularly concerning when predictions are
extended to individuals that died thousands of years ago, and for which little,
if any, phenotypic validation is possible. In this review, we discuss the
advantages of incorporating ancient genomes into studies of trait-associated
variants, the need for models that can better accommodate ancient genomes into
quantitative genetic frameworks, and the existing limits to inferences about
complex trait evolution, particularly with respect to past populations.


  This chapter reviews how molecular dynamics simulations, experimental data,
and Markov state models can synergize to map-out the mechanism of
protein-protein association and dissociation.


  Finding the set of the n items most dissimilar from each other out of a
larger population becomes increasingly difficult and computationally expensive
as either n or the population size grows large. Finding the set of the n most
dissimilar items is different than simply sorting an array of numbers because
there exists a pairwise relationship between each item and all other items in
the population. For instance, if you have a set of the most dissimilar n=4
items, one or more of the items from n=4 might not be in the set n=5. An exact
solution would have to search all possible combinations of size n in the
population, exhaustively. We present an open-source software called similarity
downselection (SDS), written in Python and freely available on GitHub. SDS
implements a heuristic algorithm for quickly finding the approximate set(s) of
the n most dissimilar items. We benchmark SDS against a Monte Carlo method,
which attempts to find the exact solution through repeated random sampling. We
show that for SDS to find the set of n most dissimilar conformers, our method
is not only orders of magnitude faster, but is also more accurate than running
the Monte Carlo for 1,000,000 iterations, each searching for set sizes n=3-7
out of a population of 50,000. We also benchmark SDS against the exact solution
for example small populations, showing SDS produces a solution close to the
exact solution in these instances.


  With recent advances in sequencing technology it has become affordable and
practical to sequence genomes to very high depth-of-coverage, allowing
researchers to discover low-frequency variants in the genome. However, due to
the errors in sequencing it is an active area of research to develop algorithms
that can separate noise from the true variants. LoFreq is a state of the art
algorithm for low-frequency variant detection but has a relatively long runtime
compared to other tools. In addition to this, the interface for running in
parallel could be simplified, allowing for multithreading as well as
distributing jobs to a cluster. In this work we describe some specific
contributions to LoFreq that remedy these issues.


  Toxicokinetic-toxicodynamic (TKTD) modelling is essential to make sense of
the time dependence of toxic effects, and to interpret and predict consequences
of time-varying exposure. These advantages have been recognised in the
regulatory arena, especially for environmental risk assessment (ERA) of
pesticides, where time-varying exposure is the norm. We critically evaluate the
link between the modelled variables in TKTD models and the observations from
laboratory ecotoxicity tests. For the endpoint reproduction, this link is far
from trivial. The relevant TKTD models for sub-lethal effects are based on
Dynamic-Energy Budget (DEB) theory, which specifies a continuous investment
flux into reproduction. In contrast, experimental tests score egg or offspring
release by the mother. The link between model and data is particularly
troublesome when a species reproduces in discrete clutches, and even more so
when eggs are incubated in the mother's brood pouch (and release of neonates is
scored in the test). This situation is quite common among aquatic invertebrates
(e.g., cladocerans, amphipods, mysids), including many popular test species. We
discuss these and other issues with reproduction data, reflect on their
potential impact on DEB-TKTD analysis, and provide preliminary recommendations
to correct them. Both modellers and users of model results need to be aware of
these complications, as ignoring them could easily lead to unnecessary failure
of DEB-TKTD models during calibration, or when validating them against
independent data for other exposure scenarios.


  Motivation: The activity of the adaptive immune system is governed by T-cells
and their specific T-cell receptors (TCR), which selectively recognize foreign
antigens. Recent advances in experimental techniques have enabled sequencing of
TCRs and their antigenic targets (epitopes), allowing to research the missing
link between TCR sequence and epitope binding specificity. Scarcity of data and
a large sequence space make this task challenging, and to date only models
limited to a small set of epitopes have achieved good performance. Here, we
establish a k-nearest-neighbor (K-NN) classifier as a strong baseline and then
propose TITAN (Tcr epITope bimodal Attention Networks), a bimodal neural
network that explicitly encodes both TCR sequences and epitopes to enable the
independent study of generalization capabilities to unseen TCRs and/or
epitopes. Results: By encoding epitopes at the atomic level with SMILES
sequences, we leverage transfer learning and data augmentation to enrich the
input data space and boost performance. TITAN achieves high performance in the
prediction of specificity of unseen TCRs (ROC-AUC 0.87 in 10-fold CV) and
surpasses the results of the current state-of-the-art (ImRex) by a large
margin. Notably, our Levenshtein-distance-based K-NN classifier also exhibits
competitive performance on unseen TCRs. While the generalization to unseen
epitopes remains challenging, we report two major breakthroughs. First, by
dissecting the attention heatmaps, we demonstrate that the sparsity of
available epitope data favors an implicit treatment of epitopes as classes.
This may be a general problem that limits unseen epitope performance for
sufficiently complex models. Second, we show that TITAN nevertheless exhibits
significantly improved performance on unseen epitopes and is capable of
focusing attention on chemically meaningful molecular structures.


  Addressed in this work is the performance of five popular algorithms, which
aim at assessing the dissemination dynamics of the COVID-19 disease on the
basis of the time series of new confirmed cases. The tests are based on
simulated data, generated by means of a deterministic compartmental
epidemiological model \cite{Matsinos2020a}, adapted herein to also include the
possibility of the loss of immunity by the group of the recovered (or
vaccinated) subjects. Assuming a simple temporal dependence of the effective
reproduction number (the exact details are of no relevance as far as the
conclusions of this work are concerned), time series of new cases were
generated in a time domain of nearly one year for the five top-ranking
countries in the cumulative number of infections by January 1, 2021. These
countries are (in descending order of infections): the United States of
America, India, Brazil, Russia, and the United Kingdom. The processing of each
simulated time series led to the establishment of relations between the input
(actual) and the reconstructed values of the effective reproduction number for
each country and algorithm, separately; this work argues that all five
algorithms underestimate the effective reproduction number when the latter
exceeds the critical value of $1$. The five algorithms were subsequently
applied to the real-life time series of new cases for the aforementioned five
countries, which also span a temporal interval of nearly one year; corrected
values of the effective reproduction number are obtained for these countries in
2020.


  We provide a characterisation of the continuous-time Markov models where the
Markov matrices from the model can be parameterised directly in terms of the
associated rate matrices (generators). That is, each Markov matrix can be
expressed as the sum of the identity matrix and a rate matrix from the model.
We show that the existence of an underlying Jordan algebra provides a
sufficient condition, which becomes necessary for (so-called) linear models. We
connect this property to the well-known uniformization procedure for
continuous-time Markov chains by demonstrating that the property is equivalent
to all Markov matrices from the model taking the same form as the corresponding
discrete time Markov matrices in the uniformized process. We apply our results
to analyse two model hierarchies practically important to phylogenetic
inference, obtained by assuming (i) time-reversibility and (ii) permutation
symmetry, respectively.


  Information on protein-protein interactions (PPIs) not only advances our
understanding of molecular biology but also provides important clues for target
selection in drug discovery and the design of PPI inhibitors. One of the
techniques used for computational prediction of PPIs is protein-protein docking
calculations, and a variety of software has been developed. However, a friendly
interface for users who are not sufficiently familiar with the command line
interface has not been developed so far. In this study, we have developed a
graphical user interface, MEGADOCK-GUI, which enables users to easily predict
PPIs and protein complex structures. In addition to the original 3-D molecular
viewer and input file preparation functions, MEGADOCK-GUI is software that can
automatically perform complete cross-docking of $M$ vs. $N$ proteins. With
MEGADOCK-GUI, various applications related to the prediction of PPIs, such as
ensemble docking that handles multiple conformations of proteins and screening
of binding partner proteins that bind to specific proteins, can now be easily
performed.


  Human Papillomaviruses (HPV) are one of the most prevalent sexually
transmitted infections (STI) and the most oncogenic viruses known to humans.
The vast majority of HPV infections clear in less than 3 years, but the
underlying mechanisms, especially the involvement of the immune response, are
still poorly known. Building on earlier work stressing the potential importance
of stochasticity in HPV infection clearance, we develop a stochastic
mathematical model of HPV dynamics with an explicit description of the
intracellular level. We show that the random partitioning of virus episomes
upon stem cell division and the occurrence of symmetric divisions greatly
affect viral persistence. These results call for more detailed within-host
studies to better understand the relative importance of stochasticity and
immunity in HPV infection clearance.


  We study a fundamental problem in computational chemistry known as molecular
conformation generation, trying to predict stable 3D structures from 2D
molecular graphs. Existing machine learning approaches usually first predict
distances between atoms and then generate a 3D structure satisfying the
distances, where noise in predicted distances may induce extra errors during 3D
coordinate generation. Inspired by the traditional force field methods for
molecular dynamics simulation, in this paper, we propose a novel approach
called ConfGF by directly estimating the gradient fields of the log density of
atomic coordinates. The estimated gradient fields allow directly generating
stable conformations via Langevin dynamics. However, the problem is very
challenging as the gradient fields are roto-translation equivariant. We notice
that estimating the gradient fields of atomic coordinates can be translated to
estimating the gradient fields of interatomic distances, and hence develop a
novel algorithm based on recent score-based generative models to effectively
estimate these gradients. Experimental results across multiple tasks show that
ConfGF outperforms previous state-of-the-art baselines by a significant margin.


  Repeated games are useful models to analyze long term interactions of living
species and complex social phenomena. Zero-determinant (ZD) strategies in
repeated games discovered by Press and Dyson in 2012 enforce a linear payoff
relationship between a focal player and the opponent. This linear relationship
can be set arbitrarily by a ZD player. Hence, a subclass of ZD strategies can
fix the opponent's expected payoff and another subclass of the strategies can
exceed the opponent for the expected payoff. Since this discovery, theories for
ZD strategies are extended to cope with various natural situations. It is
especially important to consider the theory of ZD strategies for repeated games
with a discount factor and observation errors because it allows the theory to
be applicable in the real world. Recent studies revealed their existence of ZD
strategies even in repeated games with both factors. However, the conditions
for the existence has not been sufficiently analyzed. Here, we mathematically
analyzed the conditions in repeated games with both factors. First, we derived
the thresholds of a discount factor and observation errors which ensure the
existence of Equalizer and positively correlated ZD (pcZD) strategies, which
are well-known subclasses of ZD strategies. We found that ZD strategies exist
only when a discount factor remains high as the error rates increase. Next, we
derived the conditions for the expected payoff of the opponent enforced by
Equalizer as well as the conditions for the slope and base line payoff of
linear lines enforced by pcZD. As a result, we found that, as error rates
increase or a discount factor decreases, the conditions for the linear line
that Equalizer or pcZD can enforce become strict.


  Understanding the dynamics of complex systems is a central task in many
different areas ranging from biology via epidemics to economics and
engineering. Unexpected behaviour of dynamic systems or even system failure is
sometimes difficult to comprehend. Such a data-mismatch can be caused by
endogenous model errors including misspecified interactions and inaccurate
parameter values. These are often difficult to distinguish from unmodelled
process influencing the real system like unknown inputs or faults. Localizing
the root cause of these errors or faults and reconstructing their dynamics is
only possible if the measured outputs of the system are sufficiently
informative.
  Here, we present criteria for the measurements required to localize the
position of error sources in large dynamic networks. We assume that faults or
errors occur at a limited number of positions in the network. This invariable
sparsity differs from previous sparsity definitions for inputs to dynamic
systems. We provide an exact criterion for the recovery of invariable sparse
inputs to nonlinear systems and formulate an optimization criterion for
invariable sparse input reconstruction. For linear systems we can provide exact
error bounds for this reconstruction method.


  This paper presents maplet, an open-source R package for the creation of
highly customizable, fully reproducible statistical pipelines for omics data
analysis, with a special focus on metabolomics-based methods. It builds on the
SummarizedExperiment data structure to create a centralized pipeline framework
for storing data, analysis steps, results, and visualizations. maplet's key
design feature is its modularity, which offers several advantages, such as
ensuring code quality through the individual maintenance of functions and
promoting collaborative development by removing technical barriers to code
contribution. With over 90 functions, the package includes a wide range of
functionalities, covering many widely used statistical approaches and data
visualization techniques.


  It is standard in chemistry to represent a sequence of reactions by a single
overall reaction, often called a complex reaction in contrast to an elementary
reaction. Photosynthesis $6 \text{CO}_2+6 \text{H}_2\text{O} \to \
\text{C}_6\text{H}_{12}\text{O}_6$ $+\ 6 \text{O}_2$ is an example of such
complex reaction. We introduce a mathematical operation that corresponds to
summing two chemical reactions. Specifically, we define an associative and
non-communicative operation on the product space $\mathbb{N}_0^n\times
\mathbb{N}_0^n$ (representing the reactant and the product of a chemical
reaction, respectively). The operation models the overall effect of two
reactions happening in succession, one after the other. We study the algebraic
properties of the operation and apply the results to stochastic reaction
networks, in particular to reachability of states, and to reduction of reaction
networks.


  In this paper, we estimate epidemiological data of the COVID-19 pandemic in
Hungary using only the daily number of hospitalized patients, and applying
well-known techniques from systems and control theory. We use a previously
published and validated compartmental model for the description of epidemic
spread. Exploiting the fact that an important subsystem of the model is linear,
first we compute the number of latent infected persons in time. Then an
estimate can be given for the number of people in other compartments. From
these data, it is possible to track the time dependent reproduction numbers via
a recursive least squares estimate. The credibility of the obtained results is
discussed using available data from the literature.


  We introduce ppsim, a software package for efficiently simulating population
protocols, a widely-studied subclass of chemical reaction networks (CRNs) in
which all reactions have two reactants and two products. Each step in the
dynamics involves picking a uniform random pair from a population of $n$
molecules to collide and have a (potentially null) reaction. In a recent
breakthrough, Berenbrink, Hammer, Kaaser, Meyer, Penschuck, and Tran [ESA 2020]
discovered a population protocol simulation algorithm quadratically faster than
the naive algorithm, simulating $\Theta(\sqrt{n})$ reactions in *constant* time
(independently of $n$, though the time scales with the number of species),
while preserving the *exact* stochastic dynamics.
  ppsim implements this algorithm, with a tightly optimized Cython
implementation that can exactly simulate hundreds of billions of reactions in
seconds. It dynamically switches to the CRN Gillespie algorithm for efficiency
gains when the number of applicable reactions in a configuration becomes small.
As a Python library, ppsim also includes many useful tools for data
visualization in Jupyter notebooks, allowing robust visualization of time
dynamics such as histogram plots at time snapshots and averaging repeated
trials.
  Finally, we give a framework that takes any CRN with only bimolecular (2
reactant, 2 product) or unimolecular (1 reactant, 1 product) reactions, with
arbitrary rate constants, and compiles it into a continuous-time population
protocol. This lets ppsim exactly sample from the chemical master equation
(unlike approximate heuristics such as tau-leaping or LNA), while achieving
asymptotic gains in running time. In linked Jupyter notebooks, we demonstrate
the efficacy of the tool on some protocols of interest in molecular
programming, including the approximate majority CRN and CRN models of DNA
strand displacement reactions.


  Cold atmospheric plasma (CAP) and plasma-treated liquids (PTLs) have recently
become a promising option for cancer treatment, but the underlying mechanisms
of the anti-cancer effect are still to a large extent unknown. Although
hydrogen peroxide has been recognized as the major anti-cancer agent of PTL and
may enable selectivity in a certain concentration regime, the co-existence of
nitrite can create a synergistic effect. We develop a mathematical model to
describe the key species and features of the cellular response towards PTL.
From the numerical solutions, we define a number of dependent variables, which
represent feasible measures to quantify cell susceptibility in terms of the
hydrogen peroxide membrane diffusion rate constant and the intracellular
catalase concentration. For each of these dependent variables, we investigate
the regimes of selective versus non-selective, and of synergistic versus
non-synergistic effect to evaluate their potential role as a measure of cell
susceptibility. Our results suggest that the maximal intracellular hydrogen
peroxide concentration, which in the selective regime is almost four times
greater for the most susceptible cells compared to the most resistant cells,
could be used to quantify the cell susceptibility towards exogenous hydrogen
peroxide. We believe our theoretical approach brings novelty to the field of
plasma oncology, and more broadly, to the field of redox biology, by proposing
new ways to quantify the selective and synergistic anti-cancer effect of PTL in
terms of inherent cell features.


  When an epidemic spreads into a population, it is often unpractical or
impossible to have a continuous monitoring of all subjects involved. As an
alternative, algorithmic solutions can be used to infer the state of the whole
population from a limited amount of measures. We analyze the capability of deep
neural networks to solve this challenging task. Our proposed architecture is
based on Graph Convolutional Neural Networks. As such it can reason on the
effect of the underlying social network structure, which is recognized as the
main component in the spreading of an epidemic. We test the proposed
architecture with two scenarios modeled on the CoVid-19 pandemic: a generic
homogeneous population, and a toy model of Boston metropolitan area.


  Motivation: Protein-ligand affinity prediction is an important part of
structure-based drug design. It includes molecular docking and affinity
prediction. Although molecular dynamics can predict affinity with high accuracy
at present, it is not suitable for large-scale virtual screening. The existing
affinity prediction and evaluation functions based on deep learning mostly rely
on experimentally-determined conformations. Results: We build a predictive
model of protein-ligand affinity through the ResNet neural network with added
attention mechanism. The resulting ResAtom-Score model achieves Pearson's
correlation coefficient R = 0.833 on the CASF-2016 benchmark test set. At the
same time, we evaluated the performance of a variety of existing scoring
functions in combination with ResAtom-Score in the absence of
experimentally-determined conformations. The results show that the use of
{\Delta}VinaRF20 in combination with ResAtom-Score can achieve affinity
prediction close to scoring functions in the presence of
experimentally-determined conformations. These results suggest that ResAtom
system may be used for in silico screening of small molecule ligands with
target proteins in the future. Availability: https://github.com/wyji001/ResAtom


  To generate drug molecules of desired properties with computational methods
is the holy grail in pharmaceutical research. Here we describe an AI strategy,
retro drug design, or RDD, to generate novel small molecule drugs from scratch
to meet predefined requirements, including but not limited to biological
activity against a drug target, and optimal range of physicochemical and ADMET
properties. Traditional predictive models were first trained over experimental
data for the target properties, using an atom typing based molecular descriptor
system, ATP. Monte Carlo sampling algorithm was then utilized to find the
solutions in the ATP space defined by the target properties, and the deep
learning model of Seq2Seq was employed to decode molecular structures from the
solutions. To test feasibility of the algorithm, we challenged RDD to generate
novel drugs that can activate {\mu} opioid receptor (MOR) and penetrate blood
brain barrier (BBB). Starting from vectors of random numbers, RDD generated
180,000 chemical structures, of which 78% were chemically valid. About 42,000
(31%) of the valid structures fell into the property space defined by MOR
activity and BBB permeability. Out of the 42,000 structures, only 267 chemicals
were commercially available, indicating a high extent of novelty of the
AI-generated compounds. We purchased and assayed 96 compounds, and 25 of which
were found to be MOR agonists. These compounds also have excellent BBB scores.
The results presented in this paper illustrate that RDD has potential to
revolutionize the current drug discovery process and create novel structures
with multiple desired properties, including biological functions and ADMET
properties. Availability of an AI-enabled fast track in drug discovery is
essential to cope with emergent public health threat, such as pandemic of
COVID-19.


  This paper presents the methods that have participated in the SHREC 2021
contest on retrieval and classification of protein surfaces on the basis of
their geometry and physicochemical properties. The goal of the contest is to
assess the capability of different computational approaches to identify
different conformations of the same protein, or the presence of common
sub-parts, starting from a set of molecular surfaces. We addressed two
problems: defining the similarity solely based on the surface geometry or with
the inclusion of physicochemical information, such as electrostatic potential,
amino acid hydrophobicity, and the presence of hydrogen bond donors and
acceptors. Retrieval and classification performances, with respect to the
single protein or the existence of common sub-sequences, are analysed according
to a number of information retrieval indicators.


  Chemical reaction networks (CRNs) are fundamental computational models used
to study the behavior of chemical reactions in well-mixed solutions. They have
been used extensively to model a broad range of biological systems, and are
primarily used when the more traditional model of deterministic continuous mass
action kinetics is invalid due to small molecular counts. We present a perfect
sampling algorithm to draw error-free samples from the stationary distributions
of stochastic models for coupled, linear chemical reaction networks. The state
spaces of such networks are given by all permissible combinations of molecular
counts for each chemical species, and thereby grow exponentially with the
numbers of species in the network. To avoid simulations involving large numbers
of states, we propose a subset of chemical species such that coupling of paths
started from these states guarantee coupling of paths started from all states
in the state space and we show for the well-known Reversible Michaelis-Menten
model that the subset does in fact guarantee perfect draws from the stationary
distribution of interest. We compare solutions computed in two ways with this
algorithm to those found analytically using the chemical master equation and we
compare the distribution of coupling times for the two simulation approaches.


  We introduce a minimalist outbreak forecasting model that combines
data-driven parameter estimation with variational data assimilation. By
focusing on the fundamental components of nonlinear disease transmission and
representing data in a domain where model stochasticity simplifies into a
process with independent increments, we design an approach that only requires
four core parameters to be estimated. We illustrate this novel methodology on
COVID-19 forecasts. Results include case count and deaths predictions for the
US and all of its 50 states, the District of Columbia, and Puerto Rico. The
method is computationally efficient and is not disease- or location-specific.
It may therefore be applied to other outbreaks or other countries, provided
case counts and/or deaths data are available.


  MomentClosure.jl is a Julia package providing automated derivation of the
time-evolution equations of the moments of molecule numbers for virtually any
chemical reaction network using a wide range of moment closure approximations.
It extends the capabilities of modelling stochastic biochemical systems in
Julia and can be particularly useful when exact analytic solutions of the
chemical master equation are unavailable and when Monte Carlo simulations are
computationally expensive.
  MomentClosure.jl is freely accessible under the MIT license. Source code and
documentation are available at https://github.com/augustinas1/MomentClosure.jl


  We propose a minimal model to simulate long waiting times followed by
evolutionary bursts on rugged landscapes. It combines point and inversions-like
mutations as sources of genetic variation. The inversions are intended to
simulate one of the main chromosomal rearrangements. Using the well-known
family of NK fitness landscapes, we simulate random adaptive walks, i.e.
successive mutational events constrained to incremental fitness selection. We
report the emergence of different time scales: a short-term dynamics mainly
driven by point mutations, followed by a long-term (stasis-like) waiting period
until a new mutation arises. This new mutation is an inversion which can
trigger a burst of successive point mutations, and then drives the system to
new short-term increasing-fitness period. We analyse the effect of genes
epistatic interactions on the evolutionary time scales. We suggest that the
present model mimics the process of evolutionary innovation and punctuated
equilibrium.


  Multiple instances of Zika virus epidemic have been reported around the world
in the last two decades, turning the related illness into an international
concern. In this context the use of mathematical models for epidemics is of
great importance, since they are useful tools to study the underlying outbreak
numbers and allow one to test the effectiveness of different strategies used to
combat the associated diseases. This work deals with the development and
calibration of an epidemic model to describe the 2016 outbreak of Zika virus in
Brazil. A system of 8 differential equations with 8 parameters is employed to
model the evolution of the infection through two populations. Nominal values
for the model parameters are estimated from the literature. An inverse problem
is formulated and solved by comparing the system response to real data from the
outbreak. The calibrated results presents realistic parameters and returns
reasonable descriptions, with the curve shape similar to the outbreak evolution
and peak value close to the highest number of infected people during 2016.
Considerations about the lack of data for some initial conditions are also made
through an analysis over the response behavior according to their change in
value.


  Across many domains of interaction, both natural and artificial, individuals
use past experience to shape future behaviors. The results of such learning
processes depend on what individuals wish to maximize. A natural objective is
one's own success. However, when two such "selfish" learners interact with each
other, the outcome can be detrimental to both, especially when there are
conflicts of interest. Here, we explore how a learner can align incentives with
a selfish opponent. Moreover, we consider the dynamics that arise when learning
rules themselves are subject to evolutionary pressure. By combining extensive
simulations and analytical techniques, we demonstrate that selfish learning is
unstable in most classical two-player repeated games. If evolution operates on
the level of long-run payoffs, selection instead favors learning rules that
incorporate social (other-regarding) preferences. To further corroborate these
results, we analyze data from a repeated prisoner's dilemma experiment. We find
that selfish learning is insufficient to explain human behavior when there is a
trade-off between payoff maximization and fairness.


  Free energy calculations are rapidly becoming indispensable in
structure-enabled drug discovery programs. As new methods, force fields, and
implementations are developed, assessing their expected accuracy on real-world
systems (benchmarking) becomes critical to provide users with an assessment of
the accuracy expected when these methods are applied within their domain of
applicability, and developers with a way to assess the expected impact of new
methodologies. These assessments require construction of a benchmark - a set of
well-prepared, high quality systems with corresponding experimental
measurements designed to ensure the resulting calculations provide a realistic
assessment of expected performance when these methods are deployed within their
domains of applicability. To date, the community has not yet adopted a common
standardized benchmark, and existing benchmark reports suffer from a myriad of
issues, including poor data quality, limited statistical power, and
statistically deficient analyses, all of which can conspire to produce
benchmarks that are poorly predictive of real-world performance. Here, we
address these issues by presenting guidelines for (1) curating experimental
data to develop meaningful benchmark sets, (2) preparing benchmark inputs
according to best practices to facilitate widespread adoption, and (3) analysis
of the resulting predictions to enable statistically meaningful comparisons
among methods and force fields.


  The debate on the role of school closures as a mitigation strategy against
the spread of Covid-19 is gaining relevance due to emerging variants in Europe.
According to WHO, decisions on schools "should be guided by a risk-based
approach". However, risk evaluation requires sound methods, transparent data
and careful consideration of the context at the local level. We review a recent
study by Gandini et al., on the role of school opening as a driver of the
second COVID-19 wave in Italy, which concluded that there was no connection
between school openings/closures and SARS-CoV-2 incidence. infections. This
analysis has been widely commented in Italian media as conclusive proof that
"schools are safe". However the study presents severe oversights and careless
interpretation of data.


  We construct an agent-based SEIR model to simulate COVID-19 spread at a
16000-student mostly non-residential urban university during the Fall 2021
Semester. We find that mRNA vaccine coverage above 80% makes it possible to
safely reopen to in-person instruction. If vaccine coverage is 100%, then our
model indicates that facemask use is not necessary. Our simulations with
vaccine coverage below 70% exhibit a right-skew for total infections over the
semester, which suggests that high levels of infection are not exceedingly rare
with campus social connections the main transmission route. Less effective
vaccines or incidence of new variants may require additional intervention such
as screening testing to reopen safely.


  Context - Prioritising is likely to differ depending on the species
considered for connectivity assessments, leading to a lack of consensual
decisions for territorial planning.
  Objectives - The objective was to assess the relevance of identifying
priority areas for connectivity for groups of species based on common dispersal
abilities. We aimed to assess the impact of target groups choices on predicted
priority areas.
  Method - The study was located at the Thau Lagoon territory to demonstrate
the methodological approach. Ecological niche modelling was used to quantify
species resistance and to identify suitable habitat patches. We coupled the
least-cost path methodology with circuit theory to assess species connectivity.
We classified connectivity from high to low levels and averaged the results by
dispersal groups.
  Results - We found important differences in identified priority areas between
groups with dissimilar dispersal abilities, with little overlap between highly
connected areas. We identified a gap between the level of protection of low
dispersal species and highly connected areas. We found mismatches between
existing corridors and connectivity in low dispersal species, and a greater
impact in areas of expected urban sprawl projects on favourably connected areas
for species with high dispersal capabilities.
  Conclusion - We have demonstrated that a diversity of dispersal capacity
ranges must be accounted for in order to identify ecological corridors in
programmes that aim to restore habitat connectivity at territorial levels. Our
findings are oriented to support the decisions of planning initiatives, at both
local and regional scale.


  Protein structural classification (PSC) is a supervised problem of assigning
proteins into pre-defined structural (e.g., CATH or SCOPe) classes based on the
proteins' sequence or 3D structural features. We recently proposed PSC
approaches that model protein 3D structures as protein structure networks
(PSNs) and analyze PSN-based protein features, which performed better than or
comparable to state-of-the-art sequence or other 3D structure-based approaches
in the task of PSC. However, existing PSN-based PSC approaches model the whole
3D structure of a protein as a static PSN. Because folding of a protein is a
dynamic process, where some parts of a protein fold before others, modeling the
3D structure of a protein as a dynamic PSN might further help improve the
existing PSC performance. Here, we propose for the first time a way to model 3D
structures of proteins as dynamic PSNs, with the hypothesis that this will
improve upon the current state-of-the-art PSC approaches that are based on
static PSNs (and thus upon the existing state-of-the-art sequence and other 3D
structural approaches). Indeed, we confirm this on 71 datasets spanning ~44,000
protein domains from CATH and SCOPe


  Apoptosis is a mechanism of programmed cell death in which cells engage in a
controlled demolition and prepare to be digested without damaging their
environment. In normal conditions apoptosis is repressed, until it is
irreversibly induced by an appropriate signal. In adult organisms apoptosis is
a natural way to dispose of damaged cells, and its disruption or excess is
associated with cancer and autoimmune diseases. Apoptosis is regulated by a
complex signaling network controlled by caspases, specialized enzymes that
digest essential cellular components and promote the degradation of genomic
DNA. In this work we propose an effective description of the signaling network
focused on caspase-3 as a readout of cell fate. We integrate intermediate
network interactions into a nonlinear feedback function acting on caspase-3 and
introduce the effect of pro-apoptotic stimuli and regulatory elements as a
saturating activation function. We find that the theory has a robust bistable
regime where two possible states coexist, representing survival and cell death
fates. For a broad range of parameters, strong stimuli can induce an
irreversible switch to the death fate. We use the theory to explore dynamical
stimulation conditions and determine how cell fate depends on different stimuli
patterns. This analysis reveals a critical relation between transient stimuli
intensity and duration to trigger irreversible apoptosis.


  Local clock models propose that the rate of molecular evolution is constant
within phylogenetic sub-trees. Current local clock inference procedures scale
poorly to large taxa problems, impose model misspecification, or require a
priori knowledge of the existence and location of clocks. To overcome these
challenges, we present an autocorrelated, Bayesian model of heritable clock
rate evolution that leverages heavy-tailed priors with mean zero to shrink
increments of change between branch-specific clocks. We further develop an
efficient Hamiltonian Monte Carlo sampler that exploits closed form gradient
computations to scale our model to large trees. Inference under our
shrinkage-clock exhibits an over 3-fold speed increase compared to the popular
random local clock when estimating branch-specific clock rates on a simulated
dataset. We further show our shrinkage-clock recovers known local clocks within
a rodent and mammalian phylogeny. Finally, in a problem that once appeared
computationally impractical, we investigate the heritable clock structure of
various surface glycoproteins of influenza A virus in the absence of prior
knowledge about clock placement.


  A statistical method is developed to estimate the maximum amplitude of the
base pair fluctuations in a three dimensional mesoscopic model for nucleic
acids. The base pair thermal vibrations around the helix diameter are viewed as
a Brownian motion for a particle embedded in a stable helical structure. The
probability to return to the initial position is computed, as a function of
time, by integrating over the particle paths consistent with the physical
properties of the model potential. The zero time condition for the
first-passage probability defines the constraint to select the integral cutoff
for various macroscopic helical conformations, obtained by tuning the twist,
the bending and the slide motion between adjacent base pairs along the molecule
stack. Applying the method to a short homogeneous chain at room temperature, we
obtain meaningful estimates for the maximum fluctuations in the twist
conformation with $\sim 10.5$ base pairs per helix turn, typical of double
stranded DNA helices. Untwisting the double helix, the base pair fluctuations
broaden and the integral cutoff grows. The cutoff is found to increase also in
the presence of a sliding motion which shortens the helix contour length, a
situation peculiar of dsRNA molecules.


  Predicting molecular conformations (or 3D structures) from molecular graphs
is a fundamental problem in many applications. Most existing approaches are
usually divided into two steps by first predicting the distances between atoms
and then generating a 3D structure through optimizing a distance geometry
problem. However, the distances predicted with such two-stage approaches may
not be able to consistently preserve the geometry of local atomic
neighborhoods, making the generated structures unsatisfying. In this paper, we
propose an end-to-end solution for molecular conformation prediction called
ConfVAE based on the conditional variational autoencoder framework.
Specifically, the molecular graph is first encoded in a latent space, and then
the 3D structures are generated by solving a principled bilevel optimization
program. Extensive experiments on several benchmark data sets prove the
effectiveness of our proposed approach over existing state-of-the-art
approaches. Code is available at
\url{https://github.com/MinkaiXu/ConfVAE-ICML21}.


  Head and Neck Squamous Cell Carcinoma (HNSCC) is one of cancer type that is
most distressing leading to acute pain, effecting speech and primary survival
functions such as swallowing and breathing. The morbidity and mortality of
HNSCC patients have not significantly improved even tough there has been
advancement in surgical and radiotherapy treatments. The high mortality may be
attributed to the complexity and significant changes in the clinical outcomes.
Therefore, it is important to increase the accuracy of predicting the outcome
of cancer survival. Few cancer survival prediction models of HNSCC have been
proposed so far. In this study, genomic data (whole exome sequencing) are
integrated with clinical data to improve the performance of prediction model.
The somatic mutations of every patient is processed using Multifractal
Deterended Fluctuation Analysis (MFDFA) algorithm and the parameter values of
Fractal Dimension (Dq) is included along with clinical data for cancer survival
prediction. Feature ranking proves that the new engineered feature is one of
the important feature in prediction model. In order to improve the performance
index of models, hyperparameters were also tuned in all the classifiers
considered. 10-Fold cross validation is implemented and XGBoost (98% AUROC, 94%
precision, and 93% recall) proves to be best model classifier followed by
Random Forest 93% AUROC, 93% precision, and 93% recall), Support Vector Machine
(84% AUCROC, 79% precision, and 79% recall) and Logistic Regression (80% AUROC,
77% precision, and 76% recall).


  The potential of deep learning has been recognized in the protein structure
prediction community for some time, and became indisputable after CASP13. In
CASP14, deep learning has boosted the field to unanticipated levels reaching
near-experimental accuracy. This success comes from advances transferred from
other machine learning areas, as well as methods specifically designed to deal
with protein sequences and structures, and their abstractions. Novel emerging
approaches include (i) geometric learning, i.e. learning on representations
such as graphs, 3D Voronoi tessellations, and point clouds; (ii) pre-trained
protein language models leveraging attention; (iii) equivariant architectures
preserving the symmetry of 3D space; (iv) use of large meta-genome databases;
(v) combinations of protein representations; (vi) and finally truly end-to-end
architectures, i.e. differentiable models starting from a sequence and
returning a 3D structure. Here, we provide an overview and our opinion of the
novel deep learning approaches developed in the last two years and widely used
in CASP14.


  In nanopore sequencing, electrical signal is measured as DNA molecules pass
through the sequencing pores. Translating these signals into DNA bases (base
calling) is a highly non-trivial task, and its quality has a large impact on
the sequencing accuracy. The most successful nanopore base callers to date use
convolutional neural networks (CNN) to accomplish the task.
  Convolutional layers in CNNs are typically composed of filters with constant
window size, performing best in analysis of signals with uniform speed.
However, the speed of nanopore sequencing varies greatly both within reads and
between sequencing runs. Here, we present dynamic pooling, a novel neural
network component, which addresses this problem by adaptively adjusting the
pooling ratio. To demonstrate the usefulness of dynamic pooling, we developed
two base callers: Heron and Osprey. Heron improves the accuracy beyond the
experimental high-accuracy base caller Bonito developed by Oxford Nanopore.
Osprey is a fast base caller that can compete in accuracy with Guppy
high-accuracy mode, but does not require GPU acceleration and achieves a near
real-time speed on common desktop CPUs.
  Availability: https://github.com/fmfi-compbio/osprey,
https://github.com/fmfi-compbio/heron
  Keywords: nanopore sequencing, base calling, convolutional neural networks,
pooling


  Gene promoters are the key DNA regulatory elements positioned around the
transcription start sites and are responsible for regulating gene transcription
process. Various alignment-based, signal-based and content-based approaches are
reported for the prediction of promoters. However, since all promoter sequences
do not show explicit features, the prediction performance of these techniques
is poor. Therefore, many machine learning and deep learning models have been
proposed for promoter prediction. In this work, we studied methods for vector
encoding and promoter classification using genome sequences of three distinct
higher eukaryotes viz. yeast (Saccharomyces cerevisiae), A. thaliana (plant)
and human (Homo sapiens). We compared one-hot vector encoding method with
frequency-based tokenization (FBT) for data pre-processing on 1-D Convolutional
Neural Network (CNN) model. We found that FBT gives a shorter input dimension
reducing the training time without affecting the sensitivity and specificity of
classification. We employed the deep learning techniques, mainly CNN and
recurrent neural network with Long Short Term Memory (LSTM) and random forest
(RF) classifier for promoter classification at k-mer sizes of 2, 4 and 8. We
found CNN to be superior in classification of promoters from non-promoter
sequences (binary classification) as well as species-specific classification of
promoter sequences (multiclass classification). In summary, the contribution of
this work lies in the use of synthetic shuffled negative dataset and
frequency-based tokenization for pre-processing. This study provides a
comprehensive and generic framework for classification tasks in genomic
applications and can be extended to various classification problems.


  We present a new and versatile implementation of rapid and localized
immunohistochemical staining of tissue sections. Immunohistochemistry (IHC)
allows to detect specific proteins on tissue sections and comprises a sequence
of specific biochemical reactions. For the rapid implementation of IHC, we
fabricated horizontally oriented microfluidic probes (MFP) with functionally
designed apertures to enable square and circular footprints, which we employ to
locally expose a tissue to time-optimized sequences of different biochemicals.


  As Europe is experiencing a second violent CoVid-19 storm, with the PCR-based
testing system deteriorating due to the high volumes of people to be tested
daily, there is a general reconsideration of the mathematical theories at the
basis of our contact tracing and testing approaches. Drawing upon the concept
of super spreader, we propose the use of (less sensitive) rapid tests to detect
those secondary infections that do not need the use of PCRs, thus saving the
most part of PCR tests currently used. This before the system fails.


  We study the induction and stabilization of spiral structures for the cyclic
three-species stochastic May-Leonard model with asymmetric predation rates on a
spatially inhomogeneous two-dimensional toroidal lattice using Monte Carlo
simulations. In an isolated setting, strongly asymmetric predation rates lead
to rapid extinction from coexistence of all three species to a single surviving
population. Even for weakly asymmetric predation rates, only a fraction of
ecologies in a statistical ensemble manages to maintain full three-species
coexistence. However, when the asymmetric competing system is coupled via
diffusive proliferation to a fully symmetric May-Leonard patch, the stable
spiral patterns from this region induce transient plane-wave fronts and
ultimately quasi-stationary spiral patterns in the vulnerable asymmetric
region. Thus the endangered ecological subsystem may effectively become
stabilized through immigration from even a much smaller stable region. To
describe the stabilization of spiral population structures in the asymmetric
region, we compare the increase in the robustness of these topological defects
at extreme values of the asymmetric predation rates in the spatially coupled
system with the corresponding asymmetric May-Leonard model in isolation. We
delineate the quasi-stationary nature of coexistence induced in the asymmetric
subsystem by its diffusive coupling to a symmetric May-Leonard patch, and
propose a (semi-)quantitative criterion for the spiral oscillations to be
sustained in the asymmetric region.


  Inferring parameters related to the aggregation pattern of parasites and to
their dispersal propensity is important for predicting their ecological
consequences and evolutionary potential. Nonetheless, it is notoriously
difficult to infer these parameters from wildlife parasites given the
difficulty in tracking these organisms. Molecular-based inferences constitute a
promising approach that has yet rarely been applied in the wild.Here, we
combinedseveral population genetic analyses including sibship reconstruction to
documentthe genetic structure, patterns of sibship aggregation and the
dispersal dynamics of a non-native parasite of fish, the freshwater copepod
ectoparasiteTracheliastespolycolpus. We collected parasites according to a
hierarchical sampling design,withthe sampling of all parasites from all host
individualscapturedineight sites spread along an upstream-downstream river
gradient. Individual multilocus genotypes were obtained from 14 microsatellite
markers, and used to assign parasites to full-sib families and to investigate
the genetic structure of T.polycolpus among both hosts and sampling sites. The
distribution of full-sibs obtainedamong the sampling sites was used to estimate
individual dispersal distances within families. Our results showed that T.
polycolpus sibs tend to be aggregated within sites but not withinhost
individuals. We detected important upstream-to-downstream dispersal events of
T.polycolpusbetween sites (modal distance: 25.4 km; 95% CI [22.9, 27.7]),
becoming scarcer as the geographic distance from their family core location
increases. Such a dispersal pattern likely contributes to the strong
isolation-by-distance observed at the river scale. We also detected some
downstream-to-upstream dispersal events (modal distance: 2.6 km; 95% CI
[2.2-23.3]) that likely result from movementsof infected hosts.Within each
site, the dispersal of free-living infective larvae among hosts likely
contributes to increasing genetic diversity on hosts, possibly fostering the
evolutionary potential of T. polycolpus.


  Darwinian evolution is driven by random mutations, genetic recombination
(gene shuffling) and selection that favors genotypes with high fitness. For
systems where each genotype can be represented as a bitstring of length $L$, an
overview of possible evolutionary trajectories is provided by the $L$-cube
graph with nodes labeled by genotypes and edges directed toward the genotype
with higher fitness. Peaks (sinks in the graphs) are important since a
population can get stranded at a suboptimal peak. The fitness landscape is
defined by the fitness values of all genotypes in the system. Some notion of
curvature is necessary for a more complete analysis of the landscapes,
including the effect of recombination. The shape approach uses triangulations
(shapes) induced by fitness landscapes. The main topic for this work is the
interplay between peak patterns and shapes. Because of constraints on the
shapes for $L=3$ imposed by peaks, there are in total 25 possible combinations
of peak patterns and shapes. Similar constraints exist for higher $L$.
Specifically, we show that the constraints induced by the staircase
triangulation can be formulated as a condition of {\emph{universal positive
epistasis}}, an order relation on the fitness effects of arbitrary sets of
mutations that respects the inclusion relation between the corresponding
genetic backgrounds. We apply the concept to a large protein fitness landscape
for an immunoglobulin-binding protein expressed in Streptococcal bacteria.


  SIMD||DNA is a model of DNA strand displacement allowing parallel in-memory
computation on DNA storage. We show how to simulate an arbitrary 3-symbol
space-bounded Turing machine with a SIMD||DNA program, giving a more direct and
efficient route to general-purpose information manipulation on DNA storage than
the Rule 110 simulation of [Wang, Chalk, Soloveichik, DNA 2019]. We also
develop software that can simulate SIMD||DNA programs and produce SVG figures.


  In the pharmaceutical industry, where it is common to generate many QSAR
models with large numbers of molecules and descriptors, the best QSAR methods
are those that can generate the most accurate predictions but that are also
insensitive to hyperparameters and are computationally efficient. Here we
compare Light Gradient Boosting Machine (LightGBM) to random forest,
single-task deep neural nets, and Extreme Gradient Boosting (XGBoost) on 30
in-house data sets. While any boosting algorithm has many adjustable
hyperparameters, we can define a set of standard hyperparameters at which
LightGBM makes predictions about as accurate as single-task deep neural nets,
but is a factor of 1000-fold faster than random forest and ~4-fold faster than
XGBoost in terms of total computational time for the largest models. Another
very useful feature of LightGBM is that it includes a native method for
estimating prediction intervals.


  The SARS-CoV-2 spike (S) protein facilitates viral infection, and has been
the focus of many structure determination efforts. Its flexible loop regions
are known to be involved in protein binding and may adopt multiple
conformations. This paper identifies the S protein loops and studies their
conformational variability based on the available Protein Data Bank (PDB)
structures. While most loops had essentially one stable conformation, 17 of 44
loop regions were observed to be structurally variable with multiple
substantively distinct conformations based on a cluster analysis. Loop modeling
methods were then applied to the S protein loop targets, and the prediction
accuracies discussed in relation to the characteristics of the conformational
clusters identified. Loops with multiple conformations were found to be
challenging to model based on a single structural template.


  This article considers the optimal control of the SIR model with both
transmission and treatment uncertainty. It follows the model presented in Gatto
and Schellhorn (2021). We make four significant improvements on the latter
paper. First, we prove the existence of a solution to the model. Second, our
interpretation of the control is more realistic: while in Gatto and Schellhorn
the control $\alpha$ is the proportion of the population that takes a basic
dose of treatment, so that $\alpha >1$ occurs only if some patients take more
than a basic dose, in our paper, $\alpha$ is constrained between zero and one,
and represents thus the proportion of the population undergoing treatment.
Third, we provide a complete solution for the moderate infection regime (with
constant treatment). Finally, we give a thorough interpretation of the control
in the moderate infection regime, while Gatto and Schellhorn focussed on the
interpretation of the low infection regime. Finally, we compare the efficiency
of our control to curb the COVID-19 epidemic to other types of control.


  Advances in high-throughput sequencing technologies and bioinformatics
approaches over almost the last three decades have substantially increased our
ability to explore microorganisms and their functions-including those that have
yet to be cultivated in pure isolation. Genome-resolved metagenomic approaches
have enabled linking powerful functional predictions to specific taxonomical
groups with increasing fidelity. Additionally, whole community gene expression
surveys and metabolite profiling have permitted direct surveys of
community-scale functions in specific environmental settings. These advances
have allowed for a shift in microbiome science away from descriptive studies
and towards mechanistic and predictive frameworks for designing and harnessing
microbial communities for desired beneficial outcomes. Here, we review how
modern genome-resolved metagenomic approaches have been applied to a variety of
water engineering applications from lab-scale bioreactors to full-scale
systems. We describe integrated omics analysis across engineered water systems
and the foundations for pairing these insights with modeling approaches.
Lastly, we summarize emerging omics-based technologies that we believe will be
powerful tools for water engineering applications. Overall, we provide a
framework for microbial ecologists specializing in water engineering to apply
cutting-edge omics approaches to their research questions to achieve novel
functional insights. Successful adoption of predictive frameworks in engineered
water systems could enable more economically and environmentally sustainable
bioprocesses as demand for water and energy resources increases.


  We investigate the effect of degree correlation on a
susceptible-infected-susceptible (SIS) model with a nonlinear cooperative
effect (synergy) in infectious transmissions. In a mean-field treatment of the
synergistic SIS model on a bimodal network with tunable degree correlation, we
identify a discontinuous transition that is independent of the degree
correlation strength unless the synergy is absent or extremely weak. Regardless
of synergy (absent or present), a positive and negative degree correlation in
the model reduces and raises the epidemic threshold, respectively. For networks
with a strongly positive degree correlation, the mean-field treatment predicts
the emergence of two discontinuous jumps in the steady-state infected density.
To test the mean-field treatment, we provide approximate master equations of
the present model. We quantitatively confirm that the approximate master
equations agree with not only all qualitative predictions of the mean-field
treatment but also corresponding Monte-Carlo simulations.


  Cooperation among individuals has been key to sustaining societies. However,
natural selection favors defection over cooperation. Cooperation can be favored
when the mobility of individuals allows cooperators to form a cluster (or
group). Mobility patterns of animals sometimes follow a L\'evy flight. A L\'evy
flight is a kind of random walk but it is composed of many small movements with
a few big movements. The role of L\'evy flights for cooperation has been
studied by Antonioni and Tomassini. They showed that L\'evy flights promoted
cooperation combined with conditional movements triggered by neighboring
defectors. However, the optimal condition for neighboring defectors and how the
condition changes by the intensity of L\'evy flights are still unclear. Here,
we developed an agent-based model in a square lattice where agents perform
L\'evy flights depending on the fraction of neighboring defectors. We
systematically studied the relationships among three factors for cooperation:
sensitivity to defectors, the intensity of L\'evy flights, and population
density. Results of evolutionary simulations showed that moderate sensitivity
most promoted cooperation. Then, we found that the shortest movements were best
for cooperation when the sensitivity to defectors was high. In contrast, when
the sensitivity was low, longer movements were best for cooperation. Thus,
L\'evy flights, the balance between short and long jumps, promoted cooperation
in any sensitivity, which was confirmed by evolutionary simulations. Finally,
as the population density became larger, higher sensitivity was more beneficial
for cooperation to evolve. Our study highlights that L\'evy flights are an
optimal searching strategy not only for foraging but also for constructing
cooperative relationships with others.


  Identifying early warning signs of sudden population changes and mechanisms
leading to regime shifts are highly desirable in population biology. In this
paper, a two-trophic ecosystem comprising of two species of predators,
competing for their common prey, with explicit interference competition is
considered. With proper rescaling, the model is portrayed as a singularly
perturbed system with fast prey dynamics and slow dynamics of the predators. In
a parameter regime near singular Hopf bifurcation, chaotic mixed-mode
oscillations (MMOs), featuring concatenation of small and large amplitude
oscillations are observed as long-lasting transients before the system
approaches its asymptotic state. To analyze the dynamical cause that initiates
a large amplitude oscillation in an MMO orbit, the model is reduced to a
suitable normal form near the singular-Hopf point. The normal form possesses a
separatrix surface that separates two different types of oscillations. A large
amplitude oscillation is initiated if a trajectory moves from the "inner" to
the "outer side" of this surface. A set of conditions on the normal form
variables are obtained to determine whether a trajectory would exhibit another
cycle of MMO dynamics before experiencing a regime shift (i.e. approaching its
asymptotic state). These conditions can serve as early warning signs for a
sudden population shift in an ecosystem.