  Cocaine addiction accounts for a large portion of substance use disorders and
threatens millions of lives worldwide. There is an urgent need to come up with
efficient anti-cocaine addiction drugs. Unfortunately, no medications have been
approved by the Food and Drug Administration (FDA), despite the extensive
effort in the past few decades. The main challenge is the intricate molecular
mechanisms of cocaine addiction, involving synergistic interactions among
proteins upstream and downstream of dopamine transporter (DAT) functions
impacted by cocaine. However, traditional in vivo or in vitro experiments can
not address the roles of so many proteins, highlighting the need for innovative
strategies in the field. We propose a proteome-informed machine learning/deep
learning (ML/DL) platform to discover nearly optimal anti-cocaine addiction
lead compounds. We construct and analyze proteomic protein-protein interaction
(PPI) networks for cocaine dependence to identify 141 involved drug targets and
represent over 60,000 associated drug candidates or experimental drugs in the
latent space using an autoencoder (EA) model trained from over 104 million
molecules. We build 32 ML models for cross-target analysis of these drug
candidates for side effects and repurposing potential. We further screen the
absorption, distribution, metabolism, excretion, and toxicity (ADMET)
properties of these candidates. Our platform reveals that essentially all of
the existing drug candidates, including dozens of experimental drugs, fail to
pass our cross-target and ADMET screenings. Nonetheless, we have identified two
nearly optimal leads for further optimization.


  The use of predator-prey models in theoretical ecology has a long history,
and the model equations have largely evolved since the original Lotka-Volterra
system towards more realistic descriptions of the processes of predation,
reproduction and mortality. One important aspect is the recognition of the fact
that the growth of a population can be subject to an Allee effect, where the
per capita growth rate increases with the population density. Including an
Allee effect has been shown to fundamentally change predator-prey dynamics and
strongly impact species persistence, but previous studies mostly focused on
scenarios of an Allee effect in the prey population. Here we explore a
predator-prey model with an ecologically important case of the Allee effect in
the predator population where it occurs in the numerical response of predator
without affecting its functional response. Biologically, this can result from
various scenarios such as a lack of mating partners, sperm limitation and
cooperative breeding mechanisms, among others. Unlike previous studies, we
consider here a generic mathematical formulation of the Allee effect without
specifying a concrete parameterisation of the functional form, and analyse the
possible local bifurcations in the system. Further, we explore the global
bifurcation structure of the model and its possible dynamical regimes for three
different concrete parameterisations of the Allee effect. The model possesses a
complex bifurcation structure: there can be multiple coexistence states
including two stable limit cycles. Inclusion of the Allee effect in the
predator generally has a destabilising effect on the coexistence equilibrium.
We also show that regardless of the parametrisation of the Allee effect,
enrichment of the environment will eventually result in extinction of the
predator population.


  We report the detection of Nipah virus in an infectious clone format, a
BSL4-level pathogen and CDC-designated Bioterrorism Agent, in raw RNA-Seq
sequencing reads deposited by the Wuhan Institute of Virology (WIV) produced
from five December 2019 patients infected with SARS-CoV-2. Research involving
Nipah infectious clones has never been reported to have occured at the WIV.
These patient samples have been previously reported to contain reads from
several other viruses: Influenza A, Spodoptera frugiperda rhabdovirus and
Nipah. Previous authors have interpreted the presence of these virus sequences
as indicative of co-infections of the patients in question by these pathogens
or laboratory contamination. However, our analysis shows that NiV genes are
encapsulated in synthetic vectors, which we infer was for assembly of a NiV
infectious clone. In particular, we document the finding of internal N, P-V-W-C
and L protein coding sequences as well as coverage of the G and F genes.
Furthermore, the format of Hepatitis D virus ribozyme and T7 terminator
downstream of the 5-prime end of the NiV sequence is consistent with truncation
required at the end of the genome for a full length infectious clone. This
indicates that research at WIV was being conducted on an assembled NiV
infectious clone. Contamination of patient sequencing reads by an infectious
NiV clone of the highly pathogenic Bangladesh strain could indicate a
significant breach of BSL-4 protocols. We call on WIV to explain the purpose of
this research on infectious clones of Nipah Virus, the full chronology of this
work, and to explain how and at what stage of sample preparation this
contamination occurred.


  Networks are widely used in many fields for their powerful ability to provide
vivid representations of relationships between variables. However, many of them
may be corrupted by experimental noise or inappropriate network inference
methods that inherently hamper the efficacy of network-based downstream
analysis. Consequently, it's necessary to develop systematic methods for
denoising networks, namely, improve the Signal-to-Noise Ratio (SNR) of noisy
networks. In this paper, we have explored the properties of network signal and
noise and proposed a novel network denoising framework called Network
Refinement (NR) that adjusts the edge weights by applying a nonlinear graph
operator based on a diffusion process defined by random walk on the graph.
Specifically, this unified framework consists of two closely linked approaches
named NR-F and NR-B, which improve the SNR of noisy input networks from two
different perspectives: NR-F aims at enhancing signal strength, while NR-B aims
at weakening noise strength. Users can choose from which angle to improve the
SNR of the network according to the characteristics of the network itself. We
show that NR can significantly refine the quality of many networks by several
applications on simulated networks and typical real-world biological and social
networks.


  Mass vaccination produces a reduction in virus circulation, but also
evolutive pressure towards the appearance of virus-resistant strains. We
discuss the balance between these two effects, in particular when the mass
vaccination takes place in the middle of an epidemic period.


  Understanding looping probabilities, including the particular case of
ring-closure or cyclization, of fluctuating polymers (eg DNA) is important in
many applications in molecular biology and chemistry. In a continuum limit the
configuration of a polymer is a curve in the group SE(3) of rigid body
displacements, whose energy can be modelled via the Cosserat theory of elastic
rods. Cosserat rods are a more detailed version of the classic wormlike-chain
(WLC) model, which we show to be more appropriate in short-length scale, or
stiff, regimes, where the contributions of extension and shear deformations are
not negligible and lead to noteworthy high values for the cyclization
probabilities (or J-factors). Characterizing the stochastic fluctuations about
minimizers of the energy by means of Laplace expansions in a (real) path
integral formulation, we develop efficient analytical approximations for the
two cases of full looping, in which both end-to-end relative translation and
rotation are prescribed, and of marginal looping probabilities, where only
end-to-end translation is prescribed. For isotropic Cosserat rods, certain
looping boundary value problems admit non-isolated families of critical points
of the energy due to an associated continuous symmetry. For the first time,
taking inspiration from (imaginary) path integral techniques, a quantum
mechanical probabilistic treatment of Goldstone modes in statistical rod
mechanics sheds light on J-factor computations for isotropic rods in the
semi-classical context. All the results are achieved exploiting appropriate
Jacobi fields arising from Gaussian path integrals, and show good agreement
when compared with intense Monte Carlo simulations for the target examples.


  The evolution of complex cellular life involved two major transitions: the
encapsulation of self-replicating genetic entities into cellular units and the
aggregation of individual genes into a collectively replicating genome. In this
paper, we formulate a minimal model of the evolution of proto-chromosomes
within protocells. We model a simple protocell composed of two types of genes:
a "fast gene" with an advantage for gene-level self-replication and a "slow
gene" that replicates more slowly at the gene level, but which confers an
advantage for protocell-level reproduction. Protocell-level replication
capacity depends on cellular composition of fast and slow genes. We use a
partial differential equation to describe how the composition of genes within
protocells evolves over time under within-cell and between-cell competition. We
find that the gene-level advantage of fast replicators casts a long shadow on
the multilevel dynamics of protocell evolution: no level of between-protocell
competition can produce coexistence of the fast and slow replicators when the
two genes are equally needed for protocell-level reproduction. By introducing a
"dimer replicator" consisting of a linked pair of the slow and fast genes, we
show analytically that coexistence between the two genes can be promoted in
pairwise multilevel competition between fast and dimer replicators, and provide
numerical evidence for coexistence in trimorphic competition between fast,
slow, and dimer replicators. Our results suggest that dimerization, or the
formation of a simple chromosome-like dimer replicator, can help to overcome
the shadow of lower-level selection and work in concert with deterministic
multilevel selection to allow for the coexistence of two genes that are
complementary at the protocell-level but compete at the level of individual
gene-level replication.


  Chromosomal rearrangements act as a source of genetic novelty by shuffling
DNA throughout the genome. These mutations can produce chimeric genes, induce
de novo gene formation, or alter gene expression changes for existing genes.
Here, we explore how these mutations may serve as agents of evolutionary change
as populations adapt to new environments during habitat shifts. We identify
16,480 rearrangements in mainland D. yakuba and two locally adapted populations
of D. santomea and D yakuba on Sao Tome. Three loci that are associated with
signals of strong differentiation in D. santomea lie adjacent to UV resistance
or DNA repair genes, suggesting that these rearrangements confer selective
advantages in high altitude environments with greater UV stressors. Some 55% of
these mutations are facilitated by TE insertions, and 28% are TE facilitated
ectopic recombination. In D. santomea 468 mutations are associated with strong
signals of differentiation from the mainland while in island D yakuba we
identify 383 candidates of local adaptation. A total of 49.4% of mutations
associated with signals of local adaptation also show significant changes in
transcript levels, suggesting that the adaptive value of rearrangements is
related to effects on gene expression. Together, this survey of structural
variation identifies key modes of evolutionary innovation that would be missed
in SNP-based screens. This work offers a portrait of how these mutations appear
and help organisms to survive during habitat shifts, furthering our
understanding of evolutionary processes.


  Indirect reciprocity is a key mechanism that promotes cooperation in social
