  In this study, we formulate a mathematical model incorporating age specific
transmission dynamics of COVID-19 to evaluate the role of vaccination and
treatment strategies in reducing the size of COVID-19 burden. Initially, we
establish the positivity and boundedness of the solutions of the model and
calculate the basic reproduction number. We then formulate an optimal control
problem with vaccination and treatment as control variables. Optimal
vaccination and treatment policies are analysed for different values of the
weight constant associated with the cost of vaccination and different
transmissibility levels. Findings from these suggested that the combined
strategies(vaccination and treatment) worked best in minimizing the infection
and disease induced mortality. In order to reduce COVID-19 infection and
COVID-19 induced deaths to maximum, it was observed that optimal control
strategy should be prioritized to population with age greater than 40 years.
Not much difference was found between individual strategies and combined
strategies in case of mild epidemic ($R_0 \in (0, 2)$). For higher values of
$R_0 (R_0 \in (2, 10))$ the combined strategies was found to be best in terms
of minimizing the overall infection. The infection curves varying the
efficacies of the vaccines were also analysed and it was found that higher
efficacy of the vaccine resulted in lesser number of infection and COVID
induced death.


  Discrete-time replicator map is a prototype of evolutionary selection game
dynamical models that have been very successful across disciplines in rendering
insights into the attainment of the equilibrium outcomes, like the Nash
equilibrium and the evolutionarily stable strategy. By construction, only the
fixed point solutions of the dynamics can possibly be interpreted as the
aforementioned game-theoretic solution concepts. Although more complex outcomes
like chaos are omnipresent in the nature, it is not known to which
game-theoretic solutions they correspond. Here we construct a game-theoretic
solution that is realized as the chaotic outcomes in the selection monotone
game dynamic. To this end, we invoke the idea that in a population game having
two-player--two-strategy one-shot interactions, it is the product of the
fitness and the heterogeneity (the probability of finding two individuals
playing different strategies in the infinitely large population) that is
optimized over the generations of the evolutionary process.


  Synchronization, cooperation, and chaos are ubiquitous phenomena in nature.
In a population composed of many distinct groups of individuals playing the
prisoner's dilemma game, there exists a migration dilemma: No cooperator would
migrate to a group playing the prisoner's dilemma game lest it should be
exploited by a defector; but unless the migration takes place, there is no
chance of the entire population's cooperator-fraction to increase. Employing a
randomly rewired coupled map lattice of chaotic replicator maps, modelling
replication-selection evolutionary game dynamics, we demonstrate that the
cooperators -- evolving in synchrony -- overcome the migration dilemma to
proliferate across the population when altruism is mildly incentivized making
few of the demes play the leader game.


  This paper provides a mathematical model to show that the incorrect
estimation of r, the fraction of asymptomatic COVID-19 carriers in the general
population, can account for much of the world's failure to contain the pandemic
in its early phases. The SE(A+O)R model with infectives separated into
asymptomatic and ordinary carriers, supplemented by a model of the data
generation process, is calibrated to standard datasets for several countries.
It is shown that certain fundamental parameters, notably r, are unidentifiable
with this data. A number of potential types of policy intervention are
analyzed. It is found that the lack of parameter identifiability implies that
only some, but not all, potential policy interventions can be correctly
predicted. In an example representing Italy in March 2020, a hypothetical
optimal policy of isolating confirmed cases that aims to reduce the basic
reproduction number of the outbreak to R0 = 0.8 assuming r = 10%, only achieves
R0 = 1.4 if it turns out that r = 40%.


  Game theory deals with strategic interactions among players and evolutionary
game dynamics tracks the fate of the players' populations under selection. In
this paper, we consider the replicator equation for two-player-two-strategy
games involving cooperators and defectors. We modify the equation to include
the effect of mutation and also delay that corresponds either to the delayed
information about the population state or in realizing the effect of
interaction among players. By focusing on the four exhaustive classes of
symmetrical games -- the Stag Hunt game, the Snowdrift game, the Prisoners'
Dilemma game, and the Harmony game -- we analytically and numerically analyze
delayed replicator-mutator equation to find the explicit condition for the Hopf
bifurcation bringing forth stable limit cycle. The existence of the
asymptotically stable limit cycle imply the coexistence of the cooperators and
the defectors; and in the games, where defection is a stable Nash strategy, a
stable limit cycle does provide a mechanism for evolution of cooperation. We
find that while mutation alone can never lead to oscillatory cooperation state
in two-player-two-strategy games, the delay can change the scenario. On the
other hand, there are situations when delay alone cannot lead to the Hopf
bifurcation in the absence of mutation in the selection dynamics.


  The infectious coronavirus disease 2019 (COVID-19) has become a serious
global pandemic. Different studies have shown that increasing temperature can
play a crucial role in the spread of the virus. Most of these studies were
limited to winter or moderate temperature levels and were conducted using
conventional models. However, traditional models are too simplistic to
investigate complex, non-linear relationships and suffer from some
restrictions. Therefore, we employed copula models to examine the impact of
high temperatures on virus transmission. The findings from the copula models
showed that there was a weak to moderate effect of temperature on the number of
infections and the effect almost vanished under a lockdown policy. Therefore,
this study provides new insight into the relationship between COVID-19 and
temperature, both with and without social isolation practices. Such results can
lead to improvements in our understanding of this new virus. In particular, the
results derived from the copula models examined here, unlike existing
traditional models, provide evidence that there is no substantial influence of
high temperatures on the active COVID-19 outbreak situation. In addition, the
results indicate that the transmission of COVID-19 is strongly influenced by
social isolation practices. To the best of the author knowledge, this is the
first copula model investigation applied to the COVID-19 pandemic.


  The evolution of many microbes and pathogens, including circulating viruses
such as seasonal influenza, is driven by immune pressure from the host
population. In turn, the immune systems of infected populations get updated,
chasing viruses even further away. Quantitatively understanding how these
dynamics result in observed patterns of rapid pathogen and immune adaptation is
instrumental to epidemiological and evolutionary forecasting. Here we present a
mathematical theory of co-evolution between immune systems and viruses in a
finite-dimensional antigenic space, which describes the cross-reactivity of
viral strains and immune systems primed by previous infections. We show the
emergence of an antigenic wave that is pushed forward and canalized by
cross-reactivity. We obtain analytical results for shape, speed, and angular
diffusion of the wave. In particular, we show that viral-immune co-evolution
generates a new emergent timescale, the persistence time of the wave's
direction in antigenic space, which can be much longer than the coalescence
time of the viral population. We compare these dynamics to the observed
antigenic turnover of influenza strains, and we discuss how the dimensionality
of antigenic space impacts on the predictability of the evolutionary dynamics.
Our results provide a concrete and tractable framework to describe
pathogen-host co-evolution.


  With the continued spread of coronavirus, the task of forecasting distinctive
