  Since early 2020, the world has been dealing with a raging pandemic outbreak:
COVID-19. A year later, vaccines have become accessible, but in limited
quantities, so that governments needed to devise a strategy to decide which
part of the population to prioritize when assigning the available doses, and
how to manage the interval between doses for multi-dose vaccines. In this
paper, we present an optimization framework to address the dynamic double-dose
vaccine allocation problem whereby the available vaccine doses must be
administered to different age-groups to minimize specific societal objectives.
In particular, we first identify an age-dependent
Susceptible-Exposed-Infected-Recovered (SEIR) epidemic model including an
extension capturing partially and fully vaccinated people, whereby we account
for age-dependent immunity and infectiousness levels together with disease
severity. Second, we leverage our model to frame the dynamic age-dependent
vaccine allocation problem for different societal objectives, such as the
minimization of infections or fatalities, and solve it with nonlinear
programming techniques. Finally, we carry out a numerical case study with
real-world data from The Netherlands. Our results show how different societal
objectives can significantly alter the optimal vaccine allocation strategy. For
instance, we find that minimizing the overall number of infections results in
delaying second doses, whilst to minimize fatalities it is important to fully
vaccinate the elderly first.


  We revisit here a landmark five parameter SIR-type model of [DvdD93, Sec. 4],
which is maybe the simplest example where a complete picture of all cases,
including non-trivial bistability behavior, may be obtained using simple tools.
We also generalize it by adding essential vaccination and vaccination-induced
death parameters, with the aim of revealing the role of vaccination and its
possible failure. The main result is Theorem 5, which describes the stability
behavior of our model in all possible cases.


  We have designed a computational model of a virus spread near the outbreak
threshold. Using computer simulation we studied the Susceptible - Infected -
Recovered (SIR) process where in consequence of a force of habit that is
manifested by the population mobility patterns, the recovered persons create
the spatio-temporal patterns as the barriers to a virus transmission. The
results show a spontaneous stopping of the virus spread without a need to
infect the whole population, a non-trivial random noise of daily count of
infected cases, and power laws of a cumulative count of infected cases.
Outbreak evolution strongly depends on the initial conditions thus we concluded
that the model has the features of chaotic systems that makes it difficult to
predict its behaviors.


  Recent advances in protein function prediction exploit graph-based deep
learning approaches to correlate the structural and topological features of
proteins with their molecular functions. However, proteins in vivo are not
static but dynamic molecules that alter conformation for functional purposes.
Here we apply normal mode analysis to native protein conformations and augment
protein graphs by connecting edges between dynamically correlated residue
pairs. In the multilabel function classification task, our method demonstrates
a remarkable performance gain based on this dynamics-informed representation.
The proposed graph neural network, ProDAR, increases the interpretability and
generalizability of residue-level annotations and robustly reflects structural
nuance in proteins. We elucidate the importance of dynamic information in graph
representation by comparing class activation maps for hMTH1, nitrophorin, and
SARS-CoV-2 receptor binding domain. Our model successfully learns the dynamic
fingerprints of proteins and pinpoints the residues of functional impacts, with
vast untapped potential for broad biotechnology and pharmaceutical
applications.


  Pharmaceutical and non-pharmaceutical interventions (NPIs) have been crucial
for controlling COVID-19. They are complemented by voluntary health-protective
behaviour, building a complex interplay between risk perception, behaviour, and
disease spread. We studied how voluntary health-protective behaviour and
vaccination willingness impact the long-term dynamics. We analysed how
different levels of mandatory NPIs determine how individuals use their leeway
for voluntary actions. If mandatory NPIs are too weak, COVID-19 incidence will
surge, implying high morbidity and mortality before individuals react; if they
are too strong, one expects a rebound wave once restrictions are lifted,
challenging the transition to endemicity. Conversely, moderate mandatory NPIs
give individuals time and room to adapt their level of caution, mitigating
disease spread effectively. When complemented with high vaccination rates, this
also offers a robust way to limit the impacts of the Omicron variant of
concern. Altogether, our work highlights the importance of appropriate
mandatory NPIs to maximise the impact of individual voluntary actions in
pandemic control.


  Large-scale human mobility data is a key resource in data-driven policy
making and across many scientific fields. Most recently, mobility data was
extensively used during the COVID-19 pandemic to study the effects of
governmental policies and to inform epidemic models. Large-scale mobility is
often measured using digital tools such as mobile phones. However, it remains
an open question how truthfully these digital proxies represent the actual
travel behavior of the general population. Here, we examine mobility datasets
from multiple countries and identify two fundamentally different types of bias
caused by unequal access to, and unequal usage of mobile phones. We introduce
the concept of data generation bias, a previously overlooked type of bias,
which is present when the amount of data that an individual produces influences
their representation in the dataset. We find evidence for data generation bias
in all examined datasets in that high-wealth individuals are overrepresented,
with the richest 20% contributing over 50% of all recorded trips, substantially
skewing the datasets. This inequality is consequential, as we find mobility
patterns of different wealth groups to be structurally different, where the
mobility networks of high-wealth users are denser and contain more long-range
connections. To mitigate the skew, we present a framework to debias data and
show how simple techniques can be used to increase representativeness. Using
our approach we show how biases can severely impact outcomes of dynamic
processes such as epidemic simulations, where biased data incorrectly estimates
the severity and speed of disease transmission. Overall, we show that a failure
to account for biases can have detrimental effects on the results of studies
and urge researchers and practitioners to account for data-fairness in all
future studies of human mobility.


  We consider a stochastic individual-based model of adaptive dynamics on a
finite trait graph $G=(V,E)$. The evolution is driven by a linear birth rate, a
density dependent logistic death rate an the possibility of mutations along the
(possibly directed) edges in $E$. We study the limit of small mutation rates
for a simultaneously diverging population size. Closing the gap between the
works of Bovier, Coquille and Smadi (2019) and Coquille, Kraut and Smadi
(2021), we give a precise description of transitions between evolutionary
stable conditions (ESC), where multiple mutations are needed to cross a valley
in the fitness landscape. The system shows a metastable behaviour on several
divergent time scales associated to a degree of stability. We develop the
framework of a meta graph that is constituted of ESCs and possible metastable
transitions between those. This allows for a concise description of the
multi-scale jump chain arising from concatenating several jumps. Finally, for
each of the various time scale, we prove the convergence of the population
process to a Markov jump process visiting only ESCs of sufficiently high
stability.


  Malaria is one of the deadliest infectious diseases globally, causing
hundreds of thousands of deaths each year. It disproportionately affects young
children, with two-thirds of fatalities occurring in under-fives. Individuals
acquire protection from disease through repeated exposure, and this immunity
plays a crucial role in the dynamics of malaria spread. We develop a novel
age-structured PDE malaria model, which couples vector-host epidemiological
dynamics with immunity dynamics. Our model tracks the acquisition and loss of
anti-disease immunity during transmission and its corresponding nonlinear
feedback onto the transmission parameters. We derive the basic reproduction
number ($\mathcal{R}_0$) as the threshold condition for the stability of
disease-free equilibrium; we also interpret $\mathcal{R}_0$ probabilistically
as a weighted sum of cases generated by infected individuals at different
infectious stages and different ages. We parametrize our model using
demographic and immunological data from sub-Saharan regions. Numerical
bifurcation analysis demonstrates the existence of an endemic equilibrium, and
we observe a forward bifurcation in $\mathcal{R}_0$. Our numerical simulations
reproduce the heterogeneity in the age distributions of immunity profiles and
infection status created by frequent exposure. Motivated by the recently
approved RTS,S vaccine, we also study the impact of vaccination; our results
show a reduction in severe disease among young children but a small increase in
severe malaria among older children due to lower acquired immunity from delayed
exposure.


  The first arrivals among $N$ Brownian particles is ubiquitous in the life
sciences, as it often trigger cellular processes from the molecular level. We
study here the case where stochastic particles, which represent molecules,
proteins or molecules can switch between two states inside the non-negative
real line. The switching process is modeled as a two-state Markov chain and
particles can only escape in state 1. We estimate the fastest arrival time by
solving asymptotically the Fokker-Planck equations for three different initial
distributions: Dirac-delta, uniformly distributed and long-tail decay. The
derived formulas reveal that the fastest particle avoid switching when the
switching rates are much smaller than the diffusion time scale, but switches
twice when the diffusion is state 2 is much faster than in state 1. The present
results are compared to stochastic simulations revealing the range of validity
of the derived formulas.


  ATP synthase's intrinsic molecular electrostatic potential (MESP) adds
constructively to, and hence reinforces, the chemiosmotic voltage. This ATP
synthase voltage represents a new free energy term that appears to have been
overlooked. This term is at least roughly equal in order of magnitude and
opposite in sign to the energy needed to be dissipated as a Maxwell's demon
(Landauer principle).


  We use Direct Coupling Analysis (DCA) to determine epistatic interactions
between loci of variability of the SARS-CoV-2 virus, segmenting genomes by
month of sampling. We use full-length, high-quality genomes from the GISAID
repository up to October 2021, in total over 3,500,000 genomes. We find that
DCA terms are more stable over time than correlations, but nevertheless change
over time as mutations disappear from the global population or reach fixation.
Correlations are enriched for phylogenetic effects, and in particularly
statistical dependencies at short genomic distances, while DCA brings out links
at longer genomic distance. We discuss the validity of a DCA analysis under
these conditions in terms of a transient Quasi-Linkage Equilibrium state. We
identify putative epistatic interaction mutations involving loci in Spike.


  Combinatorics of biopolymer structures, especially enumeration of various RNA