dilemmas by means of reputation. Although it has been a common practice to
represent reputations by binary values, either `good' or `bad', such a
dichotomy is a crude approximation considering the complexity of reality. In
this work, we studied norms with three different reputations, i.e., `good',
`neutral', and `bad'. Through massive supercomputing for handling more than
thirty billion possibilities, we fully identified which norms achieve
cooperation and possess evolutionary stability against behavioural mutants. By
systematically categorizing all these norms according to their behaviours, we
found similarities and dissimilarities to their binary-reputation counterpart,
the leading eight. We obtained four rules that should be satisfied by the
successful norms, and the behaviour of the leading eight can be understood as a
special case of these rules. A couple of norms that show counter-intuitive
behaviours are also presented. We believe the findings are also useful for
designing successful norms with more general reputation systems.


  Quasi-best match graphs (qBMGs) are a hereditary class of directed, properly
vertex-colored graphs. They arise naturally in mathematical phylogenetics as a
generalization of best match graphs, which formalize the notion of evolutionary
closest relatedness of genes (vertices) in multiple species (vertex colors).
They are explained by rooted trees whose leaves correspond to vertices. In
contrast to BMGs, qBMGs represent only best matches at a restricted
phylogenetic distance. We provide characterizations of qBMGs that give rise to
polynomial-time recognition algorithms and identify the BMGs as the qBMGs that
are color-sink-free. Furthermore, two-colored qBMGs are characterized as
directed graphs satisfying three simple local conditions, two of which have
appeared previously, namely bi-transitivity in the sense of Das et al. (2021)
and a hierarchy-like structure of out-neighborhoods, i.e., $N(x)\cap
N(y)\in\{N(x),N(y),\emptyset\}$ for any two vertices $x$ and $y$. Further
results characterize qBMGs that can be explained by binary phylogenetic trees.


  The evolutionary relationships among organisms have traditionally been
represented using rooted phylogenetic trees. However, due to reticulate
processes such as hybridization or lateral gene transfer, evolution cannot
always be adequately represented by a phylogenetic tree, and rooted
phylogenetic networks that describe such complex processes have been introduced
as a generalization of rooted phylogenetic trees. In fact, estimating rooted
phylogenetic networks from genomic sequence data and analyzing their structural
properties is one of the most important tasks in contemporary phylogenetics.
Over the last two decades, several subclasses of rooted phylogenetic networks
(characterized by certain structural constraints) have been introduced in the
literature, either to model specific biological phenomena or to enable
tractable mathematical and computational analyses. In the present manuscript,
we provide a thorough review of these network classes, as well as provide a
biological interpretation of the structural constraints underlying these
networks where possible. In addition, we discuss how imposing structural
constraints on the network topology can be used to address the scalability and
identifiability challenges faced in the estimation of phylogenetic networks
from empirical data.


  Evolutionary theory predicts that children fare better if they resemble their
father. However, if a man is promiscuous then his children tend to be
(unwittingly) raised within other families; such children should fare better if
they do not resemble their father. This suggests that father-son facial
similarity should decrease as a function of male promiscuity. To test this
hypothesis, the degree of similarity between each of 37 pairs of father-son
photographs was rated by 83 participants. The promiscuity score of a father was
defined as the number of sexual relationships he had been involved in. A linear
regression of father-son similarity ratings against fathers' promiscuity scores
indicated that the slope of the fitted line was significantly less than zero
(slope=-0.176, p < 0.001), and that fathers' promiscuity scores account for 75%
of the variance father-son facial similarity ratings. Additionally, a
significant negative non-parametric correlation (Spearman's rho rho=-0.886,
p=0.019) between father-son similarity ratings and fathers' promiscuity scores
was found.


  In light of the continuing emergence of new SARS-CoV-2 variants and vaccines,
we create a simulation framework for exploring possible infection trajectories
under various scenarios. The situations of primary interest involve the
interaction between three components: vaccination campaigns, non-pharmaceutical
interventions (NPIs), and the emergence of new SARS-CoV-2 variants.
Additionally, immunity waning and vaccine boosters are modeled to account for
their growing importance. New infections are generated according to a
hierarchical model in which people have a random, individual infectiousness.
The model thus includes super-spreading observed in the COVID-19 pandemic. Our
simulation functions as a dynamic compartment model in which an individual's
history of infection, vaccination, and possible reinfection all play a role in
their resistance to further infections. We present a risk measure for each
SARS-CoV-2 variant, $\rho^\V$, that accounts for the amount of resistance
within a population and show how this risk changes as the vaccination rate
increases. Furthermore, by considering different population compositions in
terms of previous infection and type of vaccination, we can learn about
variants which pose differential risk to different countries. Different control
strategies are implemented which aim to both suppress COVID-19 outbreaks when
they occur as well as relax restrictions when possible. We demonstrate that a
controller that responds to the effective reproduction number in addition to
case numbers is more efficient and effective in controlling new waves than
monitoring case numbers alone. This is of interest as the majority of the
public discussion and well-known statistics deal primarily with case numbers.


  Evolutionary histories for species that cross with one another or exchange
genetic material can be represented by leaf-labelled, directed graphs called
phylogenetic networks. A major challenge in the burgeoning area of phylogenetic
networks is to develop algorithms for building such networks by amalgamating
small networks into a single large network. The level of a phylogenetic network
is a measure of its deviation from being a tree; the higher the level of
network, the less treelike it becomes. Various algorithms have been developed
for building level-1 networks from small networks. However, level-1 networks
may not be able to capture the complexity of some data sets. In this paper, we
present a polynomial-time algorithm for constructing a rooted binary level-2
phylogenetic network from a collection of 3-leaf networks or trinets. Moreover,
we prove that the algorithm will correctly reconstruct such a network if it is
given all of the trinets in the network as input. The algorithm runs in time
$O(t\cdot n+n^4)$ with $t$ the number of input trinets and $n$ the number of
leaves. We also show that there is a fundamental obstruction to constructing
level-3 networks from trinets, and so new approaches will need to be developed
for constructing level-3 and higher level-networks.


  There is the need to represent in a standard manner all the possible
variations of a protein or peptide primary sequence, including both artefactual
and post-translational modifications of peptides and proteins. With that
overall aim, here, the Human Proteome Organization (HUPO) Proteomics Standards
Initiative (PSI) has developed a notation, called ProForma 2.0, which is a
substantial extension of the original ProForma notation, developed by the
Consortium for Top-Down Proteomics (CTDP). ProForma 2.0 aims to unify the
representation of proteoforms and peptidoforms. Therefore, this notation
supports use cases needed for bottom-up and middle/topdown proteomics
approaches and allows the encoding of highly modified proteins and peptides
using a human and machine-readable string. ProForma 2.0 covers encoding protein
modification names and accessions, cross-linking reagents including disulfides,
glycans, modifications encoded using mass shifts and/or via chemical formulas,
labile and C or N-terminal modifications, ambiguity in the modification
position and representation of atomic isotopes, among other use cases.
Notational conventions are based on public controlled vocabularies and
ontologies. Detailed information about the notation and existing
implementations are available at http://www.psidev.info/proforma and at the
corresponding GitHub repository (https://github.com/HUPO-PSI/proforma).


  Simulations with an adaptive time-dependent bias, such as metadynamics,
enable an efficient exploration of the conformational space of a system.
However, the dynamic information of the system is altered by the bias. With
infrequent metadynamics it is possible to recover the transition rate of
crossing a barrier, if the collective variables are ideal and there is no bias
deposition near the transition state. Unfortunately, for simulations of complex
molecules, these conditions are not always fulfilled. To overcome these
limitations, and inspired by single-molecule force spectroscopy, we developed a
method based on Kramers' theory for calculating the barrier-crossing rate when
a time-dependent bias is added to the system. We assess the quality of the bias
parameter by measuring how efficiently the bias accelerates the transitions
compared to ideal behavior. We present approximate analytical expressions of
the survival probability that accurately reproduce the barrier-crossing time
statistics, and enable the extraction of the unbiased transition rate even for
challenging cases, where previous methods fail.


  Populations composed of a collection of subpopulations (demes) with random
migration between them are quite common occurrences. The emergence and
sustenance of cooperation in such a population is a highly researched topic in
the evolutionary game theory. If the individuals in every deme are considered
to be either cooperators or defectors, the migration dilemma can be envisaged:
The cooperators would not want to migrate to a defector-rich deme as they fear
of facing exploitation; but without migration, cooperation can not be
established throughout the network of demes. With a view to studying the
aforementioned scenario, in this paper, we set up a theoretical model
consisting of a coupled map lattice of replicator maps based on
two-player--two-strategy games. The replicator map considered is capable of
showing a variety of evolutionary outcomes, like convergent (fixed point)
outcomes and nonconvergent (periodic and chaotic) outcomes. Furthermore, this
coupled network of the replicator maps undergoes the phenomenon of amplitude
death leading to non-oscillatory stable synchronized states. We specifically
explore the effect of (i) the nature of coupling that models migration between
the maps, (ii) the heterogenous demes (in the sense that not all the demes have
same game being played by the individuals), (iii) the degree of the network,
and (iv) the cost associated with the migration. In the course of
investigation, we are intrigued by the effectiveness of the random migration in
sustaining a uniform cooperator fraction across a population irrespective of
the details of the replicator dynamics and the interaction among the demes.


  For countries starting to receive steady supplies of vaccines against
SARS-CoV-2, the course of Covid-19 for the following months will be determined
by the emergence of new variants and successful roll-out of vaccination
campaigns. To anticipate this scenario, we used a multilayer network model
developed to forecast the transmission dynamics of Covid-19 in Costa Rica, and
to estimate the impact of the introduction of the Delta variant in the country,
under two plausible vaccination scenarios, one sustaining Costa Rica's July
2021 vaccination pace of 30,000 doses per day and with high acceptance from the
population and another with declining vaccination pace to 13,000 doses per day
and with lower acceptance. Results suggest that the introduction and gradual
dominance of the Delta variant would increase Covid-19 hospitalizations and ICU
admissions between $35\%$ and $33.25\%$, from August 2021 to December 2021,
depending on vaccine administration and acceptance. In the presence of the
Delta variant, new Covid-19 hospitalizations and ICU admission would experience
an average increase of $24.26\%$ and $27.19\%$ respectively in the same period
if the vaccination pace drops. Our results can help decision-makers better
prepare for the COVID-19 pandemic in the months to come.


  In this paper we analyze the impact of vaccinations on spread of the COVID-19
virus for different age groups. More specifically we examine the deployment of
vaccines in the Nordic countries in a comparative analysis where we analyze
factors such as healthcare stress level and severity of disease through new
infections, hospitalizations, intensive care unit (ICU) occupancy and deaths.
Moreover, we analyze the impact of the various vaccine types, vaccination rate
on the spread of the virus in each age group for Denmark, Finland, Iceland,
Norway and Sweden from the start of the vaccination period in December 2020
until the end of August 2021. We perform a three-fold analysis: i) frequency
analysis of infections and vaccine rates by age groups, ii) rolling
correlations between vaccination strategies, severity of COVID-19 and
healthcare stress level and; iii) we also employ the epidemic Renormalization
Group (eRG) framework. The latter is used to mathematically model wave
structures, as well as the impact of vaccinations on wave dynamics. We further
compare the Nordic countries with England. Our main results are the
quantification of the impact of the vaccination campaigns on age groups
epidemiological data, across countries with high vaccine uptake. The data
clearly shows that vaccines markedly reduce the number of new cases and the
risk of serious illness.


  Cancer ranks as one of the deadliest diseases worldwide. The high mortality
