  Small-angle X-ray scattering (SAXS) experiments are increasingly used to
probe RNA structure. A number of \emph{forward models} that relate measured
SAXS intensities and structural features, and that are suitable to model either
explicit-solvent effects or solute dynamics, have been proposed in the past
years. Here we introduce an approach that integrates atomistic molecular
dynamics simulations and SAXS experiments to reconstruct RNA structural
ensembles while simultaneously accounting for both RNA conformational dynamics
and explicit-solvent effects. Our protocol exploits SAXS pure-solute forward
models and enhanced sampling methods to sample an heterogenous ensemble of
structures, with no information towards the experiments provided on-the-fly.
The generated structural ensemble is then reweighted through the maximum
entropy principle so as to match reference SAXS experimental data at multiple
ionic conditions. Importantly, accurate explicit-solvent forward models are
used at this reweighting stage. We apply this framework to the
GTPase-associated center, a relevant RNA molecule involved in protein
translation, in order to elucidate its ion-dependent conformational ensembles.
We show that (a) both solvent and dynamics are crucial to reproduce
experimental SAXS data and (b) the resulting dynamical ensembles contain an
ion-dependent fraction of extended structures.


  Rapid antigen tests detect proteins at the surface of virus particles,
identifying the disease during its infectious phase. In contrast, PCR tests
detect viral genomes; they can thus diagnose COVID-19 before the infectious
phase but also react to remnants of the virus genome, even weeks after live
virus ceases to be detectable in the respiratory tract. Furthermore, the
logistics for administering the tests are different, with rapid antigen tests
being much easier to administer at-scale. In this article, we discuss the
relative advantages of the different testing procedures and summarise evidence
that shows that using antigen tests 2-3 times per week could become a powerful
tool to suppress the COVID-19 pandemic. We also discuss the results of recent
large-scale rapid antigen testing in Austrian schools. While our report on
testing predates Delta, we have updated the review with recent data on viral
loads in breakthrough infections and more information about testing efficacy,
especially in children.


  This work deals with the influence of the neighborhood in simple
rock-paper-scissors models of biodiversity. We consider the case of three
distinct species which evolve under the standard rules of mobility,
reproduction and competition. The rule of competition follows the guidance of
the rock-paper-scissors game, with the prey being annihilated, leaving an empty
site in accordance with the May-Leonard proposal for the predator and prey
competition. We use the von Neumann neighborhood, but we consider mobility
under the presence of the first, second and third neighbors in three distinct
environments, one with equal probability and the others with probability
following the power law and exponential profiles. The results are different,
but they all show that increasing the neighbourhood increases the
characteristic length of the system in an important way. We have studied other
possibilities, in particular the case where one modifies the manner a specific
species competes, unveiling the interesting result in which the strongest
individuals may constitute the less abundant population.


  Several organelles in eukaryotic cells, including mitochondria and the
endoplasmic reticulum, form interconnected tubule networks extending throughout
the cell. These tubular networks host many biochemical pathways that rely on
proteins diffusively searching through the network to encounter binding
partners or localized target regions. In this work we develop both exact
analytical methods to compute mean first passage times and efficient kinetic
Monte Carlo algorithms to simulate trajectories of particles diffusing in a
tubular network. Our approach leverages exact propagator functions for the
distribution of transition times between network nodes and allows large
simulation time steps determined by the network structure. The methodology is
applied to both synthetic planar networks and organelle network structures,
demonstrating key general features such as the heterogeneity of search times in
different network regions and the functional advantage of broadly distributing
target sites throughout the network. The proposed algorithms pave the way for
future exploration of the interrelationship between tubular network structure
and biomolecular reaction kinetics.


  Backtracking of RNA polymerase (RNAP) is an important pausing mechanism
during DNA transcription that is part of the error correction process that
enhances transcription fidelity. We model the backtracking mechanism of RNA
polymerase, which usually happens when the polymerase tries to incorporate a
mismatched nucleotide triphosphate. Previous models have made simplifying
assumptions such as neglecting the trailing polymerase behind the backtracking
polymerase or assuming that the trailing polymerase is stationary. We derive
exact analytic solutions of a stochastic model that includes locally
interacting RNAPs by explicitly showing how a trailing RNAP influences the
probability that an error is corrected or incorporated by the leading
backtracking RNAP. We also provide two related methods for computing the mean
times to error correction or incorporation given an initial local RNAP
configuration.


  Synaptic communication is based on a biological Molecular Communication (MC)
system which may serve as a blueprint for the design of synthetic MC systems.
However, the physical modeling of synaptic MC is complicated by the possible
saturation of the molecular receiver caused by the competition of
neurotransmitters (NTs) for postsynaptic receptors. Receiver saturation renders
the system behavior nonlinear in the number of released NTs and is commonly
neglected in existing analytical models. Furthermore, due to the ligands'
competition for receptors (and vice versa), the individual binding events at
the molecular receiver are in general statistically dependent and the binomial
model for the statistics of the received signal does not apply. In this work,
we propose a novel deterministic model for receptor saturation in terms of a
state-space description based on an eigenfunction expansion of Fick's diffusion
equation. The presented solution is numerically stable and computationally
efficient. Employing the proposed deterministic model, we show that saturation
at the molecular receiver reduces the peak-value of the expected received
signal and accelerates the clearance of NTs as compared to the case when
receptor occupancy is neglected. We further derive a statistical model for the
received signal in terms of the hypergeometric distribution which accounts for
the competition of NTs for receptors and the competition of receptors for NTs.
The proposed statistical model reveals how the signal statistics are shaped by
the number of released NTs, the number of receptors, and the binding kinetics
of the receptors, respectively, in the presence of competition. We show that
the impact of these parameters on the signal variance depends on the relative
numbers of NTs and receptors. The accuracy of the proposed deterministic and
statistical models is verified by particle-based computer simulations.


  Effective understanding of a disease such as cancer requires fusing multiple
sources of information captured across physical scales by multimodal data. In
this work, we propose a novel feature embedding module that derives from
canonical correlation analyses to account for intra-modality and inter-modality
correlations. Experiments on simulated and real data demonstrate how our
proposed module can learn well-correlated multi-dimensional embeddings. These
embeddings perform competitively on one-year survival classification of
TCGA-BRCA breast cancer patients, yielding average F1 scores up to 58.69% under
5-fold cross-validation.


  Natural populations often show enhanced genetic drift consistent with a