secondary structures and protein contact maps, is of significant interest for
communities of both combinatorics and computational biology. However, most of
the previous combinatorial enumeration results for these structures are
presented in terms of generating functions, and few are explicit formulas. This
paper is mainly concerned with finding explicit enumeration formulas for a
particular class of biologically relevant structures, say, saturated 2-regular
simple stacks, whose configuration is related to protein folds in the 2D
honeycomb lattice. We establish a semi-bijective algorithm that converts
saturated 2-regular simple stacks into forests of small trees, which produces a
uniform formula for saturated extended 2-regular simple stacks with any of the
six primary component types. Summarizing the six different primary component
types, we obtain a bivariate explicit formula for saturated extended 2-regular
simple stacks with $n$ vertices and $k$ arcs. As consequences, the uniform
formula can be reduced to Clote's results on $k$-saturated 2-regular simple
stacks and the optimal 2-regular simple stacks, and Guo et al.'s result on the
optimal extended 2-regular simple stacks.


  In this paper we apply the structure of genomes as second-order Markov
processes specified by the distributions of successive triplets of bases to two
bioinformatics problems: identification of outliers in genome databases and
read classification in metagenomics, using real coronavirus and adenovirus
data.


  Radiation therapy has remained as one of the main cancer treatment modalities
and a highly cost-effective single modality treatment of cancer care. Typical
regimens for fractionated external beam radiotherapy comprise a constant dose
administered on weekdays, and no radiation on weekends. However, every patient
has a tumor with distinct properties depending on intra-tumor heterogeneity,
aggressiveness, and interactive properties with other cells that may make it
more resistant or sensitive to radiation treatment. Accordingly, the concept of
personalized cancer treatment is emerging to specialize each patient treatment
case to the unique properties of the tumor. In this paper, we examine adaptive
radiation treatment strategies for heterogeneous tumors using a dynamical
system model that consists of radiation-resistant and parental cell populations
with unique interactive properties. We study different adaptive dosage
strategies for PC3 and DU145 prostate cancer cell lines. We show that stronger
doses of radiation given in longer time intervals, while keeping the overall
dosage the same, reduce final tumor volume by more than half in PC3 cell lines,
but by only five percent in DU145 cell lines. In addition, we tested an
adaptive dosing schedule by administering a stronger dosage on Friday to
compensate for the treatment-off period during the weekend, which was effective
in decreasing the final tumor volume of both cell lines. This result creates
interesting possibilities for new radiotherapy strategies at clinics that
cannot provide treatment on weekends. Finally, we propose a dosage plan
incorporating our findings.


  Breast cancer is the most common cancer in women worldwide, and discovering
the biomarkers of this disease became so vital nowadays and Cell in Cell
structure could be one of them, and it may be used as an available proxy for
tumor malignancy. (CICs) are unusual in that keep morphologically healthy cells
within another cell. They are found in various human cancers and result from
active cell-cell interaction, and it has different kinds. In this study, we
analyzed the microarray data from GEO (GSE103865) to genetically evaluate CICs'
incidence in samples obtained from breast cancer patients to understand the
relationship between the rate of CIC and the prognosis of breast cancer. The
preprocessing was performed using R software. The DAVID website was used to
analyze gene ontology (GO) and Gene and Genome (KEGG) pathways. The
protein-protein interactions (PPIs) of the obtained DEGs were assessed using
the STRING website, and hub modules in Cytoscape and cytoHubba were screened.
According to the results from analyzing the 20 hub genes, we understood that
overexpression of our Top genes is effective in focal adhesion, ECM-receptor
interaction, platelet activation and PI3K-Akt signaling pathway, which shows
that changes in these pathways could be the reason the overexpression of CICs
in breast cancer. These data and research by many others have uncovered various
genes involved in CIC formation and have started to give us an idea of why they
are formed and how they could contribute to breast cancer


  Background: Bladder cancer is the 10th most common cancer worldwide, and its
prevalence is increasing, especially in developing countries. Objective: In the
present study, we employed gene expression profiles from the GSE163209 data set
in the GEO database to identify potential molecular and genetic markers in BC
patients. Methods: The data set comprised 217 samples, with 113 stage Ta tumor
tissue samples and 104 stage T1 tumor tissue samples. The top 766 genes were
chosen. P.value<0.0001 and |logFC|=1 was used to change the cutoff criteria for
defining DEGs. Moreover, the MCODE plugin and cytoHubba plugin were employed to
produce a module and detect 20 hub genes in these DEGs. We used GO and KEGG
pathway enrichment analyses to get a better understanding of these DEGs.
Results: The KEGG pathway enrichment results indicated that the top genes were
mainly involved: Systemic lupus erythematosus, Alcoholism, and Viral
carcinogenesis. SLE activation in the renal glomeruli could explain the
connection between this disease's route and bladder cancer, and according to
our results and previous researches, heavy alcohol intake can increase the risk
of BC in males and particular populations. Conclusion: According to our hub
genes, we can consider CDK1 and NDC80 as bladder cancer biomarkers. Not much
research has been done on the effect of this gene on bladder cancer.


  Large genomic datasets are now created through numerous activities, including
recreational genealogical investigations, biomedical research, and clinical
care. At the same time, genomic data has become valuable for reuse beyond their
initial point of collection, but privacy concerns often hinder access. Over the
past several years, Beacon services have emerged to broaden accessibility to
such data. These services enable users to query for the presence of a
particular minor allele in a private dataset, information that can help care
providers determine if genomic variation is spurious or has some known clinical
indication. However, various studies have shown that even this limited access
model can leak if individuals are members in the underlying dataset. Several
approaches for mitigating this vulnerability have been proposed, but they are
limited in that they 1) typically rely on heuristics and 2) offer probabilistic
privacy guarantees, but neglect utility. In this paper, we present a novel
algorithmic framework to ensure privacy in a Beacon service setting with a
minimal number of query response flips (e.g., changing a positive response to a
negative). Specifically, we represent this problem as combinatorial
optimization in both the batch setting (where queries arrive all at once), as
well as the online setting (where queries arrive sequentially). The former
setting has been the primary focus in prior literature, whereas real Beacons
allow sequential queries, motivating the latter investigation. We present
principled algorithms in this framework with both privacy and, in some cases,
worst-case utility guarantees. Moreover, through an extensive experimental
evaluation, we show that the proposed approaches significantly outperform the
state of the art in terms of privacy and utility.


  We offer a general theoretical framework for brain and behavior that is
evolutionarily and computationally plausible. The brain in our abstract model
is a network of nodes and edges. Although it has some similarities to standard
neural network models, as we show, there are some significant differences. Both
nodes and edges in our network have weights and activation levels. They act as
probabilistic transducers that use a set of relatively simple rules to
determine how activation levels and weights are affected by input, generate
output, and affect each other. We show that these simple rules enable a
learning process that allows the network to represent increasingly complex
knowledge, and simultaneously to act as a computing device that facilitates
planning, decision-making, and the execution of behavior. By specifying the
innate (genetic) components of the network, we show how evolution could endow
the network with initial adaptive rules and goals that are then enriched
through learning. We demonstrate how the developing structure of the network
(which determines what the brain can do and how well) is critically affected by
the co-evolved coordination between the mechanisms affecting the distribution
of data input and those determining the learning parameters (used in the
programs run by nodes and edges). Finally, we consider how the model accounts
for various findings in the field of learning and decision making, how it can
address some challenging problems in mind and behavior, such as those related
to setting goals and self-control, and how it can help understand some
cognitive disorders.


  Spatial models where growth is limited to the edge of the expansions have
been instrumental to understand the population dynamics and the clone size
distribution in growing cellular populations, such as microbial colonies and
avascular tumours. A complete characterization of the coalescence process
generated by spatial growth is still lacking, limiting our ability to apply
classic population genetics inference to spatially growing populations. Here,
we start filling this gap by investigating the statistical properties of the
cell lineages generated by the two dimensional Eden model, leveraging their
physical analogy with directed polymers. Our analysis provides quantitative
estimates for population measurements that can easily be assessed via
sequencing, such as the average number of segregating sites and the clone size
distribution of a subsample of the population. Our results not only reveal
remarkable features of the genealogies generated during growth, but also
highlight new properties that can be misinterpreted as signs of selection if
non-spatial models are inappropriately applied.


  Peptide T is a synthetic octapeptide fragment, which corresponds to the
region 185-192 of the gp120 HIV coat protein and functions as a viral entry
inhibitor. In this work, a folding molecular dynamics simulation of peptide T
in a membrane-mimicking (DMSO) solution was performed with the aim of
characterizing the peptide's structural and dynamical properties. We show that
peptide T is highly flexible and dynamic. The main structural characteristics
observed were rapidly interconverting short helical stretches and turns, with a
notable preference for the formation of $\beta$-turns. The simulation also
indicated that the C-terminal part appears to be more stable than the rest of
the peptide, with the most preferred conformation for residues 5-8 being a
$\beta$-turn. In order to validate the accuracy of the simulations, we compared
our results with the experimental NMR data obtained for the T-peptide in the
same solvent. In agreement with the simulation, the NMR data indicated the
presence of a preferred structure in solution that was consistent with a
$\beta$-turn comprising the four C-terminal residues. An additional comparison
between the experimental and simulation-derived chemical shifts also showed a
reasonable agreement between experiment and simulation, further validating the
simulation-derived structural characterization of the T-peptide.


  Viruses constantly undergo mutations with genomic changes. The propagation of
variants of viruses is an interesting problem. We perform numerical simulations
of the microscopic epidemic model based on network theory for the spread of
variants. Assume that a small number of individuals infected with the variant
are added to widespread infection with the original virus. When a highly
infectious variant that is more transmissible than the original lineage is
added, the variant spreads quickly to the wide space. On the other hand, if the
infectivity is about the same as that of the original virus, the infection will
not spread. The rate of spread is not linear as a function of the infection
strength but increases non-linearly. This cannot be explained by the
compartmental model of epidemiology but can be understood in terms of the
dynamic absorbing state known from the contact process.


  Beyond identifying genetic variants, we introduce a set of Boolean relations