rate associated with cancer is partially due to the lack of reliable early
detection methods and/or inaccurate diagnostic tools such as certain protein
biomarkers. Cell-free nucleic acids (cfNA) such as circulating long non-coding
RNAs (lncRNAs) have recently been proposed as a new class of potential
biomarkers that could improve cancer diagnosis. The reported correlation
between circulating lncRNA levels and the presence of tumors has triggered a
great amount of interest among clinicians and scientists who have been actively
investigating their potentials as reliable cancer biomarkers. In this report,
we review the progress achieved (the Good) and challenges encountered (the Bad)
in the development of circulating lncRNAs as potential biomarkers for early
cancer diagnosis. We report and discuss the specificity and sensitivity issues
of blood-based lncRNAs currently considered as promising biomarkers for various
cancers such as hepatocellular carcinoma, colorectal cancer, gastric cancer and
prostate cancer. We also emphasize the potential clinical applications (the
Beauty) of circulating lncRNAs both as therapeutic targets and agents, on top
of diagnostic and prognostic capabilities. Based on different published works,
we finally provide recommendations for investigators who seek to investigate
and compare the levels of circulating lncRNAs in the blood of cancer patients
compared to healthy subjects by RT-qPCR or Next Generation Sequencing.


  A plethora of prediction models of SARS-CoV-2 pandemic were proposed in the
past. Prediction performances not only depend on the structure and features of
the model, but also on its parametrization. Official databases are often biased
due to lag in reporting of cases, changing testing policy or incompleteness of
data. Moreover, model parametrization is time-dependent e.g. due to changing
age-structures, new emerging virus variants, non-pharmaceutical interventions
and ongoing vaccination programs. To cover these aspects, we develop a
principled approach to parametrize SIR-type epidemiologic models of different
complexities by embedding the model structure as a hidden layer into a general
Input-Output Non-Linear Dynamical System (IO-NLDS). Non-explicitly modelled
impacts on the system are imposed as inputs of the system. Observable data are
coupled to hidden states of the model by appropriate data models considering
possible biases of the data. We estimate model parameters including their
time-dependence by a Bayesian knowledge synthesis process considering parameter
ranges derived from external studies as prior information. We applied this
approach on a SIR-type model and data of Germany and Saxony demonstrating good
prediction performances. By our approach, we can estimate and compare for
example the relative effectiveness of non-pharmaceutical interventions and can
provide scenarios of the future course of the epidemic under specified
conditions. Our method of parameter estimation can be translated to other data
sets, i.e. other countries and other SIR-type models even for other disease
contexts.


  We study epidemic arrival times in meta-population disease models through the
lens of front propagation into unstable states. We demonstrate that several
features of invasion fronts in the PDE context are also relevant to the network
case. We show that the susceptible-infected-recovered model on a network is
linearly determined in the sense that the arrival times in the nonlinear system
are approximated by the arrival times of the instability in the system
linearized near the disease free state. Arrival time predictions are extended
to an susceptible-exposed-infected-recovered model. We then study a recent
model of social epidemics where high order interactions of individuals lead to
faster invasion speeds. For these pushed fronts we compute corrections to the
estimated arrival time in this case. Finally, we show how inhomogeneities in
local infection rates lead to faster average arrival times.


  Monte Carlo (MC) methods are important computational tools for molecular
structure optimizations and predictions. When solvent effects are explicitly
considered, MC methods become very expensive due to the large degree of freedom
associated with the water molecules and mobile ions. Alternatively
implicit-solvent MC can largely reduce the computational cost by applying a
mean field approximation to solvent effects and meanwhile maintains the atomic
detail of the target molecule. The two most popular implicit-solvent models are
the Poisson-Boltzmann (PB) model and the Generalized Born (GB) model in a way
such that the GB model is an approximation to the PB model but is much faster
in simulation time. In this work, we develop a machine learning-based
implicit-solvent Monte Carlo (MLIMC) method by combining the advantages of both
implicit solvent models in accuracy and efficiency. Specifically, the MLIMC
method uses a fast and accurate PB-based machine learning (PBML) scheme to
compute the electrostatic solvation free energy at each step. We validate our
MLIMC method by using a benzene-water system and a protein-water system. We
show that the proposed MLIMC method has great advantages in speed and accuracy
for molecular structure optimization and prediction.


  The spread of a contagious disease clearly depends on when infected
individuals come into contact with susceptible ones. Such effects, however,
have remained largely unexplored in the study of epidemic outbreaks. In
particular, it remains unclear how the timing of contacts interacts with the
latent and infectious stages of the disease. Here, we use real-world physical
proximity data to study this interaction and find that the temporal statistics
of actual human contact patterns i) destabilize epidemic outbreaks and ii)
modulate the basic reproduction number. We explain both observations by
distinct aspects of the observed contact patterns. On the one hand, we find the
destabilization of outbreaks to be caused by temporal clustering of contacts
leading to over-dispersed offspring distributions and increased probabilities
of otherwise rare events (zero- and super-spreading). Notably, our analysis
enables us to disentangle previously elusive sources of over-dispersion in
empirical offspring distributions. On the other hand, we find the modulation of
R0 to be caused by a periodically varying contact rate. Both mechanisms are a
direct consequence of the memory in contact behavior, and we showcase a
generative process that reproduces these non-Markovian statistics. Our results
point to the importance of including non-Markovian contact timings into studies
of epidemic outbreaks.


  Over the last fifty years, geneticists have made great strides in
understanding how our species' evolutionary history gave rise to current
patterns of human genetic diversity classically summarized by Lewontin in his
1972 paper, 'The Apportionment of Human Diversity'. One evolutionary process
that requires special attention in both population genetics and statistical
genetics is admixture: gene flow between two or more previously separated
source populations to form a new admixed population. The admixture process
introduces unique patterns of genetic variation within and between populations,
which in turn influences the inference of demographic histories, identification
genetic targets of selection, and prediction of phenotypes. In this review, we
highlight recent studies and methodological advances that have leveraged
genomic signatures of admixture to gain insights into human history, natural
selection, and complex trait architecture. We also outline some challenges for
admixture population genetics, including limitations of applying methods
designed for single-ancestry populations to the study of admixed populations.


  Tree balance plays an important role in phylogenetics and other research
areas, which is why several indices to measure tree balance have been
introduced over the years. Nevertheless, a formal definition of what a balance
index actually is and what makes it a useful measure of balance (or, in other
cases, imbalance), has so far not been introduced in the literature. While the
established indices all summarize the (im)balance of a tree in a single number,
they vary in their definitions and underlying principles. It is the aim of the
present manuscript to introduce formal definitions of balance and imbalance
indices that classify desirable properties of such indices and to analyze and
categorize established indices accordingly. In this regard, we review 19
established (im)balance indices from the literature, summarize their general,
statistical and combinatorial properties (where known), prove numerous
additional results and indicate directions for future research by making
explicit open questions and gaps in the literature. We also prove that a few
tree shape statistics that have been used to measure tree balance in the
literature do not fulfill our definition of an (im)balance index, which might
indicate that their properties are not as useful for practical purposes.
Moreover, we show that five additional tree shape statistics from other
contexts actually are tree (im)balance indices according to our definition. The
manuscript is accompanied by the website \url{treebalance.wordpress.com}
containing fact sheets of the discussed indices. Moreover, we introduce the
software package \verb|treebalance| implemented in $\mathsf{R}$ that can be
used to calculate all indices discussed.


  The vulnerability of networks to targeted attacks is an issue of widespread
