  The eigenvalue spectrum of a random matrix often only depends on the first
and second moments of its elements, but not on the specific distribution from
which they are drawn. The validity of this universality principle is often
assumed without proof in applications. In this letter, we offer a pertinent
counterexample in the context of the generalised Lotka-Volterra equations.
Using dynamic mean-field theory, we derive the statistics of the interactions
between species in an evolved ecological community. We then show that the full
statistics of these interactions, beyond those of a Gaussian ensemble, are
required to correctly predict the eigenvalue spectrum and therefore stability.
Consequently, the universality principle fails in this system. Our findings
connect two previously disparate ways of modelling complex ecosystems: Robert
May's random matrix approach and the generalised Lotka-Volterra equations. We
show that the eigenvalue spectra of random matrices can be used to deduce the
stability of dynamically constructed (or `feasible') communities, but only if
the emergent non-Gaussian statistics of the interactions between species are
taken into account.


  We suggest a new hardcore Poisson-type distribution for Young diagrams with
the row lengths from some finite list. A discrete variant of the time-ordered
Mat\'{e}rn II process in 1D is employed. This approach is related to that based
on the interlacing sequences due to Kerov and others, but we restrict the
number of rows. The basic lengths are assumed comparable with the total order
of the diagram in the quasi-classical limit, which results in new methods and
new formulas. An interesting application is to random walks where the steps are
at the points satisfying the classical Poisson distribution or our
truncatedone. In the simplest case, one obtains the distribution in terms of
Bessel I-functions, which provides some probabilistic interpretation of its
many properties. An immediate application of our truncated Poisson
distributions is to modeling reinfections in epidemics.


  In the recent decade, the emergence of liquid biopsy has significantly
improved cancer monitoring and detection. Dying cells, including those
originating from tumors, shed their DNA into the bloodstream and contribute to
a pool of circulating fragments called cell-free DNA (cfDNA). Identifying the
tissue origin of these DNA fragments from their epigenetic features has
implications in various clinical contexts. Open chromatin regions (OCRs) are
important epigenetic features of DNA that reflect cell types of origin.
Profiling these features by DNase-seq, ATAC-seq, and histone ChIP-seq provides
insights into tissue-specific and disease-specific regulatory mechanisms.
Integration of genomic and epigenomic features for cancer detection by liquid
biopsy has previously been reported. However, many multimodal analyses require
large amounts of cfDNA input and/or multiple types of experiments to cover the
genomic and epigenomic aspects of a single sample which is cost and time
prohibitive. Thus, methods that capture genomic and epigenomic profiles in a
single experiment type with low input requirements are of importance.
Predicting OCRs from whole genome sequencing (WGS) data is one such approach.
Here, we applied a correlation clustering algorithm to predict OCRs. We used
local sequencing depth as input to our algorithm. Multiple processing steps
were then applied as follows: count normalization, discrete Fourier transform
conversion, graph construction, graph cut optimization by linear programming,
and clustering. To validate the proposed method, we compared the output of our
predictions (OCR vs. non-OCR) with previously validated open chromatin regions
related to human blood samples of the ATAC-db. The percentage of overlap
between them is greater than 67%.


  Learning and reasoning about 3D molecular structures with varying size is an
emerging and important challenge in machine learning and especially in drug
discovery. Equivariant Graph Neural Networks (GNNs) can simultaneously leverage
the geometric and relational detail of the problem domain and are known to
learn expressive representations through the propagation of information between
nodes leveraging higher-order representations to faithfully express the
geometry of the data, such as directionality in their intermediate layers. In
this work, we propose an equivariant GNN that operates with Cartesian
coordinates to incorporate directionality and we implement a novel attention
mechanism, acting as a content and spatial dependent filter when propagating
information between nodes. We demonstrate the efficacy of our architecture on
predicting quantum mechanical properties of small molecules and its benefit on
problems that concern macromolecular structures such as protein complexes.


  The rooted subtree prune and regraft (rSPR) distance between two rooted
binary phylogenetic trees is a well-studied measure of topological
dissimilarity that is NP-hard to compute. Here we describe an improved linear
kernel for the problem. In particular, we show that if the classical subtree
and chain reduction rules are augmented with a modified type of chain reduction
rule, the resulting trees have at most 9k-3 leaves, where k is the rSPR
distance; and that this bound is tight. The previous best-known linear kernel
had size O(28k). To achieve this improvement we introduce cyclic generators,
which can be viewed as cyclic analogues of the generators used in the
phylogenetic networks literature. As a corollary to our main result we also
give an improved weighted linear kernel for the minimum hybridization problem
on two rooted binary phylogenetic trees.


  Levels of sociality in nature vary widely. Some species are solitary; others
live in family groups; some form complex multi-family societies. Increased
levels of social interaction can allow for the spread of useful innovations and
beneficial information, but can also facilitate the spread of harmful
contagions, such as infectious diseases. It is natural to assume that these
contagion processes shape the evolution of complex social systems, but an
explicit account of the dynamics of sociality under selection pressure imposed
by contagion remains elusive.
  We consider a model for the evolution of sociality strategies in the presence
of both a beneficial and costly contagion. We study the dynamics of this model
at three timescales: using a susceptible-infectious-susceptible (SIS) model to
describe contagion spread for given sociality strategies, a replicator equation
to study the changing fractions of two different levels of sociality, and an
adaptive dynamics approach to study the long-time evolution of the population
level of sociality.
  For a wide range of assumptions about the benefits and costs of infection, we
identify a social dilemma: the evolutionarily-stable sociality strategy (ESS)
is distinct from the collective optimum -- the level of sociality that would be
best for all individuals. In particular, the ESS level of social interaction is
greater (respectively less) than the social optimum when the good contagion
spreads more (respectively less) readily than the bad contagion.
  Our results shed light on how contagion shapes the evolution of social
interaction, but reveals that evolution may not necessarily lead populations to
social structures that are good for any or all.


  We introduce cognitive-affective maps (CAMs) as a novel tool to assess
individual experiences and belief systems. CAMs were first presented by the
cognitive scientist and philosopher Paul Thagard as a graphical representation
of a mental network, visualizing attitudes, thoughts, and affective
connotations toward the topic of interest. While CAMs were originally used
primarily to visualize existing data, the recent release of the new software
tool Valence has facilitated the applicability of CAMs for empirical data
collection. In this article, we explain the concept and the theoretical
background of CAMs. We exemplify how CAMs can be applied in research practice,
including different options for analysis. We propose CAMs as a user-friendly
and versatile methodological bridge between qualitative and quantitative
research approaches and encourage incorporating the method into studies to
access and visualize human attitudes and experience.


  We present a new method for analyzing stochastic epidemic models under
minimal assumptions. The method, dubbed DSA, is based on a simple yet powerful
observation, namely that population-level mean-field trajectories described by
a system of PDE may also approximate individual-level times of infection and
recovery. This idea gives rise to a certain non-Markovian agent-based model and
provides an agent-level likelihood function for a random sample of infection
and/or recovery times. Extensive numerical analyses on both synthetic and real
epidemic data from the FMD in the United Kingdom and the COVID-19 in India show
good accuracy and confirm method's versatility in likelihood-based parameter
estimation. The accompanying software package gives prospective users a
practical tool for modeling, analyzing and interpreting epidemic data with the
help of the DSA approach.


  Proteins control many vital functions in living cells, such as cell growth
and cell division. Reliable coordination of these functions requires the
spatial and temporal organizaton of proteins inside cells, which encodes
information about the cell's geometry and the cell-cycle stage. Such protein
patterns arise from protein transport and reaction kinetics, and they can be
controlled by various guiding cues within the cell. Here, we review how protein
patterns are guided by cell size and shape, by other protein patterns that act
as templates, and by the mechanical properties of the cell. The basic
mechanisms of guided pattern formation are elucidated with reference to recent
observations in various biological model organisms. We posit that understanding
the controlled formation of protein patterns in cells will be an essential part
of understanding information processing in living systems.


  George Price introduced his famous equation to study selective and
environmental effects in discrete populations. We extend Price's framework to
the measurable and quantum cases, decomposing all evolutionary processes into
selective and environmental components. We also extend Fisher's fundamental
theorem, showing that selective change of relative fitness equals variance of
relative fitness.
  We introduce novel selective and environmental entropy functionals. Selective
entropy is non-positive, representing biological negentropy, and environmental
entropy is non-negative, representing physical entropy. Environmental entropy
further decomposes into dispersion and mixing entropies.
  We prove four novel Laws of Natural Selection, showing that selection
consistently acts to increase selection, but can be disrupted by environmental
change. We apply convex analysis to variance and entropy functionals and their
selective changes, and equilibrium processes arise to optimize these
inequalities. These laws are inspired by but distinct from the classical
Thermodynamic Laws.
  Our Zeroth Law is a refinement of Fisher's theorem, showing that variance of
relative fitness is bounded below by a quantity depending on the child-bearing
population. Our First Law shows that selective acceleration of relative fitness
is also bounded below, depending on the variance. This is a non-conservative,
selective version of the Thermodynamic First Law. Our Second Law shows that the
selective change of selective entropy and its selective acceleration are
similarly bounded by non-positive constants. This is a formal, rigorous version
of the Thermodynamic Second Law. Our Third Law shows that for a class of
equilibrium processes, selective change of environmental entropy vanishes, and
otherwise may vary in an open window around zero. This is a selective version
of the Third Law of Thermodynamics.


  Read mapping is a fundamental, yet computationally-expensive step in many
genomics applications. It is used to identify potential matches and differences
between fragments (called reads) of a sequenced genome and an already known
genome (called a reference genome). To address the computational challenges in
genome analysis, many prior works propose various approaches such as filters
that select the reads that must undergo expensive computation, efficient
heuristics, and hardware acceleration. While effective at reducing the
computation overhead, all such approaches still require the costly movement of
a large amount of data from storage to the rest of the system, which can
significantly lower the end-to-end performance of read mapping in conventional
and emerging genomics systems.
  We propose GenStore, the first in-storage processing system designed for
genome sequence analysis that greatly reduces both data movement and
computational overheads of genome sequence analysis by exploiting low-cost and
accurate in-storage filters. GenStore leverages hardware/software co-design to
address the challenges of in-storage processing, supporting reads with 1)
different read lengths and error rates, and 2) different degrees of genetic
variation. Through rigorous analysis of read mapping processes, we meticulously
design low-cost hardware accelerators and data/computation flows inside a NAND
flash-based SSD. Our evaluation using a wide range of real genomic datasets
shows that GenStore, when implemented in three modern SSDs, significantly
improves the read mapping performance of state-of-the-art software (hardware)
baselines by 2.07-6.05$\times$ (1.52-3.32$\times$) for read sets with high
similarity to the reference genome and 1.45-33.63$\times$ (2.70-19.2$\times$)
for read sets with low similarity to the reference genome.


  We introduce EINNs, a framework crafted for epidemic forecasting that builds
upon the theoretical grounds provided by mechanistic models as well as the
data-driven expressibility afforded by AI models, and their capabilities to
ingest heterogeneous information. Although neural forecasting models have been
successful in multiple tasks, predictions well-correlated with epidemic trends
and long-term predictions remain open challenges. Epidemiological ODE models
contain mechanisms that can guide us in these two tasks; however, they have
limited capability of ingesting data sources and modeling composite signals.
Thus, we propose to leverage work in physics-informed neural networks to learn
latent epidemic dynamics and transfer relevant knowledge to another neural
network which ingests multiple data sources and has more appropriate inductive
bias. In contrast with previous work, we do not assume the observability of
complete dynamics and do not need to numerically solve the ODE equations during
training. Our thorough experiments on all US states and HHS regions for
COVID-19 and influenza forecasting showcase the clear benefits of our approach
in both short-term and long-term forecasting as well as in learning the
mechanistic dynamics over other non-trivial alternatives.


  We study the qualitative behavior of a model to represent local regulation in
a metabolic network. The model is based on the end-product control structure
introduced in [A. Goelzer, F. Bekkal Brikci, I. Martin-Verstraete et al., BMC
Syst Biol 2 (2008), pp. 1--18]. In this class of regulation, the metabolite
effector is the end-product of a metabolic pathway. We suppose the input to the
pathway to switch between zero and a positive value according to the
concentration of the metabolite effector. Considering the switched system as a
differential inclusion, we prove that it converges to a globally uniformly
asymptotically stable equilibrium point, reaches the sliding mode or oscillates
around the sliding mode depending on the positive value of the input. Finally,
we show that in any case the solution of the switched system is the limit of
solutions of equation sequences with smooth or piecewise linear inputs.


  Heterogeneity in characteristics from one region (sub-population) to another,
commonly observed in complex systems, such as glasses and a collection of
cells, is hard to describe theoretically. In the context of cancer, intra-tumor
heterogeneity (ITH), characterized by cells with genetic and phenotypic
variability that co-exist within a single tumor, is often the cause of
ineffective therapy and recurrence of cancer. Multi-region sequencing (M-Seq),
obtained by sampling multiple regions of a single tumor, has vividly
demonstrated the pervasive nature of ITH, raising the need for a theory that
accounts for evolution of tumor heterogeneity. Here, we develop a statistical
mechanical theory to quantify ITH, using the Hamming distance, between genetic
mutations in distinct regions within a single tumor. An analytic expression for
ITH, expressed in terms of cell division probability ($\alpha$) and mutation
probability ($p$), is validated using cellular-automaton type simulations.
Application of the theory successfully captures ITH extracted from M-seq data
in patients with exogenous cancers (melanoma and lung). The theory, based on
punctuated evolution at the early stages of the tumor followed by neutral
evolution, is accurate provided the spatial variation in the tumor mutation
burden is not large. We show that there are substantial variations in ITH in
distinct regions of a single solid tumor, which supports the notion that
distinct subclones could co-exist. The simulations show that there are
substantial variations in the sub-populations, with the ITH increasing as the
distance between the regions increases. The analytical and simulation framework
developed here could be used in the quantitative analyses of the experimental
(M-Seq) data. More broadly, our theory is likely to be useful in analyzing
dynamic heterogeneity in complex systems such as super-cooled liquids.


  Protein-protein binding enables orderly and lawful biological
self-organization, and is therefore considered a miracle of nature.
Protein-protein binding is steered by electrostatic forces, hydrogen bonding,
van der Waals force, and hydrophobic interactions. Among these physical forces,
only the hydrophobic interactions can be considered as long-range
intermolecular attractions between proteins in intracellular and extracellular
fluid. Low-entropy regions of hydration shells around proteins drive
hydrophobic attraction among them that essentially coordinate protein-protein
docking in rotational-conformational space of mutual orientations at the
guidance stage of the binding. Here, an innovative method was developed for
identifying the low-entropy regions of hydration shells of given proteins, and
we discovered that the largest low-entropy regions of hydration shells on
proteins typically cover the binding sites. According to an analysis of
determined protein complex structures, shape matching between the largest
low-entropy hydration shell region of a protein and that of its partner at the
binding sites is revealed as a regular pattern. Protein-protein binding is thus
found to be mainly guided by hydrophobic collapse between the shape-matched
low-entropy hydration shells that is verified by bioinformatics analyses of
hundreds of structures of protein complexes. A simple algorithm is developed to
precisely predict protein binding sites.


  Proteins need to selectively interact with specific targets among a multitude
of similar molecules in the cell. But despite a firm physical understanding of
binding interactions, we lack a general theory of how proteins evolve high
specificity. Here, we present such a model that combines chemistry, mechanics
and genetics, and explains how their interplay governs the evolution of
specific protein-ligand interactions. The model shows that there are many
routes to achieving molecular discrimination - by varying degrees of
flexibility and shape/chemistry complementarity - but the key ingredient is
precision. Harder discrimination tasks require more collective and precise
coaction of structure, forces and movements. Proteins can achieve this through
correlated mutations extending far from a binding site, which fine-tune the
localized interaction with the ligand. Thus, the solution of more complicated
tasks is enabled by increasing the protein size, and proteins become more
evolvable and robust when they are larger than the bare minimum required for
discrimination. The model makes testable, specific predictions about the role
of flexibility and shape mismatch in discrimination, and how evolution can
independently tune affinity and specificity. Thus, the proposed theory of
specific binding addresses the natural question of "why are proteins so big?".
A possible answer is that molecular discrimination is often a hard task best
performed by adding more layers to the protein.


  Robust and efficient interpretation of QSAR methods is quite useful to
validate AI prediction rationales with subjective opinion (chemist or biologist
expertise), understand sophisticated chemical or biological process mechanisms,
and provide heuristic ideas for structure optimization in pharmaceutical
industry. For this purpose, we construct a multi-layer self-attention based
Graph Neural Network framework, namely Ligandformer, for predicting compound
property with interpretation. Ligandformer integrates attention maps on
compound structure from different network blocks. The integrated attention map
reflects the machine's local interest on compound structure, and indicates the
relationship between predicted compound property and its structure. This work
mainly contributes to three aspects: 1. Ligandformer directly opens the
black-box of deep learning methods, providing local prediction rationales on
chemical structures. 2. Ligandformer gives robust prediction in different
experimental rounds, overcoming the ubiquitous prediction instability of deep
learning methods. 3. Ligandformer can be generalized to predict different
chemical or biological properties with high performance. Furthermore,
Ligandformer can simultaneously output specific property score and visible
attention map on structure, which can support researchers to investigate
chemical or biological property and optimize structure efficiently. Our
framework outperforms over counterparts in terms of accuracy, robustness and
generalization, and can be applied in complex system study.


  The tryptophan (trp) operon in E. coli codes for the proteins responsible for
the synthesis of the amino acid tryptophan from chorismic acid, and has been
one of the most well-studied gene networks since its discovery in the 1960s.
The tryptophanase (tna) operon codes for proteins needed to transport and
metabolize it. Both of these have been modeled individually with delay
differential equations under the assumption of mass-action kinetics. Recent
work has provided strong evidence for bistable behavior of the tna operon. The
authors of (Orozco, 2019) identified a medium range of tryptophan in which the
system has two stable steady-states, and they reproduced these experimentally.
In this paper, we will show how a Boolean model can capture this bistability.
We will also develop and analyze a Boolean model of the trp operon. Finally, we
will combine these two to create a single Boolean model of the transport,
synthesis, and metabolism of tryptophan. In this amalgamated model, the
bistability disappears, presumably reflecting the ability of the trp operon to
produce tryptophan and drive the system toward homeostasis. All of these models
have longer attractors that we call "artifacts of synchrony", which disappear
in the asynchronous automata. This curiously matches the behavior of a recent
Boolean model of the arabinose operon in E. coli, and we discuss some
open-ended questions that arise along these lines.


  In this paper we investigate the homeostatic mechanism in two biologically
motivated models: intracellular copper regulation and self immune recognition.
The analysis is based on the notions of infinitesimal homeostasis and
near-perfect homeostasis. We introduce a theoretical background that makes it
possible to consider points of infinitesimal homeostasis that lie at the
boundary of the domain of definition of the input-output function. We show that
the two models display near-perfect homeostasis. Moreover, we show that, unlike
the examples of [M. Reed, J. Best, M. Golubitsky, I. Stewart, and H. F.
Nijhout. Analysis of homeostatic mechanisms in biochemical networks. Bull.
Math. Biol., 79(11):2534-2557, 2017], the obstruction of occurrence of
infinitesimal homeostasis in both of them is caused by the modeling assumptions
that force the point of infinitesimal homeostasis to lie at the boundary of
domain of definition of the respective input-output functions.


  The main idea in this paper is that the age associated with reported measles
cases can be used to estimate the number of undetected measles infections.
Somewhat surprisingly, even with age only to the nearest year, estimates of
underreporting can be generated at the much faster, 2 week time-scale
associated with measles transmission. I describe this idea by focusing on the
well-studied, 60 city United Kingdom data set, which covers the transition to
universal healthcare in 1948, and is, as a result, an interesting case study in
infectious disease surveillance. Finally, at the end of the paper, I comment
briefly on how the approach can be modified for application to modern contexts.


  One clear aspect of behaviour in the COVID-19 pandemic has been people's
focus on, and response to, reported or observed infection numbers in their
community. We describe a simple model of infectious disease spread in a
pandemic situation where people's behaviour is influenced by the current risk
of infection and where this behavioural response acts homeostatically to return
infection risk to a certain preferred level. This model predicts that the
reproduction rate $R$ will be centered around a median value of 1, and that a
related measure of relative change in the number of new infections will follow
the standard Cauchy distribution. Analysis of worldwide COVID-19 data shows
that the estimated reproduction rate has a median of 1, and that this measure
of relative change calculated from reported numbers of new infections closely
follows the standard Cauchy distribution at both an overall and an individual
country level.


  Mechanical properties of nucleic acids play an important role in many
biological processes which often involve physical deformations of these
molecules. At sufficiently long length scales (say above $\sim 20-30$ base
pairs) the mechanics of DNA and RNA double helices is described by a
homogeneous Twistable Wormlike Chain (TWLC), a semiflexible polymer model
characterized by twist and bending stiffnesses. At shorter scales this model
breaks down for two reasons: the elastic properties become sequence-dependent
and the mechanical deformations at distal sites gets coupled. We discuss in
this paper the origin of the latter effect using the framework of a non-local
Twistable Wormlike Chain (nlTWLC). We show, by comparing all-atom simulations
data for DNA and RNA double helices, that the non-local couplings are of very
similar nature in these two molecules: couplings between distal sites are
strong for tilt and twist degrees of freedom and weak for roll. We introduce
and analyze a simple double-stranded polymer model which clarifies the origin
of this universal distal couplings behavior. In this model, referred to as the
ladder model, a nlTWLC description emerges from the coarsening of local
(atomic) degrees of freedom into angular variables which describe the twist and
bending of the molecule. Differently from its local counterpart, the nlTWLC is
characterized by a length-scale-dependent elasticity. Our analysis predicts
that nucleic acids are mechanically softer at the scale of a few base pairs and
are asymptotically stiffer at longer length scales, a behavior which matches
experimental data.


  Geographic isolation is a central mechanism of speciation, but perfect
isolation of populations is rare. Although speciation can be hindered if gene
flow is large, intermediate levels of migration can enhance speciation by
introducing genetic novelty in the semi-isolated populations or founding small
communities of migrants. Here we consider a two island neutral model of
speciation with continuous migration and study diversity patterns as a function
of the migration probability, population size, and number of genes involved in
reproductive isolation (dubbed as genome size). For small genomes, low levels
of migration induce speciation on the islands that otherwise would not occur.
Diversity, however, drops sharply to a single species inhabiting both islands
as the migration probability increases. For large genomes, sympatric speciation
occurs even when the islands are strictly isolated. Then species richness per
island increases with the probability of migration, but the total number of
species decreases as they become cosmopolitan. For each genome size, there is
an optimal migration intensity for each population size that maximizes the
number of species. We discuss the observed modes of speciation induced by
migration and how they increase species richness in the insular system while
promoting asymmetry between the islands and hindering endemism.


  Motivation: 3D structures of proteins provide rich information for
understanding their biochemical roles. Identifying the representative protein
structures for protein sequences is essential for analysis of proteins at
proteome scale. However, there are technical difficulties in identifying the
representative structure of a given protein sequence and providing accurate
mapping of residue indices. Existing databases of mapping between structures
and sequences are usually static that are not suitable for studying proteomes
with frequent gene model revisions. They often do not provide reliable and
consistent representative structures that maximizes sequence coverage.
Furthermore, proteins isomers are usually not properly resolved. Results: To
overcome these difficulties, we have developed a computational pipeline called
SeqMapPDB to provide high-quality representative PDB structures of given
sequences. It provides mapping to structures that fully cover the sequences
when available, or to the set of partial non-overlapping structural domains
that maximally cover the query sequence. The residue indices are accurate
mapped and isomeric proteins are resolved. SeqMapPDB is efficient and can
rapidly carry out proteome-wide mapping to the selected version of reference
genomes in real-time. Furthermore, SeqMapPDB provides the flexibility of a
stand-alone pipeline for large scale mapping of in-house sequence and structure
data. Availability: Our method is available at
https://bitbucket.org/lianglabuic/seqmappdb with GNU GPL license.


  The reproductive habits of helminths are important for the study of the
dynamics of their transmission. For populations of parasites distributed by
Poisson or negative binomial models, these habits have already been studied.
However, there are other statistical models that describe these populations,
such as zero-inflated models, but where reproductive characteristics were not
analyzed. Using an arbitrary model for the parasite population, we model the
distribution of females and males per host, and from these we model the
different reproductive variables such as the mean number of fertile females,
the mean egg production, the mating probability, the mean fertilized egg
production. We show that these variables change due to the effects of a
negative density-dependence fecundity, a characteristic of helminth parasites.
We present the results obtained for some particular models.


  Natural ecosystems, in particular on the microbial scale, are inhabited by a
large number of species. The population size of each species is affected by
interactions of individuals with each other and by spatial and temporal changes
in environmental conditions, such as resource abundance. Here, we use a generic
population dynamics model to study how, and under what conditions, a periodic
temporal environmental variation can alter an ecosystem's composition and
biodiversity. We demonstrate that using time scale separation allows one to
qualitatively predict the long-term population dynamics of interacting species
in varying environments. We show that the notion of Tilman's R* rule, a
well-known principle that applies for constant environments, can be extended to
periodically varying environments if the time scale of environmental changes
(e.g., seasonal variations) is much faster than the time scale of population
growth (doubling time in bacteria). When these time scales are similar, our
analysis shows that a varying environment deters the system from reaching a
steady state, and stable coexistence between multiple species becomes possible.
Our results posit that biodiversity can in part be attributed to natural
environmental variations.


  Ecological interactions among phytoplankton occur in a moving fluid
environment. Oceanic flows can modulate the competition and coexistence between
phytoplankton populations, which in turn can affect ecosystem function and
biogeochemical cycling. We explore the impact of submesoscale velocity
gradients on phytoplankton ecology using observations, simulations, and theory.
Observations reveal that the relative abundance of Synechoccocus oligotypes
varies on 1--10 km scales at an ocean front with submesoscale velocity
gradients at the same scale. Simulations in realistic flow fields demonstrate
that regions of divergence in the horizontal flow field can substantially
modify ecological competition and dispersal on timescales of hours to days.
Regions of positive (negative) divergence provide an advantage (disadvantage)
to local populations, resulting in up to ~20% variation in community
composition in our model. We propose that submesoscale divergence is a
plausible contributor to observed taxonomic variability at oceanic fronts, and
can lead to regional variability in community composition.


  The reason of significantly higher transmissibility of SARS Covid (2019
CoV-2) compared to SARS Covid (2003 CoV) and MERS Covid (2012 MERS) can be
attributed to mutations reported in structural proteins, and the role played by
non-structural proteins (nsps) and accessory proteins (ORFs) for viral
replication, assembly, and shedding. Envelope protein E is one of the four
structural proteins of minimum length. Recent studies have confirmed critical
role played by the envelope protein in the viral life cycle including assembly
of virion exported from infected cell for its transmission. However, the
determinants of the highly complex viral - host interactions of envelope
protein, particularly with host Golgi complex, have not been adequately
characterized. CoV-2 and CoV Envelope proteins of length 75 and 76 amino acids
differ in four amino acid locations. The additional amino acid Gly (G) at
location 70 makes CoV length 76. The amino acid pair EG at location 69-70 of
CoV in place of amino acid R in location 69 of CoV-2, has been identified as a
major determining factor in the current investigation. This paper concentrates
on the design of computational model to compare the structure/function of wild
and mutants of CoV-2 with wild and mutants of CoV in the functionally important
region of the protein chain pair. We hypothesize that differences of CAML model
parameter of CoV-2 and CoV characterize the deviation in structure and function
of envelope proteins in respect of interaction of virus with host Golgi
complex; and this difference gets reflected in the difference of their
transmissibility. The hypothesis has been validated from single point
mutational study on- (i) human HBB beta-globin hemoglobin protein associated
with sickle cell anemia, (ii) mutants of envelope protein of Covid-2 infected
patients reported in recent publications.


  We consider interactions between players in groups of size $d\geq2$ with
payoffs that not only depend on the strategies used in the group but also
fluctuate at random over time. An individual can adopt either cooperation or
defection as strategy and the population is updated from one-time step to the
next by a birth-death event according to a Moran model. Assuming recurrent
symmetric mutation and payoffs with expected values, variances, and covariances
of the same small order, we derive a first-order approximation of the average
abundance of cooperation in the selection-mutation equilibrium. We show that
increasing the variance of any payoff for defection or decreasing the variance
of any payoff for cooperation increases the average abundance of cooperation.
As for the effect of the covariance between any payoff for cooperation and any
payoff for defection, we show that it depends on the number of cooperators in
the group associated with these payoffs. We study in particular the public
goods game, the stag hunt game, and the snowdrift game, all social dilemmas
based on random benefit $b$ and cost $c$ for cooperation. We show that a
decrease in the scaled variance of $b$ or $c$, or an increase in their scaled
covariance, makes it easier for weak selection to favor the abundance of
cooperation in the stag hunt game and the snowdrift game. The same conclusion
holds for the public goods game except that the covariance of $b$ has no effect
on the average abundance of $C$. On the other hand, increasing the scaled
mutation rate or the group size can enhance or lessen the condition for weak
selection to favor the abundance of $C$.


  The SARS CoV-2 virus emerged at the end of 2019 and rapidly developed several
mutated variants, specifically the Delta and Omicron, which demonstrate higher
infectivity and escalating infection cases worldwide. The dominant transmission
pathway of this virus is via human-to-human contact and aerosols, but another
possible route is through contact with surfaces contaminated with SARS-CoV-2,
often exhibiting long-period survival. Here we compare the adsorption
capacities of the S1 and S2 subunits of the spike (S) protein from the original
variant to that of the S1 subunit from the Delta and Omicron variants. The
results clearly show a significant difference in adsorption capacity between
the different variants, as well as between the S1 and S2 subunits. Overall, our
study demonstrates that while the Omicron variant is able to adsorb much more
successfully than the Delta, both variants show enhanced adsorption capacity
than the original strain. We also examined the influence of pH conditions on
the adsorption ability of the S1 subunit and found that adsorption was
strongest at pH 7.4, which is the physiological pH. The main conclusion of this
study is that there is a strong correlation between the adsorption capacity and
the infectivity of the various SARS CoV-2 variants.


  To understand the biodiversity of an ecosystem cannot be understood by solely
analyzing the pair relations of competing species. Instead, we should consider
multi-point interactions because the presence of a third party could change the
original microscopic outcome significantly. In this way an alliance may emerge
where species, who may have biased relations otherwise, can protect each other
from an external invader. Such an alliance can be formed by two, three or even
more species. By introducing a minimal model where six species compete for
space we here study how the size of an alliance determines the vitality of a
formation. We show that in the majority of parameter space the group of the
smallest size prevails and other solutions can only be observed in a limited
parameter range. These phases are separated by discontinuous phase transitions
which can only be identified by intensive numerical efforts due to serious
finite size effects and long relaxation processes.


  Repeat proteins are made with tandem copies of similar amino acid stretches
that fold into elongated architectures. Due to their symmetry, these proteins
constitute excellent model systems to investigate how evolution relates to
structure, folding and function. Here, we propose a scheme to map evolutionary
information at the sequence level to a coarse-grained model for repeat-protein
folding and use it to investigate the folding of thousands of repeat-proteins.
We model the energetics by a combination of an inverse Potts model scheme with
an explicit mechanistic model of duplications and deletions of repeats to
calculate the evolutionary parameters of the system at single residue level.
This is used to inform an Ising-like model that allows for the generation of
folding curves, apparent domain emergence and occupation of intermediate states
that are highly compatible with experimental data in specific case studies. We
analyzed the folding of thousands of natural Ankyrin-repeat proteins and found
that a multiplicity of folding mechanisms are possible. Fully cooperative
all-or-none transition are obtained for arrays with enough sequence-similar
elements and strong interactions between them, while non-cooperative
element-by-element intermittent folding arose if the elements are dissimilar
and the interactions between them are energetically weak. In between, we
characterised nucleation-propagation and multi-domain folding mechanisms.
Finally, we showed that stability and cooperativity of a repeat-array can be
quantitatively predicted from a simple energy score, paving the way for guiding
protein folding design with a co-evolutionary model.


  Climate-mediated changes in the spatiotemporal distribution of thermal stress
can destabilize animal populations and promote extinction risk. Using quantile,
spectral, and wavelet analyses of temperature projections from the latest
generation of earth system models, we show that significant regional
differences are expected to arise in the way that temperatures will increase
over time. When integrated into empirically-parameterized mathematical models
that simulate the dynamical and cumulative effects of thermal stress on the
performance of 38 global ectotherm species, the projected spatiotemporal
changes in temperature fluctuations are expected to give rise to complex
regional changes in population abundance and stability over the course of the
21st century. However, despite their idiosyncratic effects on stability,
projected temperatures universally increase extinction risk. These results show
that population changes under future climate conditions may be more extensive
and complex than the current literature suggests based on the statistical
relationship between biological performance and average temperature.


  We estimate patterns of covariation between COVID-19 vaccination rates and a
set of widely used indicators of human, social, and economic capital across 146
countries in July 2021 and February 2022. About 70% of the variability in
COVID-19 vaccination rates worldwide can be explained by differences in the
Human Development Index (HDI) and, specifically, in life expectancy at birth,
one year after the campaign debut. Trust in doctors and nurses adds predictive
value beyond the HDI, clarifying controversial discrepancies between
vaccination rates in countries with similar levels of human development and
vaccine availability. Cardiovascular disease deaths, an indicator of general
health system effectiveness, and infant measles immunization coverage, an
indicator of country-level immunization effectiveness, are also significant,
though weaker, predictors of COVID-19 vaccination success. The metrics of
economic inequality, perceived corruption, poverty, and inputs into the health
system have strong bivariate correlations with COVID-19 vaccination but no
longer remain statistically significant when controlling for the HDI. Our
analysis identified the contours of a social structure that sustains life and
is reproduced through this process. COVID-19 vaccines have proven to be part of
the Matthew effect of accumulating advantages and aggravating disadvantages
that the pandemic inflicted on societies and communities across the world. At
the same time, the remaining variability in vaccination success that cannot be
pinned down through these sets of metrics points to a considerable scope for
collective and individual agency in a time of crisis. The mobilization and
coordination in the vaccination campaigns of citizens, medical professionals,
scientists, journalists, and politicians, among others, account for at least
some of this variability in overcoming vaccine hesitancy and inequity.


  Addiction epidemiology has been an active area of mathematical research in
recent years. However, the social and mental processes involved in substance
use disorders versus contraction of a pathogenic disease have presented
challenges to advancing the epidemiological theory of substance abuse,
especially within the context of the opioids where both prescriptions and
social contagion have played a major role. In this paper, we utilize an
agent-based modeling approach on social networks to further explore these
dynamics. Using parameter estimation approaches, we compare our results to that
of the Phillips et al. SPAHR model which was previously fit to data from the
state of Tennessee. Our results show that the average path length of a social
network has a strong relationship to social contagion dynamics for drug use
initiation, while other pathways to substance use disorder should not be
constrained to social network interactions that predate the individual's drug
use.


  We study a nonhierarchical tritrophic system, whose predator-prey
interactions are described by the rock-paper-scissors game rules. In our
stochastic simulations, individuals may move strategically towards the
direction with more conspecifics to form clumps instead of moving aimlessly on
the lattice. Considering that the conditioning to move gregariously depends on
the organism's physical and cognitive abilities, we introduce a maximum
distance an individual can perceive the environment and a minimum conditioning
level to perform the gregarious movement. We investigate the pattern formation
and compute the average size of the single-species spatial domains emerging
from the grouping behaviour. The results reveal that the defence tactic reduces
the predation risk significantly, being more profitable if individuals perceive
further distances, thus creating bigger groups. Our outcomes show that the
species with more conditioned organisms dominate the cyclic spatial game,
controlling most of the territory. On the other hand, the species with fewer
individuals ready to perform aggregation strategy gives its predator the chance
to fill the more significant fraction of the grid. The spatial interactions
assumed in our numerical experiments constitute a data set that may help
biologists and data scientists understand how local interactions influence
ecosystem dynamics.


  Reaction networks are widely used models to describe biochemical processes.
Stochastic fluctuations in the counts of biological macromolecules have
amplified consequences due to their small population sizes. This makes it
necessary to favor stochastic, discrete population, continuous time models. The
stationary distributions provide snapshots of the model behavior at the
stationary regime, and as such finding their expression in terms of the model
parameters is of great interest. The aim of the present paper is to describe
when the stationary distributions of the original model, whose state space is
potentially infinite, coincide exactly with the stationary distributions of the
process truncated to finite subsets of states, up to a normalizing constant.
The finite subsets of states we identify are called copies and are inspired by
the modular topology of reaction network models. With such a choice we prove a
novel graphical characterization of the concept of complex balancing for
stochastic models of reaction networks. The results of the paper hold for the
commonly used mass-action kinetics but are not restricted to it, and are in
fact stated for more general setting.


  We introduce a stochastic household model for vector-borne diseases, in
particular as relevant to prominent vectors belonging to the Aedes genus and
hence the Zika, chikungunya, and dengue viruses. In this model, vectors remain
local to each household, while hosts mix for a proportion of their time in
their household and the remaining proportion in the population at random. This
is approximated with a two-type branching process, allowing us to efficiently
calculate a number of useful epidemiological characteristics, such as
reproductive numbers, early growth rates and household-type proportions,
offspring distributions, probabilities of a major outbreak, and
within-household final size distributions. We compare control interventions of
spraying -- reducing the number of vectors in each household -- and
social-distancing -- having individuals spend more time at home -- in terms of
these characteristics.


  Motivation: An essential step in developing computational tools for the
inference, optimization, and simulation of biochemical reaction networks is
gauging tool performance against earlier efforts using an appropriate set of
benchmarks. General strategies for the assembly of benchmark models include
collection from the literature, creation via subnetwork extraction and de novo
generation. However, with respect to biochemical reaction networks, these
approaches and their associated tools are either poorly suited to generate
models that reflect the wide range of properties found in natural biochemical
networks or to do so in numbers that enable rigorous statistical analysis.
Results: In this work we present SBbadger, a python-based software tool for the
generation of synthetic biochemical reaction or metabolic networks with
user-defined degree distributions, multiple available kinetic formalisms, and a
host of other definable properties. SBbadger thus enables the creation of
benchmark model sets that reflect properties of biological systems and generate
the kinetics and model structures typically targeted by computational analysis
and inference software. Here we detail the computational and algorithmic
workflow of SBbadger, demonstrate its performance under various settings,
provide samples outputs, and compare it to currently available biochemical
reaction network generation software.


  Modern biology and biomedicine are undergoing a big-data explosion needing
advanced computational algorithms to extract mechanistic insights on the
physiological state of living cells. We present the motivation for the Cell
Physiome: a framework and approach for creating, sharing, and using
biophysics-based computational models of single cell physiology. Using examples
in calcium signaling, bioenergetics, and endosomal trafficking, we highlight
the need for spatially detailed, biophysics-based computational models to
uncover new mechanisms underlying cell biology. We review progress and
challenges to date towards creating cell physiome models. We then introduce
bond graphs as an efficient way to create cell physiome models that integrate
chemical, mechanical, electromagnetic, and thermal processes while maintaining
mass and energy balance. Bond graphs enhance modularization and re-usability of
computational models of cells at scale. We conclude with a look forward into
steps that will help fully realize this exciting new field of mechanistic
biomedical data science.


  Studies of transient dynamics captures the time history of any dramatic
changes in the dynamics of a system. The transient dynamics is investigated
here in a classic ecological model of a bistable tri-trophic food chain. All
the species in the food chain either coexist or undergoes a partial extinction
with the loss of predator after a transient time that depends upon the initial
density of species populations. The distribution of transient time of
extinction of the predator in response to initial states shows interesting
pattern of inhomogeneity and anisotropy in the basin of predator-fee state.
Precisely, the distribution shows a multimodal character when the initial
points are near the basin boundary, and unimodal at locations far from the
border. The distribution is also anisotropic in the sense that the number of
modes depends on the direction from the basin location. We define two new
metrics, viz., homogeneity index and local isotropic index to characterize the
distinctive features of the distribution. The inhomogeneity of distribution
reduces significantly with resource enrichment, however, a similar change in
anisotropy is relatively low. We try to explain the origin of such multimodal
distributions from the dynamical system perspectives and present its ecological
implications.


  This paper investigates the effect of the structure of the contact network on
the dynamics of the epidemic outbreak. In particular, we focus on the peak
number of critically infected nodes (PCIN), determining the maximum workload of
intensive healthcare units which should be kept low. As a model and simulation
method, we develop a continuous-time Markov chain (CTMC) model and an efficient
simulation-based on Gillespie's Stochastic Simulation Algorithm (SSA). This
methods combine a realistic approximation of the stochastic process not relying
on the assumptions of mean-field models and large asymptotically large
population sizes as in differential equation models, and at the same time an
efficient way to simulate networks of moderate size. The approach is analysed
for different scenarios, based on data from the COVID-19 outbreak and
demographic data from Ukraine. From these results we extract network topology
features that need to be considered to effectively decrease the peak number of
infections. The CTMC simulation is implemented in python and integrated in a
dashboard that can be used for interactive exploration and it is made openly
available.


  Homogeneity across lineages is a common assumption in phylogenetics according
to which nucleotide substitution rates remain constant in time and do not
depend on lineages. This is a simplifying hypothesis which is often adopted to
make the process of sequence evolution more tractable. However, its validity
has been explored and put into question in several papers. On the other hand,
dealing successfully with the general case (heterogeneity across lineages) is
one of the key features of phylogenetic reconstruction methods based on
algebraic tools.
  The goal of this paper is twofold. First, we present a new weighting system
for quartets (ASAQ) based on algebraic and semi-algebraic tools, thus specially
indicated to deal with data evolving under heterogeneus rates. This method
combines the weights two previous methods by means of a test based on the
positivity of the branch length estimated with the paralinear distance. ASAQ is
statistically consistent when applied to GM data, considers rate and base
composition heterogeneity among lineages and does not assume stationarity nor
time-reversibility. Second, we test and compare the performance of several
quartet-based methods for phylogenetic tree reconstruction (namely, Quartet
Puzzling, Weight Optimization and Wilson's method) in combination with ASAQ
weights and other weights based on algebraic and semi-algebraic methods or on
the paralinear distance. These tests are applied to both simulated and real
data and support Weight Optimization with ASAQ weights as a reliable and
successful reconstruction method.


  Regulatory networks (RNs) are a well-accepted modelling formalism in
computational systems biology. The control of RNs is currently receiving a lot
of attention because it provides a computational basis for cell reprogramming
-- an attractive technology developed in regenerative medicine. By solving the
control problem, we learn which parts of a biological system should be
perturbed to stabilise the system in the desired phenotype.
  We allow the specification of the Boolean model representing a given RN to be
incomplete. To that end, we utilise the formalism of partially specified
Boolean networks which covers every possible behaviour of unspecified parts of
the system. Such an approach causes a significant state explosion. This problem
is addressed by using symbolic methods to represent both the unspecified model
parts and all possible perturbations of the system.
  Additionally, to make the control design efficient and practically
applicable, the optimal control should be minimal in terms of size. Moreover,
in a partially specified model, a control may achieve the desired stabilisation
only for a subset of the possible fully specified model instantiations. To
address these aspects, we utilise several quantitative measures. Apart from the
size of perturbation, we also examine its robustness -- a portion of
instantiations for which the control is applicable.
  We show that proposed symbolic methods solving the control problem for
partially specified BNs are efficient and scale well. We also evaluate the
robustness metrics in cases of all three studied control types. The robustness
metric tells us how big a proportion of fully defined systems the given
perturbation works. Our experiments support the hypothesis that one-step
perturbations may be less robust than temporary or permanent perturbations.
  This is a full version of a paper that is submitted to a journal.


  Coordinated responses to environmental stimuli are critical for multicellular
organisms. To overcome the obstacles of cell-to-cell heterogeneity and noisy
signaling dynamics within individual cells, cells must effectively exchange
information with peers. However, the dynamics and mechanisms of collective
information transfer driven by external signals is poorly understood. Here we
investigate the calcium dynamics of neuronal cells that form confluent
monolayers and respond to cyclic ATP stimuli in microfluidic devices. Using
Granger inference to reconstruct the underlying causal relations between the
cells, we find that the cells self-organize into spatially decentralized and
temporally stationary networks to support information transfer via gap junction
channels. The connectivity of the causal networks depend on the temporal
profile of the external stimuli, where short periods, or long periods with
small duty fractions, lead to reduced connectivity and fractured network
topology. We build a theoretical model based on communicating excitable units
that reproduces our observations. The model further predicts that connectivity
of the causal network is maximal at an optimal communication strength, which is
confirmed by the experiments. Together, our results show that information
transfer between neuronal cells is externally regulated by the temporal profile
of the stimuli, and internally regulated by cell-cell communication.


  Epidemic disease spreading is conventionally often modelled and analyzed by
means of rate and diffusion equations, following the paradigms of
well-controlled chemical reactions and diffusive dynamics in a test tube. Yet,
serious worries that this suggestive and appealing similarity might be a false
friend were already voiced by the pioneers of mathematical epidemiology. A
century later, we can draw on cross-fertilizations from network and game theory
and the emerging field of eco-evolutionary dynamics to substantiate them.
Epidemiological spreading is thereby revealed as a fundamentally heterogeneous
and erratic process that shares certain properties with more unwieldy
phenomena, such as earthquakes, hurricanes, traffic jams, and stock crashes.
They are all characterised by high tail risks that materialize very rarely but
fatally. That they arise from bursts of unlikely chains of localized random
"superspreader" events, by which micro-scale fluctuations and uncertainties may
get heftily magnified, makes their accurate prediction and control
intrinsically and notoriously hard. That epidemic disease spreading is moreover
closely intertwined with equally heterogeneous genetic drift and information
feedback adds new challenges -- and chances.


  Spatial separation is often included in models of ethnic divergence but it
has also been realised that urban subcultures can, and frequently do, emerge in
sympatry. Previous research tended to attribute this phenomenon to the human
tendency to imitate self-similar individuals and actively differentiate oneself
from individuals recognized as members of an outgroup. Application of such a
model to non-human animals has been, however, viewed as problematic. We present
a parsimonious model of subculture emergence where the algorithm of social
learning does not require the assumption of an 'imitation threshold'. All it
takes is a slight modification of Galton-Pearson's biometric model previously
used to approximate cultural inheritance. The new model includes
proportionality between the variance of inputs (cultural 'parents') and the
variance of outputs (cultural 'offspring'). In this model, assortment alone can
lead to the formation of distinct cohesive clusters of individuals
(subcultures) with a low within-group and large between-group variability even
in absence of a spatial separation or disruptive natural selection. Sympatric
emergence of arbitrary behavioural varieties preceding ecological divergence
may thus represent the norm, not the exception, in all cultural animals.


  Feature embedding methods have been proposed in literature to represent
sequences as numeric vectors to be used in some bioinformatics investigations,
such as family classification and protein structure prediction. Recent
theoretical results showed that the well-known Lyndon factorization preserves
common factors in overlapping strings. Surprisingly, the fingerprint of a
sequencing read, which is the sequence of lengths of consecutive factors in
variants of the Lyndon factorization of the read, is effective in preserving
sequence similarities, suggesting it as basis for the definition of novels
representations of sequencing reads. We propose a novel feature embedding
method for Next-Generation Sequencing (NGS) data using the notion of
fingerprint. We provide a theoretical and experimental framework to estimate
the behaviour of fingerprints and of the $k$-mers extracted from it, called
$k$-fingers, as possible feature embeddings for sequencing reads. As a case
study to assess the effectiveness of such embeddings, we use fingerprints to
represent RNA-Seq reads and to assign them to the most likely gene from which
they were originated as fragments of transcripts of the gene. We provide an
implementation of the proposed method in the tool lyn2vec, which produces
Lyndon-based feature embeddings of sequencing reads.


  Given multiple new COVID-19 variants are continuously emerging,
non-pharmaceutical interventions are still primary control strategies to curb
the further spread of coronavirus. However, implementing strict interventions
over extended periods of time is inevitably hurting the economy. With an aim to
solve this multi-objective decision-making problem, we investigate the
underlying associations between policies, mobility patterns, and virus
transmission. We further evaluate the relative performance of existing COVID-19
control measures and explore potential optimal strategies that can strike the
right balance between public health and socio-economic recovery for individual
states in the US. The results highlight the power of state of emergency
declaration and wearing face masks and emphasize the necessity of pursuing
tailor-made strategies for different states and phases of epidemiological
transmission. Our framework enables policymakers to create more refined designs
of COVID-19 strategies and can be extended to inform policy makers of any
country about best practices in pandemic response.


  Hexapeptides are widely applied as a model system for studying
amyloid-forming properties of polypeptides, including proteins. Recently, large
experimental databases have become publicly available with amyloidogenic
labels. Using these datasets for training and testing purposes, one may build
artificial intelligence (AI)-based classifiers for predicting the amyloid state
of peptides. In our previous work (Biomolecules, 11(4) 500, (2021)) we
described the Support Vector Machine (SVM)-based Budapest Amyloid Predictor
(\url{https://pitgroup.org/bap}). Here we apply the Budapest Amyloid Predictor
for discovering numerous amyloidogenic and non-amyloidogenic hexapeptide
patterns with accuracy between 80\% and 84\%, as surprising and succinct novel
rules for further understanding the amyloid state of peptides. For example, we
have shown that for any independently mutated residue (position marked by
``x''), the patterns CxFLWx, FxFLFx, or xxIVIV are predicted to be
amyloidogenic, while those of PxDxxx, xxKxEx, and xxPQxx non-amyloidogenic at
all. We note that each amyloidogenic pattern with two x's (e.g.,CxFLWx)
describes succinctly $20^2=400$ hexapeptides, while the non-amyloidogenic
patterns comprising four point mutations (e.g.,PxDxxx) gives $20^4=160,000$
hexapeptides in total. To our knowledge, no similar applications of artificial
intelligence tools or succinct amyloid patterns were described before the
present work.


  Disease outbreaks affect many ecosystems threatening species that also fight
against other natural enemies. We investigate a cyclic game system with $5$
species, whose organisms outcompete according to the rules of a generalised
spatial rock-paper-scissors game, during an epidemic. We study the effects of
behavioural movement strategies that allow individuals of one out of the
species to move towards areas with a low density of disease vectors and a high
concentration of enemies of their enemies. We perform a series of stochastic
simulations to discover the impact of self-preservation strategies in pattern
formation, calculating the species' spatial autocorrelation functions.
Considering organisms with different physical and cognitive abilities, we
compute the benefits of each movement tactic to reduce selection and infection
risks. Our findings show that the maximum profit in terms of territorial
dominance in the cyclic game is achieved if both survival movement strategies
are combined, with individuals prioritising social distancing. In the case of
an epidemic causing symptomatic illness, the drop in infection risk when
organisms identify and avoid disease vectors does not render a rise in the
species population because many refuges are disregarded, limiting the benefits
of safeguarding against natural enemies. Our results may be helpful to the
understanding of the behavioural strategies in ecosystems where organisms adapt
to face living conditions changes.


  An Allee effect occurs when the per-capita growth rate increases at low
densities. Here, we investigate the evolutionary stability of a partial
migration population with migrant population experiencing Allee effects.
Partial migration is a unique form of phenotypic diversity wherein migrant and
non-migrant individuals coexist together. It is shown that when Allee effect is
incorporated, the population undergoes a bifurcation as the fraction of
migrating population increases from zero to unity. Using an evolutionary game
theoretic approach, we prove the existence of a unique evolutionary stable
strategy (ESS). It is also shown that the ESS is the only ideal free
distribution (IFD) that arises in the context of partially migrating
population.


  In this letter, we deal with evolutionary game theoretic learning processes
for population games on networks with dynamically evolving communities.
Specifically, we propose a novel mathematical framework in which a
deterministic, continuous-time replicator equation on a community network is
coupled with a closed dynamic flow process between communities that is governed
by an environmental feedback mechanism, resulting in co-evolutionary dynamics.
Through a rigorous analysis of the system of differential equations obtained,
we characterize the equilibria of the coupled dynamical system. Moreover, for a
class of population games with two actions and symmetric rewards a Lyapunov
argument is employed to establish an evolutionary folk theorem that guarantees
convergence to the evolutionary stable states of the game. Numerical
simulations are provided to illustrate and corroborate our findings.


  Reputation-based cooperation on social networks offers a causal mechanism
between graph properties and social trust. Recent papers on the `structural
microfoundations` of the society used this insight to show how demographic
processes, such as falling fertility, urbanisation, and migration, can alter
the logic of human societies. This paper demonstrates the underlying mechanism
in a way that is accessible to scientists not specialising in networks.
Additionally, the paper shows that, when the size and degree of the network is
fixed (i.e., all graphs have the same number of agents, who all have the same
number of connections), it is the clustering coefficient that drives
differences in how cooperative social networks are.


  In this paper, we propose a parameter identification methodology of the SIRD
model, an extension of the classical SIR model, that considers the deceased as
a separate category. In addition, our model includes one parameter which is the
ratio between the real total number of infected and the number of infected that
were documented in the official statistics.
  Due to many factors, like governmental decisions, several variants
circulating, opening and closing of schools, the typical assumption that the
parameters of the model stay constant for long periods of time is not
realistic. Thus our objective is to create a method which works for short
periods of time. In this scope, we approach the estimation relying on the
previous 7 days of data and then use the identified parameters to make
predictions.
  To perform the estimation of the parameters we propose the average of an
ensemble of neural networks. Each neural network is constructed based on a
database built by solving the SIRD for 7 days, with random parameters. In this
way, the networks learn the parameters from the solution of the SIRD model.
  Lastly we use the ensemble to get estimates of the parameters from the real
data of Covid19 in Romania and then we illustrate the predictions for different
periods of time, from 10 up to 45 days, for the number of deaths. The main goal
was to apply this approach on the analysis of COVID-19 evolution in Romania,
but this was also exemplified on other countries like Hungary, Czech Republic
and Poland with similar results.
  The results are backed by a theorem which guarantees that we can recover the
parameters of the model from the reported data. We believe this methodology can
be used as a general tool for dealing with short term predictions of infectious
diseases or in other compartmental models.


  Molecular design based on generative models, such as variational autoencoders
(VAEs), has become increasingly popular in recent years due to its efficiency
for exploring high-dimensional molecular space to identify molecules with
desired properties. While the efficacy of the initial model strongly depends on
the training data, the sampling efficiency of the model for suggesting novel
molecules with enhanced properties can be further enhanced via latent space
optimization. In this paper, we propose a multi-objective latent space
optimization (LSO) method that can significantly enhance the performance of
generative molecular design (GMD). The proposed method adopts an iterative
weighted retraining approach, where the respective weights of the molecules in
the training data are determined by their Pareto efficiency. We demonstrate
that our multi-objective GMD LSO method can significantly improve the
performance of GMD for jointly optimizing multiple molecular properties.


  Contacts between people are the absolute drivers of contagious respiratory
infections. For this reason, limiting and tracking contacts is a key strategy
for the control of the COVID-19 epidemic. Digital contact tracing has been
proposed as an automated solution to scale up traditional contact tracing.
However, the required penetration of contact tracing apps within a population
to achieve a desired target in the control of the epidemic is currently under
discussion within the research community. In order to understand the effects of
digital contact tracing, several mathematical models have been proposed. In
this article, we survey the main ones and we propose a compartmental SEIR model
with which it is possible, differently from the models in the related
literature, to derive closed-form conditions regarding the control of the
epidemic as a function of the contact tracing apps penetration and the testing
efficiency. Closed-form conditions are crucial for the understandability of
models, and thus for decision makers (including digital contact tracing
designers) to correctly assess the dependencies within the epidemic. With our
model, we find that digital contact tracing alone can rarely tame an epidemic:
for unrestrained COVID-19, this would require a testing turnaround of around 1
day and app uptake above 80% of the population, which are very difficult to
achieve in practice. However, digital contact tracing can still be effective if
complemented with other mitigation strategies, such as social distancing and
mask-wearing.


  In conventional molecular communication (MC) systems, the signaling molecules
used for information transmission are stored, released, and then replenished by
a transmitter (TX). However, the replenishment of signaling molecules at the TX
is challenging in practice. Furthermore, in most envisioned MC applications,
e.g., in the medical field, it is not desirable to insert the TX into the MC
system, as this might impair natural biological processes. In this paper, we
propose the concept of media modulation based MC where the TX is placed outside
the channel and utilizes signaling molecules already present inside the system.
The signaling molecules can assume different states which can be switched by
external stimuli. Hence, in media modulation based MC, the TX modulates
information into the state of the signaling molecules. In particular, we
exploit the group of photochromic molecules, which undergo light-induced
reversible state transitions, for media modulation. We study the usage of these
molecules for information transmission in a three-dimensional duct system,
which contains an eraser, a TX, and a receiver for erasing, writing, and
reading of information via external light, respectively. We develop a
statistical model for the received signal which accounts for the distribution
of the signaling molecules in the system, the initial states of the signaling
molecules, the reliability of the state control mechanism, the randomness of
irrepressible, spontaneous state switching, and the randomness of molecule
propagation. We adopt a maximum likelihood detector and a threshold based
detector. Furthermore, we derive analytical expressions for the optimal
threshold value and the resulting bit error rate (BER), respectively. Our
results reveal that media modulation enables reliable information transmission,
validating it as a promising alternative to MC based on molecule emitting TXs.


  There are two phases of Wnt signalling in early vertebrate embryogenesis:
very early, maternal Wnt signalling promotes dorsal development, and slightly
later, zygotic Wnt signalling promotes ventral and lateral mesoderm induction.
However, recent molecular biology analysis has revealed more complexity among
the direct Wnt target genes, with at least five classes. Here in order to test
the logic and the dynamics of a new Gene Regulatory Network model suggested by
these discoveries we use mathematical modelling based on ordinary differential
equations (ODEs). Our mathematical modelling of this Gene Regulatory Network
reveals that a simplified model, with one "super-gene" for each class is
sufficient to a large extent to describe the regulatory behaviour previously
observed experimentally.


  Anti-cancer drug discoveries have been serendipitous, we sought to present
the Open Molecular Graph Learning Benchmark, named CandidateDrug4Cancer, a
challenging and realistic benchmark dataset to facilitate scalable, robust, and
reproducible graph machine learning research for anti-cancer drug discovery.
CandidateDrug4Cancer dataset encompasses multiple most-mentioned 29 targets for
cancer, covering 54869 cancer-related drug molecules which are ranged from
pre-clinical, clinical and FDA-approved. Besides building the datasets, we also
perform benchmark experiments with effective Drug Target Interaction (DTI)
prediction baselines using descriptors and expressive graph neural networks.
Experimental results suggest that CandidateDrug4Cancer presents significant
challenges for learning molecular graphs and targets in practical application,
indicating opportunities for future researches on developing candidate drugs
for treating cancers.


  We investigate a quantitative network of gene expression dynamics describing
the competence development in Bacillus subtilis. First, we introduce an
Onsager-Machlup approach to quantify the most probable transition pathway for
both excitable and bistable dynamics. Then, we apply a machine learning method
to calculate the most probable transition pathway via the Euler-Lagrangian
equation. Finally, we analyze how the noise intensity affects the transition
phenomena.


  Networks exist everywhere in nature from the physical, chemical, biological
or social worlds to the designed spheres. To explore, if there exists some
higher-order organization that can be exploited to distinguish different types
of networks, we study 4,738 protein interaction networks (PINs) belonging to 16
phyla encompassing all the three domains of life. Our method utilizes
positional information of a network's nodes by appropriately normalizing the
frequency of automorphic orbits appearing in the induced graphlets of sizes
2-5. There are some evolutionary constraints imposed on the network's topology
which shape its local architecture as well as its behavior. According to these
rules (features), each type of network occupies its respective position within
a common network space. A deep neural network was trained on differentially
expressed orbits resulting in a prediction accuracy of 85%. Our results
indicate that nature has, probably, allocated a specific band of design space
to various superfamilies of PINs.


  Phylogenetics is a classical methodology in computational biology that today
has become highly relevant for medical investigation of single-cell data, e.g.,
in the context of cancer development. The exponential size of the tree space
is, unfortunately, a substantial obstacle for Bayesian phylogenetic inference
using Markov chain Monte Carlo based methods since these rely on local
operations. And although more recent variational inference (VI) based methods
offer speed improvements, they rely on expensive auto-differentiation
operations for learning the variational parameters. We propose VaiPhy, a
remarkably fast VI based algorithm for approximate posterior inference in an
augmented tree space. VaiPhy produces marginal log-likelihood estimates on par
with the state-of-the-art methods on real data and is considerably faster since
it does not require auto-differentiation. Instead, VaiPhy combines coordinate
ascent update equations with two novel sampling schemes: (i) SLANTIS, a
proposal distribution for tree topologies in the augmented tree space, and (ii)
the JC sampler, to the best of our knowledge, the first-ever scheme for
sampling branch lengths directly from the popular Jukes-Cantor model. We
compare VaiPhy in terms of density estimation and runtime. Additionally, we
evaluate the reproducibility of the baselines. We provide our code on GitHub:
\url{https://github.com/Lagergren-Lab/VaiPhy}.


  Motivation: This paper presents libRoadRunner 2.0, an extensible,
high-performance, cross-platform, open-source software library for the
simulation and analysis of models expressed using Systems Biology Markup
Language SBML).
  Results: libRoadRunner is a self-contained library, able to run both as a
component inside other tools via its C++ and C bindings, and interactively
through its Python or Julia interface. libRoadRunner uses a custom Just-In-Time
JIT compiler built on the widely-used LLVM JIT compiler framework. It compiles
SBML-specified models directly into native machine code for a large variety of
processors, making it appropriate for solving extremely large models or
repeated runs. libRoadRunner is flexible, supporting the bulk of the SBML
specification (except for delay and nonlinear algebraic equations) and
including several SBML extensions such as composition and distributions. It
offers multiple deterministic and stochastic integrators, as well as tools for
steady-state, sensitivity, stability analysis, and structural analysis of the
stoichiometric matrix.
  Availability: libRoadRunner binary distributions are available for Mac OS X,
Linux, and Windows. The library is licensed under the Apache License Version
2.0. libRoadRunner is also available for ARM-based computers such as the
Raspberry Pi and can in principle be compiled on any system supported by
LLVM-13. http://sys-bio.github.io/roadrunner/index.html provides online
documentation, full build instructions, binaries, and a git source repository.


  Boolean networks are popular tools for the exploration of qualitative
dynamical properties of biological systems. Several dynamical interpretations
have been proposed based on the same logical structure that captures the
interactions between Boolean components. They reproduce, in different degrees,
the behaviours emerging in more quantitative models. In particular, regulatory
conflicts can prevent the standard asynchronous dynamics from reproducing some
trajectories that might be expected upon inspection of more detailed models. We
introduce and study the class of networks with linear cuts, where linear
components -- intermediates with a single regulator and a single target --
eliminate the aforementioned regulatory conflicts. The interaction graph of a
Boolean network admits a linear cut when a linear component occurs in each
cycle and in each path from components with multiple targets to components with
multiple regulators. Under this structural condition the attractors are in
one-to-one correspondence with the minimal trap spaces, and the reachability of
attractors can also be easily characterized. Linear cuts provide the base for a
new interpretation of the Boolean semantics that captures all behaviours of
multi-valued refinements with regulatory thresholds that are uniquely defined
for each interaction, and contribute a new approach for the investigation of
behaviour of logical models.


  The investigation of dynamical processes on networks has been one focus for
the study of contagion processes. It has been demonstrated that contagions can
be used to obtain information about the embedding of nodes in a Euclidean
space. Specifically, one can use the activation times of threshold contagions
to construct contagion maps as a manifold-learning approach. One drawback of
contagion maps is their high computational cost. Here, we demonstrate that a
truncation of the threshold contagions may considerably speed up the
construction of contagion maps. Finally, we show that contagion maps may be
used to find an insightful low-dimensional embedding for single-cell
RNA-sequencing data in the form of cell-similarity networks and so reveal
biological manifolds. Overall, our work makes the use of contagion maps as
manifold-learning approaches on empirical network data more viable.


  The severe acute respiratory syndrome of coronavirus 2 spread globally very
quickly, causing great concern at the international level due to the severity
of the associated respiratory disease, the so-called COVID-19. Considering Rio
de Janeiro city (Brazil) as an example, the first diagnosis of this disease
occurred in March 2020, but the exact moment when the local spread of the virus
started is uncertain as the Brazilian epidemiological surveillance system was
not widely prepared to detect suspected cases of COVID-19 at that time.
Improvements in this surveillance system occurred over the pandemic, but due to
the complex nature of the disease transmission process, specifying the exact
moment of emergence of new community contagion outbreaks is a complicated task.
This work aims to propose a general methodology to determine possible start
dates for the multiple community outbreaks of COVID-19, using for this purpose
a parametric statistical approach that combines surveillance data, nonlinear
regression, and information criteria to obtain a statistical model capable of
describing the multiple waves of contagion observed. The dynamics of COVID-19
in the city of Rio de Janeiro is taken as a case study, and the results suggest
that the original strain of the virus was already circulating in Rio de Janeiro
city as early as late February 2020, probably being massively disseminated in
the population during the carnival festivities.


  High-throughput structure-based screening of drug-like molecules has become a
common tool in biomedical research. Recently, acceleration with graphics
processing units (GPUs) has provided a large performance boost for molecular
docking programs. Both cloud and high-performance computing (HPC) resources
have been used for large screens with molecular docking programs; while NVIDIA
GPUs have dominated cloud and HPC resources, new vendors such as AMD and Intel
are now entering the field, creating the problem of software portability across
different GPUs. Ideally, software productivity could be maximized with portable
programming models that are able to maintain high performance across
architectures. While in many cases compiler directives have been used as an
easy way to offload parallel regions of a CPU-based program to a GPU
accelerator, they may also be an attractive programming model for providing
portability across different GPU vendors, in which case the porting process may
proceed in the reverse direction: from low-level, architecture-specific code to
higher-level directive-based abstractions. MiniMDock is a new mini-application
(miniapp) designed to capture the essential computational kernels found in
molecular docking calculations, such as are used in pharmaceutical drug
discovery efforts, in order to test different solutions for porting across GPU
architectures. Here we extend MiniMDock to GPU offloading with OpenMP
directives, and compare to performance of kernels using CUDA, and HIP on both
NVIDIA and AMD GPUs, as well as across different compilers, exploring
performance bottlenecks. We document this reverse-porting process, from highly
optimized device code to a higher-level version using directives, compare code
structure, and describe barriers that were overcome in this effort.


  Xenon can produce general anesthesia. Its main protein target is the
N-methyl-D-aspartate receptor, a ionotropic channel playing a pivotal role in
the function of the central nervous system.
  The molecular mechanisms allowing this noble gas to have such a specific
effect remain obscure, probably as a consequence of the lack of structural data
at the atomic level of detail. Herein, as a result of five independent
molecular dynamics simulations, three different binding sites were found for
xenon in the glycine binding domain of the N-methyl-D-aspartate receptor. The
absolute binding free energy of xenon in these sites ranges between -8 and -14
kJ/mole. However, it depends significantly upon the protein conformer chosen
for performing the calculation, suggesting that larger values could probably be
obtained, if other conformers were considered.
  These three sites are next to each other, one of them being next to the
glycine site. This could explain why the F758W and F758Y mutations can prevent
competitive inhibition by xenon without affecting glycine binding.


  Structural templates are 3D signatures representing protein functional sites,
such as ligand binding cavities, metal coordination motifs or catalytic sites.
Here we explore methods to generate template libraries and algorithms to query
structures for conserved 3D motifs. Applications of templates are discussed, as
well as some exemplar cases for examining evolutionary links in enzymes. We
also introduce the concept of using more than one template per structure to
represent flexible sites, as an approach to better understand catalysis through
snapshots captured in enzyme structures. Functional annotation from structure
is an important topic that has recently resurfaced due to the new more accurate
methods of protein structure prediction. Therefore, we anticipate that
template-based functional site detection will be a powerful tool in the task of
characterising a vast number of new protein models.


  For the last few years there has been a resurgence in the use of
phenomenological growth models for predicting the early dynamics of infectious
diseases. These models assume that time is a continuous variable whereas in the
present contribution, the discrete versions of Gompertz and Generalized
Logistic models are used for early monitoring and short-term forecasting of the
spread of an epidemic in a region. The time-continuous models are represented
mathematically by first-order differential equations while their discrete
versions are represented by first-order difference equations that involve
parameters that should be estimated prior to forecasting. The methodology for
estimating such parameters is described in detail. Real data of COVID-19
infection in Cuba is used to illustrate this methodology. The proposed
methodology was implemented for the first thirty-five days, being able to
predict with very good precision the data reported for the following twenty
days. The codes implemented to study the Gompertz model in differences are
included in an appendix with each step of the methodology identified.


  The recent development of imaging and sequencing technologies enables
systematic advances in the clinical study of lung cancer. Meanwhile, the human
mind is limited in effectively handling and fully utilizing the accumulation of
such enormous amounts of data. Machine learning-based approaches play a
critical role in integrating and analyzing these large and complex datasets,
which have extensively characterized lung cancer through the use of different
perspectives from these accrued data. In this article, we provide an overview
of machine learning-based approaches that strengthen the varying aspects of
lung cancer diagnosis and therapy, including early detection, auxiliary
diagnosis, prognosis prediction and immunotherapy practice. Moreover, we
highlight the challenges and opportunities for future applications of machine
learning in lung cancer.


  Artificial ecosystems provide an additional experimental tool to support
laboratory work, field work, and theoretical development in competitive
exclusion research. A novel application of a spatiotemporal agent based model
is presented which simulates two foraging species of different intrinsic growth
rates competing for the same replenishing resource in stable and seasonal
environments. This experimental approach provides precise control over the
parameters for the environment, the species, and individual movements. Detailed
trajectories of these non-equilibrium populations and their characteristics are
produced. Narrow zones of potential coexistence are identified within the
environmental and intrinsic growth parameter spaces. An example of commensalism
driven by the local spatial dynamics is identified. Results of these
experiments are discussed in context of modern coexistence theory and research
in movement-mediated community assembly. Constraints on possible origination
scenarios are identified.


  We demonstrate an approach to replicate and forecast the spread of the
SARS-CoV-2 (COVID-19) pandemic using the toolkit of probabilistic programming
languages (PPLs). Our goal is to study the impact of various modeling
assumptions and motivate policy interventions enacted to limit the spread of
infectious diseases. Using existing compartmental models we show how to use
inference in PPLs to obtain posterior estimates for disease parameters. We
improve popular existing models to reflect practical considerations such as the
under-reporting of the true number of COVID-19 cases and motivate the need to
model policy interventions for real-world data. We design an SEI3RD model as a
reusable template and demonstrate its flexibility in comparison to other
models. We also provide a greedy algorithm that selects the optimal series of
policy interventions that are likely to control the infected population subject
to provided constraints. We work within a simple, modular, and reproducible
framework to enable immediate cross-domain access to the state-of-the-art in
probabilistic inference with emphasis on policy interventions. We are not
epidemiologists; the sole aim of this study is to serve as an exposition of
methods, not to directly infer the real-world impact of policy-making for
COVID-19.


  One of the most contested questions about human behaviour is whether there
are inherent sex or gender differences in the formation and maintenance of
social bonds. On one hand, female and male brains are structurally almost
identical, and while there are sex differences in the endocrine system, these
are small, while much of gendered identity and behaviour is learned. On the
other hand, sex differences in some aspects of social behaviour have deep
evolutionary roots, and are widely present in non-human animals. This
observational study recorded the frequency of same-aged, adult human groups
appearing in public spaces through 2636 hours, recording group formation by
1.2mn people via 170 research assistants in 46 countries across the world. The
results show (a) a significant sex-gender difference in same-sex-same-age
frequency, in that ~50% more female-female than male-male pairs are observed in
public spaces globally, and (b) that despite regional variation, the patterns
holds up in every global region. This is the first study of sex-gender
difference in dyadic social behaviour across the world on this scale, and the
first global study that uses direct rather than internet-based observations.


  Most natural and engineered information-processing systems transmit
information via signals that vary in time. Computing the information
transmission rate or the information encoded in the temporal characteristics of
these signals, requires the mutual information between the input and output
signals as a function of time, i.e. between the input and output trajectories.
Yet, this is notoriously difficult because of the high-dimensional nature of
the trajectory space, and all existing techniques require approximations. We
present an exact Monte Carlo technique called Path Weight Sampling (PWS) that,
for the first time, makes it possible to compute the mutual information between
input and output trajectories for any stochastic system that is described by a
master equation. The principal idea is to use the master equation to evaluate
the exact conditional probability of an individual output trajectory for a
given input trajectory, and average this via Monte Carlo sampling in trajectory
space to obtain the mutual information. We present three variants of PWS, which
all generate the trajectories using the standard stochastic simulation
algorithm. While Direct PWS is a brute-force method, Rosenbluth-Rosenbluth PWS
exploits the analogy between signal trajectory sampling and polymer sampling,
and Thermodynamic Integration PWS is based on a reversible work calculation in
trajectory space. PWS also makes it possible to compute the mutual information
between input and output trajectories for systems with hidden internal states
as well as systems with feedback from output to input. Applying PWS to the
bacterial chemotaxis system, consisting of 182 coupled chemical reactions,
demonstrates not only that the scheme is highly efficient, but also that the
number of receptor clusters is much smaller than hitherto believed, while their
size is much larger.


  Despite the progress in medical data collection the actual burden of
SARS-CoV-2 remains unknown due to under-ascertainment of cases. This was
apparent in the acute phase of the pandemic and the use of reported deaths has
been pointed out as a more reliable source of information, likely less prone to
under-reporting. Since daily deaths occur from past infections weighted by
their probability of death, one may infer the total number of infections
accounting for their age distribution, using the data on reported deaths. We
adopt this framework and assume that the dynamics generating the total number
of infections can be described by a continuous time transmission model
expressed through a system of non-linear ordinary differential equations where
the transmission rate is modelled as a diffusion process allowing to reveal
both the effect of control strategies and the changes in individuals behavior.
We develop this flexible Bayesian tool in Stan and study 3 pairs of European
countries, estimating the time-varying reproduction number($R_t$) as well as
the true cumulative number of infected individuals. As we estimate the true
number of infections we offer a more accurate estimate of $R_t$. We also
provide an estimate of the daily reporting ratio and discuss the effects of
changes in mobility and testing on the inferred quantities.


  Despite their noted potential in polynomial-system solving, there are few
concrete examples of Newton-Okounkov bodies arising from applications.
Accordingly, in this paper, we introduce a new application of Newton-Okounkov
body theory to the study of chemical reaction networks, and compute several
examples. An important invariant of a chemical reaction network is its maximum
number of positive steady states, which is realized as the maximum number of
positive real roots of a parametrized polynomial system. Here, we introduce a
new upper bound on this number, namely the `Newton-Okounkov body bound' of a
chemical reaction network. Through explicit examples, we show that the
Newton-Okounkov body bound of a network gives a good upper bound on its maximum
number of positive steady states. We also compare this Newton-Okounkov body
bound to a related upper bound, namely the mixed volume of a chemical reaction
network, and find that it often achieves better bounds.