that allows for a comprehensive classification of the relations for every pair
of variants by taking all minimal alignments into account. We present an
efficient algorithm to compute these relations, including a novel way of
efficiently computing all minimal alignments within the best theoretical
complexity bounds. We show that for variants of the CFTR gene in dbSNP these
relations are common and many non-trivial. Ultimately, we present an approach
for the storing and indexing of variants in the context of a database that
enables efficient querying for all these relations.


  Humankind has spread worldwide supported by cultural and technological
knowledge, but the environmental sustainability on the human niche evolution
depends on a new human beings relationship with the biosphere. Human lifestyles
nowadays are very Antropocentric and in many ways deleterious to the other life
forms. Here we try to identify future scenarios, where the less deleterious is
the Natural-Technological Model that points the urgent need to change the
evolutionary direction of the human niche seeking the resumption of original
ecological relations. New cultural habits and novel technologies, thereby,
would reverse the current anthropogenic impacts. The middle way is the
Bio-Anthropogenic Model that predicts the success of the emerging ecosystems
and the symbiotic relationship of humans and anthropogenic-favored species,
hybrids, aliens and genetically modified organisms. For such, we must also
change our way of life and adopt new conscious ways of consumption aiming at
the socio-environmental good. Lastly, the Wear Out Model, which depends only on
maintaining current patterns of human expansion. The lack of investments on new
technologies and new cultural habits, added to the current patterns of human
niche evolution that are based on the massive exploitation of world resources,
will lead to a fearsome scenario with a precarious global health, biodiversity
losses and food scarcity. This theoretical models indicates some pathways and
can help us to choose a better future.


  Natural changes have been played a vital part within the consider of
irresistible illnesses uncommonly coronavirus. It is important to think about
the effect of irregular unsettling influence on epidemic models. In this work,
first we present a SIR coronavirus mathematical model and consider that the
contact rate is perturbed by L\'evy noise. To protect and control diseases, the
L\'evy process should be used. Then, we discuss the dynamics of the
deterministic model and for the global positive solution of the stochastic
model, present the existence and uniqueness. Further, explore a few conditions
for the termination and persistence of the infection. Then derive the basic
reproduction number that determines the extinction and the persistence of the
disease (infection). When the noise is huge or little, the numerical results
appear that the COVID19 vanishes from the individuals on the chance that the
reproduction number is less than one; though control the epidemic illnesses in
the event that reproduction number is greater than one. Finally, for the
demonstration of this phenomenon numerical simulations are presented.


  Molecular motors belonging to the kinesin and myosin super family hydrolyze
ATP by cycling through a sequence of chemical states. These cytoplasmic motors
are dimers made up of two linked identical monomeric globular proteins. Fueled
by the free energy generated by ATP hydrolysis, the motors walk on polar tracks
(microtubule or filamentous actin) processively, which means that only one head
detaches and executes a mechanical step while the other stays bound to the
track. Thus, the one motor head must regulate chemical state of the other,
referred to as "gating", a concept that is not fully understood. Inspired by
experiments, showing that only a fraction of the energy from ATP hydrolysis is
used to advance the kinesin motors against load, we demonstrate that additional
energy is used for coordinating the chemical cycles of the two heads in the
dimer - a feature that characterizes gating. To this end, we develop a general
framework based on information theory and stochastic thermodynamics, and
establish that gating could be quantified in terms of information flow between
the motor heads. Applications of the theory to kinesin-1 and Myosin V show that
information flow occurs, with positive cooperativity, at external resistive
loads that are less than a critical value, $F_c$. When force exceeds $F_c$,
effective information flow ceases. Interestingly, $F_c$, which is independent
of the input energy generated through ATP hydrolysis, coincides with force at
which the probability of backward steps starts to increase. Our findings
suggest that transport efficiency is optimal only at forces less than $F_c$,
which implies that these motors must operate at low loads under $\textit{in
vivo}$ conditions.


  Motivation: Cutting the cost of DNA sequencing technology led to a quantum
leap in the availability of genomic data. While sharing genomic data across
researchers is an essential driver of advances in health and biomedical
research, the sharing process is often infeasible due to data privacy concerns.
Differential privacy is one of the rigorous mechanisms utilized to facilitate
the sharing of aggregate statistics from genomic datasets without disclosing
any private individual-level data. However, differential privacy can still
divulge sensitive information about the dataset participants due to the
correlation between dataset tuples. Results: Here, we propose GenShare model
built upon Laplace-perturbation-mechanism-based DP to introduce a
privacy-preserving query-answering sharing model for statistical genomic
datasets that include dependency due to the inherent correlations between
genomes of individuals (i.e., family ties). We demonstrate our privacy
improvement over the state-of-the-art approaches for a range of practical
queries including cohort discovery, minor allele frequency, and chi^2
association tests. With a fine-grained analysis of sensitivity in the Laplace
perturbation mechanism and considering joint distributions, GenShare results
near-achieve the formal privacy guarantees permitted by the theory of
differential privacy as the queries that computed over independent tuples (only
up to 6% differences). GenShare ensures that query results are as accurate as
theoretically guaranteed by differential privacy. For empowering the advances
in different scientific and medical research areas, GenShare presents a path
toward an interactive genomic data sharing system when the datasets include
participants with familial relationships.


  Invariants for complicated objects such as those arising in phylogenetics,
whether they are invariants as matrices, polynomials, or other mathematical
structures, are important tools for distinguishing and working with such
objects. In this paper, we generalize a complete polynomial invariant on trees
to a class of phylogenetic networks called separable networks, which will
include orchard networks. Networks are becoming increasingly important for
their ability to represent reticulation events, such as hybridization, in
evolutionary history. We provide a function from the space of internally
multi-labelled phylogenetic networks, a more generic graph structure than
phylogenetic networks where the reticulations are also labelled, to a
polynomial ring. We prove that the separability condition allows us to
characterize, via the polynomial, the phylogenetic networks with the same
number of leaves and same number of reticulations by considering their
internally labelled versions. While the invariant for trees is a polynomial in
Z[x_1,..., x_n,y] where n is the number of leaves, the invariant for internally
multi-labelled phylogenetic networks is an element of Z[x_1,...,
x_n,lambda_1,...,lambda_r,y], where r is the number of reticulations in the
network. When the networks are considered without leaf labels the number of
variables reduces to r+2.


  Motivated by recent findings of enhanced species survival when fragmented
habitats are reconnected through narrow strips of land [S. Pimm, and C. N.
Jenkins, Am. Sci. {\bf 107}(3), 162 (2019).], we study the effect of a corridor
connecting two square regions on the survival time of three models exhibiting
extinction/survival phase transitions: the basic contact process (CP), the
diffusive contact process, and the two-species symbiotic contact process
(2SCP). We find that connecting the regions generally increases the survival
time for $\lambda\ge \lambda_c$, where $\lambda$ is the reproduction rate and
$\lambda_c$ its critical value. The enhancement of the survival time increases
with $\lambda$, and is largest in the 2SCP.


  The COVID-19 pandemic has placed unprecedented demands on local public health
units in Ontario, Canada, one of which was the need for in-house
epidemiological-modelling capabilities. To address this need, Ontario Tech
University and the Durham Region Health Department developed a native Windows
desktop app that performs epidemiological modelling of infectious diseases. The
app is an implementation of a multi-stratified compartmental epidemiological
model that can accommodate multiple virus variants and levels of vaccination,
as well as public health measures such as physical distancing, contact tracing
followed by quarantine, and testing followed by isolation. This article
presents the epidemiological model and epidemiological-simulation results
obtained using the developed app. The simulations investigate the effects of
different factors on COVID-19 transmission in Durham Region, including
vaccination coverage, vaccine effectiveness, waning of vaccine-induced
immunity, advent of the Omicron variant and effect of COVID-19 booster vaccines
in reducing the number of infections and severe cases. Results indicate that,
for the Delta variant, natural immunity, in addition to vaccination-induced
immunity, is necessary to achieve herd immunity and that waning of
vaccine-induced immunity lengthens the time necessary to reach herd immunity.
In the absence of additional public health measures, a wave driven by the
Omicron variant is predicted to pose significant public health challenges with
infections predicted to peak in approximately two to three months, depending on
the rate of administration of booster doses.


  A phylogenetic network is a simplex (or 1-component tree-child) network if
the child of every reticulation node is a network leaf. Simplex networks are a
superclass of phylogenetic trees and a subclass of tree-child networks.
Generalizing the Sackin index to phylogenetic networks, we prove that the
expected Sackin index of a random simplex network is asymptotically
$\Omega(n^{7/4})$ in the uniform model.


  This report is based on the work in (1). We first review definitions and
notation developped there and provide derivations for the exact mathematical
description of an SIS epidemic on a hypergraph. We then generalise the work in
(1) to a new class of models that encompass SIS and SIR models. The exact
differential equations are derived for the expected values of the population of
each state. Focusing on Bi-uniform hypergraphs, we make suitable approximations
obtain numerical solutions to those equations. These are compared with
stochastic simulations of the model for various systems.


  Cocaine addiction is a psychosocial disorder induced by the chronic use of
cocaine and causes a large of number deaths around the world. Despite many
decades' effort, no drugs have been approved by the Food and Drug
Administration (FDA) for the treatment of cocaine dependence. Cocaine
dependence is neurological and involves many interacting proteins in the
interactome. Among them, dopamine transporter (DAT), serotonin transporter
(SERT), and norepinephrine transporter (NET) are three major targets. Each of
these targets has a large protein-protein interaction (PPI) network which must
be considered in the anti-cocaine addiction drug discovery. This work presents
DAT, SERT, and NET interactome network-informed machine learning/deep learning
(ML/DL) studies of cocaine addiction. We collect and analyze 61 protein targets
out 460 proteins in the DAT, SERT, and NET PPI networks that have sufficient
existing inhibitor datasets. Utilizing autoencoder and other ML algorithms, we
build ML/DL models for these targets with 115,407 inhibitors to predict drug
repurposing potentials and possible side effects. We further screen their
absorption, distribution, metabolism, and excretion, and toxicity (ADMET)
properties to search for nearly optimal leads for anti-cocaine addiction. Our
approach sets up a systematic protocol for artificial intelligence (AI)-based
anti-cocaine addiction lead discovery.


  Phages are one of the most present groups of organisms in the biosphere.