interest for policymakers, military strategists, network engineers and systems
biologists alike. Current approaches to circumvent targeted attacks seek to
increase the robustness of a network by changing the network structure in one
way or the other, leading to a higher size of the largest connected component
for a given fraction of nodes removed. In this work, we propose a strategy in
which there is a pre-existing, dormant spare capacity already built into the
network for an identified vulnerable node, such that the traffic of the
disrupted node can be diverted to another pre-existing node/set of nodes in the
network. Using our algorithm, the increase in robustness of canonical
scale-free networks was nearly 16-fold. We also analysed real-world networks
using our algorithm, where the mean increase in robustness was nearly
five-fold. To our knowledge, these numbers are significantly higher than those
hitherto reported in literature. The normalized cost of this spare capacity and
its effect on the operational parameters of the network have also been
discussed. Instances of spare capacity in biological networks, termed as
distributed robustness, are also presented.


  Since its selection as the method of the year in 2013, single-cell
technologies have become mature enough to provide answers to complex research
questions. With the growth of single-cell profiling technologies, there has
also been a significant increase in data collected from single-cell profilings,
resulting in computational challenges to process these massive and complicated
datasets. To address these challenges, deep learning (DL) is positioning as a
competitive alternative for single-cell analyses besides the traditional
machine learning approaches. Here we present a processing pipeline of
single-cell RNA-seq data, survey a total of 25 DL algorithms and their
applicability for a specific step in the processing pipeline. Specifically, we
establish a unified mathematical representation of all variational autoencoder,
autoencoder, and generative adversarial network models, compare the training
strategies and loss functions for these models, and relate the loss functions
of these models to specific objectives of the data processing step. Such
presentation will allow readers to choose suitable algorithms for their
particular objective at each step in the pipeline. We envision that this survey
will serve as an important information portal for learning the application of
DL for scRNA-seq analysis and inspire innovative use of DL to address a broader
range of new challenges in emerging multi-omics and spatial single-cell
sequencing.


  Seasonal influenza kills hundreds of thousands every year, with multiple
constantly-changing strains in circulation at any given time. A high mutation
rate enables the influenza virus to evade recognition by the human immune
system, including immunity acquired through past infection and vaccination.
Here, we capture the genetic similarity of influenza strains and their
evolutionary dynamics with genotype networks. We show that the genotype
networks of influenza A (H3N2) hemagglutinin are characterized by heavy-tailed
distributions of module sizes and connectivity, suggesting critical-like
behavior. We argue that: (i) genotype networks are driven by mutation and host
immunity to explore a subspace of networks predictable in structure, and (ii)
genotype networks provide an underlying structure necessary to capture the rich
dynamics of multistrain epidemic models. In particular, inclusion of
strain-transcending immunity in epidemic models is dependent upon the structure
of an underlying genotype network. This interplay suggests a self-organized
criticality where the epidemic dynamics of influenza locates critical-like
regions of its genotype network. We conclude that this interplay between
disease dynamics and network structure might be key for future network analysis
of pathogen evolution and realistic multistrain epidemic models.


  We develop a Bayesian inference framework to quantify uncertainties in
epidemiological models. We use SEIJR and SIJR models involving populations of
susceptible, exposed, infective, diagnosed, dead and recovered individuals to
infer from covid-19 data rate constants, as well as their variations in
response to lockdown measures. To account for confinement, we distinguish two
susceptible populations at different risk: confined and unconfined. We show
that transmission and recovery rates within them vary in response to facts. A
key unknown to predict the evolution of the epidemic is the fraction of the
population affected by the virus, including asymptomatic subjects. Our study
tracks its time evolution with quantified uncertainty from available official
data from the onset of the epidemic, limited, however, by the data quality. We
exemplify the technique with data from Spain, country in which late drastic
lockdowns were enforced for months. In late actions and in the absence of other
measures, spread is delayed but not stopped unless a large enough fraction of
the population is confined until the asymptomatic population is depleted. To
some extent, confinement could be replaced by strong distancing through masks
in adequate circumstances.


  Besides maintaining health precautions, vaccination has been the only
prevention from SARS-CoV-2, though no clinically proved 100% effective vaccine
has been developed till date. At this stage, to withhold the debris of this
pandemic, experts need to know the impact of the vaccine efficacy rate's
threshold and how long this pandemic may extent with vaccines that have
different efficacy rates. In this article, a mathematical model study has been
done on the importance of vaccination and vaccine efficiency rate during an
ongoing pandemic. We simulated a five compartment mathematical model to analyze
the pandemic scenario in both California, and whole U.S. We considered four
vaccines, Pfizer, Moderna, AstraZeneca, and Johnson and Johnson, which are
being used rigorously to control the COVID-19 pandemic, in addition with two
special cases: a vaccine with 100% efficacy rate and no vaccine under use. Both
the infection and death rates are very high in California. Our model suggests
that the pandemic situation in California will be under control in the last
quartile of the year 2023 if frequent vaccination is continued with the Pfizer
vaccine. During this time, six waves will happen from the beginning of the
immunization where the case fatality and recovery rates will be 1.697% and
98.30%, respectively. However, according to the considered model, this period
might be extended to the mid of 2024 when vaccines with lower efficacy rates
are used. The more effective a vaccine, the less people suffer from this malign
infection. Although specific groups of people get prioritized initially, mass
vaccination is needed to control the spread of the disease.


  Last time, an attempt was made to analyze the first wave of COVID-19 in Japan
from February to May 2020 by the Avrami equation. This time, all infection
waves that occurred repeatedly from February 2020 to September 2021 were
simulated by developing the last work. The entire waveform was basically
simulated by the superposition of five major waves, and the least-squares
method was applied to each wave to determine the parameters of the equation.
Despite the simulation using the single equation assuming that the parameters
were constants, the accuracy of the simulated cumulative infection D was as
good as 95%CI / D < 2.5% after the second wave when D exceeded 10,000. Since
the simulation was highly accurate, it was effective in predicting the number
of infections in the near future and in considering the cause of changes in the
number of infections observed as the detailed structure of the waveform.


  Drug resistance and metastasis - the major complications in cancer - both
entail adaptation of cancer cells to stress, whether a drug or a lethal new
environment. Intriguingly, these adaptive processes share similar features that
cannot be explained by a pure Darwinian scheme, including dormancy, increased
heterogeneity, and stress-induced plasticity. Here, we propose that learning
theory offers a framework to explain these features and may shed light on these
two intricate processes. In this framework, learning is performed at the single
cell level, by stress-driven exploratory trial-and-error. Such a process is not
contingent on pre-existing pathways but on a random search for a state that
diminishes the stress. We review underlying mechanisms that may support this
search, and show by using a learning model that such exploratory adaptation is
feasible in a high dimensional system as the cell. At the population level, we
view the tissue as a network of exploring agents that communicate and restrain
cancer formation in health. In this view, disease results from the breakdown of
homeostasis between cellular exploratory drive and tissue homeostasis.


  Networks of gene regulation govern morphogenesis, determine cell identity and
regulate cell function. But we have little understanding, at the local level,
of which logics are biologically preferred or even permitted. To solve this
puzzle, we studied the consequences of a fundamental aspect of gene regulatory
networks: genes and transcription factors talk to each other but not
themselves. Remarkably, this bipartite structure severely restricts the number
of logical dependencies that a gene can have on other genes. We developed a
theory for the number of permitted logics for different regulatory building
blocks of genes and transcription factors. We tested our predictions against a
simulation of the 19 simplest building blocks, and found complete agreement.
The restricted range of biological logics is a key insight into how information
is processed at the genetic level. It constraints global network function and
makes it easier to reverse engineer regulatory networks from observed behavior.


  In controlling transmission of COVID-19, the effectiveness of border
quarantine strategies is a key concern for jurisdictions in which the local
prevalence of disease and immunity is low. In settings like this such as China,
Australia, and New Zealand, rare outbreak events can lead to escalating
epidemics and trigger the imposition of large scale lockdown policies. Here, we
examine to what degree vaccination status of incoming arrivals and the
quarantine workforce can allow relaxation of quarantine requirements. To do so,
we develop and apply a detailed model of COVID-19 disease progression and
transmission taking into account nuanced timing factors. Key among these are
disease incubation periods and the progression of infection detectability
during incubation. Using the disease characteristics associated with the
ancestral lineage of SARS-CoV-2 to benchmark the level of acceptable risk, we
examine the performance of the border quarantine system for vaccinated
arrivals. We examine disease transmission and vaccine efficacy parameters over
a wide range, covering plausible values for the Delta variant currently
circulating globally. Our results indicate a threshold in outbreak potential as
a function of vaccine efficacy, with the time until an outbreak increasing by
up to two orders of magnitude as vaccine efficacy against transmission
increases from 70% to 90%. For parameters corresponding to the Delta variant,
vaccination is able to maintain the capacity of quarantine systems to reduce
case importation and outbreak risk, by counteracting the pathogen's increased
infectiousness. To prevent outbreaks, heightened vaccination in border
quarantine systems must be combined with mass vaccination. The ultimate success
of these programs will depend sensitively on the efficacy of vaccines against
viral transmission.


  We present an information-theoretic formalism to study signal transduction in
four architectural variants of a model two-step cascade with increasing input
population. Our results categorize these four types into two classes depending
upon the effect played out by activation and repression on mutual information,
net synergy, and signal-to-noise ratio. Within the Gaussian framework and using
the linear noise approximation, we derive the analytic expressions for these
metrics to establish their underlying relationships in terms of the biochemical
parameters. We also verify our approximations through stochastic simulations.


  Being a dual purpose enzyme, the DNA polymerase is responsible for elongation