strong skew in their offspring number distribution. The skew arises because the
variability of family sizes is either inherently strong or amplified by
population expansions, leading to so-called `jackpot' events. The resulting
allele frequency fluctuations are large and, therefore, challenge standard
models of population genetics, which assume sufficiently narrow offspring
distributions. While the neutral dynamics backward in time can be readily
analyzed using coalescent approaches, we still know little about the effect of
broad offspring distributions on the dynamics forward in time, especially with
selection. Here, we employ an exact asymptotic analysis combined with a scaling
hypothesis to demonstrate that over-dispersed frequency trajectories emerge
from the competition of conventional forces, such as selection or mutations,
with an emerging time-dependent sampling bias against the minor allele. The
sampling bias arises from the characteristic time-dependence of the largest
sampled family size within each allelic type. Using this insight, we establish
simple scaling relations for allele frequency fluctuations, fixation
probabilities, extinction times, and the site frequency spectra that arise when
offspring numbers are distributed according to a power law $~n^{-(1+\alpha)}$.
To demonstrate that this coarse-grained model captures a wide variety of
non-equilibrium dynamics, we validate our results in traveling waves, where the
phenomenon of 'gene surfing' can produce any exponent $1<\alpha <2$. We argue
that the concept of a dynamic sampling bias is useful generally to develop both
intuition and statistical tests for the unusual dynamics of populations with
skewed offspring distributions, which can confound commonly used tests for
selection or demographic history.


  We develop a theory for the susceptible-infected-susceptible (SIS) epidemic
model on networks that incorporate both network structure and dynamic
correlations. This theory can account for the multistage onset of the epidemic
phase in scale-free networks. This phenomenon is characterized by multiple
peaks in the susceptibility as a function of the infection rate. It can be
explained by that, even under the global epidemic threshold, a hub can sustain
the epidemics for an extended period. Moreover, our approach improves
theoretical calculations of prevalence close to the threshold in heterogeneous
networks and also can predict the average risk of infection for neighbors of
nodes with different degree and state on uncorrelated static networks.


  In this work, using a detailed dataset furnished by National Health
Authorities concerning the Province of Pavia (Lombardy, Italy), we propose to
determine the essential features of the ongoing COVID-19 pandemic in term of
contact dynamics. Our contribution is devoted to provide a possible planning of
the needs of medical infrastructures in the Pavia Province and to suggest
different scenarios about the vaccination campaign which possibly help in
reducing the fatalities and/or reducing the number of infected in the
population. The proposed research combines a new mathematical description of
the spread of an infectious diseases which takes into account both age and
average daily social contacts with a detailed analysis of the dataset of all
traced infected individuals in the Province of Pavia. These information are
used to develop a data-driven model in which calibration and feeding of the
model are extensively used. The epidemiological evolution is obtained by
relying on an approach based on statistical mechanics. This leads to study the
evolution over time of a system of probability distributions characterizing the
age and social contacts of the population. One of the main outcomes shows that,
as expected, the spread of the disease is closely related to the mean number of
contacts of individuals. The model permits to forecast thanks to an uncertainty
quantification approach and in the short time horizon, the average number and
the confidence bands of expected hospitalized classified by age and to test
different options for an effective vaccination campaign with age-decreasing
priority.


  The development and authorization of COVID-19 vaccines has provided the
clearest path forward to eliminate community spread hence end the ongoing
SARS-CoV-2 pandemic. However, the limited pace at which the vaccine can be
administered motivates the question, to what extent must we continue to adhere
to social intervention measures such as mask wearing and social distancing? To
address this question, we develop a mathematical model of COVID-19 spread
incorporating both vaccine dynamics and socio-epidemiological parameters. We
use this model to study two important measures of disease control and
eradication, the effective reproductive number $R_t$ and the peak intensive
care unit (ICU) caseload, over three key parameters: social measure adherence,
vaccination rate, and vaccination coverage. Our results suggest that, due to
the slow pace of vaccine administration, social measures must be maintained by
a large proportion of the population until a sufficient proportion of the
population becomes vaccinated for the pandemic to be eradicated. By contrast,
with reduced adherence to social measures, hospital ICU cases will greatly
exceed capacity, resulting in increased avoidable loss of life. These findings
highlight the complex interplays involved between vaccination and social
protective measures, and indicate the practical importance of continuing with
extent social measures while vaccines are scaled up to allow the development of
the herd immunity needed to end or control SARS-CoV-2 sustainably.


  Graph network science is becoming increasingly popular, notably in big-data
perspective where understanding individual entities for individual functional
roles is complex and time consuming. It is likely when a set of genes are
regulated by a set of genetic variants, the genes set is recruited for a common
or related functional purpose. Grouping and extracting communities from network
of associations becomes critical to understand system complexity, thus
prioritizing genes for dis-ease and functional associations. Workload is
reduced when studying entities one at a time. For this, we present GraphBreak,
a suite of tools for community detection application, such as for gene
co-expression, protein interaction, regulation network, etc.Although developed
for use case of eQTLs regulatory genomic net-work community study -- results
shown with our analysis with sample eQTL data. Graphbreak can be deployed for
other studies if input data has been fed in requisite format, including but not
limited to gene co-expression networks, protein-protein interaction network,
signaling pathway and metabolic network. Graph-Break showed critical use case
value in its downstream analysis for disease association of communities
detected. If all independent steps of community detection and analysis are a
step-by-step sub-part of the algorithm, GraphBreak can be considered a new
algorithm for community based functional characterization. Combination of
various algorithmic implementation modules into a single script for this
purpose illustrates GraphBreak novelty. Compared to other similar tools, with
GraphBreak we can better detect communities with over-representation of its
member genes for statistical association with diseases, therefore target genes
which can be prioritized for drug-positioning or drug-re-positioning as the
case be.


  We introduce a Bayesian sequential data assimilation method for COVID-19
forecasting. It is assumed that suitable transmission, epidemic and observation
models are available and previously validated and the transmission and epidemic
models are coded into a dynamical system. The observation model depends on the
dynamical system state variables and parameters, and is cast as a likelihood
function. We elicit prior distributions of the effective population size, the
dynamical system initial conditions and infectious contact rate, and use Markov
Chain Monte Carlo sampling to make inference and prediction of quantities of
interest (QoI) at the onset of the epidemic outbreak. The forecast is
sequentially updated over a sliding window of epidemic records as new data
becomes available. Prior distributions for the state variables at the new
forecasting time are assembled using the dynamical system, calibrated for the
previous forecast. Moreover, changes in the contact rate and effective
population size are naturally introduced through auto-regressive models on the
corresponding parameters. We show our forecasting method's performance using a
SEIR type model and COVID-19 data from several Mexican localities.


  Mass vaccination offers a promising exit strategy for the COVID-19 pandemic.
However, as vaccination progresses, demands to lift restrictions increase,
despite most of the population remaining susceptible. Using our age-stratified
SEIRD-ICU compartmental model and curated epidemiological and vaccination data,
we quantified the rate (relative to vaccination progress) at which countries
can lift non-pharmaceutical interventions without overwhelming their healthcare
systems. We analyzed scenarios ranging from immediately lifting restrictions
(accepting high mortality and morbidity) to reducing case numbers to a level
where test-trace-and-isolate (TTI) programs efficiently compensate for local
spreading events. In general, the age-dependent vaccination roll-out implies a
transient decrease of more than ten years in the average age of ICU patients
and deceased. The pace of vaccination determines the speed of lifting
restrictions; Taking the European Union (EU) as an example case, all considered
scenarios allow for steadily increasing contacts starting in May 2021 and
relaxing most restrictions by autumn 2021. Throughout summer 2021, only mild
contact restrictions will remain necessary. However, only high vaccine uptake
can prevent further severe waves. Across EU countries, seroprevalence impacts
the long-term success of vaccination campaigns more strongly than age
demographics. In addition, we highlight the need for preventive measures to
reduce contagion in school settings throughout the year 2021, where children
might be drivers of contagion because of them remaining susceptible...


  The wide spread of coronavirus disease 2019 (COVID-19) has declared a global
health emergency. As one of the most important targets for antibody and drug
developments, Spike RBD-ACE2 interface has received extensive attention. Here,
using molecular dynamics simulations, we explicitly evaluated the binding
energetic features of the RBD-ACE2 complex of both SARS-CoV and SARS-CoV-2 to
find the key residues. Although the overall ACE2-binding mode of the SARS-CoV-2
RBD is nearly identical to that of the SARS-CoV RBD, the difference in binding
affinity is as large as -16.35 kcal/mol. Energy decomposition analyses
identified three binding patches in the SARS-CoV-2 RBD and eleven key residues
(Phe486, Tyr505, Asn501, Tyr489, Gln493, Leu455 and etc) which are believed to
be the main targets for drug development. The dominating forces are from van
der Waals attractions and dehydration of these residues. It is also worth
mention that we found seven mutational sites (Lys417, Leu455, Ala475, Gly476,
Glu484, Gln498 and Val503) on SARS-CoV-2 which unexpectedly weakened the
RBD-ACE2 binding. Very interestingly, the most repulsive residue at the
RBD-ACE2 interface (E484), is found to be mutated in the latest UK variant,
B1.1.7, cause complete virus neutralization escapes from highly neutralizing
COVID-19 convalescent plasma. Our present results indicate that at least from
the energetic point of view such E484 mutation may have beneficial effects on
ACE2 binding. The present study provides a systematical understanding, from the
energetic point of view, of the binding features of SARS-CoV-2 RBD with ACE2
acceptor. We hope that the present findings of three binding patches, key
attracting residues and unexpected mutational sites can provide insights to the
design of SARS-CoV-2 drugs and identification of cross-active antibodies.


  Best match graphs (BMGs) are vertex-colored digraphs that naturally arise in
mathematical phylogenetics to formalize the notion of evolutionary closest
genes w.r.t. an a priori unknown phylogenetic tree. BMGs are explained by
unique least resolved trees. We prove that the property of a rooted,
leaf-colored tree to be least resolved for some BMG is preserved by the
contraction of inner edges. For the special case of two-colored BMGs, this
leads to a characterization of the least resolved trees (LRTs) of
binary-explainable trees and a simple, polynomial-time algorithm for the
minimum cardinality completion of the arc set of a BMG to reach a BMG that can
be explained by a binary tree.


  This study presents a new risk-averse multi-stage stochastic
epidemics-ventilator-logistics compartmental model to address the resource
allocation challenges of mitigating COVID-19. This epidemiological logistics
model involves the uncertainty of untested asymptomatic infections and
incorporates short-term human migration. Disease transmission is also
forecasted through a new formulation of transmission rates that evolve over
space and time with respect to various non-pharmaceutical interventions, such
as wearing masks, social distancing, and lockdown. The proposed multi-stage
stochastic model overviews different scenarios on the number of asymptomatic
individuals while optimizing the distribution of resources, such as
ventilators, to minimize the total expected number of newly infected and
deceased people. The Conditional Value at Risk (CVaR) is also incorporated into
the multi-stage mean-risk model to allow for a trade-off between the weighted
expected loss due to the outbreak and the expected risks associated with
experiencing disastrous pandemic scenarios. We apply our multi-stage mean-risk
epidemics-ventilator-logistics model to the case of controlling the COVID-19 in
highly-impacted counties of New York and New Jersey. We calibrate, validate,
and test our model using actual infection, population, and migration data. The
results indicate that short-term migration influences the transmission of the
disease significantly. The optimal number of ventilators allocated to each
region depends on various factors, including the number of initial infections,
disease transmission rates, initial ICU capacity, the population of a
geographical location, and the availability of ventilator supply. Our
data-driven modeling framework can be adapted to study the disease transmission
dynamics and logistics of other similar epidemics and pandemics.


  A novel approach to simulate simple protein-ligand systems at large time- and
length-scales is to couple Markov state models (MSMs) of molecular kinetics
with particle-based reaction-diffusion (RD) simulations, MSM/RD. Currently,
MSM/RD lacks a mathematical framework to derive coupling schemes; is limited to
isotropic ligands in a single conformational state, and is lacking a
multi-particle extensions. In this work, we address these needs by developing a
general MSM/RD framework by coarse-graining molecular dynamics into hybrid
switching diffusion processes. Given enough data to parametrize the model, it
is capable of modeling protein-protein interactions over large time- and
length-scales, and it can be extended to handle multiple molecules. We derive
the MSM/RD framework, and we implement and verify it for two protein-protein
benchmark systems and one multiparticle implementation to model the formation
of pentameric ring molecules. To enable reproducibility, we have published our
code in the MSM/RD software package.


  Background: To prevent future outbreaks of COVID-19, Australia is pursuing a
mass-vaccination approach in which a targeted group of the population
comprising healthcare workers, aged-care residents and other individuals at
increased risk of exposure will receive a highly effective priority vaccine.
The rest of the population will instead have access to a less effective
vaccine.
  Methods: We apply a large-scale agent-based model of COVID-19 in Australia to
investigate the possible implications of this hybrid approach to
mass-vaccination. The model is calibrated to recent epidemiological and
demographic data available in Australia, and accounts for several components of
vaccine efficacy.
  Findings: Within a feasible range of vaccine efficacy values, our model
supports the assertion that complete herd immunity due to vaccination is not
likely in the Australian context. For realistic scenarios in which herd
immunity is not achieved, we simulate the effects of mass-vaccination on
epidemic growth rate, and investigate the requirements of lockdown measures
applied to curb subsequent outbreaks. In our simulations, Australia's
vaccination strategy can feasibly reduce required lockdown intensity and
initial epidemic growth rate by 43% and 52%, respectively. The severity of
epidemics, as measured by the peak number of daily new cases, decreases by up
to two orders of magnitude under plausible mass-vaccination and lockdown
strategies.
  Interpretation: The study presents a strong argument for a large-scale
vaccination campaign in Australia, which would substantially reduce both the
intensity of future outbreaks and the stringency of non-pharmaceutical
interventions required for their suppression.


  Open-source biodiversity databases contain a large amount of species
occurrence records, but these are often spatially biased, which affects the
reliability of species distribution models based on these records. Sample bias
correction techniques include data filtering at the cost of record numbers or
require considerable additional sampling effort. However, independent data are
rarely available and assessment of the correction technique must rely on
performance metrics computed with subsets of the only available (biased) data,
which may be misleading. Here we assess the extent to which an acknowledged
sample bias correction technique is likely to improve models' ability to
predict species distributions in the absence of independent data. We assessed
the variation in model predictions induced by the correction and model
stochasticity. We present an index of the effect of correction relative to
model stochasticity, the Relative Overlap Index (ROI). We tested whether the
ROI better represented the effect of correction than classic performance
metrics and absolute overlap metrics using 64 vertebrate species and 21 virtual
species with a generated sample bias. When based on absolute overlaps and
cross-validation performance metrics, we found no effect of correction, except
for cAUC. When considering its effect relative to model stochasticity, the
effect of correction depended on the site and the species. Virtual species
enabled us to verify that the correction actually improved distribution
predictions and the biological relevance of the selected variables at the sites
with a clear gradient of sample bias, and when species distribution predictors
are not correlated with sample bias patterns.


  Best match graphs (BMGs) are a class of colored digraphs that naturally
appear in mathematical phylogenetics and can be approximated with the help of
similarity measures between gene sequences, albeit not without errors. The
corresponding graph editing problem can be used as a means of error correction.
Since the arc set modification problems for BMGs are NP-complete, efficient
heuristics are needed if BMGs are to be used for the practical analysis of
biological sequence data. Since BMGs have a characterization in terms of
consistency of a certain set of rooted triples, we consider heuristics that
operate on triple sets. As an alternative, we show that there is a close
connection to a set partitioning problem that leads to a class of top-down
recursive algorithms that are similar to Aho's supertree algorithm and give
rise to BMG editing algorithms that are consistent in the sense that they leave
BMGs invariant. Extensive benchmarking shows that community detection
algorithms for the partitioning steps perform best for BMG editing.


  This article presents an overview of future truly personal communications,
ranging from networking inside the human body to the exchange of data with
external wireless devices in the surrounding environment. At the nano- and
micro-scales, communications can be realized with the aid of molecular
mechanisms, Forster resonance energy transfer phenomenon, electromagnetic or
ultrasound waves. At a larger scale, in the domain of Body Area Networks, a
wide range of communication mechanisms is available, including smart-textiles,
inductive- and body-couplings, ultrasounds, optical and wireless radio
transmissions, a number of mature technologies existing already. The main goal
of this article is to identify the potential mechanisms that can be exploited
to provide interfaces in between nano- and micro-scale systems and Body Area
Networks. These interfaces have to bridge the existing gap between the two
worlds, in order to allow for truly personal communication systems to become a
reality. The extraordinary applications of such systems are also discussed, as
they are strong drivers of the research in this area.


  CoV2019 has evolved to be much more dangerous than CoV2003. Experiments
suggest that structural rearrangements dramatically enhance CoV2019 activity.
We identify a new first stage of infection that precedes structural
rearrangements by using biomolecular evolutionary theory to identify sequence
differences enhancing viral attachment rates. We find a small cluster of four
single mutations which show that CoV-2 has a new feature that promotes much
stronger viral attachment and enhances contagiousness. The extremely dangerous
dynamics of human coronavirus infection is a dramatic example of evolutionary
approach of self-organized networks to criticality. It may favor a very
successful vaccine. The identified mutations can be used to test the present
theory experimentally. The theory also works well for the newer strains and
explains their increased contagiousness.


  Proteins play crucial roles in every cellular process by interacting with
each other, with nucleic acids, metabolites, and other molecules. The resulting
assemblies can be very large and intricate and pose challenges to experimental
methods. In the current era of integrative modeling, it is often only by a
combination of various experimental techniques and computations that 3D models
of those molecular machines can be obtained. Among the various computational
approaches available, molecular docking is often the method of choice when it
comes to predicting 3D structures of complexes. Docking can generate
particularly accurate models when taking into account the available information
on the complex of interest. We review here the use of experimental and
bioinformatics data in protein-protein docking, describing recent software
developments and highlighting applications for the modeling of antibody-antigen
complexes and membrane protein complexes, and the use of evolutionary and shape
information.


  It has been a century since the species-area relationship (SAR) was first
proposed as a power law to explain how species richness scales with area. There
have been many attempts to explain the origin of this predominant form. Apart
from the power law, numerous empirical studies also report a semi-log form of
the SAR, but very few have addressed its incidence. In this work, we test
whether these relationships could emerge from the assembly of large random
communities on island-like systems. We reformulate the generalized
Lotka-Volterra model by introducing an area parameter that determines the
species richness of the assembled communities. Our analysis demonstrates that
the two most widely reported relationship forms can emerge due to differences
in immigration rates and skewness towards weak interactions. We particularly
highlight the incidence of the semi-log SAR for low immigration rates from a
source pool, which is consistent with several previous empirical studies. The
two SAR forms might show good fits to data over a large span of areas but a
power-law overestimates species richness on smaller islands in remote
archipelagoes.


  Heparin has been found to have antiviral activity against SARS-CoV-2. Here,
by means of sliding window docking, molecular dynamics simulations and
biochemical assays, we investigate the binding mode of heparin to the virus
spike glycoprotein and the molecular basis for its antiviral activity. The
simulations show that heparin binds at long, mostly positively charged patches
on the spike, thereby masking the basic residues of the receptor binding domain
and of the S1/S2 site. Experiments corroborated the simulation results by
showing that heparin inhibits the cleavage of spike by furin by binding to the
basic S1/S2 site. Our results indicate that heparin exerts its antiviral
activity by both direct and allosteric mechanisms. Furthermore, the simulations
provide insights into how heparan sulfate proteoglycans on the host cell can
facilitate viral infection. Our results will aid the rational optimization of
heparin derivatives for SARS-CoV-2 antiviral therapy.


  In a geographically distributed population, assortative clustering plays an
important role in evolution by modifying local environments. To examine its
effects in a linear habitat, we consider a one-dimensional grid of cells, where
each cell is either empty or occupied by an organism whose replication strategy
is genetically inherited to offspring. The strategy determines whether to have
offspring in surrounding cells, as a function of the neighborhood
configuration. If more than one offspring compete for a cell, then they can be
all exterminated due to the cost of conflict depending on environmental
conditions. We find that the system is more densely populated in an unfavorable
environment than in a favorable one because only the latter has to pay the cost
of conflict. This observation agrees reasonably well with a mean-field analysis
which takes assortative clustering of strategies into consideration. Our
finding suggests a possibility of intrinsic nonlinearity between environmental
conditions and population density when an evolutionary process is involved.


  We study dynamics of a variant of the SIR the model, where we assume that
individuals respond differently to dynamics of the epidemic. Their
heterogeneous response is modeled by the Preisach hysteresis operator. The
degree of heterogeneity of the response is measured by the variance $\sigma^2$
of the corresponding distribution (the Preisach density function). The proposed
model has a continuum of endemic equilibrium states characterized by different
proportions of susceptible, infected and recovered populations. We consider how
the limit point of the epidemic trajectory and the infection peak along this
trajectory depend on $\sigma$. The heterogeneous model is compared to the the
switched model with an ideally uniform response corresponding to $\sigma=0$.


  Recently, the SARS-CoV-2 variants from the United Kingdom (UK), South Africa,
and Brazil have received much attention for their increased infectivity,
potentially high virulence, and possible threats to existing vaccines and
antibody therapies. The question remains if there are other more infectious
variants transmitted around the world. We carry out a large-scale study of
252,874 SARS-CoV-2 genome isolates from patients to identify many other rapidly
growing mutations on the spike (S) protein receptor-binding domain (RDB). We
reveal that 88 out of 95 significant mutations that were observed more than 10
times strengthen the binding between the RBD and the host
angiotensin-converting enzyme 2 (ACE2), indicating the virus evolves toward
more infectious variants. In particular, we discover new fast-growing RBD
mutations N439K, L452R, S477N, S477R, and N501T that also enhance the RBD and
ACE2 binding. We further unveil that mutation N501Y involved in United Kingdom
(UK), South Africa, and Brazil variants may moderately weaken the binding
between the RBD and many known antibodies, while mutations E484K and K417N
found in South Africa and Brazilian variants can potentially disrupt the
binding between the RDB and many known antibodies. Among three newly identified
fast-growing RBD mutations, L452R, which is now known as part of the California
variant B.1.427, and N501T are able to effectively weaken the binding of many
known antibodies with the RBD. Finally, we hypothesize that RBD mutations that
can simultaneously make SARS-CoV-2 more infectious and disrupt the existing
antibodies, called vaccine escape mutations, will pose an imminent threat to
the current crop of vaccines. A list of most likely vaccine escape mutations is
given, including N501Y, L452R, E484K, N501T, S494P, and K417N.


  Bacterial cells navigate around their environment by directing their movement
along chemical gradients. This process, known as chemotaxis, can promote the
rapid expansion of bacterial populations into previously unoccupied
territories. However, despite numerous experimental and theoretical studies on
this classical topic, chemotaxis-driven population expansion is not understood
in quantitative terms. Building on recent experimental progress, we here
present a detailed analytical study that provides a quantitative understanding
of how chemotaxis and cell growth lead to rapid and stable expansion of
bacterial populations. We provide analytical relations that accurately describe
the dependence of the expansion speed and density profile of the expanding
population on important molecular, cellular, and environmental parameters. In
particular, expansion speeds can be boosted by orders of magnitude when the
environmental availability of chemicals relative to the cellular limits of
chemical sensing is high. As analytical understanding of such complex
spatiotemporal dynamic processes is rare, the results derived here provide a
mathematical framework for further investigations of the different roles
chemotaxis plays in diverse ecological contexts.


  We develop a landscape-flux framework to investigate observed frequency
distributions of vegetation and the stability of these ecological systems under
fluctuations. The frequency distributions can characterize the
population-potential landscape related to the stability of ecological states.
We illustrate the practical utility of this approach by analyzing a
forest-savanna model. Savanna, and Forest states coexist under certain
conditions, consistent with past theoretical work and empirical observations.
However, a new Grassland state, unseen in the corresponding deterministic
model, emerges as an alternative quasi-stable state under fluctuations,
providing a novel theoretical basis for the appearance of widespread grasslands
in some empirical analyses. The ecological dynamics are determined by both the
population-potential landscape gradient and the steady-state probability flux.
The flux quantifies the net input/output to the ecological system and therefore
the degree of nonequilibriumness. Landscape and flux together determine the
transitions between stable states characterized by dominant paths and switching
rates. The intrinsic potential landscape admits a Lyapunov function, which
provides a quantitative measure of global stability. We find that the average
flux, entropy production rate, and free energy have significant changes near
bifurcations under both finite and zero fluctuation. These may provide both
dynamical and thermodynamic origins of the bifurcations. We identified the
variances in observed frequency time traces, fluctuations and time
irreversibility as kinematic measures for bifurcations. This new framework
opens the way to characterize ecological systems globally, to uncover how they
change among states, and to quantify the emergence of new quasi-stable states
under stochastic fluctuations.


  Liquid-liquid phase separation (LLPS) is currently of great interest in cell
biology. LLPS is an example of what is called an emergent phenomenon -- an idea
that comes from condensed-matter physics. Emergent phenomena have the
characteristic feature of having a switch-like response. I show that the Hill
equation of biochemistry can be used as a simple model of strongly cooperative,
switch-like, behaviour. One result is that a switch-like response requires
relatively few molecules, even ten gives a strongly switch-like response. Thus
if a biological function enabled by LLPS relies on LLPS to provide a
switch-like response to a stimulus, then condensates large enough to be visible
in optical microscopy are not needed.


  This paper is concerned with a novel method allowing communication between
FRET nanonetworks and nerve cells. It is focused on two system components:
fluorophores and channelrhodopsins which serve as transmitters and receivers,
respectively. Channelrhodopsins are used here also as a FRET signal-to-voltage
converter. The trade-off between throughput and bit error rate is also
investigated.


  With the advent of high-throughput sequencing (HTS) in molecular biology and
medicine, the need for scalable statistical solutions for modeling complex
biological systems has become of critical importance. The increasing number of
platforms and possible experimental scenarios raised the problem of integrating
large amounts of new heterogeneous data and current knowledge, to test novel
hypotheses and improve our comprehension of physiological processes and
diseases. Although network theory provided a framework to represent biological
systems and study their hidden properties, different algorithms still offer low
reproducibility and robustness, dependence on user-defined setup, and poor
interpretability. Here we discuss the R package SEMgraph, combining network
analysis and causal inference within the framework of structural equation
modeling (SEM). It provides a fully automated toolkit, managing complex
biological systems as multivariate networks, ensuring robustness and
reproducibility through data-driven evaluation of model architecture and
perturbation, that is readily interpretable in terms of causal effects among
system components. In addition, SEMgraph offers several functions for perturbed
path finding, model reduction, and parallelization options for the analysis of
large interaction networks.


  Initial hopes of quickly eradicating the COVID-19 pandemic proved futile, and
the goal shifted to controlling the peak of the infection, so as to minimize
the load on healthcare systems. To that end, public health authorities
intervened aggressively to institute social distancing, lock-down policies, and
other Non-Pharmaceutical Interventions (NPIs). Given the high social,
educational, psychological, and economic costs of NPIs, authorities tune them,
alternatively tightening up or relaxing rules, with the result that, in effect,
a relatively flat infection rate results. For example, during the summer in
parts of the United States, daily infection numbers dropped to a plateau. This
paper approaches NPI tuning as a control-theoretic problem, starting from a
simple dynamic model for social distancing based on the classical SIR epidemics
model. Using a singular-perturbation approach, the plateau becomes a
Quasi-Steady-State (QSS) of a reduced two-dimensional SIR model regulated by
adaptive dynamic feedback. It is shown that the QSS can be assigned and it is
globally asymptotically stable. Interestingly, the dynamic model for social
distancing can be interpreted as a nonlinear integral controller. Problems of
data fitting and parameter identifiability are also studied for this model. The
paper also discusses how this simple model allows for meaningful study of the
effect of population size, vaccinations, and the emergence of second waves.


  We characterize absorption-time distributions for birth-death Markov chains
with an absorbing boundary. For "extinction-prone" chains (which drift on
average toward the absorbing state) the asymptotic distribution is Gaussian,
Gumbel, or belongs to a family of skewed distributions. The latter two cases
arise when the dynamics slow down dramatically near the boundary. Several
models of evolution, epidemics, and chemical reactions fall into these classes;
in each case we establish new results for the absorption-time distribution.
Applications to African sleeping sickness are discussed.


  The spread of COVID-19 caused by the recently discovered SARS-CoV-2 virus has
become a worldwide problem with devastating consequences. To slow down the
spread of the pandemic, mass quarantines have been implemented globally,
provoking further social and economic disruptions. Here, we implement a
comprehensive contact tracing network analysis to find an optimized quarantine
protocol to dismantle the chain of transmission of coronavirus with minimal
disruptions to society. We track billions of anonymized GPS human mobility
datapoints from a compilation of hundreds of mobile apps deployed in Latin
America to monitor the evolution of the contact network of disease transmission
before and after the confinements. As a consequence of the lockdowns, people's
mobility across the region decreases by $\sim$53\%, which results in a drastic
disintegration of the transmission network by $\sim$90\%. However, this
disintegration did not halt the spreading of the disease. Our analysis
indicates that superspreading k-core structures persist in the transmission
network to prolong the pandemic. Once the k-cores are identified, an optimized
strategy to break the chain of transmission is to quarantine a minimal number
of 'weak links' with high betweenness centrality connecting the large k-cores.
As countries built contact tracing apps to fight the pandemic, our results
could turn into a valuable resource to help deploy quarantine protocols with
minimized disruptions.


  The subset selection problem of linear algebra is applied to identify
independent patterns of COVID-19 evolution within Brazil. The data consist of a
set of mortality curves in states of Brazil. A subset of the most independent
curves is selected by using a functional version of the QR matrix decomposition
technique with column pivoting. The selected subset is used next as a basis to
represent the remaining curves filtering out any data redundancy. For each
independent curve, an associated epidemiological region of influence is
defined. The results show two main independent curves with a similar two-peak
pattern and a 50-day shift between the patterns. Two main epidemiological
regions are next identified: one encompassing most of the country from the
center and northeast states to the south, an another one containing the
Amazonian region at the northwest.


  We present a numerical implementation for a multilayer network used to model
the transmission of Covid-19 or other diseases with a similar transmission
mechanism. The model incorporates different contact types between individuals
(household, social contacts, and strangers), which allows flexibility compared
to standard SIR type models. The algorithm described in this paper is a
simplification of the model used to give public health authorities an
additional tool for the decision-making process in Costa Rica, by simulating
extensive possible scenarios and projections.


  Inside a cell, heterotypic proteins assemble in inhomogeneous, crowded
systems where the abundance of these proteins vary with cell types. While some
protein complexes form putative structures that can be visualized with imaging,
there are far more protein complexes that are yet to be solved because of their
dynamic associations with one another. Yet, it is possible to infer these
protein complexes through a physical model. However, it is often not clear to
physicists what kind of data from biology is necessary for such a modeling
endeavor. Here, we aim to model these clusters of coarse-grained protein
assemblies from multiple subunits through the constraints of interactions among
the subunits and the chemical potential of each subunit. We obtained the
constraints on the interactions among subunits from the known protein
structures. We inferred the chemical potential, that dictates the particle
number distribution of each protein subunit, from the knowledge of protein
abundance from experimental data. Guided by the maximum entropy principle, we
formulate an inverse statistical mechanical method to infer the distribution of
particle numbers from the data of protein abundance as chemical potentials for
a grand canonical multi-component mixture. Using grand canonical Monte Carlo
simulations, we captured a distribution of high-order clusters in a protein
complex of Succinate Dehydrogenase (SDH) with four known subunits. The
complexity of hierarchical clusters varies with the relative protein abundance
of each subunit in distinctive cell types such as lung, heart, and brain. When
the crowding content increases, we observed that crowding stabilizes emergent
clusters that do not exist in dilute conditions. We, therefore, proposed a
testable hypothesis that the hierarchical complexity of protein clusters on a
molecular scale is a plausible biomarker of predicting the phenotypes of a
cell.


  Two stochastic models are proposed to describe the evolution of the COVID-19
pandemic. In the first model the population is partitioned into four
compartments: susceptible $S$, infected $I$, removed $R$ and dead people $D$.
In order to have a cross validation, a deterministic version of such a model is
also devised which is represented by a system of ordinary differential
equations with delays. In the second stochastic model two further compartments
are added: the class $A$ of asymptomatic individuals and the class $L$ of
isolated infected people. Effects such as social distancing measures are easily
included and the consequences are analyzed. Numerical solutions are obtained
with Monte Carlo simulations. Quantitative predictions are provided which can
be useful for the evaluation of political measures, e.g. the obtained results
suggest that strategies based on herd immunity are too risky.


  The COVID-19 infection cases have surged globally, causing devastations to
both the society and economy. A key factor contributing to the sustained
spreading is the presence of a large number of asymptomatic or hidden
spreaders, who mix among the susceptible population without being detected or
quarantined. Here we propose an effective non-pharmacological intervention
method of detecting the asymptomatic spreaders in contact-tracing networks, and
validated it on the empirical COVID-19 spreading network in Singapore. We find
that using pure physical spreading equations, the hidden spreaders of COVID-19
can be identified with remarkable accuracy. Specifically, based on the unique
characteristics of COVID-19 spreading dynamics, we propose a computational
framework capturing the transition probabilities among different infectious
states in a network, and extend it to an efficient algorithm to identify
asymptotic individuals. Our simulation results indicate that a screening method
using our prediction outperforms machine learning algorithms, e.g. graph neural
networks, that are designed as baselines in this work, as well as random
screening of infection's closest contacts widely used by China in its early
outbreak. Furthermore, our method provides high precision even with incomplete
information of the contract-tracing networks. Our work can be of critical
importance to the non-pharmacological interventions of COVID-19, especially
with increasing adoptions of contact tracing measures using various new
technologies. Beyond COVID-19, our framework can be useful for other epidemic
diseases that also feature asymptomatic spreading


  The statistical analysis of enzyme kinetic reactions usually involves models
of the response functions which are well defined on the basis of
Michaelis-Menten type equations. The error structure however is often without
good reason assumed as additive Gaussian noise. This simple assumption may lead
to undesired properties of the analysis, particularly when simulations are
involved and consequently negative simulated reaction rates may occur. In this
study we investigate the effect of assuming multiplicative lognormal errors
instead. While there is typically little impact on the estimates, the
experimental designs and their efficiencies are decisively affected,
particularly when it comes to model discrimination problems.


  Models of disease spreading are critical for predicting infection growth in a
population and evaluating public health policies. However, standard models
typically represent the dynamics of disease transmission between individuals
using macroscopic parameters that do not accurately represent person-to-person
variability. To address this issue, we present a dynamic network model that
provides a straightforward way to incorporate both disease transmission
dynamics at the individual scale as well as the full spatiotemporal history of
infection at the population scale. We find that disease spreads through a
social network as a traveling wave of infection, followed by a traveling wave
of recovery, with the onset and dynamics of spreading determined by the
interplay between disease transmission and recovery. We use these insights to
develop a scaling theory that predicts the dynamics of infection for diverse
diseases and populations. Furthermore, we show how spatial heterogeneities in
susceptibility to infection can either exacerbate or quell the spread of
disease, depending on its infectivity. Ultimately, our dynamic network approach
provides a simple way to model disease spreading that unifies previous findings
and can be generalized to diverse diseases, containment strategies, seasonal
conditions, and community structures.


  Transient dynamics can often differ drastically from the asymptotic dynamics
of systems. In this paper we provide a unifying framework for analysing
transient dynamics in marine metapopulations, from the choice of norms to the
addition of stage structure. We use the $\ell_1$ norm, because of its
biological interpretation, to extend the transient metrics of reactivity and
attenuation to marine metapopulations, and use examples to compare these
metrics under the more commonly used $\ell_2$ norm. We then connect the
reactivity and attenuation of marine metapopulations to the source-sink
distribution of habitat patches and demonstrate how to meaningfully measure
reactivity when metapopulations are stage-structured.


  The discovery of disease subtypes is an essential step for developing
precision medicine, and disease subtyping via omics data has become a popular
approach. While promising, subtypes obtained from existing approaches are not
necessarily associated with clinical outcomes. With the rich clinical data
along with the omics data in modern epidemiology cohorts, it is urgent to
develop an outcome-guided clustering algorithm to fully integrate the
phenotypic data with the high-dimensional omics data. Hence, we extended a
sparse K-means method to an outcome-guided sparse K-means (GuidedSparseKmeans)
method. An unified objective function was proposed, which was comprised of (i)
weighted K-means to perform sample clusterings; (ii) lasso regularizations to
perform gene selection from the high-dimensional omics data; (iii)
incorporation of a phenotypic variable from the clinical dataset to facilitate
biologically meaningful clustering results. By iteratively optimizing the
objective function, we will simultaneously obtain a phenotype-related sample
clustering results and gene selection results. We demonstrated the superior
performance of the GuidedSparseKmeans by comparing with existing clustering
methods in simulations and applications of high-dimensional transcriptomic data
of breast cancer and Alzheimer's disease. Our algorithm has been implemented
into an R package, which is publicly available on GitHub
(https://github.com/LingsongMeng/GuidedSparseKmeans).


  The risk for severe illness and mortality from COVID-19 significantly
increases with age. As a result, age-stratified modeling for COVID-19 dynamics
is the key to study how to reduce hospitalizations and mortality from COVID-19.
By taking advantage of network theory, we develop an age-stratified epidemic
model for COVID-19 in complex contact networks. Specifically, we present an
extension of standard SEIR (susceptible-exposed-infectious-removed)
compartmental model, called age-stratified SEAHIR
(susceptible-exposedasymptomatic-hospitalized-infectious-removed) model, to
capture the spread of COVID-19 over multitype random networks with general
degree distributions. We derive several key epidemiological metrics and then
propose an age-stratified vaccination strategy to decrease the mortality and
hospitalizations. Through extensive study, we discover that the outcome of
vaccination prioritization depends on the reproduction number R0. Specifically,
the elderly should be prioritized only when R0 is relatively high. If ongoing
intervention policies, such as universal masking, could suppress R0 at a
relatively low level, prioritizing the high-transmission age group (i.e.,
adults aged 20-39) is most effective to reduce both mortality and
hospitalizations. These conclusions provide useful recommendations for
age-based vaccination prioritization for COVID-19.


  As camera pixel arrays have grown larger and faster, and optical microscopy
techniques ever more refined, there has been an explosion in the quantity of
data acquired during routine light microcopy. At the single-molecule level,
analysis involves multiple steps and can rapidly become computationally
expensive, in some cases intractable on office workstations. Complex bespoke
software can present high activation barriers to entry for new users. Here, we
redevelop our quantitative single-molecule analysis routines into an optimized
and extensible Python program, with GUI and command-line implementations to
facilitate use on local machines and remote clusters, by beginners and advanced
users alike. We demonstrate that its performance is on par with previous MATLAB
implementations but runs an order of magnitude faster. We tested it against
challenge data and demonstrate its performance is comparable to
state-of-the-art analysis platforms. We show the code can extract fluorescence
intensity values for single reporter dye molecules and, using these, estimate
molecular stoichiometries and cellular copy numbers of fluorescently-labeled
biomolecules. It can evaluate 2D diffusion coefficients for the
characteristically short single-particle tracking data. To facilitate
benchmarking we include data simulation routines to compare different analysis
programs. Finally, we show that it works with 2-color data and enables
colocalization analysis based on overlap integration, to infer interactions
between differently labelled biomolecules. By making this freely available we
aim to make complex light microscopy single-molecule analysis more
democratized.


  In this article we propose a novel method to estimate the frequency
distribution of linguistic variables while controlling for statistical
non-independence due to shared ancestry. Unlike previous approaches, our
technique uses all available data, from language families large and small as
well as from isolates, while controlling for different degrees of relatedness
on a continuous scale estimated from the data. Our approach involves three
steps: First, distributions of phylogenies are inferred from lexical data.
Second, these phylogenies are used as part of a statistical model to
statistically estimate transition rates between parameter states. Finally, the
long-term equilibrium of the resulting Markov process is computed. As a case
study, we investigate a series of potential word-order correlations across the
languages of the world.


  Searching for novel molecules with desired chemical properties is crucial in
drug discovery. Existing work focuses on developing neural models to generate
either molecular sequences or chemical graphs. However, it remains a big
challenge to find novel and diverse compounds satisfying several properties. In
this paper, we propose MARS, a method for multi-objective drug molecule
discovery. MARS is based on the idea of generating the chemical candidates by
iteratively editing fragments of molecular graphs. To search for high-quality
candidates, it employs Markov chain Monte Carlo sampling (MCMC) on molecules
with an annealing scheme and an adaptive proposal. To further improve sample
efficiency, MARS uses a graph neural network (GNN) to represent and select
candidate edits, where the GNN is trained on-the-fly with samples from MCMC.
Experiments show that MARS achieves state-of-the-art performance in various
multi-objective settings where molecular bio-activity, drug-likeness, and
synthesizability are considered. Remarkably, in the most challenging setting
where all four objectives are simultaneously optimized, our approach
outperforms previous methods significantly in comprehensive evaluations. The
code is available at https://github.com/yutxie/mars.


  Separating environmental effects from those of interspecific interactions on
species distributions has always been a central objective of community ecology.
Despite years of effort in analysing patterns of species co-occurrences and the
developments of sophisticated tools, we are still unable to address this major
objective. A key reason is that the wealth of ecological knowledge is not
sufficiently harnessed in current statistical models, notably the knowledge on
interspecific interactions. Here, we develop ELGRIN, a statistical model that
simultaneously combines knowledge on interspecific interactions (i.e., the
metanetwork), environmental data and species occurrences to tease apart their
relative effects on species distributions. Instead of focusing on single
effects of pairwise species interactions, which have little sense in complex
communities, ELGRIN contrasts the overall effect of species interactions to
that of the environment. Using various simulated and empirical data, we
demonstrate the suitability of ELGRIN to address the objectives for various
types of interspecific interactions like mutualism, competition and trophic
interactions. We then apply the model on vertebrate trophic networks in the
European Alps to map the effect of biotic interactions on species
distributions.Data on ecological networks are everyday increasing and we
believe the time is ripe to mobilize these data to better understand
biodiversity patterns. ELGRIN provides this opportunity to unravel how
interspecific interactions actually influence species distributions.


  Sequencing by Emergence (SEQE) is a new single-molecule nucleic acid
(DNA/RNA) sequencing technology that estimates sequence as an emergent property
of the binding and localization of a repertoire of short oligonucleotide
probes. SEQE promises to deliver accurate, ultra-long, haplotype-phased reads
at the whole genome-scale for very low cost within 10 minutes. The data SEQE
generates requires entirely new inference techniques. In this paper we
introduce a probabilistic model of the SEQE measurement process and an
algorithm that estimates sequence by solving a convex relaxation of the
corresponding maximum likelihood problem. We demonstrate the effectiveness of
our algorithm on a variety of simulated datasets.


  The CRISPR/Cas9 system acts as the prokaryotic immune system and has
important applications in gene editing. The protein Cas9 is one of its crucial
components. The role of Cas9 is to search for specific target sequences on the
DNA and cleave them. In this Letter, we introduce a model of facilitated
diffusion for Cas9 and fit its parameters to single-molecule experiments. Our
model confirms that Cas9 search for targets by sliding, but shows that its
sliding length is rather short. We then investigate how Cas9 explores a long
stretch of DNA containing randomly placed targets. We solve this problem by
mapping it into the theory of Anderson localization in condensed matter
physics. Our theoretical approach rationalizes experimental evidences on the
distribution of Cas9 molecules along the DNA.


  Protein-protein interactions are the basis of many important physiological
processes and are currently promising, yet difficult, targets for drug
discovery. In this context, inhibitor of apoptosis proteins (IAPs)-mediated
interactions are pivotal for cancer cell survival; the interaction of the BIR1
domain of cIAP2 with TRAF2 was shown to lead the recruitment of cIAPs to the
TNF receptor, promoting the activation of the NF-\kappa B survival pathway. In
this work, using a combined in silico-in vitro approach, we identified a
drug-like molecule, NF023, able to disrupt cIAP2 interaction with TRAF2. We
demonstrated in vitro its ability to interfere with the assembly of the
cIAP2-BIR1/TRAF2 complex and performed a thorough characterization of the
compound's mode of action through 248 parallel unbiased molecular dynamics
simulations of 300 ns (totaling almost 75 {\mu}s of all-atom sampling), which
identified multiple binding modes to the BIR1 domain of cIAP2 via clustering
and ensemble docking. NF023 is, thus, a promising protein-protein interaction
disruptor, representing a starting point to develop modulators of NF-\kappa
B-mediated cell survival in cancer. This study represents a model procedure
that shows the use of large-scale molecular dynamics methods to typify
promiscuous interactors.


  Recent experiments showed that multiple copies of the molecular machine RNA
polymerase (RNAP) can efficiently synthesize mRNA collectively in the active
state of the promoter. However, environmentally-induced promoter repression
results in long-distance antagonistic interactions that drastically reduce the
speed of RNAPs and cause a quick arrest of mRNA synthesis. The mechanism
underlying this transition between cooperative and antagonistic dynamics
remains poorly understood. In this Letter, we introduce a continuum
deterministic model for the translocation of RNAPs, where the speed of an RNAP
is coupled to the local DNA supercoiling as well as the density of RNAPs on the
gene. We assume that torsional stress experienced by individual RNAPs is
exacerbated by high RNAP density on the gene and that transcription factors act
as physical barriers to the diffusion of DNA supercoils. We show that this
minimal model exhibits two transcription modes mediated by the torsional
stress: a fluid mode when the promoter is active and a torsionally stressed
mode when the promoter is repressed, in quantitative agreement with
experimentally observed dynamics of co-transcribing RNAPs. Our work provides an
important step towards understanding the collective dynamics of molecular
machines involved in gene expression.


  Cooperative interactions pervade in a broad range of many-body populations,
such as ecological communities, social organizations, and economic webs. We
investigate the dynamics of a population of two equivalent species A and B that
are driven by cooperative and symmetric interactions between these species. For
an isolated population, we determine the probability to reach fixation, where
only one species remains, as a function of the initial concentrations of the
two species, as well as the time to reach fixation. The latter scales
exponentially with the population size. When members of each species migrate
into the population at rate $\lambda$ and replace a randomly selected
individual, surprisingly rich dynamics ensues. Ostensibly, the population
reaches a steady state, but the steady-state population distribution undergoes
a unimodal to trimodal transition as the migration rate decreases below a
critical value $\lambda_c$. In the low-migration regime, $\lambda<\lambda_c$,
the steady state is not truly steady, but instead strongly fluctuates between
near-fixation states, where the population consists of mostly A's or of mostly
B's. The characteristic time scale of these fluctuations diverges as
$\lambda^{-1}$. Thus in spite of the cooperative interaction, a typical
snapshot of the population will contain almost all A's or almost all B's.


  The chemical master equation (CME), which describes the discrete and
stochastic molecule number dynamics associated with biological processes like
transcription, is difficult to solve analytically. It is particularly hard to
solve for models involving bursting/gene switching, a biological feature that
tends to produce heavy-tailed single cell RNA counts distributions. In this
paper, we present a novel method for computing exact and analytic solutions to
the CME in such cases, and use these results to explore approximate solutions
valid in different parameter regimes, and to compute observables of interest.
Our method leverages tools inspired by quantum mechanics, including ladder
operators and Feynman-like diagrams, and establishes close formal parallels
between the dynamics of bursty transcription, and the dynamics of bosons
interacting with a single fermion. We focus on two problems: (i) the chemical
birth-death process coupled to a switching gene/the telegraph model, and (ii) a
model of transcription and multistep splicing involving a switching gene and an
arbitrary number of downstream splicing steps. We work out many special cases,
and exhaustively explore the special functionology associated with these
problems. This is Part I in a two-part series of papers; in Part II, we explore
an alternative solution approach that is more useful for numerically solving
these problems, and apply it to parameter inference on simulated RNA counts
data.


  We propose a novel model to explain the mechanisms underlying dominance
hierarchical structures. Guided by a predetermined social convention, agents
with limited cognitive abilities optimize their strategies in a Hawk-Dove game.
We find that several commonly observed hierarchical structures in nature such
as linear hierarchy and despotism, emerge as the total fitness-maximizing
social structures given different levels of cognitive abilities.


  Here I investigate some mathematical aspects of the maximum entropy theory of
ecology (METE). In particular I address the geometrical structure of METE
endowed by information geometry. As novel results, the macrostate entropy is
calculated analytically by the Legendre transformation of the log-normalizer in
METE. This result allows for the calculation of the metric terms in the
information geometry arising from METE and, by consequence, the covariance
matrix between METE variables.


  We introduce a one-parameter family of polymatrix replicators defined in a
three-dimensional cube and study its bifurcations. For a given interval of
parameters, this family exhibits suspended horseshoes and persistent strange
attractors. The proof relies on the existence of a homoclinic cycle to the
interior equilibrium. We also describe the phenomenological steps responsible
for the transition from regular to chaotic dynamics in our system (route to
chaos).


  We present an epidemiological compartment model, SAIR(S), that explicitly
captures the dynamics of asymptomatic infected individuals in an epidemic
spread process. We first present a group model and then discuss networked
versions. We provide an investigation of equilibria and stability properties
for these models, and present simulation results illustrating the effects of
asymptomatic-infected individuals on the spread of the disease. We also discuss
local isolation effects on the epidemic dynamics in terms of the networked
models. Finally, we provide initial parameter estimation results based on
simple least-squares approaches and local test-site data.\par Keywords:
Epidemic dynamics, networks, data-informed modeling, stability analysis,
parameter estimation


  The spatial segregation of species is fundamental to ecosystem formation and
stability. Behavioural strategies may determine where species are located and
how their interactions change the local environment arrangement. In response to
stimuli in the environment, individuals may move in a specific direction
instead of walking randomly. This behaviour can be innate or learned from
experience, and allow the individuals to conquer or the maintain territory,
foraging or taking refuge. We study a generalisation of the spatial
rock-paper-scissors model where individuals of one out of the species may
perform directional movement tactics. Running a series of stochastic
simulations, we investigate the effects of the behavioural tactics on the
spatial pattern formation and the maintenance of the species diversity. We also
explore a more realistic scenario, where not all individuals are conditioned to
perform the behavioural strategy or have different levels of neighbourhood
perception. Our outcomes show that self-preservation behaviour is more
profitable in terms of territorial dominance, with the best result being
achieved when all individuals are conditioned and have a long-range vicinity
perception. On the other hand, invading is more advantageous if part of
individuals is conditioned and if they have short-range neighbourhood
perception. Finally, our findings reveal that the self-defence strategy is the
least jeopardising to biodiversity which can help biologists to understand
population dynamics in a setting where individuals may move strategically.


  Seed formation is part of the reproductive cycle, leading to the accumulation
of resistance stages that can withstand harsh environmental conditions for long
periods of time. At the community level, multiple species with such
long-lasting life stages can be more likely to coexist. While the implications
of this process for biodiversity have been studied in terrestrial plants, seed
banks are usually neglected in phytoplankton multispecies dynamic models, in
spite of widespread empirical evidence for such seed banks. In this study, we
build a metacommunity model of interacting phytoplankton species, including a
resting stage supplying the seed bank. The model is parameterized with
empirically-driven growth rate functions and field-based interaction estimates,
which include both facilitative and competitive interactions. Exchanges between
compartments (coastal pelagic cells, coastal resting cells on the seabed, and
open ocean pelagic cells) are controlled by hydrodynamical parameters to which
the sensitivity of the model is assessed. We consider two models, i.e., with
and without a saturating effect of the interactions on the growth rates. Our
results are consistent between models, and show that a seed bank allows to
maintain all species in the community over 30 years. Indeed, a fraction of the
species are vulnerable to extinction at specific times within the year, but
this process is buffered by their survival in their resting stage. We thus
highlight the potential role of the seed bank in the recurrent re-invasion of
the coastal community, and of coastal environments in re-seeding oceanic
regions. Moreover, the seed bank enables populations to tolerate stronger
interactions within the community as well as more severe changes to the
environment, such as those predicted in a climate change context. Our study
therefore shows how resting stages may help phytoplanktonic diversity
maintenance.


  In this article, we review the different models and methods to understand and
estimate correlated evolution within the same genome, individual or species. We
describe correlated evolution among traits, among genetic components and
finally between traits and genetics.


  This article considers a mechanism to explain the emergence and evolution of
social cooperation. Selfish individuals tend to benefit themselves, which makes
it hard for the maintenance of cooperation between unrelated individuals. We
propose and validate that a smart group strategy can effectively facilitate the
evolution of cooperation, provided cooperators spatially align whilst
cooperating at a new level of alliance. The general evolutionary model
presented here shows that a non-kin selection effect is a possible cause for
cooperation between unrelated individuals and highlights that non-kin selection
may be a hallmark of biological evolution.


  Domain motions involved in the function of proteins can often be well
described as a combination of motions along a handfull of low-frequency modes,
that is, with the values of a few normal coordinates. This means that, when the
functional motion of a protein is unknown, it should prove possible to predict
it, since it amounts to guess a few values. However, without the help of
additional experimental data, using normal coordinates for generating accurate
conformers far away from the initial one is not so straightforward. To do so, a
new approach is proposed: instead of building conformers directly with the
values of a subset of normal coordinates, they are built in two steps, the
conformer built with normal coordinates being just used for defining a set of
distance constraints, the final conformer being built so as to match them. Note
that this approach amounts to transform the problem of generating accurate
protein conformers using normal coordinates into a better known one: the
distance-geometry problem, which is herein solved with the help of the ROSETTA
software. In the present study, this approach allowed to rebuild accurately six
large amplitude conformational changes, using at most six low-frequency normal
coordinates. As a consequence of the low-dimensionality of the corresponding
subspace, random exploration also proved enough for generating low-energy
conformers close to the known end-point of the conformational change of the LAO
binding protein, lysozyme T4 and adenylate kinase.


  FoF1-ATP synthases are ubiquitous membrane-bound, rotary motor enzymes that
can catalyze ATP synthesis and hydrolysis. Their enzyme kinetics are controlled
by internal subunit rotation, by substrate and product concentrations, by
mechanical inhibitory mechanisms, but also by the electrochemical potential of
protons across the membrane. Single-molecule F\"orster resonance energy
transfer (smFRET) has been used to detect subunit rotation within FoF1-ATP
synthases embedded in freely diffusing liposomes. We now report that kinetic
monitoring of functional rotation can be prolonged from milliseconds to seconds
by utilizing an Anti-Brownian electrokinetic trap (ABEL trap). These extended
observation times allowed us to observe fluctuating rates of functional
rotation for individual FoF1-liposomes in solution. Broad distributions of
ATP-dependent catalytic rates were revealed. The buildup of an electrochemical
potential of protons was confirmed to limit the maximum rate of ATP hydrolysis.
In the presence of ionophores or uncouplers, the fastest subunit rotation
speeds measured in single reconstituted FoF1-ATP synthases were 180 full rounds
per second. This was much faster than measured by biochemical ensemble
averaging, but not as fast as the maximum rotational speed reported previously
for isolated single F1 complexes uncoupled from the membrane-embedded Fo
complex. Further application of ABEL trap measurements should help resolve the
mechanistic causes of such fluctuating rates of subunit rotation.


  We investigate the application of the line-graph operator to one-dimensional
spin models with periodic boundary conditions. The spins (or interactions) in
the original spin structure become the interactions (or spins) in the resulting
spin structure. We identify conditions which ensure that each new spin
structure is stable, that is, its spin configuration minimises its internal
energy. Then, making a correspondence between spin configurations and binary
sequences, we propose a model of information growth and evolution based on the
line-graph operator. Since this operator can generate frustrations in newly
formed spin chains, in the proposed model such frustrations are immediately
removed. Also, in some cases, the previously frustrated chains are allowed to
recombine into new stable chains. As a result, we obtain a population of spin
chains whose dynamics is studied using Monte Carlo simulations. Lastly, we
discuss potential applications to areas of research such as combinatorics and
theoretical biology.


  In a series of recent papers and in a book, this author put forward a
mathematical model capable of embracing the search for extra-terrestrial
intelligence (SETI), Darwinian Evolution and Human History into a single,
unified statistical picture, concisely called Evo-SETI. The relevant
mathematical tools are: (1) Geometric Brownian motion (GBM), the stochastic
process representing evolution as the stochastic increase of the number of
species living on Earth over the last 3.5 billion years. This GBM is well known
in the mathematics of finances (Black-Sholes models). (2) The probability
distributions known as b-lognormals, i.e. lognormals starting at a certain
positive instant b>0 rather than at the origin. In the framework of Darwinian
Evolution, the resulting mathematical construction was shown to be what
evolutionary biologists call Cladistics. (3) The (Shannon) entropy of such
b-lognormals is then seen to represent the 'degree of progress' reached by each
living organism or by each big set of living organisms, like historic human
civilizations. (4) All these results also match with SETI in that the
statistical Drake equation (generalization of the ordinary Drake equation to
encompass statistics) leads just to the lognormal distribution as the
probability distribution for the number of extra-terrestrial civilizations
existing in the Galaxy. (5) The well-known 'Molecular Clock of Evolution',
namely the 'constant rate of Evolution at the molecular level' as shown by
Kimura's Neutral Theory of Molecular Evolution, identifies with growth rate of
the entropy of our Evo-SETI model. (6) Furthermore, we apply our Evo-SETI model
to lognormal stochastic processes other than GBMs. For instance, we provide two
models for the mass extinctions that occurred in the past. (7) Finally, we show
that the Markov & Korotayev model for Darwinian Evolution identifies with an
Evo-SETI model.


  Savanna ecosystems are shaped by the frequency and intensity of regular
fires. We model savannas via an ordinary differential equation (ODE) encoding a
one-sided inhibitory Lotka-Volterra interaction between trees and grass. By
applying fire as a discrete disturbance, we create an impulsive dynamical
system that allows us to identify the impact of variation in fire frequency and
intensity. The model exhibits three different bistability regimes: between
savanna and grassland; two savanna states; and savanna and woodland. The
impulsive model reveals rich bifurcation structures in response to changes in
fire intensity and frequency -- structures that are largely invisible to
analogous ODE models with continuous fire. In addition, by using the amount of
grass as an example of a socially-valued function of the system state, we
examine the resilience of the social value to different disturbance regimes. We
find that large transitions ("tipping") in the valued quantity can be triggered
by small changes in disturbance regime.


  Ecological systems are studied using many different approaches and
mathematical tools. One approach, based on the Jacobian of Lotka-Volterra type
models, has been a staple of mathematical ecology for years, leading to many
ideas such as on questions of system stability. Instability in such methods is
determined by the presence of an eigenvalue of the community matrix lying in
the right half plane. The coefficients of the characteristic polynomial derived
from community matrices contain information related to the specific matrix
elements that play a greater destabilising role. Yet the destabilising
circuits, or cycles, constructed by multiplying these elements together, form
only a subset of all the feedback loops comprising a given system. This paper
looks at the destabilising feedback loops in predator-prey, mutualistic and
competitive systems in terms of sets of the matrix elements to explore how sign
structure affects how the elements contribute to instability. This leads to
quite rich combinatorial structure among the destabilising cycle sets as set
size grows within the coefficients of the characteristic polynomial.


  It has long been hypothesized that operating close to the critical state is
beneficial for natural and artificial systems. We test this hypothesis by
evolving foraging agents controlled by neural networks that can change the
system's dynamical regime throughout evolution. Surprisingly, we find that all
populations, regardless of their initial regime, evolve to be subcritical in
simple tasks and even strongly subcritical populations can reach comparable
performance. We hypothesize that the moderately subcritical regime combines the
benefits of generalizability and adaptability brought by closeness to
criticality with the stability of the dynamics characteristic for subcritical
systems. By a resilience analysis, we find that initially critical agents
maintain their fitness level even under environmental changes and degrade
slowly with increasing perturbation strength. On the other hand, subcritical
agents originally evolved to the same fitness, were often rendered utterly
inadequate and degraded faster. We conclude that although the subcritical
regime is preferable for a simple task, the optimal deviation from criticality
depends on the task difficulty: for harder tasks, agents evolve closer to
criticality. Furthermore, subcritical populations cannot find the path to
decrease their distance to criticality. In summary, our study suggests that
initializing models near criticality is important to find an optimal and
flexible solution.


  Immune checkpoint therapy is one of the most promising immunotherapeutic
methods that are likely able to give rise to durable treatment response for
various cancer types. Despite much progress in the past decade, there are still
critical open questions with particular regards to quantifying and predicting
the efficacy of treatment and potential optimal regimens for combining
different immune-checkpoint blockades. To shed light on this issue, here we
develop clinically-relevant, dynamical systems models of cancer immunotherapy
with a focus on the immune checkpoint PD-1/PD-L1 blockades. Our model allows
the acquisition of adaptive immune resistance in the absence of treatment,
whereas immune checkpoint blockades can reverse such resistance and boost
anti-tumor activities of effector cells. Our numerical analysis predicts that
anti-PD-1 agents are commonly less effective than anti-PD-L1 agents for a wide
range of model parameters. We also observe that combination treatment of
anti-PD-1 and anti-PD-L1 blockades leads to a desirable synergistic effect. Our
modeling framework lays the ground for future data-driven analysis on
combination therapeutics of immune-checkpoint treatment regimes and thorough
investigation of optimized treatment on a patient-by-patient basis.


  Most of the ongoing projects aimed at the development of specific therapies
and vaccines against COVID-19 use the SARS-CoV-2 spike (S) protein as the main
target [1-3]. The binding of the spike protein with the ACE2 receptor (ACE2) of
the host cell constitutes the first and key step for virus entry. During this
process, the receptor binding domain (RBD) of the S protein plays an essential
role, since it contains the receptor binding motif (RBM), responsible for the
docking to the receptor. So far, mostly biochemical methods are being tested in
order to prevent binding of the virus to ACE2 [4]. Here we show, with the help
of atomistic simulations, that external electric fields of easily achievable
and moderate strengths can dramatically destabilise the S protein, inducing
long-lasting structural damage. One striking field-induced conformational
change occurs at the level of the recognition loop L3 of the RBD where two
parallel beta sheets, believed to be responsible for a high affinity to ACE2
[5], undergo a change into an unstructured coil, which exhibits almost no
binding possibilities to the ACE2 receptor (Figure 1a). Remarkably, while the
structural flexibility of S allows the virus to improve its probability of
entering the cell, it is also the origin of the surprising vulnerability of S
upon application of electric fields of strengths at least two orders of
magnitude smaller than those required for damaging most proteins. Our findings
suggest the existence of a clean physical method to weaken the SARS-CoV-2 virus
without further biochemical processing. Moreover, the effect could be used for
infection prevention purposes and also to develop technologies for in-vitro
structural manipulation of S. Since the method is largely unspecific, it can be
suitable for application to mutations in S, to other proteins of SARS-CoV-2 and
in general to membrane proteins of other virus types.


  Respiratory syncytial virus (RSV) is a leading cause of acute lower
respiratory tract infection worldwide, resulting in approximately sixty
thousand annual hospitalizations of <5-year-olds in the United States alone and
three million annual hospitalizations globally. The development of over 40
vaccines and immunoprophylactic interventions targeting RSV has the potential
to significantly reduce the disease burden from RSV infection in the near
future. In the context of RSV, a highly contagious pathogen, dynamic
transmission models (DTMs) are valuable tools in the evaluation and comparison
of the effectiveness of different interventions. This review, the first of its
kind for RSV DTMs, provides a valuable foundation for future modelling efforts
and highlights important gaps in our understanding of RSV epidemics.
Specifically, we have searched the literature using Web of Science, Scopus,
Embase, and PubMed to identify all published manuscripts reporting the
development of DTMs focused on the population transmission of RSV. We reviewed
the resulting studies and summarized the structure, parameterization, and
results of the models developed therein. We anticipate that future RSV DTMs,
combined with cost-effectiveness evaluations, will play a significant role in
shaping decision making in the development and implementation of intervention
programs.


  While a common trend in disease modeling is to develop models of increasing
complexity, it was recently pointed out that outbreaks appear remarkably simple
when viewed in the incidence vs. cumulative cases (ICC) plane. This article
details the theory behind this phenomenon by analyzing the stochastic SIR
(Susceptible, Infected, Recovered) model in the cumulative cases domain. We
prove that the Markov chain associated with this model reduces, in the ICC
plane, to a pure birth chain for the cumulative number of cases, whose limit
leads to an independent increments Gaussian process that fluctuates about a
deterministic ICC curve. We calculate the associated variance and quantify the
additional variability due to estimating incidence over a finite period of
time. We also illustrate the universality brought forth by the ICC concept on
real-world data for Influenza A and for the COVID-19 outbreak in Arizona.


  The photo-triggered unbinding of the intrinsically disordered S-peptide from
the RNase S complex is studied with the help of transient IR spectroscopy,
covering a wide range of time scales from 100 ps to 10 ms. To that end, an
azobenzene moiety has been linked to the S-peptide in a way that its helicity
is disrupted by light, thereby initiating its complete unbinding. The full
sequence of events is observed, starting from unfolding of the helical
structure of the S-peptide on a 20 ns timescale while still being in the
binding pocket of the S-protein, S-peptide unbinding after 300 microseconds,
and the structural response of the S-protein after 3 ms. With regard to the
S-peptide dynamics, the binding mechanism can be classified as an induced fit,
while the structural response of the S-protein is better described as
conformational selection.


  The war of attrition in game theory is a model of a stand-off situation
between two opponents where the winner is determined by its persistence. We
model a stand-off between a predator and a prey when the prey is hiding and the
predator is waiting for the prey to come out from its refuge, or when the two
are locked in a situation of mutual threat of injury or even death. The
stand-off is resolved when the predator gives up or when the prey tries to
escape. Instead of using the asymmetric war of attrition, we embed the
stand-off as an integral part of the predator-prey model of Rosenzweig and
MacArthur derived from first principles. We apply this model to study the
coevolution of the giving-up rates of the prey and the predator, using the
adaptive dynamics approach. We find that the long term evolutionary process
leads to three qualitatively different scenarios: the predator gives up
immediately, while the prey never gives up; the predator never gives up, while
the prey adopts any giving-up rate greater than or equal to a given positive
threshold value; the predator goes extinct. We observe that some results are
the same as for the asymmetric war of attrition, but others are quite
different.


  Each gene has its own evolutionary history which can substantially differ
from the evolutionary histories of other genes. For example, some individual
genes or operons can be affected by specific horizontal gene transfer and
recombination events. Thus, the evolutionary history of each gene should be
represented by its own phylogenetic tree which may display different
evolutionary patterns from the species tree that accounts for the main patterns
of vertical descent. The output of traditional consensus tree or supertree
inference methods is a unique consensus tree or supertree. We describe a new
efficient method for inferring multiple alternative consensus trees and
supertrees to best represent the most important evolutionary patterns of a
given set of gene phylogenies. We show how an adapted version of the popular
k-means clustering algorithm, based on some interesting properties of the
Robinson and Foulds distance, can be used to partition a given set of trees
into one (for homogeneous data) or multiple (for heterogeneous data) cluster(s)
of trees. Moreover, we adapt the popular Cali\'nski-Harabasz, Silhouette, Ball
and Hall, and Gap cluster validity indices to tree clustering with k-means. A
special attention is given to the relevant but very challenging problem of
inferring alternative supertrees. The use of the Euclidean property of the
objective function of the method makes it faster than the existing tree
clustering techniques, and thus perfectly suitable for analyzing large
evolutionary datasets. We apply the new method to discover alternative
supertrees characterizing the main patterns of evolution of SARS-CoV-2 and
genetically related betacoronaviruses.


  We propose a new modeling framework for evaluating the risk of disease
transmission during a pandemic in small-scale settings driven by stochasticity
in the arrival and service processes, i.e., congestion-prone confined-space
service facilities, such as grocery stores. We propose a novel metric
($R_0^{\mathsf{sys}}$) inspired by $R_0$, the "basic reproduction number"
concept from epidemiology, which measures the transmissibility of infectious
diseases. We derive the $R_0^{\mathsf{sys}}$ metric for various queueing models
of service facilities by leveraging a novel queueing-theoretic notion: sojourn
time overlaps. We showcase how our metric can be used to explore the efficacy
of a variety of interventions aimed at curbing the spread of disease inside
service facilities. Specifically, we focus on some prevalent interventions
employed during the COVID-19 pandemic: limiting the occupancy of service
facilities, protecting high-risk customers (via prioritization or designated
time windows), and increasing the service speed (or limiting patronage
duration). We discuss a variety of directions for adapting our transmission
model to incorporate some more nuanced features of disease transmission,
including heterogeneity in the population immunity level, varying levels of
mask usage, and spatial considerations in disease transmission.


  For a fixed set $X$ containing $n$ taxon labels, an ordered pair consisting
of a gene tree topology $G$ and a species tree $S$ bijectively labeled with the
labels of $X$ possesses a set of coalescent histories -- mappings from the set
of internal nodes of $G$ to the set of edges of $S$ describing possible lists
of edges in $S$ on which the coalescences in $G$ take place. Enumerations of
coalescent histories for gene trees and species trees have produced suggestive
results regarding the pairs $(G,S)$ that, for a fixed $n$, have the largest
number of coalescent histories. We define a class of 2-cherry binary tree
topologies that we term $p$-pseudocaterpillars, examining coalescent histories
for non-matching pairs $(G,S)$, in the case in which $S$ has a caterpillar
shape and $G$ has a $p$-pseudocaterpillar shape. Using a construction that
associates coalescent histories for $(G,S)$ with a class of "roadblocked"
monotonic paths, we identify the $p$-pseudocaterpillar labeled gene tree
topology that, for a fixed caterpillar labeled species tree topology, gives
rise to the largest number of coalescent histories. The shape that maximizes
the number of coalescent histories places the "second" cherry of the
$p$-pseudocaterpillar equidistantly from the root of the "first" cherry and
from the tree root. A symmetry in the numbers of coalescent histories for
$p$-pseudocaterpillar gene trees and caterpillar species trees is seen to exist
around the maximizing value of the parameter $p$. The results provide insight
into the factors that influence the number of coalescent histories possible for
a given gene tree and species tree.