Their identification continues and their taxonomies are divergent. However, due
to their evolution mode and the complexity of their species ecosystem, their
classification is not complete. Here, we present a new approach to the phages
classification that combines the methods of horizontal gene transfer detection
and ancestral sequence reconstruction.


  Phylogenetic comparative methods are new in our field and are shrouded, for
most linguists, in at least a little mystery. Yet the path that led to their
discovery in comparative biology is so similar to the methodological history of
balanced sampling, that it is only an accident of history that they were not
discovered by a typologist. Here we clarify the essential logic behind
phylogenetic comparative methods and their fundamental relatedness to a deep
intellectual tradition focussed on sampling. Then we introduce concepts,
methods and tools which will enable typologists to use these methods in
everyday typological research. The key commonality of phylogenetic comparative
methods and balanced sampling is that they attempt to deal with statistical
non-independence due to genealogy. Whereas sampling can never achieve
independence and requires most comparative data to be discarded, phylogenetic
comparative methods achieve independence while retaining and using all data. We
discuss the essential notions of phylogenetic signal; uncertainty about trees;
typological averages and proportions that are sensitive to genealogy;
comparison across language families; and the effects of areality. Extensive
supplementary materials illustrate computational tools for practical analysis
and we illustrate the methods discussed with a typological case study of the
laminal contrast in Pama-Nyungan.


  Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative
agent of the current global COVID-19 pandemic, in which millions of lives have
been lost. Understanding the zoonotic evolution of the coronavirus may provide
insights for developing effective vaccines, monitoring the transmission trends,
and preventing new zoonotic infections. Homopolymeric nucleotide repeats (HP),
the most simple tandem repeats, are a ubiquitous feature of eukaryotic genomes.
Yet the HP distributions and roles in coronavirus genome evolution are poorly
investigated. In this study, we characterize the HP distributions and trends in
the genomes of bat and human coronaviruses and SARS-CoV-2 variants. The results
show that the SARS-CoV-2 genome is abundant in HPs, and has augmented HP
contents during evolution. Especially, the disparity of HP poly-(A/T) and
ploy-(C/G) of coronaviruses increases during the evolution in human hosts. The
disparity of HP poly-(A/T) and ploy-(C/G) is correlated to host adaptation and
the virulence level of the coronaviruses. Therefore, we propose that the HP
disparity can be a quantitative measure for the zoonotic evolution levels of
coronaviruses. Peculiarly, the HP disparity measure infers that SARS-CoV-2
Omicron variants have a high disparity of HP poly-(A/T) and ploy-(C/G),
suggesting a high adaption to the human hosts.


  A molecule of dimension 125nm has caused around 479 Million human infections
(80M for the USA) & 6.1 Million human deaths (977,000 for the USA) worldwide
and slashed the global economy by US$ 8.5 Trillion over two years. The only
other events in recent history that caused comparative human life loss through
direct usage (either by (wo)man or nature, respectively) of structure-property
relations of 'nano-structures' (either (wo)man-made or nature, respectively)
were nuclear bomb attacks of Japanese cities by the USA during World War II and
1918 Flu Pandemic. This molecule is SARS-CoV-2, which causes a disease known as
COVID-19. The high liability cost of the pandemic had incentivized various
private, government, and academic entities to work towards finding a cure for
these & emerging diseases. As result, multiple vaccine candidates are
discovered to avoid the infection in first place. But so far, there has been no
success in finding fully effective therapeutics candidates. In this paper, we
attempted to provide multiple therapy candidates based upon a sophisticated
multi-scale in-silico framework. We have used the following robust framework to
screen the ligands; Step-I: high throughput docking, Step-II: molecular
dynamics, Step-III: density functional theory analysis. In total, we have
analyzed 2.2 Million unique protein binding site/ligand combinations. The
proteins were selected based on recent experimental studies. Step-I had
filtered that number down to 10 ligands/protein based on molecular docking
binding energy, further screening down to 2 ligands/protein based on
drug-likeness analysis. Additionally, these two ligands/proteins were
investigated in Step-II with a molecular dynamic based RMSD analysis. It
finally suggested three ligands (ZINC1176619532, ZINC517580540, ZINC952855827)
attacking different binding sites of the protein(7BV2), which were further
analyzed in Step III.


  Carbon isotope labeling method is a standard metabolic engineering tool for
flux quantification in living cells. To cope with the high dimensionality of
isotope labeling systems, diverse algorithms have been developed to reduce the
number of variables or operations in metabolic flux analysis (MFA), but lacks
generalizability to non-stationary metabolic conditions. In this study, we
present a stochastic simulation algorithm (SSA) derived from the chemical
master equation of the isotope labeling system. This algorithm allows to
compute the time evolution of isotopomer concentrations in non-stationary
conditions, with the valuable property that computational time does not scale
with the number of isotopomers. The efficiency and limitations of the algorithm
is benchmarked for the forward and inverse problems of 13C-DMFA in the pentose
phosphate pathways. Overall, SSA constitute an alternative class to
deterministic approaches for metabolic flux analysis that is well adapted to
comprehensive dataset including parallel labeling experiments, and whose
limitations associated to the sampling size can be overcome by using Monte
Carlo sampling approaches.


  The effect of bias on hypothesis formation is characterized for an automated
data-driven projection pursuit neural network to extract and select features
for binary classification of data streams. This intelligent exploratory process
partitions a complete vector state space into disjoint subspaces to create
working hypotheses quantified by similarities and differences observed between
two groups of labeled data streams. Data streams are typically time sequenced,
and may exhibit complex spatio-temporal patterns. For example, given atomic
trajectories from molecular dynamics simulation, the machine's task is to
quantify dynamical mechanisms that promote function by comparing protein
mutants, some known to function while others are nonfunctional. Utilizing
synthetic two-dimensional molecules that mimic the dynamics of functional and
nonfunctional proteins, biases are identified and controlled in both the
machine learning model and selected training data under different contexts. The
refinement of a working hypothesis converges to a statistically robust
multivariate perception of the data based on a context-dependent perspective.
Including diverse perspectives during data exploration enhances
interpretability of the multivariate characterization of similarities and
differences.


  Prokaryotic viruses, which infect bacteria and archaea, are key players in
microbial communities. Predicting the hosts of prokaryotic viruses helps
decipher the dynamic relationship between microbes. Experimental methods for
host prediction cannot keep pace with the fast accumulation of sequenced
phages. Thus, there is a need for computational host prediction. Despite some
promising results, computational host prediction remains a challenge because of
the limited known interactions and the sheer amount of sequenced phages by
high-throughput sequencing technologies. The state-of-the-art methods can only
achieve 43\% accuracy at the species level. In this work, we formulate host
prediction as link prediction in a knowledge graph that integrates multiple
protein and DNA-based sequence features. Our implementation named CHERRY can be
applied to predict hosts for newly discovered viruses and to identify viruses
infecting targeted bacteria. We demonstrated the utility of CHERRY for both
applications and compared its performance with 11 popular host prediction
methods. To our best knowledge, CHERRY has the highest accuracy in identifying
virus-prokaryote interactions. It outperforms all the existing methods at the
species level with an accuracy increase of 37\%. In addition, CHERRY's
performance on short contigs is more stable than other tools.


  Microbial biomass carbon (MBC), a crucial soil labile carbon fraction, is the
most active component of the soil organic carbon (SOC) that regulates
bio-geochemical processes in terrestrial ecosystems. Some studies in the
literature ignore the effect of microbial population growth on carbon
decomposition rates. In reality, we might expect that the decomposition rate
should be related to the population of microbes in the soil and have a positive
relationship with the size of the microbial biomass pool. In this study, we
explore the effect of microbial population growth on the accuracy of modelling
soil carbon sequestration by developing and comparing two soil carbon models
that consider a carrying capacity and limit to the growth of the microbial
pool. We apply our models to three datasets, two small and one large datasets,
and we select the best model in terms of having the best predictive performance
through two model selection methods. Through this analysis we reveal that
commonly used complex soil carbon models can over-fit in the presence of both
small and large time-series datasets, and our simpler model can produce more
accurate predictions. We conclude that considering the microbial population
growth in a soil carbon model improves the accuracy of a model in the presence
of a large dataset.


  Reaction networks (RNs) comprise a set $X$ of species and a set $\mathscr{R}$
of reactions $Y\to Y'$, each converting a multiset of educts $Y\subseteq X$
into a multiset $Y'\subseteq X$ of products. RNs are equivalent to directed
hypergraphs. However, not all RNs necessarily admit a chemical interpretation.
Instead, they might contradict fundamental principles of physics such as the
conservation of energy and mass or the reversibility of chemical reactions. The
consequences of these necessary conditions for the stoichiometric matrix
$\mathbf{S} \in \mathbb{R}^{X\times\mathscr{R}}$ have been discussed
extensively in the literature. Here, we provide sufficient conditions for
$\mathbf{S}$ that guarantee the interpretation of RNs in terms of balanced sum
formulas and structural formulas, respectively.
  Chemically plausible RNs allow neither a perpetuum mobile, i.e., a "futile