of the newly formed DNA strand as well as cleaving the erroneous growth in case
of a misincorporation. The efficiency of replication depends on the
coordination of the polymerization and exonuclease activity of DNA polymerase.
Here we propose and analyze a minimal kinetic model of DNA replication and
determine exact expressions for the velocity of elongation and the accuracy of
replication. We first analyze the case without exonuclease activity. In that
case, accuracy is determined by a kinetic competition between stepping and
unbinding, with discrimination between correct and incorrect nucleotides in
both transitions. We then include exonuclease activity and ask how different
modes of additional discrimination in the exonuclease pathway can improve the
accuracy while limiting the detrimental effect of exonuclase on the speed of
replication. In this way, we ask how the kinetic parameters of the model have
to be set to coordinate the two activities of the enzyme for high accuracy and
high speed. The analysis also shows that the design of a replication system
does not universally have to follow the speed-accuracy trade-off rule, although
it does in the biologically realized parameter range. The accuracy of the
process is mainly controlled by the crucial role of stepping after erroneous
incorporation, which has impact on both polymerase and exonuclease activities
of DNA polymerase.


  We analyze the harvesting and stocking of a population that is affected by
random and seasonal environmental fluctuations. The main novelty comes from
having three layers of environmental fluctuations. The first layer is due to
the environment switching at random times between different environmental
states. This is similar to having sudden environmental changes or catastrophes.
The second layer is due to seasonal variation, where there is a significant
change in the dynamics between seasons. Finally, the third layer is due to the
constant presence of environmental stochasticity -- between the seasonal or
random regime switches, the species is affected by fluctuations which can be
modelled by white noise. This framework is more realistic because it can
capture both significant random and deterministic environmental shifts as well
as small and frequent fluctuations in abiotic factors. Our framework also
allows for the price or cost of harvesting to change deterministically and
stochastically, something that is more realistic from an economic point of
view.
  The combined effects of seasonal and random fluctuations make it impossible
to find the optimal harvesting-stocking strategy analytically. We get around
this roadblock by developing rigorous numerical approximations and proving that
they converge to the optimal harvesting-stocking strategy. We apply our methods
to multiple population models and explore how prices, or costs, and
environmental fluctuations influence the optimal harvesting-stocking strategy.
We show that in many situations the optimal way of harvesting and stocking is
not of threshold type.


  Given the ongoing Covid-19 pandemic, it is of interest to understand how the
infections spread as the combined result of measures taken by central planners
(governments) and individual behavior. In this work, the spread of Covid-19 is
modelled as a differentiable game between the planner and population with
appropriate disease spread dynamical equations. We first characterise the
equilibrium dynamics of only the population with modifed
Susceptible-Infected-Recovered (SIR) equations to highlight the qualitative
nature of the equilbrium. Using this result, we formulate the joint equilibrium
exposure profile between the planner and population. Additionally, as in case
of Covid-19, the role of asymptomatic carriers, inadequacies in testing,
contact tracing and quarantining can lead to a significant underestimate of the
true infected numbers as compared to just the detected numbers. Therefore, it
is vital to model the true infected numbers within the context of choices made
by individuals within the population. To incorporate this, we extend our
framework by modifying the dynamics to include additional sub-compartments of
`undetected infected' and `detected infected' in the disease dynamics. The
individuals make their own estimates of the total infected from the detected
numbers and base their strategies on those estimates. We show that these
considerations lead to a retarded optimal control problem for the players. We
present some simulation results based on these results to demonstrate how
population behavior, planner control, detection rates and trust in the reported
numbers play a key role in how the disease spreads.


  In the subtropical gyres, phytoplankton rely on eddies for transporting
nutrients from depth to the euphotic zone. But, what controls the rate of
nutrient supply for new production? We show that vertical nutrient flux depends
both on the vertical motion within the eddying flow and varies nonlinearly with
the growth rate of the phytoplankton itself. Flux is maximized when the growth
rate matches the inverse of the decorrelation timescale for vertical motion.
Using a three-dimensional ocean model and a linear nutrient uptake model, we
find that phytoplankton productivity is maximized for a growth rate of 1/3
day$^{-1}$, which corresponds to the timescale of submesoscale dynamics.
Variability in the frequency of vertical motion across different physical
features of the flow favors phytoplankton production with different growth
rates. Such a growth-transport feedback can generate diversity in the
phytoplankton community structure at submesoscales and higher net productivity
in the presence of community diversity.


  Identifying mechanisms able to sustain costly cooperation among
self-interested agents is a central problem across social and biological
sciences. One possible solution is peer punishment: when agents have an
opportunity to sanction defectors, classical behavioral experiments suggest
that cooperation can take root. Overlooked from standard experimental designs,
however, is the fact that real-world human punishment -- the administration of
justice -- is intrinsically noisy. Here we show that stochastic punishment
falls short of sustaining cooperation in the repeated public good game. As
punishment noise increases, we find that contributions decrease and punishment
efforts intensify, resulting in a $45\%$ drop in gains compared to a noiseless
control. Moreover, we observe that uncertainty causes a rise in antisocial
punishment, a mutually harmful behavior previously associated with societies
with a weak rule of law. Our approach brings to light challenges to cooperation
that cannot be explained by economic rationality and strengthens the case for
further investigations of the effect of noise -- and not just bias -- on human
behavior.


  With the unprecedented global vaccination campaign against SARS-CoV-2
attention has now turned to the potential impact of this large-scale
intervention on the evolution of the virus. In this perspective we summarize
what is currently known about evolution in the context of vaccination from
research on other pathogen species, with an eye towards the future evolution of
SARS-CoV-2.


  Identifying the genes and mutations that drive the emergence of tumors is a
major step to improve understanding of cancer and identify new directions for
disease diagnosis and treatment. Despite the large volume of genomics data, the
precise detection of driver mutations and their carrying genes, known as cancer
driver genes, from the millions of possible somatic mutations remains a
challenge. Computational methods play an increasingly important role in
identifying genomic patterns associated with cancer drivers and developing
models to predict driver events. Machine learning (ML) has been the engine
behind many of these efforts and provides excellent opportunities for tackling
remaining gaps in the field. Thus, this survey aims to perform a comprehensive
analysis of ML-based computational approaches to identify cancer driver
mutations and genes, providing an integrated, panoramic view of the broad data
and algorithmic landscape within this scientific problem. We discuss how the
interactions among data types and ML algorithms have been explored in previous
solutions and outline current analytical limitations that deserve further
attention from the scientific community. We hope that by helping readers become
more familiar with significant developments in the field brought by ML, we may
inspire new researchers to address open problems and advance our knowledge
towards cancer driver discovery.


  Mathematical models of infectious diseases exhibit robust dynamics such as
stable endemic or a disease-free equilibrium, or convergence of the solutions
to periodic epidemic waves. The present work shows that the accuracy of such
dynamics can be significantly improved by incorporating both local and global
dynamics of the infection in disease models. To demonstrate improved
accuracies, we extended a standard Susceptible-Infected-Recovered (SIR) model
by incorporating global dynamics of the COVID-19 pandemic. The extended SIR
model assumes three possibilities for the susceptible individuals traveling
outside of their community: They can return to the community without any
exposure to the infection, they can be exposed and develop symptoms after
returning to the community, or they can be tested positive during the trip and
remain quarantined until fully recovered. To examine the predictive accuracies
of the extended SIR model, we studied the prevalence of the COVID-19 infection
in Kansas City, Missouri influenced by the COVID-19 global pandemic. Using a
two-step model-fitting algorithm, the extended SIR model was parameterized
using the Kansas City, Missouri COVID-19 data during March to October 2020. The
extended SIR model significantly outperformed the standard SIR model and
revealed oscillatory behaviors with an increasing trend of infected
individuals. In conclusion, the analytics and predictive accuracies of disease
models can be significantly improved by incorporating the global dynamics of
the infection in the models.


  This chapter describes the application of constrained geometric simulations
for prediction of antibody structural dynamics. We utilize constrained
geometric simulations method FRODAN, which is a low computational complexity
alternative to Molecular Dynamics (MD) simulations that can rapidly explore
flexible motions in protein structures. FRODAN is highly suited for
conformational dynamics analysis of large proteins, complexes, intrinsically
disordered proteins and dynamics that occurs on longer biologically relevant
time scales which are normally inaccessible to classical MD simulations. This
approach predicts protein dynamics at an all-atom scale while retaining
realistic covalent bonding, maintaining dihedral angles in energetically good
conformations while avoiding steric clashes in addition to performing other
geometric and stereochemical criteria checks. In this chapter, we apply FRODAN
to showcase its applicability for probing functionally relevant dynamics of
IgG2a, including large amplitude domain-domain motions and motions of
complementarity determining region (CDR) loops. As was suggested in previous
experimental studies, our simulations show that antibodies can explore a large
range of conformational space.


  One of the hallmarks of cancer cells is their exceptional ability to migrate
within the extracellular matrix (ECM) for gaining access to the circulatory
system, a critical step of cancer metastasis. RhoA, a small GTPase, is known to
be a key molecular switch that toggles between actomyosin contractility and
lamellipodial protrusion during cell migration. Current understanding of RhoA
activity in cell migration has been largely derived from studies of cells
plated on a two-dimensional (2D) substrate using a FRET biosensor. There has
been increasing evidence that cells behave differently in a more
physiologically relevant three-dimensional (3D) environment, however, studies
of RhoA activities in 3D have been hindered by low signal-to-noise ratio in
fluorescence imaging. In this paper, we present a machine learning-assisted
FRET technique to follow the spatiotemporal dynamics of RhoA activities of
single breast tumor cells (MDA-MB-231) migrating in a 3D as well as a 2D
environment using a RhoA biosensor. We found that RhoA activity is more
polarized along the long axis of the cell for single cells migrating on 2D
fibronectin-coated glass versus those embedded in 3D collagen matrices. In
particular, RhoA activities of cells in 2D exhibit a distinct front-to-back and
back-to-front movement during migration in contrast to those in 3D. Finally,
regardless of dimensionality, RhoA polarization is found to be correlated with
cell shape.


  Genomic prediction can be a powerful tool to achieve greater rates of genetic