COVID-19 growth curves in different cities, which remain inadequately explained
by standard epidemiological models, is critical for medical supply and
treatment. Predictions must take into account non-pharmaceutical interventions
to slow the spread of coronavirus, including stay-at-home orders, social
distancing, quarantine and compulsory mask-wearing, leading to reductions in
intra-city mobility and viral transmission. Moreover, recent work associating
coronavirus with human mobility and detailed movement data suggest the need to
consider urban mobility in disease forecasts. Here we show that by
incorporating intra-city mobility and policy adoption into a novel
metapopulation SEIR model, we can accurately predict complex COVID-19 growth
patterns in U.S. cities ($R^2$ = 0.990). Estimated mobility change due to
policy interventions is consistent with empirical observation from Apple
Mobility Trends Reports (Pearson's R = 0.872), suggesting the utility of
model-based predictions where data are limited. Our model also reproduces urban
"superspreading", where a few neighborhoods account for most secondary
infections across urban space, arising from uneven neighborhood populations and
heightened intra-city churn in popular neighborhoods. Therefore, our model can
facilitate location-aware mobility reduction policy that more effectively
mitigates disease transmission at similar social cost. Finally, we demonstrate
our model can serve as a fine-grained analytic and simulation framework that
informs the design of rational non-pharmaceutical interventions policies.


  We propose novel chemical reaction networks to translate levels of
concentration into unique DNA strand species, which we call concentration
translators. Our design of the concentration translators is based on
combination of two chemical reaction networks, consensus network and conversion
network with any number of chemical species. We give geometric analysis of the
proposed CRNs from the viewpoint of nonlinear dynamical systems and show that
the CRNs can actually operate as translator. Our concentration translators
exploit DNA strand displacement (DSD) reaction, which is known for a universal
reaction that can implement arbitrary chemical reaction networks. We
demonstrate two specific types of concentration translators (translator A and
B) with different switching behavior and biochemical cost and compared their
characteristics computationally. The proposed concentration translators have an
advantage of being able to readout the concentration of targeted nucleic acid
strand without any fluorescence-based techniques. These characteristics can be
tailored according to requirements from applications, including dynamic range,
sensitivity and implementation cost.


  Human mobility, contact patterns, and their interplay are key aspects of our
social behavior that shape the spread of infectious diseases across different
regions. In the light of new evidence and data sets about these two elements,
epidemic models should be refined to incorporate both the heterogeneity of
human contacts and the complexity of mobility patterns. Here, we propose a
theoretical framework that allows accommodating these two aspects in the form
of a set of Markovian equations. We validate these equations with extensive
mechanistic simulations and derive analytically the epidemic threshold. The
expression of this critical value allows us to evaluate its dependence on the
specific demographic distribution, the structure of mobility flows, and the
heterogeneity of contact patterns, thus shedding light on the microscopic
mechanisms responsible for the epidemic detriment driven by recurrent mobility
patterns reported in the literature.


  A central goal of synthetic biology is the design of molecular controllers
that can manipulate the dynamics of intracellular networks in a stable and
accurate manner. To address the fact that detailed knowledge about
intracellular networks is unavailable, integral-feedback controllers (IFCs)
have been put forward for controlling molecular abundances. These controllers
can maintain accuracy in spite of the uncertainties in the controlled networks.
However, this desirable feature is achieved only if stability is also
maintained. In this paper, we show that molecular IFCs can suffer from a
hazardous instability called negative-equilibrium catastrophe (NEC), whereby
all nonnegative equilibria vanish under the action of the controllers, and some
of the molecular abundances blow up. We show that unimolecular IFCs do not
exist due to a NEC. We then derive a family of bimolecular IFCs that are
safeguarded against NECs when uncertain unimolecular networks, with any number
of molecular species, are controlled. However, when IFCs are applied on
uncertain bimolecular (and hence most intracellular) networks, we show that
preventing NECs generally becomes an intractable problem as the number of
interacting molecular species increases.


  There are many sources of data giving information about the number of
SARS-CoV-2 infections in the population, but all have major drawbacks,
including biases and delayed reporting. For example, the number of confirmed
cases largely underestimates the number of infections, deaths lag infections
substantially, while test positivity rates tend to greatly overestimate
prevalence. Representative random prevalence surveys, the only putatively
unbiased source, are sparse in time and space, and the results come with a big
delay. Reliable estimates of population prevalence are necessary for
understanding the spread of the virus and the effects of mitigation strategies.
We develop a simple Bayesian framework to estimate viral prevalence by
combining the main available data sources. It is based on a discrete-time SIR
model with time-varying reproductive parameter. Our model includes likelihood
components that incorporate data of deaths due to the virus, confirmed cases,
and the number of tests administered on each day. We anchor our inference with
data from random sample testing surveys in Indiana and Ohio. We use the results
from these two states to calibrate the model on positive test counts and
proceed to estimate the infection fatality rate and the number of new
infections on each day in each state in the USA. We estimate the extent to
which reported COVID cases have underestimated true infection counts, which was
large, especially in the first months of the pandemic. We explore the
implications of our results for progress towards herd immunity.


  Corona Virus Disease 2019 (COVID-19), due to its extremely high infectivity,
has been spreading rapidly around the world and bringing huge influence to
socioeconomic development as well as people's daily life. Taking for example
the virus transmission that may occur after college students return to school,
we analyze the quantitative influence of the key factors on the virus spread,
including crowd density and self-protection. One Campus Virus Infection and
Control Simulation model (CVICS) of the novel coronavirus is proposed in this
paper, fully considering the characteristics of repeated contact and strong
mobility of crowd in the closed environment. Specifically, we build an
agent-based infection model, introduce the mean field theory to calculate the
probability of virus transmission, and micro-simulate the daily prevalence of
infection among individuals. The experimental results show that the proposed
model in this paper efficiently simulate how the virus spread in the dense
crowd in frequent contact under closed environment. Furthermore, preventive and
control measures such as self-protection, crowd decentralization and isolation
during the epidemic can effectively delay the arrival of infection peak and
reduce the prevalence, and finally lower the risk of COVID-19 transmission
after the students return to school.


  Based on the information communicated in press releases, and finally
published towards the end of 2020 by Pfizer, Moderna and AstraZeneca, we have
built up a simple Bayesian model, in which the main quantity of interest plays
the role of {\em vaccine efficacy} (`$\epsilon$'). The resulting Bayesian
Network is processed by a Markov Chain Monte Carlo (MCMC), implemented in JAGS
interfaced to R via rjags. As outcome, we get several probability density
functions (pdf's) of $\epsilon$, each conditioned on the data provided by the
three pharma companies. The result is rather stable against large variations of
the number of people participating in the trials and it is `somehow' in good
agreement with the results provided by the companies, in the sense that their
values correspond to the most probable value (`mode') of the pdf's resulting
from MCMC, thus reassuring us about the validity of our simple model. However
we maintain that the number to be reported as `vaccine efficacy' should be the
mean of the distribution, rather than the mode, as it was already very clear to
Laplace about 250 years ago (its `rule of succession' follows from the simplest
problem of the kind). This is particularly important in the case in which the
number of successes equals the numbers of trials, as it happens with the
efficacy against `severe forms' of infection, claimed by Moderna to be 100%.
The implication of the various uncertainties on the predicted number of
vaccinated infectees is also shown, using both MCMC and approximated formulae.


  In evolutionary game theory, it is customary to be partial to the dynamical
models possessing fixed points so that they may be understood as the attainment
of evolutionary stability, and hence, Nash equilibrium. Any show of periodic or
chaotic solution is many a time perceived as a shortcoming of the corresponding
game dynamic because (Nash) equilibrium play is supposed to be robust and
persistent behaviour, and any other behaviour in nature is deemed transient.
Consequently, there is a lack of attempt to connect the non-fixed point
solutions with the game theoretic concepts. Here we provide a way to render
game theoretic meaning to periodic solutions. To this end, we consider a
replicator map that models Darwinian selection mechanism in unstructured
infinite-sized population whose individuals reproduce asexually forming
non-overlapping generations. This is one of the simplest evolutionary game
dynamic that exhibits periodic solutions giving way to chaotic solutions (as
parameters related to reproductive fitness change) and also obeys the folk
theorems connecting fixed point solutions with Nash equilibrium. Interestingly,
we find that a modified Darwinian fitness -- termed heterogeneity payoff -- in
the corresponding population game must be put forward as (conventional) fitness
times the probability that two arbitrarily chosen individuals of the population
adopt two different strategies. The evolutionary dynamics proceeds as if the
individuals optimize the heterogeneity payoff to reach an evolutionarily stable
orbit, should it exist. We rigorously prove that a locally asymptotically
stable period orbit must be heterogeneity stable orbit -- a generalization of
evolutionarily stable state.


  We present a differential equations model in which contagious disease
transmission is affected by contagious fear of the disease and contagious fear
of the control, in this case vaccine. The three contagions are coupled. The two
fears evolve and interact in ways that shape distancing behavior, vaccine
uptake, and their relaxation. These behavioral dynamics in turn can amplify or
suppress disease transmission, which feeds back to affect behavior. The model
reveals several coupled contagion mechanisms for multiple epidemic waves.
Methodologically, the paper advances infectious disease modeling by including
human behavioral adaptation, drawing on the neuroscience of fear learning,
extinction, and transmission.


  The dynamics of epidemic spreading is often reduced to the single control
parameter $R_0$, whose value, above or below unity, determines the state of the
contagion. If, however, the pathogen evolves as it spreads, $R_0$ may change
over time, potentially leading to a mutation-driven spread, in which an
initially sub-pandemic pathogen undergoes a breakthrough mutation. To predict
the boundaries of this pandemic phase, we introduce here a modeling framework
to couple the network spreading patterns with the intra-host evolutionary
dynamics. For many pathogens these two processes, intra- and inter-host, are
driven by different selection forces. And yet here we show that even in the
extreme case when these two forces are mutually independent, mutations can
still fundamentally alter the pandemic phase-diagram, whose transitions are now
shaped, not just by $R_0$, but also by the balance between the epidemic and the
evolutionary timescales. If mutations are too slow, the pathogen prevalence
decays prior to the appearance of a critical mutation. On the other hand, if
mutations are too rapid, the pathogen evolution becomes volatile and, once
again, it fails to spread. Between these two extremes, however, we identify a
broad range of conditions in which an initially sub-pandemic pathogen can break
through to gain widespread prevalence.


  Non-pharmaceutical interventions have been critical in the fight against the
COVID-19 pandemic. However, these sanitary measures have been partially lifted
due to socioeconomic factors causing a worrisome rebound of the epidemic in
several countries. In this work, we assess the effectiveness of the mitigation
implemented to constrain the spread of SARS-CoV-2 in the Mexican territory
during 2020. We also investigate to what extent the initial deployment of the
vaccine will help to mitigate the pandemic and reduce the need for social
distancing and other mobility restrictions. Our modeling approach is based on a
simple mechanistic Kermack-McKendrick-type model. To quantify the effect of
NPIs, we perform a monthly Bayesian inference using officially published data.
The results suggest that in the absence of the sanitary measures, the
cumulative number of infections, hospitalizations, and deaths would have been
at least twice the official number. Moreover, for low vaccine coverage levels,
relaxing NPIs may dramatically increase the disease burden; therefore, safety
measures are of critical importance at the early stages of vaccination. The
simulations also suggest that it may be more desirable to employ a vaccine with
low efficacy but reach a high coverage than a vaccine with high effectiveness
but low coverage levels. This supports the hypothesis that single doses to more
individuals will be more effective than two doses for every person.


  DNA replication in all organisms must overcome nucleoprotein blocks to
complete genome duplication. Accessory replicative helicases in Escherichia
coli, Rep and UvrD, help replication machinery overcome blocks by removing
incoming nucleoprotein complexes or aiding the re-initiation of replication.
Mechanistic details of Rep function have emerged from recent live cell studies,
however, the activities of UvrD in vivo remain unclear. Here, by integrating
biochemical analysis and super-resolved single-molecule fluorescence
microscopy, we discovered that UvrD self-associates into a tetramer and, unlike
Rep, is not recruited to a specific replisome protein despite being found at
approximately 80% of replication forks. By deleting rep and DNA repair factors
mutS and uvrA, perturbing transcription by mutating RNA polymerase, and
antibiotic inhibition; we show that the presence of UvrD at the fork is
dependent on its activity. This is likely mediated by the very high frequency
of replication blocks due to DNA bound proteins, including RNA polymerase, and
DNA damage. UvrD is recruited to sites of nucleoprotein blocks via distinctly
different mechanisms to Rep and therefore plays a more important and
complementary role than previously realised in ensuring successful DNA
replication.


  Currently, due to the COVID-19 pandemic the public life in most European
countries stopped almost completely due to measures against the spread of the
virus. Efforts to limit the number of new infections are threatened by the
advent of new variants of the SARS-COV-2 virus, most prominent the B.1.1.7
strain with higher infectivity. In this article we consider a basic two-strain
SIR model to explain the spread of those variants in Germany on small time
scales. For a linearized version of the model we calculate relevant variables
like the time of minimal infections or the dynamics of the share of variants
analytically. These analytical approximations and numerical simulations are in
a good agreement to data reported by the Robert Koch Institute (RKI) in
Germany.


  We show that in a situation where individuals have a choice between a costly
institute and a free institute to perform a collective action task, the
existence of a participation cost promotes cooperation in the costly institute.
Despite paying for a participation cost, costly cooperators, who join the
costly institute and cooperate, can out-perform defectors, who predominantly
join a free institute. This, not only promotes cooperation in the costly
institute but also facilitates the evolution of cooperation in the free
institute. A costly institute out-performs a free institute when the
profitability of the collective action is low. On the other hand, a free
institute performs better when the collective action's profitability is high.
Furthermore, we show that in a structured population, when individuals have a
choice between different institutes, a mutualistic relation between cooperators
with different institute preferences emerges and helps the evolution of
cooperation.


  Simulating soft matter systems such as the cytoskeleton can enable deep
understanding of experimentally observed phenomena. One challenge of modeling
such systems is realistic description of the steric repulsion between nearby
polymers. Previous models of the polymeric excluded volume interaction have the
deficit of being non-analytic, being computationally expensive, or allowing
polymers to erroneously cross each other. A recent solution to these issues,
implemented in the MEDYAN simulation platform, uses analytical expressions
obtained from integrating an interaction kernel along the lengths of two
polymer segments to describe their repulsion. Here, we extend this model by
re-deriving it for lower-dimensional geometrical configurations, deriving
similar expressions using a steeper interaction kernel, comparing it to other
commonly used potentials, and showing how to parameterize these models. We also
generalize this new integrated style of potential by introducing a segmental
Lennard-Jones potential, which enables modelling both attractive and repulsive
interactions in semi-flexible polymer networks. These results can be further
generalized to facilitate the development of effective interaction potentials
for other finite elements in simulations of soft-matter systems.


  Can a regulated, legal market for wildlife products protect species
threatened by poaching? It is one of the most controversial ideas in
biodiversity conservation. Perhaps the most convincing reason for legalizing
wildlife trade is that trade revenue could fund the protection and conservation
of poached species. In this paper, we examine the possible poacher-population
dynamic consequences of legal trade funding conservation. The model consists of
a manager scavenging carcasses for wildlife products, who then sells the
products, and directs a portion of the revenue towards funding anti-poaching
law enforcement. Through a global analysis of the model, we derive the critical
proportion of product the manager must scavenge, and the critical proportion of
trade revenue the manager must allocate towards increased enforcement, in order
for legal trade to lead to abundant long-term wildlife populations. We
illustrate how the model could inform management with parameter values derived
from the African elephant literature, under a hypothetical scenario where a
manager scavenges elephant carcasses to sell ivory. We find that there is a
large region of parameter space where populations go extinct under legal trade,
unless a significant portion of trade revenue is directed towards protecting
populations from poaching. The model is general and therefore can be used as a
starting point for exploring the consequences of funding many conservation
programs using wildlife trade revenue.


  Sexual selection is a fundamental aspect of evolution for all eukaryotic
organisms with mating types. This paper suggests intersexual selection is best
viewed as a mechanism to compensate for the unavoidable dynamics of coevolution
between sexes that emerge with isogamy. Using the NK model of fitness
landscapes, the conditions under which allosomes emerge are first explored.
This extends previous work on the evolution of sex where the fitness landscape
smoothing of a rudimentary form of the Baldwin effect is suggested as the
underlying cause. The NKCS model of coevolution is then used to show how
varying fitness landscape size, ruggedness, and connectedness can vary the
conditions under which a very simple sexual selection mechanism proves
beneficial. This is found to be the case whether one or both sexes exploit
sexual selection.


  A general framework for obtaining exact transition rate matrices for
stochastic systems on networks is presented and applied to many well-known
compartmental models of epidemiology. The state of the population is described
as a vector in the tensor product space of $N$ individual probability vector
spaces, whose dimension equals the number of compartments of the
epidemiological model $n_c$. The transition rate matrix for the
$n_c^N$-dimensional Markov chain is obtained by taking suitable linear
combinations of tensor products of $n_c$-dimensional matrices. The resulting
transition rate matrix is a sum over bilocal linear operators, which gives
insight in the microscopic dynamics of the system. The more familiar and
non-linear node-based mean-field approximations are recovered by restricting
the exact models to uncorrelated (separable) states. We show how the exact
transition rate matrix for the susceptible-infected (SI) model can be used to
find analytic solutions for SI outbreaks on trees and the cycle graph for
finite $N$.


  Existing compartmental models in epidemiology are limited in terms of
optimizing the resource allocation to control an epidemic outbreak under
disease growth uncertainty. In this study, we address this core limitation by
presenting a multi-stage stochastic programming compartmental model, which
integrates the uncertain disease progression and resource allocation to control
an infectious disease outbreak. The proposed multi-stage stochastic program
involves various disease growth scenarios and optimizes the distribution of
treatment centers and resources while minimizing the total expected number of
new infections and funerals. We define two new equity metrics, namely infection
and capacity equity, and explicitly consider equity for allocating treatment
funds and facilities over multiple time stages. We also study the multi-stage
value of the stochastic solution (VSS), which demonstrates the superiority of
the proposed stochastic programming model over its deterministic counterpart.
We apply the proposed formulation to control the Ebola Virus Disease (EVD) in
Guinea, Sierra Leone, and Liberia of West Africa to determine the optimal and
fair resource-allocation strategies. Our model balances the proportion of
infections over all regions, even without including the infection equity or
prevalence equity constraints. Model results also show that allocating
treatment resources proportional to population is sub-optimal, and enforcing
such a resource allocation policy might adversely impact the total number of
infections and deaths, and thus resulting in a high cost that we have to pay
for the fairness. Our multi-stage stochastic epidemic-logistics model is
practical and can be adapted to control other infectious diseases in
meta-populations and dynamically evolving situations.


  The COVID-19 pandemic has had a great impact in most countries worldwide, in
particular impacting Chile greatly to become one of the worst hit countries in
the world, despite its low population compared to other affected countries
around the world. In this study we report the evolution of the COVID-19
pandemic and the spread of the SARS-CoV-2 virus in Chile since the first
positive case was announced on March 3, 2020, and until November 30, 2020. We
provide a detailed description of the data provided by the Chilean Ministry of
Science per administrative region on the number of new cases, tests per capita,
and deaths in the country.


  The site frequency spectrum (SFS) is a popular summary statistic of genomic
data. While the SFS of a constant-sized population undergoing neutral mutations
has been extensively studied in population genetics, the rapidly growing amount
of cancer genomic data has attracted interest in the spectrum of an
exponentially growing population. Recent theoretical results have generally
dealt with special or limiting cases, such as considering only cells with an
infinite line of descent, assuming deterministic tumor growth, or taking
large-time or large-population limits. In this work, we derive exact
expressions for the expected SFS of a cell population that evolves according to
a stochastic branching process, first for cells with an infinite line of
descent and then for the total population, evaluated either at a fixed time
(fixed-time spectrum) or at the stochastic time at which the population reaches
a certain size (fixed-size spectrum). We find that while the rate of mutation
scales the SFS of the total population linearly, the rates of cell birth and
cell death change the shape of the spectrum at the small-frequency end,
inducing a transition between a $1/j^2$ power-law spectrum and a $1/j$ spectrum
as cell viability decreases. We show that this insight can in principle be used
to estimate the ratio between the rate of cell death and cell birth, as well as
the mutation rate, using the site frequency spectrum alone. Although the
discussion is framed in terms of tumor dynamics, our results apply to any
exponentially growing population of individuals undergoing neutral mutations.


  The interaction between individuals in biological populations, dilute
components of chemical systems, or particles transported by turbulent flows
depends critically on their contact statistics. This work clarifies those
statistics under the simplifying assumptions that the underlying motions
approximate a Brownian random walk and that the particles can be treated as
noninteracting. We measure the contact-interval (also called the waiting-time
or inter-arrival-time), contact-count, and contact-duration distributions in
populations of individuals undergoing noninteracting continuous-space-time
random walks on a periodic two-dimensional plane (a torus), as functions of the
population number density, walker radius, and random-walk step size. The
contact-interval is exponentially distributed for times longer than the
mean-free-collision time but not for times shorter than that, and the contact
duration distribution is strongly peak at the ballistic-crossing time for
head-on collisions, when the ballistic-crossing time is short compared to the
mean step duration. While successive contacts between individuals are
independent, the probability of repeat contact decreases with time after a
previous contact. This leads to a negative duration dependence of the
waiting-time interval and over dispersion of the contact-count probability
density function for all time intervals. The paper demonstrates that for
populations of small particles (walker radius small compared to the
mean-separation or random-walk step size) the mean-free-collision interval, the
ballistic-crossing time, and the random-walk-step duration can be used to
construct temporal scalings which allow for common waiting-time, contact-count,
and contact-duration distributions across different populations. Semi-analytic
approximations for both the waiting-time and contact-duration distributions are
also presented.


  Background: By the beginning of December 2020, some vaccines against COVID-19
already presented efficacy and security, which qualify them to be used in mass
vaccination campaigns. Thus, setting up strategies of vaccination became
crucial to control the COVID19 pandemic. Methods: We use daily COVID-19 reports
from Chicago and NYC from 01-Mar2020 to 28- Nov-2020 to estimate the parameters
of an SEIR-like epidemiological model that accounts for different severity
levels. To achieve data adherent predictions, we let the model parameters to be
time-dependent. The model is used to forecast different vaccination scenarios,
where the campaign starts at different dates, from 01-Oct-2020 to 01-Apr-2021.
To generate realistic scenarios, disease control strategies are implemented
whenever the number of predicted daily hospitalizations reaches a preset
threshold. Results: The model reproduces the empirical data with remarkable
accuracy. Delaying the vaccination severely affects the mortality,
hospitalization, and recovery projections. In Chicago, the disease spread was
under control, reducing the mortality increment as the start of the vaccination
was postponed. In NYC, the number of cases was increasing, thus, the estimated
model predicted a much larger impact, despite the implementation of contention
measures. The earlier the vaccination campaign begins, the larger is its
potential impact in reducing the COVID-19 cases, as well as in the
hospitalizations and deaths. Moreover, the rate at which cases,
hospitalizations and deaths increase with the delay in the vaccination
beginning strongly depends on the shape of the incidence of infection in each
city.


  At the end of 2020, policy responses to the SARS-CoV-2 outbreak have been
shaken by the emergence of virus variants, impacting public health and policy
measures worldwide. The emergence of these strains suspected to be more
contagious, more severe, or even resistant to antibodies and vaccines, seem to
have taken by surprise health services and policymakers, struggling to adapt to
the new variants constraints. Anticipating the emergence of these mutations to
plan ahead adequate policies, and understanding how human behaviors may affect
the evolution of viruses by coevolution, are key challenges. In this article,
we propose coevolution with genetic algorithms (GAs) as a credible approach to
model this relationship, highlighting its implications, potential and
challenges. Because of their qualities of exploration of large spaces of
possible solutions, capacity to generate novelty, and natural genetic focus,
GAs are relevant for this issue. We present a dual GA model in which both
viruses aiming for survival and policy measures aiming at minimising infection
rates in the population, competitively evolve. This artificial coevolution
system may offer us a laboratory to "debug" our current policy measures,
identify the weaknesses of our current strategies, and anticipate the evolution
of the virus to plan ahead relevant policies. It also constitutes a decisive
opportunity to develop new genetic algorithms capable of simulating much more
complex objects. We highlight some structural innovations for GAs for that
virus evolution context that may carry promising developments in evolutionary
computation, artificial life and AI.


  The effective reproduction number, R(t), is a central point in the study of
infectious diseases. It establishes in an explicit way the extent of an
epidemic spread process in a population. The current estimation methods for the
time evolution of R(t), using incidence data, rely on the generation interval
distribution, g(\tau), which is usually obtained from empirical data or already
known distributions from the literature. However, there are systems, especially
highly heterogeneous ones, in which there is a lack of data and an adequate
methodology to obtain g(\tau). In this work, we use mathematical models to
bridge this gap. We present a general methodology for obtaining an explicit
expression of the reproduction numbers and the generation interval
distributions provided by an arbitrary compartmental model. Additionally, we
present the appropriate expressions to evaluate those reproduction numbers
using incidence data. To highlight the relevance of such methodology, we apply
it to the spread of Covid-19 in municipalities of the state of Rio de janeiro,
Brazil. Using two meta-population models, we estimate the reproduction numbers
and the contributions of each municipality in the generation of cases in all
others. Our results point out the importance of mathematical modelling to
provide epidemiological meaning of the available data.


  Virtual screening is a computational technique widely used for identifying
small molecules which are most likely to bind to a protein target. Here, we
performed a molecular docking study to propose potential candidates to prevent
the RBD/ACE2 attachment. These candidates are sixteen different flavonoids
present in the peppermint leaf. Results showed that Luteolin
7-O-neohesperidoside is the peppermint flavonoid with a higher binding affinity
regarding the RBD/ACE2 complex (about -9.18 Kcal/mol). On the other hand,
Sakuranetin presented the lowest affinity (about -6.38 Kcal/mol). Binding
affinities of the other peppermint flavonoids ranged from -6.44 Kcal/mol up to
-9.05 Kcal/mol. The binding site surface analysis showed pocket-like regions on
the RBD/ACE2 complex that yield several interactions (mostly hydrogen bonds)
between the flavonoid and the amino acid residues of the proteins. This study
can open channels for the understanding of the roles of flavonoids against
COVID-19 infection.


  Human pluripotent stem cells (hPSCs) have promising clinical applications in
regenerative medicine, drug-discovery and personalised medicine due to their
potential to differentiate into all cell types, a property know as
pluripotency. A deeper understanding of how pluripotency is regulated is
required to assist in controlling pluripotency and differentiation trajectories
experimentally. Mathematical modelling provides a non-invasive tool through
which to explore, characterise and replicate the regulation of pluripotency and
the consequences on cell fate. Here we use experimental data of the expression
of the pluripotency transcription factor OCT4 in a growing hPSC colony to
develop and evaluate mathematical models for temporal pluripotency regulation.
We consider fractional Brownian motion and the stochastic logistic equation and
explore the effects of both additive and multiplicative noise. We illustrate
the use of time-dependent carrying capacities and the introduction of Allee
effects to the stochastic logistic equation to describe cell differentiation.
This mathematical framework for describing intra-cellular OCT4 regulation can
be extended to other transcription factors and developed into sophisticated
predictive models.


  We investigate the effect of mobile polymer brushes on proteins embedded in
biological membranes by employing both Asakura-Oosawa type of theoretical model
and coarse-grained molecular dynamics simulations. The brush polymer-induced
depletion attraction between proteins changes non-monotonically with the size
of brush. The depletion interaction, which is determined by the ratio of
protein size to the grafting distance between brush polymers, increases
linearly with brush size as long as the polymer brush height is shorter than
the protein size. When the brush height exceeds the protein size, however, the
depletion attraction among proteins is slightly reduced. We also explore the
possibility of brush polymer-induced assembly of a large protein cluster, which
can be related to one of many molecular mechanisms underlying recent
experimental observations of integrin nanocluster formation and signaling.


  We investigated daily COVID-19 cases and deaths in the 337 lower tier local
authority regions in England and Wales to better understand how the disease
propagated over a 15-month period. Population density scaling models revealed
residual variance and skewness to be sensitive indicators of the dynamics of
propagation. Lockdowns and schools reopening triggered increased variance
indicative of outbreaks with local impact and country scale heterogeneity.
University reopening and December holidays triggered reduced variance
indicative of country scale homogenisation which reached a minimum in
mid-January 2021. Homogeneous propagation was associated with better
correspondence with normally distributed residuals while heterogeneous
propagation was more consistent with skewed models. Skewness varied from
strongly negative to strongly positive revealing an unappreciated feature of
community propagation. Hot spots and super-spreading events are well understood
descriptors of regional disease dynamics that would be expected to be
associated with positively skewed distributions. Positively skewed behaviour
was observed; however, negative skewness indicative of cold-spots and
super-isolation dominated for approximately 8 months during the period of
study. In contrast, death metrics showed near constant behaviour in scaling,
variance, and skewness metrics over the full period with rural regions
preferentially affected, an observation consistent with regional age
demographics in England and Wales. Regional positions relative to density
scaling laws were remarkably persistent after the first 5-9 days of the
available data set. The determinants of this persistent behaviour probably
precede the pandemic and remain unchanged.


  We prove a large-deviation principle (LDP) for the sample paths of jump
Markov processes in the small noise limit when, possibly, all the jump rates
vanish uniformly, but slowly enough, in a region of the state space. We further
discuss the optimality of our assumptions on the decay of the jump rates. As a
direct application of this work we relax the assumptions needed for the
application of LDPs to, e.g., Chemical Reaction Network dynamics, where
vanishing reaction rates arise naturally particularly the context of mass
action kinetics.


  Understanding the consequences of the combined effects of multiple
stressors-including stress from man-made chemicals is important for
conservation management, the ecological risk assessment of chemicals, and many
other ecological applications. Our current ability to predict and analyse the
joint effects of multiple stressors is insufficient to make the prospective
risk assessment of chemicals more ecologically relevant because we lack a full
understanding of how organisms respond to stress factors alone and in
combination. Here, we describe a Dynamic Energy Budget (DEB) based
bioenergetics model that predicts the potential effects of single or multiple
natural and chemical stressors on life history traits. We demonstrate the
plausibility of the model using a meta-analysis of 128 existing studies on
freshwater invertebrates. We then validate our model by comparing its
predictions for a combination of three stressors (i.e. chemical, temperature,
and food availability) with new, independent experimental data on life history
traits in the daphnid Ceriodaphnia dubia. We found that the model predictions
are in agreement with observed growth curves and reproductive traits. To the
best of our knowledge, this is the first time that the combined effects of
three stress factors on life history traits observed in laboratory studies have
been predicted successfully in invertebrates. We suggest that a re-analysis of
existing studies on multiple stressors within the modelling framework outlined
here will provide a robust null model for identifying stressor interactions,
and expect that a better understanding of the underlying mechanisms will arise
from these new analyses. Bioenergetics modelling could be applied more broadly
to support environmental management decision making.


  This survey analyses the role of data-driven methodologies for pandemic
modelling and control. We provide a roadmap from the access to epidemiological
data sources to the control of epidemic phenomena. We review the available
methodologies and discuss the challenges in the development of data-driven
strategies to combat the spreading of infectious diseases. Our aim is to bring
together several different disciplines required to provide a holistic approach
to epidemic analysis, such as data science, epidemiology, and
systems-and-control theory. A 3M-analysis is presented, whose three pillars
are: Monitoring, Modelling and Managing. The focus is on the potential of
data-driven schemes to address three different challenges raised by a pandemic:
(i) monitoring the epidemic evolution and assessing the effectiveness of the
adopted countermeasures; (ii) modelling and forecasting the spread of the
epidemic; (iii) making timely decisions to manage, mitigate and suppress the
contagion. For each step of this roadmap, we review consolidated theoretical
approaches (including data-driven methodologies that have been shown to be
successful in other contexts) and discuss their application to past or present
epidemics, such as Covid-19, as well as their potential application to future
epidemics.


  Covid-19 is particularly mild with children, and its severity escalates with
age. Several theories have been proposed to explain these facts. In particular,
it was proposed that the lower expression of the viral receptor ACE2 in
children protects them from severe Covid. However, other works suggested an
inverse relationship between ACE2 expression and disease severity. Here we try
to reconcile seemingly contradicting observations noting that ACE2 is not
monotonically related with age but it reaches a maximum at a young age that
depends on the cell type and then decreases. This pattern is consistent with
most existing data from humans and rodents and it is expected to be more marked
for ACE2 cell protein than for mRNA because of the increase with age of the
protease TACE/ADAM17 that sheds ACE2 from the cell membrane to the serum. The
negative relation between ACE2 level and Covid-19 severity at old age is not
paradoxical but it is consistent with a mathematical model of virus propagation
that predicts that higher viral receptor does not necessarily favour virus
propagation and it can even slow it down. More importantly, ACE2 is known to
protect organs from chronic and acute inflammation, which are worsened by low
ACE2 levels. Here we propose that ACE2 contributes essentially to reverse the
inflammatory process by downregulating the pro-inflammatory peptides of the
angiotensin and bradykinin system, and that failure to revert the inflammation
triggered by SARS-COV-2 may underlie both severe CoViD-19 infection and its
many post-infection manifestations, including the multi-inflammatory syndrome
of children (MIS-C). Within this view, lower severity in children despite lower
ACE2 expression may be consistent with their higher expression of the
alternative angiotensin II receptor ATR2 and in general of the
anti-inflammatory arm of the Renin-Angiotensin System (RAS) at young age.


  Modeling the distribution of high dimensional data by a latent tree graphical
model is a prevalent approach in multiple scientific domains. A common task is
to infer the underlying tree structure, given only observations of its terminal
nodes. Many algorithms for tree recovery are computationally intensive, which
limits their applicability to trees of moderate size. For large trees, a common
approach, termed divide-and-conquer, is to recover the tree structure in two
steps. First, recover the structure separately of multiple, possibly random
subsets of the terminal nodes. Second, merge the resulting subtrees to form a
full tree. Here, we develop Spectral Top-Down Recovery (STDR), a deterministic
divide-and-conquer approach to infer large latent tree models. Unlike previous
methods, STDR partitions the terminal nodes in a non random way, based on the
Fiedler vector of a suitable Laplacian matrix related to the observed nodes. We
prove that under certain conditions, this partitioning is consistent with the
tree structure. This, in turn, leads to a significantly simpler merging
procedure of the small subtrees. We prove that STDR is statistically consistent
and bound the number of samples required to accurately recover the tree with
high probability. Using simulated data from several common tree models in
phylogenetics, we demonstrate that STDR has a significant advantage in terms of
runtime, with improved or similar accuracy.


  Background: Structural Variations, SVs, in a genome can be linked to a
disease or characteristic phenotype. The variations come in many types and it
is a challenge, not only determining the variations accurately, but also
conducting the downstream statistical and analytical procedure. Method:
Structural variations, SVs, with size 1 base-pair to 1000s of base-pairs with
their precise breakpoints and single-nucleotide polymorphisms, SNPs, were
determined for members of a family. The genome was assembled using optimal
metrics of ABySS and SOAPdenovo assembly tools using paired-end DNA sequence.
Results: An interesting discovery was the mitochondrial DNA could have paternal
leakage of inheritance or that the mutations could be high from maternal
inheritance. It is also discovered that the mitochondrial DNA is less prone to
SVs re-arrangements than SNPs, which propose better standards for determining
ancestry and divergence between races and species over a long-time frame. Sex
determination of an individual is found to be strongly confirmed using calls of
nucleotide bases of SVs to the Y chromosome, more strongly determined than
SNPs. We note that in general there is a larger variance -and thus the standard
deviation, in the sum of SVs nucleotide compared to sum of SNPs of an
individual when compared to reference sequence, and thus SVs serve as a
stronger means to characterize an individual for a given trait or phenotype or
to determine sex. The SVs and SNPs in HLA loci would also serve as a medical
transformation method for determining the success of an organ transplant for a
patient, and predisposition to diseases apriori.


  Univariate and multivariate methods for association of the genom-ic
variations with the end-or-endo phenotype have been widely used for genome wide
association studies. In addition to encoding the SNPs, we advocate usage of
clustering as a novel method to encode the structural variations, SVs, in
genomes, such as the deletions and insertions polymorphism (DIPs), Copy Number
Variations (CNVs), translocation, inversion, etc., that can be used as an
independent fea-ture variable value for downstream computation by artificial
intelli-gence methods to predict the endo-or-end phenotype. We introduce a
clustering based encoding scheme for structural variations and om-ics based
analysis. We conducted a complete all genomic variants association with the
phenotype using deep learning and other ma-chine learning techniques, though
other methods such as genetic al-gorithm can also be applied. Applying this
encoding of SVs and one-hot encoding of SNPs on GTEx V7 pilot DNA variation
dataset, we were able to get high accuracy using various methods of DMWAS, and
particularly found logistic regression to work the best for death due to
heart-attack (MHHRTATT) phenotype. The genom-ic variants acting as feature sets
were then arranged in descending order of power of impact on the disease or
trait phenotype, which we call optimization and that also uses top univariate
association into account. Variant Id P1_M_061510_3_402_P at chromosome 3 &
position 192063195 was found to be most highly associated to MHHRTATT. We
present here the top ten optimized genomic va-riant feature set for the
MHHRTATT phenotypic cause of death.


  This paper proposes a very simple deterministic mathematical model, which, by
using a power-law, is a \emph{non-integer power model} (or \emph{fractional
power model (FPM)}). Such a model, in non-integer power of time, namely $t^m$
up to constants, enables representing at each day, with a good precision, the
totality of the contaminated individuals. Despite being enriched with knowledge
through an internal structure based on a geometric sequence "with variable
ratio", the model (in its non-integer representation) has only three
parameters, among which the non-integer power, $m$, that determines on its own,
according to its value, an aggravation or an improvement of the viral
spreading. Its simplicity comes from the power-law, $t^m$, which simply
expresses the singular dynamics of the operator of non-integer differentiation
or integration, of high parametric compactness, that governs diffusion
phenomena and, as shown in this paper, the spreading phenomena by
contamination. The proposed model is indeed validated with the official data of
Ministry of Health on the COVID-19 spreading. Used in prediction, it well
enables justifying the choice of a lockdown, without which the spreading would
have highly worsened. The comparison of this model in $t^m$ with two known
models having the same number of parameters, well shows that its
representativity of the real data is better or more general. Finally, in a more
fundamental context and particularly in terms of complexity and simplicity, a
self-filtering action enables showing the compatibility between the
\emph{internal complexity} that the internal structure and its stochastic
behavior present, and the \emph{global simplicity} that the model in $t^m$
offers in a deterministic manner: it is true that the non-integer power of a
power-law is well a marker of complexity.


  Vaccination against COVID-19 with the recently approved mRNA vaccines
BNT162b2 (BioNTech/Pfizer) and mRNA-1273 (Moderna) is currently underway in a
large number of countries. However, high incidence rates and rapidly spreading
SARS-CoV-2 variants are concerning. In combination with acute supply deficits
in Europe in early 2021, the question arises of whether stretching the vaccine,
for instance by delaying the second dose, can make a significant contribution
to preventing deaths, despite associated risks such as lower vaccine efficacy,
the potential emergence of escape mutants, enhancement, waning immunity,
reduced social acceptance of off-label vaccination, and liability shifts. A
quantitative epidemiological assessment of risks and benefits of non-standard
vaccination protocols remains elusive. To clarify the situation and to provide
a quantitative epidemiological foundation we develop a stochastic
epidemiological model that integrates specific vaccine rollout protocols into a
risk-group structured infectious disease dynamical model. Using the situation
and conditions in Germany as a reference system, we show that delaying the
second vaccine dose is expected to prevent deaths in the four to five digit
range, should the incidence resurge. We show that this considerable public
health benefit relies on the fact that both mRNA vaccines provide substantial
protection against severe COVID-19 and death beginning 12 to 14 days after the
first dose. The benefits of protocol change are attenuated should vaccine
compliance decrease substantially. To quantify the impact of protocol change on
vaccination adherence we performed a large-scale online survey. We find that,
in Germany, changing vaccination protocols may lead to small reductions in
vaccination intention. In sum, we therefore expect the benefits of a strategy
change to remain substantial and stable.


  The presence of metamorphism in the protein's native state is not yet fully
understood. In an attempt to throw light on this issue here we present an
assessment, in terms of the amide hydrogen exchange protection factor, that
aims to determine the likely existence of structural fluctuations in the
native-state consistent with both the upper bound marginal stability of
proteins and the metamorphism presence. The preliminary results enable us to
conclude that the native-state metamorphism is, indeed, more probable than
thought.


  The beautiful structures of single and multi-domain proteins are clearly
ordered in some fashion but cannot be readily classified using group theory
methods that are successfully used to describe periodic crystals. For this
reason, protein structures are considered to be aperiodic, and may have evolved
this way for functional purposes, especially in instances that require a
combination of softness and rigidity within the same molecule. By analyzing the
solved protein structures, we show that orientational symmetry is broken in the
aperiodic arrangement of the secondary structural elements (SSEs), which we
deduce by calculating the nematic order parameter, $P_{2}$. We find that the
folded structures are nematic droplets with a broad distribution of $P_{2}$. We
argue that non-zero values of $P_{2}$, leads to an arrangement of the SSEs that
can resist mechanical forces, which is a requirement for allosteric proteins.
Such proteins, which resist mechanical forces in some regions while being
flexible in others, transmit signals from one region of the protein to another
(action at a distance) in response to binding of ligands (oxygen, ATP or other
small molecules).


  How cooperation evolves and manifests itself in the thermodynamic or infinite
player limit of social dilemma games is a matter of intense speculation.
Various analytical methods have been proposed to analyze the thermodynamic
limit of social dilemmas. In this work, we compare two analytical methods,
i.e., Darwinian evolution and Nash equilibrium mapping, with a numerical
agent-based approach. For completeness, we also give results for another
analytical method, Hamiltonian dynamics. In contrast to Hamiltonian dynamics,
which involves the maximization of payoffs of all individuals, in Darwinian
evolution, the payoff of a single player is maximized with respect to its
interaction with the nearest neighbor. While the Hamiltonian dynamics method
utterly fails as compared to Nash equilibrium mapping, the Darwinian evolution
method gives a false positive for game magnetization -- the net difference
between the fraction of cooperators and defectors -- when payoffs obey the
condition a + d = b + c, wherein a,d represents the diagonal elements and b,c
the off-diagonal elements in a symmetric social dilemma game payoff matrix.
When either a + d =/= b + c or when one looks at the average payoff per player,
the Darwinian evolution method fails, much like the Hamiltonian dynamics
approach. On the other hand, the Nash equilibrium mapping and numerical
agent-based method agree well for both game magnetization and average payoff
per player for the social dilemmas in question, i.e., the Hawk-Dove game and
the Public goods game. This paper thus brings to light the inconsistency of the
Darwinian evolution method vis-a-vis both Nash equilibrium mapping and a
numerical agent-based approach.


  The site frequency spectrum describes variation among a set of n DNA
sequences. Its i'th entry (i=1,2,...,n-1) is the number of nucleotide sites at
which the mutant allele is present in i copies. Under selective neutrality,
random mating, and constant population size, the expected value of the spectrum
is well known but somewhat puzzling. Each additional sequence added to a sample
adds an entry to the end of the expected spectrum but does not affect existing
entries. This note reviews the reasons for this behavior.


  Preliminary epidemiologic, phylogenetic and clinical findings suggest that
several novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)
variants have increased transmissibility and decreased efficacy of several
existing vaccines. Four mutations in the receptor-binding domain (RBD) of the
spike protein that are reported to contribute to increased transmission.
Understanding physical mechanism responsible for the affinity enhancement
between the SARS-CoV-2 variants and ACE2 is the "urgent challenge" for
developing blockers, vaccines and therapeutic antibodies against the
coronavirus disease 2019 (COVID-19) pandemic. Based on a
hydrophobic-interaction-based protein docking mechanism, this study reveals
that the mutation N501Y obviously increased the hydrophobic attraction and
decrease hydrophilic repulsion between the RBD and ACE2 that most likely caused
the transmissibility increment of the variants. By analyzing the
mutation-induced hydrophobic surface changes in the attraction and repulsion at
the binding site of the complexes of the SARS-CoV-2 variants and antibodies, we
found out that all the mutations of N501Y, E484K, K417N and L452R can
selectively decrease or increase their binding affinity with some antibodies.


  In this study, we have developed a parallel version of the random time
simulation algorithm.
  Firstly, we gave a rigorous basis of the random time description of the
stochastic process of chemical reaction network time evolution.
  And then we reviewed the random time simulation algorithm and gave the
implementations for the parallel version of next reaction random time
algorithm.
  The discussion of computational complexity suggested a factor of $M$ (which
is the connection number of the network) folds time consuming reduction for
random time simulation algorithm as compared to other exact stochastic
simulation algorithms, such as the Gillespie algorithm.
  For large-scale system, such like the protein-protein interaction network,
$M$ is on order of $10^8$.
  We further demonstrate the power of random time simulation with a GPGPU
parallel implementation which achieved roughly 100 folds acceleration as
compared with CPU implementations.
  Therefore the stochastic simulation method we developed here can be of great
application value for simulating time evolution process of large-scale network.


  eQuilibrator (equilibrator.weizmann.ac.il) is a calculator for biochemical
equilibrium constants and Gibbs free energies, originally designed as a
web-based interface. While the website now counts ${\sim}1000$ distinct monthly
users, its design could not accommodate larger compound databases and it lacked
an application programming interface (API) for integration in other tools
developed by the systems biology community. Here, we report a new python-based
package for eQuilibrator, that comes with many new features such as a 50-fold
larger compound database, the ability to add novel compound structures,
improvements in speed and memory use, and correction for Mg2+ ion
concentrations. Moreover, it adds the ability to compute the covariance matrix
of the uncertainty between estimates, for which we show the advantages and
describe the application in metabolic modeling. We foresee that these
improvements will make thermodynamic modeling more accessible and facilitate
the integration of eQuilibrator into other software platforms.


  When faced with an imminent risk of predation, many animals react to escape
consumption. Antipredator strategies are performed by individuals acting as a
group to intimidate predators and minimize the damage when attacked. We study
the antipredator prey response in spatial tritrophic systems with cyclic
species dominance using the rock-paper-scissors game. The impact of the
antipredator behavior is local, with the predation probability reducing
exponentially with the number of preys in the predator's neighborhood. In
contrast to the standard Lotka-Volterra implementation of the
rock-paper-scissors model, where no spiral waves appear, our outcomes show that
the antipredator behavior leads to spiral patterns from random initial
conditions. The results show that the predation risk decreases exponentially
with the level of antipredator strength. Finally, we investigate the
coexistence probability and verify that antipredator behavior may jeopardize
biodiversity for high mobility. Our findings may help biologists to understand
ecosystems formed by species whose individuals behave strategically to resist
predation.


  COVID-19 has devastated the entire global community. Vaccines present an
opportunity to mitigate the pandemic; however, the effect of vaccination
coupled with the behavioral response of the population is not well understood.
We propose a model that incorporates two important dynamically varying
population behaviors: level of caution and sense of safety. Level of caution
increases with the number of infectious cases, while an increasing sense of
safety with increased vaccination lowers precautionary behaviors. To the best
of our knowledge, this is the first model that can effectively reproduce the
complete time history of COVID-19 infections for various regions of the United
States and provides relatable measures of dynamic changes in the population
behavior and disease transmission rates. We propose a parameter d_I as a direct
measure of a population's caution against an infectious disease, that can be
obtained from the ongoing new infectious cases. The model provides a method for
quantitative measure of critical infectious disease attributes for a population
including highest disease transmission rate, effective disease transmission
rate, and disease related precautionary behavior. We predict future COVID-19
pandemic trends in the United States accounting for vaccine rollout and
behavioral response. Although a high rate of vaccination is critical to quickly
end the pandemic, we find that a return towards pre-pandemic social behavior
due to increased sense of safety during vaccine deployment, can cause an
alarming surge in infections. Our results indicate that at the current rate of
vaccination, the new infection cases for COVID-19 in the United States will
approach zero by the end of August 2021. The model can be used for predicting
future epidemic and pandemic dynamics before and during vaccination.


  The Neolithic Revolution began c. 10000 years ago and is characterised by the
ultimate, near complete transition from hunting and gathering to agricultural
food production on land. The Neolithic Revolution is thought to have been
catalysed by a combination of local population pressure, cultural diffusion,
property rights and climate change. We undertake a thought experiment that
examines trends in these key hypothesised catalysts and patters of today to
explore whether society could be on a path towards another paradigm shift in
food production: away from hunting of wild fish towards a transition to mostly
fish farming. We find similar environmental and cultural pressures have driven
the rapid rise of aquaculture, during a period that has now been coined the
Blue Revolution, providing impetus for such a transition in coming decades to
centuries. We also highlight the interacting and often mutually reinforcing
impacts of 1)technological and scientific advancement, 2)environmental
awareness and collective action and 3)globalisation and trade influencing the
trajectory and momentum of the Blue Revolution. We present two qualitative
narratives that broadly fall within two future trajectories: 1)a ubiquitous
aquaculture transition and 20commercial aquaculture and fisheries coexistence.
This scenarios approach aims to encourage logical, forward thinking, and
innovative solutions to complex systems dynamics. Scenario-based thought
experiments are useful to explore large scale questions, increase the
accessibility to a wider readership and ideally catalyse discussion around
proactive governance mechanisms. We argue the future is not fixed and society
now has greater foresight and capacity to choose the workable balance between
fisheries sand aquaculture that supports economic, environmental, cultural and
social objectives through combined planning, policies and management.


  In the wake of the 2020 COVID-19 epidemic, much work has been performed on