cycle" of reactions with non-vanishing energy production, nor the creation or
annihilation of mass. Such RNs are said to be thermodynamically sound and
conservative. For finite RNs, both conditions can be expressed equivalently as
properties of $\mathbf{S}$. The first condition is vacuous for reversible
networks, but it excludes irreversible futile cycles and - in a stricter sense
- futile cycles that even contain an irreversible reaction. The second
condition is equivalent to the existence of a strictly positive reaction
invariant. Furthermore, it is sufficient for the existence of a realization in
terms of sum formulas, obeying conservation of "atoms". In particular, these
realizations can be chosen such that any two species have distinct sum
formulas, unless $\mathbf{S}$ implies that they are "obligatory isomers". In
terms of structural formulas, every compound is a labeled multigraph, in
essence a Lewis formula, and reactions comprise only a rearrangement of bonds
such that the total bond order is preserved.


  Biomolecular condensates, physically underpinned to a significant extent by
liquid-liquid phase separation (LLPS), are now widely recognized by numerous
experimental studies to be of fundamental biological, biomedical, and
biophysical importance. In the face of experimental discoveries, analytical
formulations emerged as a powerful yet tractable tool in recent theoretical
investigations of the role of LLPS in the assembly and dissociation of these
condensates. The pertinent LLPS often involves, though not exclusively,
intrinsically disordered proteins engaging in multivalent interactions that are
governed by their amino acid sequences. For researchers interested in applying
these theoretical methods, here we provide a practical guide to a set of
computational techniques devised for extracting sequence-dependent LLPS
properties from analytical formulations. The numerical procedures covered
include those for the determinination of spinodal and binodal phase boundaries
from a general free energy function with examples based on the random phase
approximation in polymer theory, construction of tie lines for
multiple-component LLPS, and field-theoretic simulation of multiple-chain
heteropolymeric systems using complex Langevin dynamics. Since a more accurate
physical picture often requires comparing analytical theory against
explicit-chain model predictions, a commonly utilized methodology for
coarse-grained molecular dynamics simulations of sequence-specific LLPS is also
briefly outlined.


  The Moran process is a classic stochastic process that models invasion
dynamics on graphs. A single "mutant" (e.g., a new opinion, strain, social
trait etc.) invades a population of residents spread over the nodes of a graph.
The mutant fitness advantage $\delta\geq 0$ determines how aggressively mutants
propagate to their neighbors. The quantity of interest is the fixation
probability, i.e., the probability that the initial mutant eventually takes
over the whole population. However, in realistic settings, the invading mutant
has an advantage only in certain locations. E.g., a bacterial mutation allowing
for lactose metabolism only confers an advantage on places where dairy products
are present. In this paper we introduce the positional Moran process, a natural
generalization in which the mutant fitness advantage is only realized on
specific nodes called active nodes. The associated optimization problem is
fixation maximization: given a budget $k$, choose a set of $k$ active nodes
that maximize the fixation probability of the invading mutant. We show that the
problem is NP-hard, while the optimization function is not submodular, thus
indicating strong computational hardness. Then we focus on two natural limits.
In the limit of $\delta\to\infty$ (strong selection), although the problem
remains NP-hard, the optimization function becomes submodular and thus admits a
constant-factor approximation using a simple greedy algorithm. In the limit of
$\delta\to 0$ (weak selection), we show that in $O(m^\omega)$ time we can
obtain a tight approximation, where $m$ is the number of edges and $\omega$ is
the matrix-multiplication exponent. Finally, we present an experimental
evaluation of the new algorithms together with some proposed heuristics.


  Certain defense mechanisms of phages against the immune system of their
bacterial host rely on cooperation of phages. Motivated by this example we
analyse invasion probabilities of cooperative parasites in host populations
that are moderately structured. More precisely we assume that hosts are
arranged on the vertices of a configuration model and that offspring of
parasites move to nearest neighbours sites to infect new hosts. We consider
parasites that generate many offspring at reproduction, but do this (usually)
only when infecting a host simultaneously. In this regime we identify and
analyse the spatial scale of the population structure at which invasion of
parasites turns from being an unlikely to an highly probable event.


  COVID-19 pandemic, is still unknown and is an important open question. There
are speculations that bats are a possible origin. Likewise, there are many
closely related (corona-) viruses, such as SARS, which was found to be
transmitted through civets. The study of the different hosts which can be
potential carriers and transmitters of deadly viruses to humans is crucial to
understanding, mitigating and preventing current and future pandemics. In
coronaviruses, the surface (S) protein, or spike protein, is an important part
of determining host specificity since it is the point of contact between the
virus and the host cell membrane. In this paper, we classify the hosts of over
five thousand coronaviruses from their spike protein sequences, segregating
them into clusters of distinct hosts among avians, bats, camels, swines, humans
and weasels, to name a few. We propose a feature embedding based on the
well-known position-weight matrix (PWM), which we call PWM2Vec, and use to
generate feature vectors from the spike protein sequences of these
coronaviruses. While our embedding is inspired by the success of PWMs in
biological applications such as determining protein function, or identifying
transcription factor binding sites, we are the first (to the best of our
knowledge) to use PWMs in the context of host classification from viral
sequences to generate a fixed-length feature vector representation. The results
on the real world data show that in using PWM2Vec, we are able to perform
comparably well as compared to baseline models. We also measure the importance
of different amino acids using information gain to show the amino acids which
are important for predicting the host of a given coronavirus.


  Deep autoencoders are often extended with a supervised or adversarial loss to
learn latent representations with desirable properties, such as greater
predictivity of labels and outcomes or fairness with respects to a sensitive
variable. Despite the ubiquity of supervised and adversarial deep latent factor
models, these methods should demonstrate improvement over simpler linear
approaches to be preferred in practice. This necessitates a reproducible linear
analog that still adheres to an augmenting supervised or adversarial objective.
We address this methodological gap by presenting methods that augment the
principal component analysis (PCA) objective with either a supervised or an
adversarial objective and provide analytic and reproducible solutions. We
implement these methods in an open-source Python package, AugmentedPCA, that
can produce excellent real-world baselines. We demonstrate the utility of these
factor models on an open-source, RNA-seq cancer gene expression dataset,
showing that augmenting with a supervised objective results in improved
downstream classification performance, produces principal components with
greater class fidelity, and facilitates identification of genes aligned with
the principal axes of data variance with implications to development of
specific types of cancer.


  Genetic data are now routinely used to study the history of population size,
subdivision, and gene flow. A variety of formal statistical methods is
available for testing hypotheses and fitting models to data. Yet it is often
unclear which hypotheses are worth testing, which models worth fitting. There
is a need for less formal methods that can be used in exploratory analysis of
genetic data. One approach to this problem uses *nucleotide site patterns*,
which provide a simple summary of the pattern in genetic data. This article
shows how to use them in exploratory data analysis.


  In this paper we model the spreading of the SARS-CoV-2 in Mexico by
introducing a new stochastic approximation constructed from first principles,
structured on the basis of a Latent-Infectious- (Recovered or Deceased)
(LI(RD)) compartmental approximation, where the number of new infected
individuals caused by a single infectious individual per unit time (a day), is
a random variable of a Poisson distribution and whose parameter is modulated
through a weight-like time-dependent function. The weight function serves to
introduce a time dependence to the average number of new infections and as we
will show, this information can be extracted from empirical data, giving to the
model self-consistency and provides a tool to study information about periodic
patterns encoded in the epidemiological dynamics


  Effective strategies of vaccine prioritization are essential to mitigate the
impacts of severe infectious diseases. We investigate the role of infection
fatality ratio (IFR) and social contact matrices on vaccination prioritization
using a compartmental epidemic model fueled by real-world data of different
diseases and countries. Our study confirms that massive and early vaccination
is extremely effective to reduce the disease fatality if the contagion is
mitigated, but the effectiveness is increasingly reduced as vaccination
beginning delays in an uncontrolled epidemiological scenario. The optimal and
least effective prioritization strategies depend non-linearly on
epidemiological variables. Regions of the epidemiological parameter space, in
which prioritizing the most vulnerable population is more effective than the
most contagious individuals, depend strongly on the IFR age profile being, for
example, substantially broader for COVID-19 in comparison with seasonal
influenza. Demographics and social contact matrices deform the phase diagrams
but do not alter their qualitative shapes.


  Some genes can promote or repress their own expressions, which is called
autoregulation. Although gene regulation is a central topic in biology,
autoregulation is much less studied. In general, it is extremely difficult to
determine the existence of autoregulation with direct biochemical approaches.
Nevertheless, some papers have observed that certain types of autoregulations
are linked to noise levels in gene expression. We generalize these results by
two propositions on discrete-state continuous-time Markov chains. These two
propositions form a simple but robust method to infer the existence of
autoregulation from gene expression data. This method only needs to compare the
mean and variance of the gene expression level. Compared to other methods for
inferring autoregulation, our method only requires non-interventional one-time
data, and does not need to estimate parameters. Besides, our method has few
restrictions on the model. We apply this method to four groups of experimental
data and find some genes that might have autoregulation. Some inferred
autoregulations have been verified by experiments or other theoretical works.


  In natural ecosystems, species often compete with multiple species to survive
and acquire resources. These species can be characterized by their different
growth profile and life history strategies, namely r and K selected species
(Williams et.al., The journal of wildlife management, 2013). In this paper, we
examine the role of different life-history strategies on the stability of the
competitive community along with network complexity. We obtain that the
r-selected species enhance community stability and resilience; on the other
hand, K-selected species show the opposite trend. Network connectance shows a
consistent trend of increasing stability; on the other side, the number of
species shows a context-dependent effect on stability. Stability decreases with
species number when the proportion of r-selected species is low or moderate and
shows precisely the opposite trend when the ratio of r-selected species is
high. We also show that our result is robust irrespective of network structure.


  Colorectal cancer (CRC) is the third most frequent malignant disease in the