gain for quantitative traits if thoroughly integrated into a breeding strategy.
In rice as in other crops, the interest in genomic prediction is very strong
with a number of studies addressing multiple aspects of its use, ranging from
the more conceptual to the more practical. In this chapter, we review the
literature on rice (Oryza sativa) and summarize important considerations for
the integration of genomic prediction in breeding programs. The irrigated
breeding program at the International Rice Research Institute is used as a
concrete example on which we provide data and R scripts to reproduce the
analysis but also to highlight practical challenges regarding the use of
predictions. The adage: "To someone with a hammer, everything looks like a
nail" describes a common psychological pitfall that sometimes plagues the
integration and application of new technologies to a discipline. We have
designed this chapter to help rice breeders avoid that pitfall and appreciate
the benefits and limitations of applying genomic prediction, as it is not
always the best approach nor the first step to increasing the rate of genetic
gain in every context.


  The well-known issue of reconstructing regulatory networks from gene
expression measurements has been somewhat disrupted by the emergence and rapid
development of single-cell data. Indeed, the traditional way of seeing a gene
regulatory network as a deterministic system affected by small noise is being
challenged by the highly stochastic, bursty nature of gene expression revealed
at single-cell level. In previous work, we described a promising strategy in
which network inference is seen as a calibration procedure for a mechanistic
model driven by transcriptional bursting: this model inherently captures the
typical variability of single-cell data without requiring ad hoc external
noise, unlike ordinary or even stochastic differential equations often used in
this context. The resulting algorithm, based on approximate resolution of the
related master equation using a self-consistent field, was derived in detail
but only applied as a proof of concept to simulated two-gene networks. Here we
derive a simplified version of the algorithm and apply it, in more relevant
situations, to both simulated and real single-cell RNA-Seq data. We point out
three interesting features of this approach: it is computationally tractable
with realistic numbers of cells and genes, it provides inferred networks with
biological interpretability, and the underlying mechanistic model allows
testable predictions to be made. A practical implementation of the inference
procedure, together with an efficient stochastic simulation algorithm for the
model, is available as a Python package.


  The tragedy of the commons (TOC) is an unfortunate situation where a shared
resource is exhausted due to uncontrolled exploitation by the selfish
individuals of a population. Recently, the paradigmatic replicator equation has
been used in conjunction with a phenomenological equation for the state of the
shared resource to gain insight into the influence of the games on the TOC. The
replicator equation, by construction, models a fixed infinite population
undergoing microevolution. Thus, it is unable to capture any effect of the
population growth and the carrying capacity of the population although the TOC
is expected to be dependent on the size of the population. Therefore, in this
paper, we present a mathematical framework that incorporates the density
dependent payoffs and the logistic growth of the population in the
eco-evolutionary dynamics modelling the game-resource feedback. We discover a
bistability in the dynamics: a finite carrying capacity can either avert or
cause the TOC depending on the initial states of the resource and the initial
fraction of cooperators. In fact, depending on the type of strategic
game-theoretic interaction, a finite carrying capacity can either avert or
cause the TOC when it is exactly the opposite for the corresponding case with
infinite carrying capacity.


  In synaptic molecular communication, the activation of postsynaptic receptors
by neurotransmitters (NTs) is governed by a stochastic reaction-diffusion
process and, hence, inherently random. It is currently not fully understood how
this randomness impacts downstream signaling in the target cell and,
ultimately, neural computation and learning. The statistical characterization
of the reaction-diffusion process is difficult because the reversible
bi-molecular reaction of NTs and receptors renders the system nonlinear.
Consequently, existing models for the receptor occupancy in the synaptic cleft
rely on simplifying assumptions and approximations which limit their practical
applicability. In this work, we propose a novel statistical model for the
reaction-diffusion process governing synaptic signal transmission in terms of
the chemical master equation (CME). We show how to compute the CME efficiently
and verify the accuracy of the obtained results with stochastic particle-based
computer simulations (PBSs). Furthermore, we compare the proposed model to two
benchmark models proposed in the literature and show that it provides more
accurate results when compared to PBSs. Finally, the proposed model is used to
study the impact of the system parameters on the statistical dependence between
binding events of NTs and receptors. In summary, the proposed model provides a
step forward towards a complete statistical characterization of synaptic signal
transmission.


  Theoretical models show that in a non-constant environment two predator
species feeding on one and the same prey may coexist because the two species
occupy different temporal niches: the one with the longer handling time has the
advantage when prey is rare so that holding on to the same catch is the better
option, while the species with the shorter handling time has the advantage when
prey is common and easy to catch. In this paper we address the question whether
a predator species with a handling time that is not fixed but decreases with
prey density could be selective superior regardless of whether the prey is rare
or common, as such predator would be able to occupy both temporal niches all by
itself.
  To that end we study the Rosenzweig-MacArthur model with a modified Holling
type II functional response with a density dependent handling time and a
handling time dependent conversion factor. We find that the population dynamics
tend to be richer than that of the standard model with fixed handling times
because of the possibility of multiple positive equilibria and positive
attractors. Increasing the strength of the density dependence eventually
stabilises the population.
  Using the framework of adaptive dynamics, we study the evolution of the
strength of the density dependence. We find that a predator with even a weakly
density dependent handling time can invade both monomorphic and dimorphic
populations of predators with fixed handling times. Eventually the strength of
the density dependence of the handling time evolves to a level where population
cycles are lost, with that the possibility of predator coexistence as well.


  Zoonotic diseases spread through pathogens-infected animal carriers. In the
case of Ebola Virus Disease (EVD), evidence supports that the main carriers are
fruit bats and non-human primates. Further, EVD spread is a multi-factorial
problem that depends on sociodemographic and economic (SDE) factors. Here we
inquire into this phenomenon and aim at determining, quantitatively, the Ebola
spillover infection exposure map and try to link it to SDE factors. To that
end, we designed and conducted a survey in Sierra Leone and implement a
pipeline to analyze data using regression and machine learning techniques. Our
methodology is able (1) to identify the features that are best predictors of an
individual's tendency to partake in behaviors that can expose them to Ebola
infection, (2) to develop a predictive model about the spillover risk
statistics that can be calibrated for different regions and future times, and
(3) to compute a spillover exposure map for Sierra Leone. Our results and
conclusions are relevant to identify the regions in Sierra Leone at risk of EVD
spillover and, consequently, to design and implement policies for an effective
deployment of resources (e.g., drug supplies) and other preventative measures
(e.g., educational campaigns).


  This paper develops an agent-based disease spread model on a contact network
in an effort to guide efforts at surveillance testing in small to moderate
facilities such as nursing homes and meat-packing plants. The model employs
Monte Carlo simulations of viral spread sample paths in the contact network.
The original motivation was to detect COVID-19 outbreaks quickly in such
facilities, but the model can be applied to any communicable disease. In
particular, the model provides guidance on how many test to administer each day
and on the importance of the testing order among staff or workers.


  In this survey we report some recent results in the mathematical modeling of
epidemic phenomena through the use of kinetic equations. We initially consider
models of interaction between agents in which social characteristics play a key
role in the spread of an epidemic, such as the age of individuals, the number
of social contacts, and their economic wealth. Subsequently, for such models,
we discuss the possibility of containing the epidemic through an appropriate
optimal control formulation based on the policy maker's perception of the
progress of the epidemic. The role of uncertainty in the data is also discussed
and addressed. Finally, the kinetic modeling is extended to spatially dependent
settings using multiscale transport models that can characterize the impact of
movement dynamics on epidemic advancement on both one-dimensional networks and
realistic two-dimensional geographic settings.


  The properties of competition models where all individuals are identical are
relatively well-understood; however, juveniles and adults can experience or
generate competition differently. We study here less well-known structured
competition models in discrete time that allow multiple life history parameters
to depend on adult or juvenile population densities. A numerical study with
Ricker density-dependence suggested that when competition coefficients acting
on juvenile survival and fertility reflect opposite competitive hierarchies,
stage structure could foster coexistence. We revisit and expand those results.
First, through a Beverton-Holt two-species juvenile-adult model, we confirm
that these findings do not depend on the specifics of density-dependence or
life cycles, and obtain analytical expressions explaining how this coexistence
emerging from stage structure can occur. Second, we show using a
community-level sensitivity analysis that such emergent coexistence is robust
to perturbations of parameter values. Finally, we ask whether these results
extend from two to many species, using simulations. We show that they do not,
as coexistence emerging from stage structure is only seen for very similar
life-history parameters. Such emergent coexistence is therefore not likely to
be a key mechanism of coexistence in very diverse ecosystems, although it may
contribute to explaining coexistence of certain pairs of intensely competing
species.


  Compartmental models are popular in the mathematics of epidemiology for their
simplicity and wide range of applications. Although they are typically solved
as initial value problems for a system of ordinary differential equations, the
observed data is typically akin of a boundary value type problem: we observe
some of the dependent variables at given times, but we do not know the initial
conditions. In this paper, we reformulate the classical
Susceptible-Infectious-Recovered system in terms of the number of detected
positive infected cases at different times, we then prove the existence and
uniqueness of a solution to the derived boundary value problem and then present
a numerical algorithm to approximate the solution.


  Background: Single-cell RNA sequencing (scRNA-seq) yields valuable insights