the development of mathematical models for the simulation of the epidemic, and
of disease models generally. Most works follow the susceptible-infected-removed
(SIR) compartmental framework, modeling the epidemic with a system of ordinary
differential equations. Alternative formulations using a partial differential
equation (PDE) to incorporate both spatial and temporal resolution have also
been introduced, with their numerical results showing potentially powerful
descriptive and predictive capacity. In the present work, we introduce a new
variation to such models by using delay differential equations (DDEs). The
dynamics of many infectious diseases, including COVID-19, exhibit delays due to
incubation periods and related phenomena. Accordingly, DDE models allow for a
natural representation of the problem dynamics, in addition to offering
advantages in terms of computational time and modeling, as they eliminate the
need for additional, difficult-to-estimate, compartments (such as exposed
individuals) to incorporate time delays. Here, we introduce a DDE epidemic
model in both an ordinary- and partial differential equation framework. We
present a series of mathematical results assessing the stability of the
formulation. We then perform several numerical experiments, validating both the
mathematical results and establishing model's ability to reproduce measured
data on realistic problems.


  Institutions can provide incentives to increase cooperation behaviour in a
population where this behaviour is infrequent. This process is costly, and it
is thus important to optimize the overall spending. This problem can be
mathematically formulated as a multi-objective optimization problem where one
wishes to minimize the cost of providing incentives while ensuring a desired
level of cooperation within the population. In this paper, we provide a
rigorous analysis for this problem. We study cooperation dilemmas in both the
pairwise (the Donation game) and multi-player (the Public Goods game) settings.
We prove the regularity of the (total incentive) cost function, characterize
its asymptotic limits (infinite population, weak selection and large selection)
and show exactly when reward or punishment is more efficient. We prove that the
cost function exhibits a phase transition phenomena when the intensity of
selection varies. We calculate the critical threshold in regards to the phase
transition and study the optimization problem when the intensity of selection
is under and above the critical value. It allows us to provide an exact
calculation for the optimal cost of incentive, for a given intensity of
selection. Finally, we provide numerical simulations to demonstrate the
analytical results. Overall, our analysis provides for the first time a
selection-dependent calculation of the optimal cost of institutional incentives
(for both reward and punishment) that guarantees a minimum amount of
cooperation. It is of crucial importance for real-world applications of
institutional incentives since intensity of selection is specific to a given
population and the underlying game payoff structure.


  Detecting cancers at early stages can dramatically reduce mortality rates.