world. In some countries with established screening programs, its incidence and
mortality have decreased, and survival has improved. AIMS: To obtain reliable
data about the epidemiology of CRC in Chile, we analyzed the trends in the last
ten years and the influence of observable factors on survival, including
explicit guarantees in CRC treatment access (GES program).
  Publicly available data published by the Health Ministry and National
Institute of Statistics were used. Data were obtained from registries of
mortality and hospital discharges, making follow-up of the individuals
possible. Crude and age-standardized incidence and mortality rates were
calculated, and individual survival was studied by constructing Kaplan-Meier
curves. Finally, a Cox statistical model was established to estimate the impact
of the observable factors.


  Nowcasting and forecasting of epidemic spreading rely on incidence series of
reported cases to derive the fundamental epidemiological parameters for a given
pathogen. Two relevant drawbacks for predictions are the unknown fractions of
undocumented cases and levels of nonpharmacological interventions, which span
highly heterogeneously across different places and times. We describe a simple
data-driven approach using a compartmental model including asymptomatic and
presymptomatic contagions that allows to estimate both the level of
undocumented infections and the value of effective reproductive number R t from
time series of reported cases, deaths, and epidemiological parameters. The
method was applied to epidemic series for COVID-19 across different
municipalities in Brazil allowing to estimate the heterogeneity level of
under-reporting across different places. The reproductive number derived within
the current framework is little sensitive to both diagnosis and infection rates
during the asymptomatic states. The methods described here can be extended to
more general cases if data is available and adapted to other epidemiological
approaches and surveillance data.


  We developed a simple mathematical model to describe criminality and the
justice system composed of the police investigation and court trial. The model
assessed two features of organized crime -- the cost-benefit analysis done by
the crime-susceptible to commit a crime and the whistleblowing of the law
offenders. The model was formulated considering the mass action law commonly
used in the disease propagation modelings, which can shed light on the model's
analysis. The crime-susceptible individuals analyze two opposing forces --
committing crime influenced by the law offenders not caught by police neither
imprisonment by the court trial (benefit of enjoying the corruption incoming),
and the refraction to commit crime influenced by those caught by police or
condemned by a court (cost of incarceration). Moreover, we assessed the dilemma
for those captured by police investigation to participate in the rewarding
whistleblowing program. The model was applied to analyze the "car wash
operation" against corruption in Brazil. The model analysis showed that the
cost-benefit analysis of crime-susceptible individuals whether the act of
bribery is worth or not determined the basic crime reproduction number
(threshold); however, the rewarding whistleblowing policies improved the combat
to corruption arising a sub-threshold. Some adopted mechanisms to control the
Covid-19 pandemic shed light on understanding the "car wash peration" and
threatens to the fight against corruption. Appropriate coverage of corruption
by media, enhancement of laws against white-collar crimes, well-functioning
police investigation and court trial, and the rewarding whistleblowing policies
inhibited and decreased the corruption.


  A wide variety of cultural practices take the form of "tacit" knowledge,
where the rules and principles are neither obvious to an observer nor known
explicitly by the practitioners. This poses a problem for cultural evolution:
if beginners cannot simply imitate experts, and experts cannot simply say or
demonstrate what they are doing, how can tacit knowledge pass from generation
to generation? We present a domain-general model of "tacit teaching", that
shows how high-fidelity transmission of tacit knowledge is possible. It applies
in cases where the underlying features of the practice are subject to
interacting and competing constraints, as is expected both in embodied and in
social practices. Our model makes predictions for key features of the teaching
process. It predicts a tell-tale distribution of teaching outcomes: some
students will be nearly perfect performers while others receiving the same
instruction will be disastrously bad. This differs from most mainstream
cultural evolution models centered on high-fidelity transmission with minimal
copying errors, which lead to a much narrower distribution where students are
mostly equally mediocre. The model also predicts generic features of the
cultural evolution of tacit knowledge. The evolution of tacit knowledge is
expected to be bursty, with long periods of stability interspersed with brief
periods of dramatic change, and where tacit knowledge, once lost, becomes
essentially impossible to recover.


  The capability of current force fields to reproduce RNA structural dynamics
is limited. Several methods have been developed to take advantage of
experimental data in order to enforce agreement with experiments. We herein
extend an existing framework, which allows arbitrarily chosen force-field
correction terms to be fitted by quantification of the discrepancy between
observables back-calculated from simulation and corresponding experiments. We
apply a robust regularization protocol to avoid overfitting, and additionally
introduce and compare a number of different regularization strategies, namely
L1-, L2-, Kish Size-, Relative Kish Size- and Relative Entropy-penalties. The
training set includes a GACC tetramer as well as more challenging systems,
namely gcGAGAgc and gcUUCGgc RNA tetraloops. Specific intramolecular hydrogen
bonds in the AMBER RNA force field are corrected with automatically determined
parameters that we call gHBfix$_{opt}$. A validation involving a separate
simulation of a system present in the training set (gcUUCGgc) and new systems
not seen during training (CAAU and UUUU tetramers) displays improvements
regarding native population of the tetraloop as well as good agreement with
NMR-experiments for tetramers when using the new parameters. Then we simulate
folded RNAs (a kink-turn and L1 stalk rRNA) including hydrogen bond types not
sufficiently present in the training set. This allows a final modification of
the parameter set which is named gHBfix21 and is suggested to be applicable to
a wider range of RNA systems.


  The relationship between interactions, flexibility and disorder in proteins
has been explored from many angles: folding upon binding, flexibility of the
core relative to the periphery, entropy changes, etc. In this work, we provide
statistical evidence for the involvement of highly mobile and disordered
regions in complex assembly. We ordered the entire set of X-ray
crystallographic structures in the PDB into hierarchies of progressive
interactions involving identical or very similar protein chains, yielding 40205
hierarchies of complexes with increasing numbers of partners. We then examine
them as proxies for the assembly pathways. Using this database, we show that
upon oligomerisation, new interfaces tend to be observed at residues that were
characterised as softly disordered (flexible, amorphous or missing residues) in
the complexes preceding them in the hierarchy. We also rule out the possibility
that this correlation is just a surface effect by restricting the analysis to
residues on the surface. Interestingly, we find that the location of soft
disordered residues in the sequence changes as the number of partners
increases. Our results show that there is a general mechanism for protein
assembly that involves soft disorder and modulates the way protein complexes
are assembled. This work highlights the difficulty of predicting the structure
of large protein complexes from sequence and emphasises the importance of
linking predictors of disorder to the next generation of predictors of complex
structure. Finally, we investigate the relationship between the AF2's
confidence metric pLDDT for structure prediction in unbound vs. bound
structures, and soft disorder. We show a strong correlation between AF2 low
confidence residues and the union of all regions of soft disorder observed in
the hierarchy. This paves the way for using the pLDDT metric as a proxy for
predicting assembly paths.


  During recent years, the study of long transients has been expanded in
ecological theory to account for shifts in long-term behavior of ecological
systems. These long transients may lead to regime shifts between alternative
states that resemble the dynamics of alternative stable states for a prolonged
period of time. One dynamic that potentially leads to long transients is the
group defense of a resource in a consumer-resource interaction. Furthermore,
time lags in the population caused by discrete reproductive pulses have the
potential to produce long transients, either independently or in conjunction to
the transients caused by the group defense. In this work, we analyze the
potential for long transients in a model for a consumer-resource system in
which the resource exhibits group defense and reproduces in discrete
reproductive pulses. We develop this discrete-time model by discretizing a
pulse differential equation. This system exhibits crawl-by transients near the
extinction and carrying capacity states of resource. In addition, we identify a
transcritical bifurcation in our system, under which a ghost limit cycle
appears. These transients resemble stable states for a prolonged transient time
period. We estimate the transient time of our system from these transients
using perturbation theory. This work advances an understanding of how systems
shift between alternate states and their duration of staying in a given regime
and what ecological dynamics may lead to long transients.


  Recent theory shows that extortioners taking advantage of the
zero-determinant (ZD) strategy can unilaterally claim an unfair share of the
payoffs in the Iterated Prisoner's Dilemma. It is thus suggested that against a
fixed extortioner, any adapting co-player should be subdued with full
cooperation as their best response. In contrast, recent experiments demonstrate
that human players often choose not to accede to extortion out of concern for
fairness, actually causing extortioners to suffer more loss than themselves. In
light of this, here we reveal fair-minded strategies that are unbending to
extortion such that any payoff-maximizing extortioner ultimately will concede
in their own interest by offering a fair split in head-to-head matches. We find
and characterize multiple general classes of such unbending strategies,
including generous zero-determinant strategies and Win-Stay, Lose-Shift as
particular examples. When against fixed unbending players, extortioners are
forced with consequentially increasing losses whenever intending to demand more
unfair share. Our analysis also pivots to the importance of payoff structure in
determining the superiority of zero-determinant strategies and in particular
their extortion ability. We show that an extortionate ZD player can be even
outperformed by, for example, Win-Stay Lose-Shift, if the total payoff of
unilateral cooperation is smaller than that of mutual defection. Unbending
strategies can be used to outlearn evolutionary extortioners and catalyze the
evolution of Tit-for-Tat-like strategies out of ZD players. Our work has
implications for promoting fairness and resisting extortion so as to uphold a
just and cooperative society.


  This paper describes a novel approach to modeling homphily, i.e. the tendency
of nodes that share (or differ in) certain attributes to be linked; we consider
dynamic networks in which nodes can be added over time but not removed. Our
application is to HIV genetic linkage analysis that has been used to
investigate HIV transmission dynamics. In this setting, two HIV sequences from
different persons with HIV (PWH) are said to be linked if the genetic distance
between these sequences is less than a given threshold. Such linkage suggests
that that the nodes representing the two infected PWH, are close to each other
in a transmission network; such proximity would imply that either one of the
infected people directly transmitted the virus to the other or indirectly
transmitted it through a small number of intermediaries. These viral genetic
linkage networks are dynamic in the sense that, over time, a group or cluster
of genetically linked viral sequences may increase in size as new people are
infected by those in the cluster either directly or through intermediaries. Our
approach makes use of a logistic model to describe homophily with regard to
demographic and behavioral characteristics that is we investigate whether
similarities (or differences) between PWH in these characteristics impacts the
probability that their sequences are be linked. Such analyses provide
information about HIV transmission dynamics within a population.


  Reaction mechanisms are often presented as sequences of elementary steps,