about gene expression and gives critical information about complex tissue
cellular composition. In the analysis of single-cell RNA sequencing, the
annotations of cell subtypes are often done manually, which is time-consuming
and irreproducible. Garnett is a cell-type annotation software based the on
elastic net method. Besides cell-type annotation, supervised machine learning
methods can also be applied to predict other cell phenotypes from genomic data.
Despite the popularity of such applications, there is no existing study to
systematically investigate the performance of those supervised algorithms in
various sizes of scRNA-seq data sets.
  Methods and Results: This study evaluates 13 popular supervised machine
learning algorithms to classify cell phenotypes, using published real and
simulated data sets with diverse cell sizes. The benchmark contained two parts.
In the first part, we used real data sets to assess the popular supervised
algorithms' computing speed and cell phenotype classification performance. The
classification performances were evaluated using AUC statistics, F1-score,
precision, recall, and false-positive rate. In the second part, we evaluated
gene selection performance using published simulated data sets with a known
list of real genes.
  Conclusion: The study outcomes showed that ElasticNet with interactions
performed best in small and medium data sets. NB was another appropriate method
for medium data sets. In large data sets, XGB works excellent. Ensemble
algorithms were not significantly superior to individual machine learning
methods. Adding interactions to ElasticNet can help, and the improvement was
significant in small data sets.


  We develop methods for investigating protein drift-diffusion dynamics in
heterogeneous cell membranes and the roles played by geometry, diffusion,
chemical kinetics, and phase separation. Our hybrid stochastic numerical
methods combine discrete particle descriptions with continuum-level models for
tracking the individual protein drift-diffusion dynamics when coupled to
continuum fields. We show how our approaches can be used to investigate
phenomena motivated by protein kinetics within dendritic spines. The spine
geometry is hypothesized to play an important biological role regulating
synaptic strength, protein kinetics, and self-assembly of clusters. We perform
simulation studies for model spine geometries varying the neck size to
investigate how phase-separation and protein organization is influenced by
different shapes. We also show how our methods can be used to study the roles
of geometry in reaction-diffusion systems including Turing instabilities. Our
methods provide general approaches for investigating protein kinetics and
drift-diffusion dynamics within curved membrane structures.


  This paper is the report of the problem proposed for the !Optimizer 2021
competition, and the solutions of the gold medalist team, i.e., the Panda team.
The competition was held in two stages, the research and development stage and
a two-week contest stage, consisting of five rounds, and seven teams succeeded
in finishing both stages to the end. In this joint report of the winner team
Panda and the problem design committee coordinated by Mojtaba Tefagh, we first
explain each of the five rounds and then provide the solutions proposed by our
team (Panda) to fulfill the required tasks in the fastest and most accurate
way. Afterward, some preprocessing and data manipulating ideas used to enhance
the algorithms would be presented. All codes are written in the Julia language,
which showed a better performance than Python on optimization problems in our
comparisons during the R&D stage, and are publicly available in the Github
repository: https://github.com/Optimizer-Competition-Panda


  Increased Atmospheric CO2 to over 400 ppm has prompted global climate
irregularities. Reducing the released CO2 from biotechnological processes could
remediate these phenomena. In this study, we sought to find a solution to
reduce the amount of CO2 in the process of growth and reproduction by
preventing the conversion of formic acid into CO2. The (bio)chemical conversion
of formic acid to CO2 is a key reaction. Therefore, we compared the growth of
BL21, being a subfamily of K12, alongside two strains in which two different
genes related to the formate metabolism were deleted, in complex and simple
media. Experimental results were entirely consistent with metabolic
predictions. Subsequently, the knockout bacteria grew more efficiently than
BL21. Interestingly, TsFDH, a formate dehydrogenase with the tendency of
converting CO2 to formate, increased the growth of all strains compared with
cells without the TsFDH. Most mutants grew in a simple medium containing
glycerol, which showed that glycerol is the preferred carbon source compared to
glucose for the growth of E. coli. These results explain the reasons for the
inconsistency of predictions in previous metabolic models that declared
glycerol as a suitable carbon source for the growth of E. coli but failed to
achieve it in practice. To conduct a more mechanistic evaluation of our
observations, RNA sequencing data analysis was conducted on an E. coli RNA-seq
dataset. The gene expression correlation outcome revealed the increased
expression levels of several genes related to protein biosynthesis and glycerol
degradation as a possible explanation of our observations.


  Recent studies revealed an important interplay between the detailed structure
of fibration symmetric circuits and the functionality of biological and
non-biological networks within which they have be identified. The presence of
these circuits in complex networks are directed related to the phenomenon of
cluster synchronization, which produces patterns of synchronized group of
nodes. Here we present a fast, and memory efficient, algorithm to identify
fibration symmetries over information-processing networks. This algorithm is
specially suitable for large and sparse networks since it has runtime of
complexity $O(M\log N)$ and requires $O(M+N)$ of memory resources, where $N$
and $M$ are the number of nodes and edges in the network, respectively. We
propose a modification on the so-called refinement paradigm to identify
circuits symmetrical to information flow (i.e., fibers) by finding the coarsest
refinement partition over the network. Finally, we show that the presented
algorithm provides an optimal procedure for identifying fibers, overcoming the
current approaches used in the literature.


  Despite of the fast development of highly effective vaccines to control the
current COVID$-$19 pandemic, the unequal distribution and availability of these
vaccines worldwide and the number of people infected in the world lead to the
continuous emergence of SARS-CoV-2 (Severe Acute Respiratory Syndrome
coronavirus 2) variants of concern. It is likely that real-time genomic
surveillance will be continuously needed as an unceasing monitoring tool,
necessary to follow the spillover of the disease spread and the evolution of
the virus. In this context, new genomic variants of SARS-CoV-2 that may emerge
as a response to selective pressure, including variants refractory to current
vaccines, makes genomic surveillance programs tools of utmost importance. Here
propose a statistical model for the estimation of the relative frequencies of
SARS-CoV-2 variants in pooled samples. This model is built by considering a
previously defined selection of genomic polymorphisms that characterize
SARS-CoV-2 variants. The methods described here support both raw sequencing
reads for polymorphisms-based markers calling and predefined markers in the VCF
format. Results obtained by using simulated data show that our method is quite
effective in recovering the correct variant proportions. Further, results
obtained by considering longitudinal data from wastewater samples of two
locations in Switzerland agree well with those describing the epidemiological
evolution of COVID-19 variants in clinical samples of these locations. Our
results show that the described method can be a valuable tool for tracking the
proportions of SARS-CoV-2 variants.


  We consider the dynamics of two coupled three-species population patches,
incorporating the Allee Effect, focussing on the onset of extreme events in the
coupled system. First we show that the interplay between coupling and the Allee
effect may change the nature of the dynamics, with regular periodic dynamics
becoming chaotic in a range of Allee parameters and coupling strengths.
Further, the growth in the vegetation population displays an explosive blow-up
beyond a critical value of coupling strength and Allee parameter. Most
interestingly, we observe that beyond a threshold of coupling strength and
Allee parameter, the population densities of all three species exhibit non-zero
probability of yielding extreme events. The emergence of extreme events in the
predator populations in the patches is the most prevalent, and the probability
of obtaining large deviations in the predator populations is not affected
significantly by either the coupling strength or the Allee effect. In the
absence of the Allee effect the prey population in the coupled system exhibits
no extreme events for low coupling strengths, but yields a sharp increase in
extreme events after a critical strength of coupling. The vegetation population
in the patches display a small finite probability of extreme events for strong
enough coupling, only in the presence of Allee effect. Lastly we consider the
influence of additive noise on the continued prevalence of extreme events. Very
significantly, we find that noise suppresses the unbounded vegetation growth
that was induced by a combination of Allee effect and coupling. Further, we
demonstrate that noise mitigates extreme events in all three populations, and
beyond a noise level we do not observe any extreme events in the system at all.
This finding has important bearing on the potential observability of extreme
events in natural and laboratory systems.


  This paper proposes an iterative algorithm for choosing gradient directions
use to reconstruct white matter fibers in the brain. The present study is not
focusing on data acquisition where scanning is performed. The Adaptive Gradient
Directions (AGD) approach is extended to refine the position and area of the
grid, resulting in an admissible reduction in angular error. We begin with the
gradient directions distributed uniformly inside a grid of bigger size and with
larger spacing between the points. Both (size of the grid and spacing between
the points) reduce iteratively. The proposed algorithm ensures that the actual
position of fiber comes inside the grid at each iteration, unlike as in the AGD
approach. As a result, the solution tends to actual orientation in each
iteration followed by better estimation of fibers. The proposed algorithm is
validated by associating it with mixture of Gaussian diffusion and mixture of
non-central Wishart distribution models. The proposed approach significantly
reduce the angular error for multiple computer-generated experiments on
synthetic simulations and real data. Moreover, we have also performed
simulations with fibers not residing in the XY-plane. For this set-up also, the
proposed work outperforms, giving lesser angular error with both the models.
Synthetic simulations have been performed with Rician distributed (R-D) noise
of standard deviation ranging from 0.02-0.1. This work helps in better
understanding of the anatomy of the brain using the MRI signal data.


  This paper presents a disease-severity-structured epidemic model with
treatment necessary only to severe infective individuals to discuss the effect
of the treatment capacity on the disease transmission. It is shown that a
backward bifurcation occurs in the basic reproduction number R0, where a stable
endemic equilibrium co-exists with a stable disease-free equilibrium when R0
<1, if the capacity is relatively small. This epidemiological implication is
that, when there is not enough capacity for treatment, the requirement R0 <1 is
not sufficient for effective disease control and disease outbreak can happen to
a high endemic level even though R0 <1.


  With winter coming in the northern hemisphere, disadvantageous seasonality of