Therefore, practical cancer screening at the population level is needed. Here,
we develop a comprehensive detection system to classify all common cancer
types. By integrating artificial intelligence deep learning neural network and
noncoding RNA biomarkers selected from massive data, our system can accurately
detect cancer vs healthy object with 96.3% of AUC of ROC (Area Under Curve of a
Receiver Operating Characteristic curve). Intriguinely, with no more than 6
biomarkers, our approach can easily discriminate any individual cancer type vs
normal with 99% to 100% AUC. Furthermore, a comprehensive marker panel can
simultaneously multi-classify all common cancers with a stable 78% of accuracy
at heterological cancerous tissues and conditions. This provides a valuable
framework for large scale cancer screening. The AI models and plots of results
were available in https://combai.org/ai/cancerdetection/


  We propose a minimal off-lattice model of living organisms where just a very
few dynamical rules of growth are assumed. The stable coexistence of many
clusters is detected when we replace the global restriction rule by a locally
applied one. A rich variety of evolving patterns is revealed where players
movement has a decisive role on the evolutionary outcome. For example,
intensive individual mobility may jeopardize the survival of the population,
but if we increase players movement further then it can save the population.
Notably, the collective drive of population members is capable to compensate
the negative consequence of intensive movement and keeps the system alive. When
the drive becomes biased then the resulting unidirectional flow alters the
stable pattern and produce a stripe-like state instead of the previously
observed hexagonal arrangement of clusters. Interestingly, the rotation of
stripes can be flipped if the individual movement exceeds a threshold value.


  In the recent COVID-19 pandemic, a wide range of epidemiological modelling