such as codified by arrow pushing. We propose an approach for representing such
mechanisms using graph transformation. In this framework, each elementary step
is a rule for modifying a molecular graph and a mechanism is a sequence of such
rules. To generate a compact representation of a multi-step reaction, we
compose the rules of individual steps into a composite rule, providing a
rigorous and fully automated approach to coarse-graining. While the composite
rule retains the graphical conditions necessary for the execution of a
mechanism, it also records information about transient changes not visible by
comparing educts and products. By projecting the rule onto a single "overlay
graph", we generalize Fujita's idea of an Imaginary Transition Structure from
elementary reactions to composite reactions. The utility of the overlay graph
construct is exemplified in the context of enzyme-catalyzed reactions. In a
first application, we exploit mechanistic information in the Mechanism and
Catalytic Site Atlas to construct overlay graphs of hydrolase reactions listed
in the database. These graphs point at a spectrum of catalytic entanglement of
enzyme and substrate, de-emphasizing the notion of a singular catalyst in favor
of a collection of catalytic sites that can be distributed across enzyme and
substrate. In a second application, we deploy composite rules to search the
Rhea database for reactions of known or unknown mechanism that are, in
principle, compatible with the mechanisms implied by the composite rules. We
believe this work adds to the utility of graph-transformation formalisms in
representing and reasoning about chemistry in an automated yet insightful
fashion.


  Scientists world-wide are putting together massive efforts to understand how
the biodiversity that we see on Earth evolved from single-cell organisms at the
origin of life and this diversification process is represented through the Tree
of Life. Low sampling rates and high heterogeneity in the rate of evolution
across sites and lineages produce a phenomenon denoted "long branch attraction"
(LBA) in which long non-sister lineages are estimated to be sisters regardless
of their true evolutionary relationship. LBA has been a pervasive problem in
phylogenetic inference affecting different types of methodologies from
distance-based to likelihood-based. Here, we present a novel neural network
model that outperforms standard phylogenetic methods and other neural network
implementations under LBA settings. Furthermore, unlike existing neural network
models, our model naturally accounts for the tree isomorphisms via permutation
invariant functions which ultimately result in lower memory and allows the
seamless extension to larger trees.


  One of the most effective strategies to mitigate the global spreading of a
pandemic (e.g., COVID-19) is to shut down international airports. From a
network theory perspective, this is since international airports and flights,
essentially playing the roles of bridge nodes and bridge links between
countries as individual communities, dominate the epidemic spreading
characteristics in the whole multi-community system. Among all epidemic
characteristics, the peak fraction of infected, $I_{\max}$, is a decisive
factor in evaluating an epidemic strategy given limited capacity of medical
resources, but is seldom considered in multi-community models. In this paper,
we study a general two-community system interconnected by a fraction $r$ of
bridge nodes and its dynamic properties, especially $I_{\max}$, under the
evolution of the Susceptible-Infected-Recovered (SIR) model. Comparing the
characteristic time scales of different parts of the system allows us to
analytically derive the asymptotic behavior of $I_{\max}$ with $r$, as
$r\rightarrow 0$, which follows different power-law relations in each regime of
the phase diagram. We also detect crossovers when $I_{\max}$ changes from one
power law to another, crossing different power-law regimes as driven by $r$.
Our results enable a better prediction of the effectiveness of strategies
acting on bridge nodes, denoted by the power-law exponent $\epsilon_I$ as in
$I_{\max}\propto r^{1/\epsilon_I}$.


  The analysis of thousands of time series in different languages reveals that
word usage presents oscillations with a prevalence of 16-year cycles, mounted
on slowly varying trends. These components carry different information: while
similar oscillatory patterns gather semantically related words, similar trends
group together keywords representative of cultural and historical periods. We
interpreted the regular oscillations as cycles of interest and saturation,
whose behavior could be captured using a simple mathematical model. Driving the
model with the empirical trends, we were able to explain word frequency traces
across multiple languages throughout the last three centuries. Our results
suggest that word frequency usage is poised at dynamical criticality, close to
a Hopf bifurcation which signals the emergence of oscillatory dynamics.
Crucially, our model explains the oscillatory synchronization observed within
groups of words and provides an interpretation of this phenomenon in terms of
the cultural context driving collective cognition. These findings contribute to
unravel how our use of language is shaped by the interplay between human
cognition and sociocultural forces.


  In a non-sustainable, "over-populated" world, what might the use of
nanotechnology-based targeted, autonomous weapons mean for the future of
humanity? In order to gain some insights, we make a simplified game-theoretical
thought experiment. We consider a population where agents play the public goods
game, and where in parallel an epidemic unfolds. Agents that are infected
defectors are killed with a certain probability and replaced by susceptible
cooperators. We show that such "nanowars", even if aiming to promote good
behavior and planetary health, fail not only to promote cooperation, but they
also significantly increase the probability of repetitive epidemic waves. In
fact, newborn cooperators turn out to be easy targets for defectors in their
neighborhood. Therefore, counterintuitively, the discussed intervention may
even have the opposite effect as desired, promoting defection. We also find a
critical threshold for the death rate of infected defectors, beyond which
resurgent epidemic waves become a certainty. In conclusion, we urgently call
for international regulation of nanotechnology and autonomous weapons.


  Motivation: Bacteriophages are viruses infecting bacteria. Being key players
in microbial communities, they can regulate the composition/function of
microbiome by infecting their bacterial hosts and mediating gene transfer.
Recently, metagenomic sequencing, which can sequence all genetic materials from
various microbiome, has become a popular means for new phage discovery.
However, accurate and comprehensive detection of phages from the metagenomic
data remains difficult. High diversity/abundance, and limited reference genomes
pose major challenges for recruiting phage fragments from metagenomic data.
Existing alignment-based or learning-based models have either low recall or
precision on metagenomic data. Results: In this work, we adopt the
state-of-the-art language model, Transformer, to conduct contextual embedding
for phage contigs. By constructing a protein-cluster vocabulary, we can feed
both the protein composition and the proteins' positions from each contig into
the Transformer. The Transformer can learn the protein organization and
associations using the self-attention mechanism and predicts the label for test
contigs. We rigorously tested our developed tool named PhaMer on multiple
datasets with increasing difficulty, including quality RefSeq genomes, short
contigs, simulated metagenomic data, mock metagenomic data, and the public
IMG/VR dataset. All the experimental results show that PhaMer outperforms the
state-of-the-art tools. In the real metagenomic data experiment, PhaMer
improves the F1-score of phage detection by 27\%.


  The role of transportation vehicles, pig movement between farms, proximity to
infected premises, and feed deliveries has not been fully considered in the
dissemination dynamics of porcine epidemic diarrhea virus (PEDV). This has
limited efforts for disease control and elimination restricting the development
of risk-based resource allocation to the most relevant modes of PEDV
dissemination. Here, we modeled nine modes of between-farm transmission
pathways including farm-to-farm proximity (local transmission), contact network
of pig farm movements between sites, four different contact networks of
transportation vehicles (vehicles that transport pigs from farm-to-farm, pigs
to markets, feed distribution and crew), the volume of animal by-products
within feed diets (e.g. animal fat and meat and bone meal) to reproduce PEDV
transmission dynamics. The model was calibrated in space and time with weekly
PEDV outbreaks. We investigated the model performance to identify outbreak
locations and the contribution of each route in the dissemination of PEDV. The
model estimated that 42.7% of the infections in sow farms were related to
vehicles transporting feed, 34.5% of infected nurseries were associated with
vehicles transporting pigs to farms, and for both farm types, pig movements or
local transmission were the next most relevant routes. On the other hand,
finishers were most often (31.4%) infected via local transmission, followed by
the vehicles transporting feed and pigs to farm networks. Feed ingredients did
not significantly improve model calibration metrics. The proposed modeling
framework provides an evaluation of PEDV transmission dynamics, ranking the
most important routes of PEDV dissemination and granting the swine industry
valuable information to focus efforts and resources on the most important
transmission routes.


  Accurate knowledge of RNA hybridization is essential for understanding RNA
structure and function. Here we mechanically unzip and rezip a 2-kbp RNA
hairpin and derive the 10 nearest-neighbor base pair (NNBP) RNA free energies
in sodium and magnesium with 0.1 kcal/mol precision using optical tweezers.
Notably, force-distance curves (FDCs) exhibit strong irreversible effects with
hysteresis and several intermediates, precluding the extraction of the NNBP
energies with currently available methods. The combination of a suitable RNA
synthesis with a tailored pulling protocol allowed us to obtain the fully
reversible FDCs necessary to derive the NNBP energies. We demonstrate the
equivalence of sodium and magnesium free-energy salt corrections at the level
of individual NNBP. To characterize the irreversibility of the
unzipping-rezipping process, we introduce a barrier energy landscape of the
stem-loop structures forming along the complementary strands, which compete
against the formation of the native hairpin. This landscape correlates with the
hysteresis observed along the FDCs. RNA sequence analysis shows that base
stacking and base pairing stabilize the stem-loops that kinetically trap the
long-lived intermediates observed in the FDC. Stem-loops formation appears as a
general mechanism to explain a wide range of behaviors observed in RNA folding.


  Rapid growth of genetic databases means huge savings from improvements in