SARS-CoV-2 requires high immunity levels in the population or increasing
non-pharmaceutical interventions (NPIs), compared to summer. Otherwise
intensive care units (ICUs) might fill up. However, compliance with mandatory
NPIs, vaccine uptake, and individual protective measures depend on individuals'
opinions and behavior. Opinions, in turn, depend on information, e.g., about
vaccine safety or current infection levels. Therefore, understanding how
information about the pandemic affects its spread through the modulation of
voluntary protection-seeking behaviors is crucial for better preparedness this
winter and for future crises.


  The RecA protein and RecBCD complex are key bacterial components for the
maintenance and repair of DNA, RecBCD a helicase-nuclease that uses homologous
recombination to resolve double-stranded DNA breaks and also facilitating
decoration of single-stranded DNA with RecA to form RecA filaments, a vital
step in the double-stranded break DNA repair pathway. However, questions remain
about the mechanistic roles of RecA and RecBCD in live cells. Here, we use
millisecond super-resolved fluorescence microscopy to pinpoint the spatial
localization of fluorescent reporters of RecA and the RecB at physiological
levels of expression in individual live Escherichia coli cells. By introducing
the DNA crosslinker mitomycin C, we induce DNA damage and quantify the
resulting changes in stoichiometry, copy number and molecular mobilities of
RecA and RecB. We find that both proteins accumulate in molecular hotspots to
effect repair, resulting in RecA filamental stoichiometries equivalent to
several hundred molecules that act largely in RecA tetramers before DNA damage,
but switch to approximately hexameric subunits when mature filaments are
formed. Unexpectedly, we find that the physiologically predominant form of RecB
is a dimer.


  Covid-19 is one of the most dreaded pandemics/epidemics in the world
threatening the human population. The dynamics of this pandemic is quite
complicated and prediction of pandemic states often fails. In this work, we
study and correlate the SIR epidemiological model with the ongoing pandemic and
found that pandemic dynamics and states are quite sensitively dependent on
model parameters. The analysis of the exact parametric solution of the
deterministic SIR model shows that the fixed points ($g^*$) depend on the SIR
parameters, where, if $g^*\langle\frac{\gamma}{\beta x_0}$ then $g^*$ is stable
and the pandemic can be controlled, whereas, if $g^*\rangle\frac{\gamma}{\beta
x_0}$ then $g^*$ is unstable indicating active pandemic state, and $g^*=0$
corresponds to an endemic state. The dynamics show asymptotic stability
bifurcation. The analytical solution of the stochastic SIR model allows to
predict endemic time $t_E\propto\frac{1}{\gamma}$ and mean infected population
becomes constant for small values of $\gamma$. The noise in the system even can
modulate pandemic dynamics. Hence, we estimated the values of these parameters
during lockdown periods using the ABC SSA algorithm and found strong dependence
on control strategy, namely, lockdown and follows power-law behaviors. The
numerical results show strong agreement with the real data and have the
positive effect of social distancing or lockdown in the pandemic control.


  The viability of a new drug molecule is a time and resource intensive task
that makes computer-aided assessments a vital approach to rapid drug discovery.
Here we develop a machine learning algorithm, iMiner, that generates novel
inhibitor molecules for target proteins by combining deep reinforcement
learning with real-time 3D molecular docking using AutoDock Vina, thereby
simultaneously creating chemical novelty while constraining molecules for shape
and molecular compatibility with target active sites. Moreover, through the use
of various types of reward functions, we can generate new molecules that are
chemically similar to a target ligand, which can be grown from known protein
bound fragments, as well as to create molecules that enforce interactions with
target residues in the protein active site. The iMiner algorithm is embedded in
a composite workflow that filters out Pan-assay interference compounds,
Lipinski rule violations, and poor synthetic accessibility, with options for
cross-validation against other docking scoring functions and automation of a
molecular dynamics simulation to measure pose stability. Because our approach
only relies on the structure of the target protein, iMiner can be easily
adapted for future development of other inhibitors or small molecule
therapeutics of any target protein.


  Recent decades have seen a rise in the use of physics methods to study
different societal phenomena. This development has been due to physicists
venturing outside of their traditional domains of interest, but also due to
scientists from other disciplines taking from physics the methods that have
proven so successful throughout the 19th and the 20th century. Here we dub this
field 'social physics' and pay our respect to intellectual mavericks who
nurtured it to maturity. We do so by reviewing the current state of the art.
Starting with a set of topics that are at the heart of modern human societies,
we review research dedicated to urban development and traffic, the functioning
of financial markets, cooperation as the basis for our evolutionary success,
the structure of social networks, and the integration of intelligent machines
into these networks. We then shift our attention to a set of topics that
explore potential threats to society. These include criminal behaviour,
large-scale migrations, epidemics, environmental challenges, and climate
change. We end the coverage of each topic with promising directions for future
research. Based on this, we conclude that the future for social physics is
bright. Physicists studying societal phenomena are no longer a curiosity, but
rather a force to be reckoned with. Notwithstanding, it remains of the utmost
importance that we continue to foster constructive dialogue and mutual respect
at the interfaces of different scientific disciplines.


  Cooperative behaviour has been extensively studied as a choice between
cooperation and defection. However, the possibility to not participate is also
frequently available. This type of problem can be studied through the optional
public goods game. The introduction of the "Loner" strategy, allows players to
withdraw from the game, which leads to a cooperator-defector-loner cycle. While
prosocial punishment can help increase cooperation, anti-social punishment --
where defectors punish cooperators -- causes its downfall in both experimental
and theoretical studies.
  In this paper, we introduce social norms that allow agents to condition their
behaviour to the reputation of their peers. We benchmark this both with respect
to the standard optional public goods game and to the variant where all types
of punishment are allowed. We find that a social norm imposing a more moderate
reputational penalty for opting out than for defecting, increases cooperation.
When, besides reputation, punishment is also possible, the two mechanisms work
synergically under all social norms that do not assign to loners a strictly
worse reputation than to defectors. Under this latter setup, the high levels of
cooperation are sustained by conditional strategies, which largely reduce the
use of pro-social punishment and almost completely eliminate anti-social
punishment.


  Computational methods for predicting the interface contacts between proteins
come highly sought after for drug discovery as they can significantly advance
the accuracy of alternative approaches, such as protein-protein docking,
protein function analysis tools, and other computational methods for protein
bioinformatics. In this work, we present the Geometric Transformer, a novel
geometry-evolving graph transformer for rotation and translation-invariant
protein interface contact prediction, packaged within DeepInteract, an
end-to-end prediction pipeline. DeepInteract predicts partner-specific protein
interface contacts (i.e., inter-protein residue-residue contacts) given the 3D
tertiary structures of two proteins as input. In rigorous benchmarks,
DeepInteract, on challenging protein complex targets from the 13th and 14th
CASP-CAPRI experiments as well as Docking Benchmark 5, achieves 14% and 1.1%
top L/5 precision (L: length of a protein unit in a complex), respectively. In
doing so, DeepInteract, with the Geometric Transformer as its graph-based
backbone, outperforms existing methods for interface contact prediction in
addition to other graph-based neural network backbones compatible with
DeepInteract, thereby validating the effectiveness of the Geometric Transformer
for learning rich relational-geometric features for downstream tasks on 3D
protein structures.


  Linus Pauling contributions span structural biology, chemistry in its
broadest definition, quantum mechanical theory, valence bond theory, and even
nuclear physics. A principal tool developed and used by Pauling is Xray, and
electron, diffraction. One possible extension of the Pauling oeuvre could be
the marriage of crystallography and quantum mechanics. Such an effort dates
back to the sixties and has now flourished into an entire subfield termed
quantum crystallography. Quantum crystallography could be achieved through the
application of Clinton equations to yield N-representable density matrices
consistent with experimental data. The implementation of the Clinton equations
is qualitatively different for small and for large systems. For a small system,
quantum mechanics is extracted from Xray data while for a large system, the
quantum mechanics is injected into the system. In both cases, Nrepresentability
is imposed by the use of the Clinton equations.


  We study two fringe subtree counting statistics, the number of cherries and
that of pitchforks for Ford's $\alpha$ model, a one-parameter family of random
phylogenetic tree models that includes the uniform and the Yule models, two
tree models commonly used in phylogenetics. Based on a nonuniform version of
the extended P\'olya urn models in which negative entries are permitted for
their replacement matrices, we obtain the strong law of large numbers and the
central limit theorem for the joint distribution of these two count statistics
for the Ford model. Furthermore, we derive a recursive formula for computing
the exact joint distribution of these two statistics. This leads to exact
formulas for their means and higher order asymptotic expansions of their second
moments, which allows us to identify a critical parameter value for the
correlation between these two statistics. That is, when $n$ is sufficiently
large, they are negatively correlated for $0\le \alpha \le 1/2$ and positively
correlated for $1/2<\alpha<1$.


  The Rock-Paper-Scissors (RPS) model successfully reproduces some of the main
features of simple cyclic predator-prey systems with interspecific competition
observed in nature. Still, lattice-based simulations of the spatial stochastic
RPS model are known to give rise to significantly different results, depending
on whether the three state Lotka-Volterra or the four state May-Leonard
formulation is employed. This is true independently of the values of the model
parameters and of the use of either a von Neumann or a Moore neighborhood. With
the objective of reducing the impact of the use of a discrete lattice, in this
paper we introduce a simple modification to the standard spatial stochastic RPS
model in which the range of the search of the nearest neighbor may be extended
up to a maximum euclidean radius $R$. We show that, with this adjustment, the
Lotka-Volterra and May-Leonard formulations can be designed to produce similar
results, both in terms of dynamical properties and spatial features, by means
of an appropriate parameter choice. In particular, we show that this modified
spatial stochastic RPS model naturally leads to the emergence of spiral
patterns in both its three and four state formulations.