approaches have been used to predict the effective reproduction number, R(t),
and other COVID-19 related measures such as the daily rate of exponential
growth, r(t). These candidate models use different modelling approaches or
differing assumptions about spatial or age mixing, and some capture genuine
uncertainty in scientific understanding of disease dynamics. Combining
estimates using appropriate statistical methodology from multiple candidate
models is important to better understand the variation of these outcome
measures to help inform decision making. In this paper, we combine these
estimates for specific UK nations and regions using random effects meta
analyses techniques, utilising the restricted maximum likelihood (REML) method
to estimate the heterogeneity variance parameter, and two approaches to
calculate the confidence interval for the combined estimate: the standard
Wald-type intervals; and the Knapp and Hartung (KNHA) method. As estimates in
this setting are derived using model predictions, each with varying degrees of
uncertainty, equal weighting is favoured over the more standard
inverse-variance weighting in order avoid potential up-weighting of models
providing estimates with lower levels of uncertainty that are not fully
accounting for inherent uncertainties. Utilising these meta-analysis techniques
has allowed for statistically robust combined estimates to be calculated for
key COVID-19 outcome measures. This in turn allows timely and informed decision
making based on all of the available information.


  In 2020 and 2021, the spread of COVID-19 was globally addressed by imposing
restrictions on the distance of individual travel. Recent literature has
uncovered a clear pattern in human mobility that underlies the complexity of
urban mobility: $r \cdot f$, the product of distance traveled $r$ and frequency
of return $f$ per user to a given location, is invariant across space. This
paper asks whether the invariant $r\cdot f$ also serves as a driver for
epidemic spread, so that the risk associated with human movement can be modeled
by a unifying variable $r\cdot f$. We use two large-scale datasets of
individual human mobility to show that there is in fact a simple relation
between $r$ and $f$ and both speed and spatial dispersion of disease spread.
This discovery could assist in modeling spread of disease and inform travel
policies in future epidemics -- based not only on travel distance $r$ but also
on frequency of return $f$.


  During Covid-19 outbreaks, school closures are employed as part of