their data compression, what requires better inexpensive statistical models.
This article proposes automatized optimizations e.g. of Markov-like models,
especially context binning and model clustering. While it is popular to just
remove low bits of the context, proposed context binning automatically
optimizes such reduction as tabled: state=bin[context] determining probability
distribution, this way extracting nearly all useful information also from very
large contexts, into a relatively small number of states. The second proposed
approach: model clustering uses k-means clustering in space of general
statistical models, allowing to optimize a few models (as cluster centroids) to
be chosen e.g. separately for each read. There are also briefly discussed some
adaptivity techniques to include data non-stationarity.


  Genome sequencing is the basis for many modern biological and medicinal
studies. With recent technological advances, metagenomics has become a problem
of interest. This problem entails the analysis and reconstruction of multiple
DNA sequences from different sources. Shotgun genome sequencing works by
breaking up long DNA sequences into shorter segments called reads. Given this
collection of reads, one would like to reconstruct the original collection of
DNA sequences. For experimental design in metagenomics, it is important to
understand how the minimal read length necessary for reliable reconstruction
depends on the number and characteristics of the genomes involved. Utilizing
simple probabilistic models for each DNA sequence, we analyze the
identifiability of collections of M genomes of length N in an asymptotic regime
in which N tends to infinity and M may grow with N. Our first main result
provides a threshold in terms of M and N so that if the read length exceeds the
threshold, then a simple greedy algorithm successfully reconstructs the full
collection of genomes with probability tending to one. Our second main result
establishes a lower threshold in terms of M and N such that if the read length
is shorter than the threshold, then reconstruction of the full collection of
genomes is impossible with probability tending to one.


  Deep sequencing has become one of the most popular tools for transcriptome
profiling in biomedical studies. While an abundance of computational methods
exists for "normalizing" sequencing data to remove unwanted between-sample
variations due to experimental handling, there is no consensus on which
normalization is the most suitable for a given data set. To address this
problem, we developed "DANA" - an approach for assessing the performance of
normalization methods for microRNA sequencing data based on biology-motivated
and data-driven metrics. Our approach takes advantage of well-known biological
features of microRNAs for their expression pattern and chromosomal clustering
to simultaneously assess (1) how effectively normalization removes handling
artifacts, and (2) how aptly normalization preserves biological signals. With
DANA, we confirm that the performance of eight commonly used normalization
methods vary widely across different data sets and provide guidance for
selecting a suitable method for the data at hand. Hence, it should be adopted
as a routine preprocessing step (preceding normalization) for microRNA
sequencing data analysis. DANA is implemented in R and publicly available at
https://github.com/LXQin/DANA.


  Pattern detection and string matching are fundamental problems in computer
science and the accelerated expansion of bioinformatics and computational
biology have made them a core topic for both disciplines. The SARS-CoV-2
pandemic has made such problems more demanding with hundreds or thousands of
new genome variants discovered every week, because of constant mutations, and
there is a desperate need for fast and accurate analyses. The requirement for
computational tools for genomic analyses, such as sequence alignment, is very
important, although, in most cases the resources and computational power
required are enormous. The presented Multiple Genome Analytics Framework
combines data structures and algorithms, specifically built for text mining and
pattern detection, that can help to efficiently address several computational
biology and bioinformatics problems concurrently with minimal resources. A
single execution of advanced algorithms, with space and time complexity
O(nlogn), is enough to acquire knowledge on all repeated patterns that exist in
multiple genome sequences and this information can be used from other
meta-algorithms for further meta-analyses. The potential of the proposed
framework is demonstrated with the analysis of more than 300,000 SARS-CoV-2
genome sequences and the detection of all repeated patterns with length up to
60 nucleotides in these sequences. These results have been used to provide
answers to questions such as common patterns among all variants, sequence
alignment, palindromes and tandem repeats detection, different organism genome
comparisons, polymerase chain reaction primers detection, etc.


  We study a stochastic individual-based model of interacting plant and
pollinator species through a bipartite graph: each species is a node of the
graph, an edge representing interactions between a pair of species. The
dynamics of the system depends on the between- and within-species interactions:
pollination by insects increases plant reproduction rate but has a cost which
can increase plant death rate, depending on pollinators density. Pollinators
reproduction is increased by the resources harvested on plants. Each species is
characterized by a trait corresponding to its degree of generalism. This trait
determines the structure of the interactions graph and the quantity of
resources exchanged between species. Our model includes in particular nested or
modular networks. Deterministic approximations of the stochastic measure-valued
process by systems of ordinary differential equations or integro-differential
equations are established and studied, when the population is large or when the
graph is dense and can be replaced with a graphon. The long-time behaviors of
these limits are studied and central limit theorems are established to quantify
the difference between the discrete stochastic individual-based model and the
deterministic approximations. Finally, studying the continuous limits of the
interaction network and the resulting PDEs, we show that nested
plant-pollinator communities are expected to collapse towards a coexistence
between a single pair of species of plants and pollinators.


  Objectives: To enhance monitoring of high-burden foodborne pathogens, there
is opportunity to combine pangenome data with network analysis.
  Methods: Salmonella enterica subspecies Enterica serovar Enteritidis isolates
were referred to the New South Wales (NSW) Enteric Reference Laboratory between
August 2015 and December 2019 (1033 isolates in total), inclusive of a
confirmed outbreak. All isolates underwent whole genome sequencing. Distances
between genomes were quantified by in silico MLVA as well as core SNPs, which
informed construction of undirected networks. Prevalence-centrality spaces were
generated from the undirected networks. Components on the undirected SNP
network were considered alongside a phylogenetic tree representation.
  Results: Outbreak isolates were identifiable as distinct components on the
MLVA and SNP networks. The MLVA network based centrality/prevalence space did
not delineate the outbreak, whereas the outbreak was clearly delineated in the
SNP network based centrality/prevalence space. Components on the undirected SNP
network showed a high concordance to the SNP clusters based on phylogenetic
analysis.
  Conclusions: Bacterial whole genome data in network based analysis can
improve the resolution of population analysis. High concordance of network
components and SNP clusters is promising for rapid population analyses of
foodborne Salmonella spp. due to the low overhead of network analysis.


  Range expansion is a universal process in biological systems, and therefore
plays a part in biological evolution. Using a quantitative individual-based
method based on the stochastic process, we identify that enhancing the inherent
self-proliferation advantage of cooperators relative to defectors is a more
effective channel to promote the evolution of cooperation in range expansion
than weakening the benefit acquisition of defectors from cooperators. With this
self-proliferation advantage, cooperators can rapidly colonize virgin space and
establish spatial segregation more readily, which acts like a protective shield
to further promote the evolution of cooperation in return. We also show that
lower cell density and migration rate have a positive effect on the competition
of cooperators with defectors. Biological evolution is based on competition
between individuals and should therefore favor selfish behaviors. However, we
observe a counterintuitive phenomenon that the evolution of a population is
impeded by the fitness-enhancing chemotactic movement of individuals. This
highlights a conflict between the interests of the individual and the
population. The short-sighted selfish behavior of individuals may not be that
favored in the competition between populations. Such information provides
important implications for the handling of cooperation.


  In synaptic molecular communication (MC), the activation of postsynaptic
receptors by neurotransmitters (NTs) is governed by a stochastic
reaction-diffusion process. This randomness of synaptic MC contributes to the
randomness of the electrochemical downstream signal in the postsynaptic cell,
called postsynaptic membrane potential (PSP). Since the randomness of the PSP
is relevant for neural computation and learning, characterizing the statistics
of the PSP is critical. However, the statistical characterization of the
synaptic reaction-diffusion process is difficult because the reversible
bi-molecular reaction of NTs with receptors renders the system nonlinear.
Consequently, there is currently no model available which characterizes the
impact of the statistics of postsynaptic receptor activation on the PSP. In
this work, we propose a novel statistical model for the synaptic
reaction-diffusion process in terms of the chemical master equation (CME). We
further propose a novel numerical method which allows to compute the CME
efficiently and we use this method to characterize the statistics of the PSP.
Finally, we present results from stochastic particle-based computer simulations
which validate the proposed models. We show that the biophysical parameters
governing synaptic transmission shape the autocovariance of the receptor
activation and, ultimately, the statistics of the PSP. Our results suggest that
the processing of the synaptic signal by the postsynaptic cell effectively
mitigates synaptic noise while the statistical characteristics of the synaptic
signal are preserved. The results presented in this paper contribute to a
better understanding of the impact of the randomness of synaptic signal
transmission on neuronal information processing.


  The widespread, and in many countries unprecedented, use of
non-pharmaceutical interventions (NPIs) during the COVID-19 pandemic has
highlighted the need for mathematical models which can estimate the impact of
these measures while accounting for the highly heterogeneous risk profile of
COVID-19. Models accounting either for age structure or the household structure
necessary to explicitly model many NPIs are commonly used in infectious disease
modelling, but models incorporating both levels of structure present
substantial computational and mathematical challenges due to their high
dimensionality. Here we present a modelling framework for the spread of an
epidemic that includes explicit representation of age structure and household
structure. Our model is formulated in terms of tractable systems of ordinary
differential equations for which we provide an open-source Python
implementation. Such tractability leads to significant benefits for model
calibration, exhaustive evaluation of possible parameter values, and
interpretability of results. We demonstrate the flexibility of our model
through four policy case studies, where we quantify the likely benefits of the
following measures which were either considered or implemented in the UK during
the current COVID-19 pandemic: control of within- and between-household mixing
through NPIs; formation of support bubbles during lockdown periods;
out-of-household isolation (OOHI); and temporary relaxation of NPIs during
holiday periods. Our ordinary differential equation formulation and associated
analysis demonstrate that multiple dimensions of risk stratification and social
structure can be incorporated into infectious disease models without
sacrificing mathematical tractability. This model and its software
implementation expand the range of tools available to infectious disease policy
analysts.