governments' non-pharmaceutical interventions around the world to reduce the
number of contacts and keep the reproduction number below 1. Yet, prolonged
school closures have profound negative impact on the future opportunities of
pupils, particularly from disadvantaged backgrounds, as well as additional
economic and social impacts by preventing their parents from returning to work.
Data on Covid-19 in children are sparse and policy frameworks are evolving
quickly. We compare a set of potential policies to accompany the reopening of
schools by means of an agent-based simulation tool. The policies and scenarios
we model reflect the public discussion and government guidelines in early March
2021 in England before the planned nationwide reopening of schools on the 8th
of March. A point of particular interest is the potential contribution of a
more wide-spread use of screening tests based on lateral flow devices. We
compare policies both with respect to their potential to contain new outbreaks
of Covid-19 in schools and the proportion of schooldays lost due to isolation
of pupils. We find that regular asymptomatic screening of the whole school as
an addition to a policy built around isolation of symptomatic pupils and their
closest contacts is beneficial across a wide range of scenarios, including when
screening tests with relatively low test sensitivity are used. Multiple
screening tests per week bring only small additional benefits in some
scenarios. These findings remain valid when test compliance is not enforced
although the effectiveness of outbreak control is reduced.


  The present paper introduces a data-driven framework for describing the
time-varying nature of an SIRD model in the context of COVID-19. By embedding a
rolling regression in a mixed integer bilevel nonlinear programming problem,
our aim is to provide the research community with a model that reproduces
accurately the observed changes in the number of infected, recovered, and death
cases, while providing information about the time dependency of the parameters
that govern the SIRD model. We propose this optimization model and a genetic
algorithm to tackle its solution. Moreover, we test this algorithm with 2020
COVID-19 data from the state of Minnesota and found that our results are
consistent both qualitatively and quantitatively, thus proving that the
framework proposed is an effective an flexible tool to describe the dynamics of
a pandemic.


  Niche and neutral theory are two prevailing, yet much debated, ideas in
ecology proposed to explain the patterns of biodiversity. Whereas niche theory
emphasizes selective differences between species and interspecific interactions
in shaping the community, neutral theory supposes functional equivalence
between species and points to stochasticity as the primary driver of ecological
dynamics. In this work, we draw a bridge between these two opposing theories.
Starting from a Lotka-Volterra (LV) model with demographic noise and random
symmetric interactions, we analytically derive the stationary population
statistics and species abundance distribution (SAD). Using these results, we
demonstrate that the model can exhibit three classes of SADs commonly found in
niche and neutral theories and found conditions that allow an ecosystem to
transition between these various regimes. Thus, we reconcile how neutral-like
statistics may arise from a diverse community with niche differentiation.


  Facing a global epidemic of new infectious diseases such as COVID-19,
non-pharmaceutical interventions (NPIs), which reduce transmission rates
without medical actions, are being implemented around the world to mitigate
spreads. One of the problems in assessing the effects of NPIs is that different
NPIs have been implemented at different times based on the situation of each
country; therefore, few assumptions can be shared about how the introduction of
policies affects the patient population. Mathematical models can contribute to
further understanding these phenomena by obtaining analytical solutions as well
as numerical simulations. In this study, an NPI was introduced into the SIR
model for a conceptual study of infectious diseases under the condition that
the transmission rate was reduced to a fixed value only once within a finite
time duration, and its effect was analyzed numerically and theoretically. It
was analytically shown that the maximum fraction of infected individuals and
the final size could be larger if the intervention starts too early. The
analytical results also suggested that more individuals may be infected at the
peak of the second wave with a stronger intervention. This study provides
quantitative relationship between the strength of a one-shot intervention and
the reduction in the number of patients with no approximation. This suggests
the importance of the strength and time of NPIs, although detailed studies are
necessary for the implementation of NPIs in complicated real-world environments
as the model used in this study is based on various simplifications.


  The near-infinite chemical diversity of natural and artificial macromolecules
arises from the vast range of possible component monomers, linkages, and
polymers topologies. This enormous variety contributes to the ubiquity and
indispensability of macromolecules but hinders the development of general
machine learning methods with macromolecules as input. To address this, we
developed GLAMOUR, a framework for chemistry-informed graph representation of
macromolecules that enables quantifying structural similarity, and
interpretable supervised learning for macromolecules.


  The spread of infectious epidemics is often accelerated by super-spreader
events. Understanding their effect is important, particularly in the context of
standard epidemiological models, which require estimates for parameters such as
$R_0$. In this letter, we show that the effective value of $R_0$ in
super-spreader situations is significantly large, of the order of hundreds,
suggesting a delta-function-like behavior during the event. Use of a well-mixed
room model supports these findings. They elucidate infection kinetic modeling
in enclosed environments, which differ from the standard SIR model, and provide
expressions for $R_0$ in terms of physical and operational parameters. The
overall impact of super-spreader events can be significant, depending on the
state of the epidemic and how the infections generated by the event
subsequently spread in the community.


  A mathematical model of COVID-19 with minimal compartments is developed. The
model is simple enough to fit data on confirmed cases, estimate the hidden
infection figure and incorporate the effect of vaccination. With the effect of
the new variant being considered, the model fits well with Philippine data of
confirmed cases. With the fitted parameters, results show that a second peak,
strikingly similar in magnitude to the first, is to be expected if the current
transmission rate is increased by about 5\% and already a complete disaster if
it will go up to 20\%. The model is then used to estimate time to achieve herd
immunity considering natural immunity from being infected and vaccination plan
of the country. Results show that if the plan of 200,000 vaccinated people per
day is achieved then herd immunity can be attained by March 2022. But if only
half of that is achieved the time could go up to February 2023.


  Generative models emerge as promising candidates for novel sequence-data
driven approaches to protein design, and for the extraction of structural and
functional information about proteins deeply hidden in rapidly growing sequence
databases. Here we propose simple autoregressive models as highly accurate but
computationally efficient generative sequence models. We show that they perform
similarly to existing approaches based on Boltzmann machines or deep generative
models, but at a substantially lower computational cost (by a factor between
$10^2$ and $10^3$). Furthermore, the simple structure of our models has
distinctive mathematical advantages, which translate into an improved
applicability in sequence generation and evaluation. Within these models, we
can easily estimate both the probability of a given sequence, and, using the
model's entropy, the size of the functional sequence space related to a
specific protein family. In the example of response regulators, we find a huge
number of ca. $10^{68}$ possible sequences, which nevertheless constitute only
the astronomically small fraction $10^{-80}$ of all amino-acid sequences of the
same length. These findings illustrate the potential and the difficulty in
exploring sequence space via generative sequence models.


  The number of elements (N) and types (S) sampled from an ecological system
are among the most powerful constraints on observations of abundance,
distribution, and diversity. Together, N and S determine sets of possible forms
(i.e., feasible sets) for the species abundance distribution (SAD). There are
three approaches to the description of the null SAD (= the average feasible
SAD). The first approach is based on the random uniform sampling of
surjections. I calculate the probability of a given SAD, given N and S, under
this approach (Eq. 4). The second approach is based on the random sampling of
compositions. I calculate the probability of a given SAD, given N and S, under
this approach (Eq. 8). The third approach is based on the random uniform
sampling of partitions. I review the approach, which was developed by Locey &
White (2013), and provide some asymptotic results useful for ecologists. The
center of a feasible set is a null expectation, which should deviate enough
from the alternative models before invoking for biological or ecological
mechanisms underlying the SAD. Here, I integrate the feasible set approach with
the typical framework of inference in ecology (goodness-of-fit, null hypothesis
testing, model comparison, null modelling). I describe how to perform numerical
simulations to describe expectations under different approaches to the feasible
set. I develop objective or fitness functions to allow the estimation of the
most likely SAD using numerical optimization. I provide tools to compare null
expectations based on the feasible set approach with the observations, in the
context of model comparison and null hypothesis testing.


  Mutations sometimes increase contagiousness for evolving pathogens. During an
epidemic, scientists use viral genome data to infer a shared evolutionary
history and connect this history to geographic spread. We propose a model that
directly relates a pathogen's evolution to its spatial contagion dynamics --
effectively combining the two epidemiological paradigms of phylogenetic
inference and self-exciting process modeling -- and apply this
\emph{phylogenetic Hawkes process} to a Bayesian analysis of 23,422 viral cases
from the 2014-2016 Ebola outbreak in West Africa. The proposed model is able to
detect individual viruses with significantly elevated rates of spatiotemporal
propagation for a subset of 1,610 samples that provide genome data. Finally, to
facilitate model application in big data settings, we develop massively
parallel implementations for the gradient and Hessian of the log-likelihood and
apply our high performance computing framework within an adaptively
preconditioned Hamiltonian Monte Carlo routine.


  Elucidating the mechanisms that lead to cooperation is still one of the main
scientific challenges of current times, as many common cooperative scenarios
remain elusive and at odds with Darwin's natural selection theory. Here, we
study evolutionary games on populations that are structured beyond pairwise
interactions. Specifically, we introduce a general evolutionary approach that
allows studying situations in which indirect interactions via a neighbor other
than the direct pairwise connection (or via a group of neighbors), impact the
strategy of the focal player. To this end, we consider simplicial graphs that
encode two- and three-body interactions, which enables to study competition
between all possible pairs of social dilemmas and to scrutinize the role of
three-body interactions in all the observed phenomenology. We simultaneously
investigate how social dilemma with different Nash equilibria compete in
simplicial structures and how such a competition is modulated by the unbalance
of 2- and 1-simplices, which in its turn reflects the relative prevalence of
pairwise or group interactions among the players. We report a number of results
that: (i) support that higher-order games allow for non-dominant strategists to
emerge and coexist with dominant ones, a scenario that can't be explained by
any pairwise schemes, no matter the network of contacts; (ii) characterize a
novel transition from dominant defection to dominant cooperation as a function
of the simplicial structure of the population; and (iii) demonstrate that
2-simplex interactions are a source of strategy diversity, i.e. increasing the
relative prevalence of group interactions always promotes diverse strategic
identities of individuals. Our study constitutes a step forward in the quest
for understanding the roots of cooperation and the mechanisms that sustain it
in real-world and social environments.


  Antibody therapeutics has been extensively studied in drug discovery and
development within the past decades. One increasingly popular focus in the
antibody discovery pipeline is the optimization step for therapeutic leads.
Both traditional methods and in silico approaches aim to generate candidates
with high binding affinity against specific target antigens. Traditional in
vitro approaches use hybridoma or phage display for candidate selection, and
surface plasmon resonance (SPR) for evaluation, while in silico computational
approaches aim to reduce the high cost and improve efficiency by incorporating
mathematical algorithms and computational processing power in the design
process. In the present study, we investigated different graph-based designs
for depicting antibody-antigen interactions in terms of antibody affinity
prediction using deep learning techniques. While other in silico computations
require experimentally determined crystal structures, our study took interest
in the capability of sequence-based models for in silico antibody maturation.
Our preliminary studies achieved satisfying prediction accuracy on binding
affinities comparing to conventional approaches and other deep learning
approaches. To further study the antibody-antigen binding specificity, and to
simulate the optimization process in real-world scenario, we introduced
pairwise prediction strategy. We performed analysis based on both baseline and
pairwise prediction results. The resulting prediction and efficiency prove the
feasibility and computational efficiency of sequence-based method to be adapted
as a scalable industry practice.


  We develop a modeling framework for bioluminescence light found in the deep
sea near neutrino telescopes by combining a hydrodynamic model with a
stochastic one. The bioluminescence is caused by organisms when exposed to a
non-constant water flow, such as past the neutrino telescopes. We model the
flow using the incompressible Navier-Stokes equations for Reynolds numbers
between 4000 and 23000. The discretization relies on a finite element method
which includes upwind-stabilization for the velocity field. On top of the flow
model, we simulate a population of random microscopic organisms. Their movement
and emission are stochastic processes which we model using Monte Carlo methods.
We observe unique time-series for the photon counts depending on the flow
velocity and detector specifications. This opens up the possibility of
categorizing organisms using neutrino detectors. We show that the average
light-yield and pulse shapes require precise flow modeling, while the emission
timing is chaotic. From this we construct a fast modeling scheme, requiring
only a subset of computationally expensive flow and population modeling.


  A single gene can encode for different protein versions through a process
called alternative splicing. Since proteins play major roles in cellular
functions, aberrant splicing profiles can result in a variety of diseases,
including cancers. Alternative splicing is determined by the gene's primary
sequence and other regulatory factors such as RNA-binding protein levels. With
these as input, we formulate the prediction of RNA splicing as a regression
task and build a new training dataset (CAPD) to benchmark learned models. We
propose discrete compositional energy network (DCEN) which leverages the
hierarchical relationships between splice sites, junctions and transcripts to
approach this task. In the case of alternative splicing prediction, DCEN models
mRNA transcript probabilities through its constituent splice junctions' energy
values. These transcript probabilities are subsequently mapped to relative
abundance values of key nucleotides and trained with ground-truth experimental
measurements. Through our experiments on CAPD, we show that DCEN outperforms
baselines and ablation variants.


  The widespread of Coronavirus has led to a worldwide pandemic with a high
mortality rate. Currently, the knowledge accumulated from different studies
about this virus is very limited. Leveraging a wide-range of biological
knowledge, such as gene ontology and protein-protein interaction (PPI) networks
from other closely related species presents a vital approach to infer the
molecular impact of a new species. In this paper, we propose the transferred
multi-relational embedding model Bio-JOIE to capture the knowledge of gene
ontology and PPI networks, which demonstrates superb capability in modeling the
SARS-CoV-2-human protein interactions. Bio-JOIE jointly trains two model
components. The knowledge model encodes the relational facts from the protein
and GO domains into separated embedding spaces, using a hierarchy-aware
encoding technique employed for the GO terms. On top of that, the transfer
model learns a non-linear transformation to transfer the knowledge of PPIs and
gene ontology annotations across their embedding spaces. By leveraging only
structured knowledge, Bio-JOIE significantly outperforms existing
state-of-the-art methods in PPI type prediction on multiple species.
Furthermore, we also demonstrate the potential of leveraging the learned
representations on clustering proteins with enzymatic function into enzyme
commission families. Finally, we show that Bio-JOIE can accurately identify
PPIs between the SARS-CoV-2 proteins and human proteins, providing valuable
insights for advancing research on this new disease.