  Contact tracing and quarantine are well established non-pharmaceutical
epidemic control tools. The paper aims to clarify the impact of these measures
in COVID-19 epidemic. A new deterministic model is introduced (SEIRQ:
susceptible, exposed, infectious, removed, quarantined) with Q compartment
capturing individuals and releasing them with delay. We obtain a simple rule
defining the reproduction number $\mathcal{R}$ in terms of quarantine
parameters, ratio of diagnosed cases and transmission parameters. The model is
applied to the epidemic in Poland in March - April 2020, when social distancing
measures were in place. We investigate 3 scenarios corresponding to different
ratios of diagnosed cases. Our results show that depending on the scenario
contact tracing could have prevented from 50\% to over 90\% of cases. The
effects of quarantine are limited by fraction of undiagnosed cases. Taking into
account the transmission intensity in Poland prior to introduction of social
restrictions it is unlikely that the control of the epidemic could be achieved
without any social distancing measures.


  The ongoing Coronavirus Disease 2019 (COVID-19) pandemic threatens the health
of humans and causes great economic losses. Predictive modelling and
forecasting the epidemic trends are essential for developing countermeasures to
mitigate this pandemic. We develop a network model, where each node represents
an individual and the edges represent contacts between individuals where the
infection can spread. The individuals are classified based on the number of
contacts they have each day (their node degrees) and their infection status.
The transmission network model was respectively fitted to the reported data for
the COVID-19 epidemic in Wuhan (China), Toronto (Canada), and the Italian
Republic using a Markov Chain Monte Carlo (MCMC) optimization algorithm. Our
model fits all three regions well with narrow confidence intervals and could be
adapted to simulate other megacities or regions. The model projections on the
role of containment strategies can help inform public health authorities to
plan control measures.


  The focus of this article is on the dynamics of a new susceptible-infected
model which consists of a susceptible group ($S$) and two different infectious
groups ($I_1$ and $I_2$). Once infected, an individual becomes a member of one
of these infectious groups which have different clinical forms of infection. In
addition, during the progress of the illness, an infected individual in group
$I_1$ may pass to the infectious group $I_2$ which has a higher mortality rate.
In this study, positiveness of the solutions for the model is proved. Stability
analysis of species extinction, $I_1$-free equilibrium and endemic equilibrium
as well as disease-free equilibrium is studied. Relation between the basic
reproduction number of the disease and the basic reproduction number of each
infectious stage is examined. The model is investigated under a specific
condition and its exact solution is obtained.


  Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infectivity is a
major concern in coronavirus disease 2019 (COVID-19) prevention and economic
reopening. However, rigorous determination of SARS-COV-2 infectivity is
essentially impossible owing to its continuous evolution with over 13752 single
nucleotide polymorphisms (SNP) variants in six different subtypes. We develop
an advanced machine learning algorithm based on the algebraic topology to
quantitatively evaluate the binding affinity changes of SARS-CoV-2 spike
glycoprotein (S protein) and host angiotensin-converting enzyme 2 (ACE2)
receptor following the mutations. Based on mutation-induced binding affinity
changes, we reveal that five out of six SARS-CoV-2 subtypes have become either
moderately or slightly more infectious, while one subtype has weakened its
infectivity. We find that SARS-CoV-2 is slightly more infectious than SARS-CoV
according to computed S protein-ACE2 binding affinity changes. Based on a
systematic evaluation of all possible 3686 future mutations on the S protein
receptor-binding domain (RBD), we show that most likely future mutations will
make SARS-CoV-2 more infectious. Combining sequence alignment, probability
analysis, and binding affinity calculation, we predict that a few residues on
the receptor-binding motif (RBM), i.e., 452, 489, 500, 501, and 505, have very
high chances to mutate into significantly more infectious COVID-19 strains.


  The "New York State on Pause" executive order came into effect on March 22
with the goal of ensuring adequate social distancing to alleviate the spread of
COVID-19. Pause will remain effective in New York City in some form until early
June. We use a compartmentalized model to study the effects of testing capacity
and social distancing measures on the evolution of the pandemic in the
"post-Pause" period in the City. We find that testing capacity must increase
dramatically if it is to counterbalance even relatively small relaxations in
social distancing measures in the immediate post-Pause period. In particular,
if the City performs 20,000 tests per day and relaxes the social distancing
measures to the pre-Pause norms, then the total number of deaths by the end of
September can reach 250,000. By keeping the social distancing measures to
somewhere halfway between the pre- and in-Pause norms and performing 100,000
tests per day, the total number of deaths by the end of September can be kept
at around 27,000. Going back to the pre-Pause social distancing norms quickly
must be accompanied by an exorbitant testing capacity, if one is to suppress
excessive deaths. If the City is to go back to the "pre-Pause" social
distancing norms in the immediate post-Pause period and keep the total number
of deaths by the end of September at around 35,000, then it should be
performing 500,000 tests per day. Our findings have important implications on
the magnitude of the testing capacity the City needs as it relaxes the social
distancing measures to reopen its economy.


  Homeostasis of protein concentrations in cells is crucial for their proper
functioning, and this requires concentrations (at their steady-state levels) to
be stable to fluctuations. Since gene expression is regulated by proteins such
as transcription factors (TFs), the full set of proteins within the cell
constitutes a large system of interacting components. Here, we explore factors
affecting the stability of this system by coupling the dynamics of mRNAs and
protein concentrations in a growing cell. We find that it is possible for
protein concentrations to become unstable if the regulation strengths or system
size becomes too large, and that other global structural features of the
networks can dramatically enhance the stability of the system. In particular,
given the same number of proteins, TFs, number of interactions, and regulation
strengths, a network that resembles a bipartite graph with a lower fraction of
interactions that target TFs has a higher chance of being stable. By scrambling
the $\textit{E. coli.}$ transcription network, we find that the randomized
network with the same number of regulatory interactions is much more likely to
be unstable than the real network. These findings suggest that constraints
imposed by system stability could have played a role in shaping the existing
regulatory network during the evolutionary process. We also find that contrary
to what one might expect from random matrix theory and what has been argued in
the literature, the degradation rate of mRNA does not affect whether the system
is stable.


  A quantitative COVID-19 model that incorporates hidden asymptomatic patients
is developed, and an analytic solution in parametric form is given. The model
incorporates the impact of lockdown and resulting spatial migration of
population due to announcement of lockdown. A method is presented for
estimating the model parameters from real-world data. It is shown that increase
of infections slows down and herd immunity is achieved when symptomatic
patients are 4-6\% of the population for the European countries we studied,
when the total infected fraction is between 50-56 \%. Finally, a method for
estimating the number of asymptomatic patients, who have been the key hidden
link in the spread of the infections, is presented.


  We adapt a SEIRD differential model with asymptomatic population and Covid
deaths, which we call SEAIRD, to simulate the evolution of COVID-19, and add a
control function affecting both the diffusion of the virus and GDP, featuring
all direct and indirect containment policies; to model feasibility, the control
is assumed to be a piece-wise linear function satisfying additional
constraints. We describe the joint dynamics of infection and the economy and
discuss the trade-off between production and fatalities. In particular, we
carefully study the conditions for the existence of the optimal policy response
and its uniqueness. Uniqueness crucially depends on the marginal rate of
substitution between the statistical value of a human life and GDP; we show an
example with a phase transition: above a certain threshold, there is a unique
optimal containment policy; below the threshold, it is optimal to abstain from
any containment; and at the threshold itself there are two optimal policies. We
then explore and evaluate various profiles of various control policies
dependent on a small number of parameters.


  Movement is fundamental to life, shaping population dynamics, biodiversity
patterns, and ecosystem structure. Recent advances in tracking technology have
enabled fundamental questions about movement to be tackled, leading to the
development of the movement ecology framework (MEF), considered a milestone in
the field [1]. The MEF introduced an integrative theory of organismal movement,
linking internal state, motion capacity and navigation capacity to external
factors. Here, a decade later, we investigated the current state of research in
the field. Using a text mining approach on >8000 peer-reviewed papers in
movement ecology, we explored the main research topics, evaluated the impact of
the MEF, and assessed changes in the use of technological devices, software and
statistical methods. The number of publications has increased considerably and
there have been major technological changes in the past decade (i.e.~increased
use of GPS devices, accelerometers and video cameras, and a convergence towards
R), yet we found that research focuses on the same questions, specifically, on
the effect of environmental factors on movement and behavior. In practice, it
appears that movement ecology research does not reflect the MEF. We call on
researchers to transform the field from technology-driven to embrace
interdisciplinary collaboration, in order to reveal key processes underlying
movement (e.g.~navigation), as well as evolutionary, physiological and
life-history consequences of particular strategies.


  Condensation of hundreds of mega-base-pair-long human chromosomes in a small
nuclear volume is a spectacular phenomenon. This process is driven by the
formation of chromosome loops. ATP consuming motor, condensin, interacts with
chromatin segments to actively extrude loops. Motivated by real-time imaging of
loop extrusion (LE) and measurements using magnetic tweezer experiments, we
created an analytically solvable model, predicting the LE velocity and step
size distribution as a function of external load. The theory fits the
experimental data quantitatively, and suggests that condensin must undergo a
large conformational change, induced by ATP binding, bringing distant parts of
the motor to proximity. Simulations using a simple model confirm that the motor
transitions between an open to closed state in order to extrude loops by a
scrunching mechanism, similar to that proposed in DNA bubble formation during
bacterial transcription. Changes in the orientation of the motor domains are
transmitted over $\sim$ 50 nm, connecting the motor head and the hinge, thus
providing an allosteric basis for LE.


  The ongoing global pandemic of infection disease COVID-19 caused by the 2019
novel coronavirus (SARS-COV-2, formerly 2019-nCoV) presents critical threats to
public health and the economy since it was identified in China, December 2019.
The genome of SARS-CoV-2 had been sequenced and structurally annotated, yet
little is known of the intrinsic organization and evolution of the genome. To
this end, we present a mathematical method for the genomic spectrum, a kind of
barcode, of SARS-CoV-2 and common human coronaviruses. The genomic spectrum is
constructed according to the periodic distributions of nucleotides, and
therefore reflects the unique characteristics of the genome. The results
demonstrate that coronavirus SARS-CoV-2 exhibits dinucleotide TT islands in the
non-structural proteins 3, 4, 5, and 6. Further analysis of the dinucleotide
regions suggests that the dinucleotide repeats are increased during evolution
and may confer the evolutionary fitness of the virus. The special dinucleotide
regions in the SARS-CoV-2 genome identified in this study may become diagnostic
and pharmaceutical targets in monitoring and curing the COVID-19 disease.


  Traditionally, the identification of parameters in the formulation and
solution of inverse problems considers that models, variables and mathematical
parameters are free of uncertainties. This aspect simplifies the estimation
process, but does not consider the influence of relatively small changes in the
design variables in terms of the objective function. In this work, the SIDR
(Susceptible, Infected, Dead and Recovered) model is used to simulate the
dynamic behavior of the novel coronavirus disease (COVID-19), and its
parameters are estimated by formulating a robust inverse problem, that is,
considering the sensitivity of design variables. For this purpose, a robust
multi-objective optimization problem is formulated, considering the
minimization of uncertainties associated to the estimation process and the
maximization of the robustness parameter. To solve this problem, the
Multi-objective Stochastic Fractal Search algorithm is associated with the
Effective Mean concept for the evaluation of robustness. The results obtained
considering real data of the epidemic in China demonstrate that the evaluation
of the sensitivity of the design variables can provide more reliable results.


  We derive both the finite and infinite population spatial replicator dynamics
as the fluid limit of a stochastic cellular automaton. The infinite population
spatial replicator is identical to the model used by Vickers and our derivation
justifies the addition of a diffusion to the replicator. The finite population
form generalizes the results by Durett and Levin on finite spatial replicator
games. We study the differences in the two equations as they pertain to a
one-dimensional rock-paper-scissors game.


  Molecular testing is rapidly becoming integral to the global tuberculosis
(TB) control effort. Uncommon mechanisms of resistance can escape detection by
these platforms and lead to the development of Multi-Drug Resistant (MDR)
strains. This article is a systematic review of published articles that
reported isoniazid (INH) resistance-conferring mutations between September-2013
and December-2019. The aims were to catalogue mutations associated with INH
resistance, estimate their global prevalence and co-occurrence, and their
utility in molecular diagnostics. The genes commonly associated with INH
resistance, katG, inhA, fabG1, and the intergenic region oxyR-ahpC were
considered in this review. In total, 52 articles were included describing 5,632
INHR clinical isolates from 31 countries. The three most frequently mutated
loci continue to be katG315 (4,100), inhA-15 (786), and inhA-8 (105). However,
the diagnostic value of inhA-8 is far lower than previously thought, only
appearing in 25 (0.4%) INHR isolates that lacked a mutation at the first two
loci. Importantly, of the four katG loci recommended by the previous systematic
review for diagnostics, only katG315 was observed in our INHR isolates. This
indicates continued evolution and regional differences in INH resistance. We
have identified 58 loci (common to both systematic reviews) in three genomic
regions as a reliable basis for molecular diagnostics. We also report 49 new
loci associated with INH resistance. Including all observed mutations provides
a cumulative sensitivity of 85.1%. The most disconcerting is the remaining
14.9% of isolates that harbor an unknown mechanism of resistance, will escape
molecular detection, and likely convert to MDR-TB, further complicating
treatment. Integrating the information cataloged in this and other similar
studies into current diagnostic tools is essential for combating the emergence
of MDR-TB.


  The World Health Organisation (WHO) has very strongly recommended testing and
isolation as a strategy for controlling the ongoing COVID-19 pandemic. The goal
of this paper is to quantify the effects of detection and isolation in formal
models of epidemics of varying complexity. A key parameter of such models is
the basic reproduction ratio. We show that an effective detection and isolation
strategy leads to a reduction of the basic reproduction ratio and can even lead
to this ratio becoming lower than one.


  The pandemic of COVID-19 has caused severe public health consequences around
the world. Many interventions of COVID-19 have been implemented. It is of great
public health and societal importance to evaluate the effects of interventions
in the pandemic of COVID-19. In this paper, with help of synthetic control
method, regression discontinuity and a Susceptible-Infected and infectious
without isolation-Hospitalized in isolation-Removed (SIHR) model, we evaluate
the horizontal and longitudinal effects of stringent interventions implemented
in Wenzhou, a representative urban city of China, where stringent interventions
were enforced to curb its own epidemic situation with rapidly increasing newly
confirmed cases. We found that there were statistically significant treatment
effects of those stringent interventions which reduced the cumulative confirmed
cases of COVID-19. Those reduction effects would increase over time. Also, if
the stringent interventions were delayed by 2 days or mild interventions were
implemented instead, the expected number of cumulative confirmed cases would
have been nearly 2 times or 5 times of the actual number. The effects of
stringent interventions are significant in mitigating the epidemic situation of
COVID-19. The slower the interventions were implemented, the more severe the
epidemic would have been, and the stronger the interventions would have been
required.


  We develop a multiple-events model and exploit within and between country
variation in the timing, type and level of intensity of various public policies
to study their dynamic effects on the daily incidence of COVID-19 and on
population mobility patterns across 135 countries. We remove concurrent policy
bias by taking into account the contemporaneous presence of multiple
interventions. The main result of the paper is that cancelling public events
and imposing restrictions on private gatherings followed by school closures
have quantitatively the most pronounced effects on reducing the daily incidence
of COVID-19. They are followed by workplace as well as stay-at-home
requirements, whose statistical significance and levels of effect are not as
pronounced. Instead, we find no effects for international travel controls,
public transport closures and restrictions on movements across cities and
regions. We establish that these findings are mediated by their effect on
population mobility patterns in a manner consistent with time-use and
epidemiological factors.


  Any epidemiological compartmental model with constant population is shown to
be a Hamiltonian dynamical system in which the total population plays the role
of the Hamiltonian function. Moreover, some particular cases within this large
class of models are shown to be bi-Hamiltonian. New interacting compartmental
models among different populations, which are endowed with a Hamiltonian
structure, are introduced. The Poisson structures underlying the Hamiltonian
description of all these dynamical systems are explicitly presented, and their
associated Casimir functions are shown to provide an efficient tool in order to
find exact analytical solutions for epidemiological models, such as the ones
describing the dynamics of the COVID-19 pandemic.


  Many countries are now investigating what the effects of summer vacation
might be on the COVID-19 pandemic. Here one particular such question is
addressed: what will happen if large numbers of metropolitan people visit a
less populated province during the summer vacation? By means of a simple
epidemic model, allowing for both short and long-term visitors to the province,
it is studied which features are most influential in determining if such summer
movements will result in large number of infections among the province
population. The method is applied to the island of Gotland off the South East
coast of Sweden. It is shown that the amount of mixing between the metropolitan
and province groups and the fraction of metropolitan people being infectious
upon arrival are most influential. Consequently, minimizing events gathering
both the province and metropolitan groups and/or reducing the number of
short-term visitors could substantially decrease spreading, as could measures
to lower the fraction initially infectious upon arrival.


  Presently, there are no approved drugs or vaccines to treat COVID-19 which
has spread to over 200 countries and is responsible for over 3,65,000 deaths
worldwide. Recent studies have shown that two human proteases, TMPRSS2 and
cathepsin L, play a key role in host cell entry of SARS-CoV-2. Importantly,
inhibitors of these proteases were shown to block SARS-CoV-2 infection. Here,
we perform virtual screening of 14010 phytochemicals produced by Indian
medicinal plants to identify natural product inhibitors of TMPRSS2 and
cathepsin L. We built a homology model of TMPRSS2 as an experimentally
determined structure is not available. AutoDock Vina was used to perform
molecular docking of phytochemicals against TMPRSS2 model structure and
cathepsin L crystal structure. Potential phytochemical inhibitors were filtered
by comparing their docked binding energies with those of known inhibitors of
TMPRSS2 and cathepsin L. Further, the ligand binding site residues and
non-covalent protein-ligand interactions were used as an additional filter to
identify phytochemical inhibitors that either bind to or form interactions with
residues important for the specificity of the target proteases. We have
identified 96 inhibitors of TMPRSS2 and 9 inhibitors of cathepsin L among
phytochemicals of Indian medicinal plants. The top inhibitors of TMPRSS2 are
Edgeworoside C, Adlumidine and Qingdainone, and of cathepsin L is Ararobinol.
Interestingly, several herbal sources of identified phytochemical inhibitors
have antiviral or anti-inflammatory use in traditional medicine. Further in
vitro and in vivo testing is needed before clinical trials of the promising
phytochemical inhibitors identified here.


  We introduce a deterministic model that partitions the total population into
the susceptible, infected, quarantined, and those traced after exposure, the
recovered and the deceased. We hypothesize 'accessible population for
transmission of the disease' to be a small fraction of the total population,
for instance when interventions are in force. This hypothesis, together with
the structure of the set of coupled nonlinear ordinary differential equations
for the populations, allows us to decouple the equations into just two
equations. This further reduces to a logistic type of equation for the total
infected population. The equation can be solved analytically and therefore
allows for a clear interpretation of the growth and inhibiting factors in terms
of the parameters in the full model. The validity of the 'accessible
population' hypothesis and the efficacy of the reduced logistic model is
demonstrated by the ease of fitting the United Kingdom data for the cumulative
infected and daily new infected cases. The model can also be used to forecast
further progression of the disease. In an effort to find optimized parameter
values compatible with the United Kingdom coronavirus data, we first determine
the relative importance of the various transition rates participating in the
original model. Using this we show that the original model equations provide a
very good fit with the United Kingdom data for the cumulative number of
infections and the daily new cases. The fact that the model calculated daily
new cases exhibits a turning point, suggests the beginning of a slow-down in
the spread of infections. However, since the rate of slowing down beyond the
turning point is small, the cumulative number of infections is likely to
saturate to about $3.52 \times 10^5$ around late July, provided the lock-down
conditions continue to prevail.


  Due to the intrinsic complexity and nonlinearity of chemical reactions,
direct applications of traditional machine learning algorithms may face with
many difficulties. In this study, through two concrete examples with biological
background, we illustrate how the key ideas of multiscale modeling can help to
reduce the computational cost of machine learning a lot, as well as how machine
learning algorithms perform model reduction automatically in a time-scale
separated system. Our study highlights the necessity and effectiveness of an
integration of machine learning algorithms and multiscale modeling during the
study of chemical reactions.


  In the present article a mathematical model is proposed to research the
current coronavirus outbreak (SARS-CoV-2) in Paraguay, describing the multiple
transmission paths in the infection dynamics, following the susceptible,
exposed, infectious, hospitalized, with lost immunity and recovered
individuals, as well as the role of the virus in the environment and deaths
from COVID-19 or other reasons. In order to reflect the impact of the control
measures adopted by the government and the population, the model employs
variable transmission rates that change with the epidemiological status and
environmental conditions. The model is validated, and its application is
demonstrated with data publicly available.


  Disease spread is a complex phenomenon requiring an interdisciplinary
approach. Covid-19 exhibited a global spatial spread in a very short time frame
resulting in a global pandemic. Data of web search effort in Greece on Covid-19
as a topic for one year on a weekly temporal scale were analyzed using
governmental intervention measures such a s school closures, movement
restrictions, national and international travelling restrictions, stay at home
requirements, mask requirements, financial support measures, and
epidemiological variables such as new cases and new deaths as potential
explanatory covariates. The relationship between web search effort on Covid-19
and the 16 in total explanatory covariates was analyzed with machine learning.
Web search in time was compared with the corresponding epidemiological
situation, expressed by the Rt at the same week. Results indicated that the
trained model exhibited a fit of R2 = 91% between the actual and predicted web
search effort. The top five variables for predicting web search effort were new
deaths, the opening of international borders to non-Greek nationals, new cases,
testing policy, and restrictions in internal movements. Web search peaked
during the same weeks that the Rt was peaking although new deaths or new cases
were not peaking during those dates, and Rt rarely is reported in public media.
As both web search effort and Rt peaked during 1-15 August 2020, the peak of
the tourist season, the implications of this are discussed.


  Learning accurate drug representation is essential for tasks such as
computational drug repositioning and prediction of drug side-effects. A drug
hierarchy is a valuable source that encodes human knowledge of drug relations
in a tree-like structure where drugs that act on the same organs, treat the
same disease, or bind to the same biological target are grouped together.
However, its utility in learning drug representations has not yet been
explored, and currently described drug representations cannot place novel
molecules in a drug hierarchy. Here, we develop a semi-supervised drug
embedding that incorporates two sources of information: (1) underlying chemical
grammar that is inferred from molecular structures of drugs and drug-like
molecules (unsupervised), and (2) hierarchical relations that are encoded in an
expert-crafted hierarchy of approved drugs (supervised). We use the Variational
Auto-Encoder (VAE) framework to encode the chemical structures of molecules and
use the knowledge-based drug-drug similarity to induce the clustering of drugs
in hyperbolic space. The hyperbolic space is amenable for encoding hierarchical
concepts. Both quantitative and qualitative results support that the learned
drug embedding can accurately reproduce the chemical structure and induce the
hierarchical relations among drugs. Furthermore, our approach can infer the
pharmacological properties of novel molecules by retrieving similar drugs from
the embedding space. We demonstrate that the learned drug embedding can be used
to find new uses for existing drugs and to discover side-effects. We show that
it significantly outperforms baselines in both tasks.


  The rapid spread of COVID-19 in the United States has imposed a major threat
to public health, the real economy, and human well-being. With the absence of
effective vaccines, the preventive actions of social distancing and travel
reduction are recognized as essential non-pharmacologic approaches to control
the spread of COVID-19. Prior studies demonstrated that human movement and
mobility drove the spatiotemporal distribution of COVID-19 in China. Little is
known, however, about the patterns and effects of co-location reduction on
cross-county transmission risk of COVID-19. This study utilizes Facebook
co-location data for all counties in the United States from March to early May
2020. The analysis examines the synchronicity and time lag between travel
reduction and pandemic growth trajectory to evaluate the efficacy of social
distancing in ceasing the population co-location probabilities, and
subsequently the growth in weekly new cases. The results show that the
mitigation effects of co-location reduction appear in the growth of weekly new
cases with one week of delay. Furthermore, significant segregation is found
among different county groups which are categorized based on numbers of cases.
The results suggest that within-group co-location probabilities remain stable,
and social distancing policies primarily resulted in reduced cross-group
co-location probabilities (due to travel reduction from counties with large
number of cases to counties with low numbers of cases). These findings could
have important practical implications for local governments to inform their
intervention measures for monitoring and reducing the spread of COVID-19, as
well as for adoption in future pandemics. Public policy, economic forecasting,
and epidemic modeling need to account for population co-location patterns in
evaluating transmission risk of COVID-19 across counties.


  The Coronavirus Disease 2019 (COVID-19) pandemic has infected over 10 million
people globally with a relatively high mortality rate. There are many
therapeutics undergoing clinical trials, but there is no effective vaccine or
therapy for treatment thus far. After affected by the Severe Acute Respiratory
Syndrome Coronavirus 2 (SARS-CoV-2), molecular signaling of host cells plays
critical roles during the life cycle of SARS-CoV-2. Thus, it is significant to
identify the involved molecular signaling pathways within the host cells, and
drugs targeting these molecular signaling pathways could be potentially
effective for COVID-19 treatment. In this study, we aimed to identify these
potential molecular signaling pathways, and repurpose existing drugs as a
potentially effective treatment of COVID-19 to facilitate the therapeutic
discovery, based on the transcriptional response of host cells. We first
identified dysfunctional signaling pathways associated with the infection
caused SARS-CoV-2 in human lung epithelial cells through analysis of the
altered gene expression profiles. In addition to the signaling pathway
analysis, the activated gene ontologies (GOs) and super gene ontologies were
identified. Signaling pathways and GOs such as MAPK, JNK, STAT, ERK, JAK-STAT,
IRF7-NFkB signaling, and MYD88/CXCR6 immune signaling were particularly
identified. Based on the identified signaling pathways and GOs, a set of
potentially effective drugs were repurposed by integrating the drug-target and
reverse gene expression data resources. The dexamethasone was top-ranked in the
prediction, which was the first reported drug to be able to significantly
reduce the death rate of COVID-19 patients receiving respiratory support. The
results can be helpful to understand the associated molecular signaling
pathways within host cells, and facilitate the discovery of effective drugs for
COVID-19 treatment.


  The coronavirus disease 2019 (COVID-19) pandemic is changing and impacting
lives on a global scale. In this paper, we introduce a mean-field game model in
controlling the propagation of epidemics on a spatial domain. The control
variable, the spatial velocity, is first introduced for the classical disease
models, such as the SIR model. For this proposed model, we provide fast
numerical algorithms based on proximal primal-dual methods. Numerical
experiments demonstrate that the proposed model illustrates how to separate
infected patients in a spatial domain effectively.


  We consider a model of cultural evolution for a strategy selection in a
population of individuals who interact in a game theoretic framework. The
evolution combines individual learning of the environment (population strategy
profile), reproduction, proportional to the success of the acquired knowledge,
and social transmission of the knowledge to the next generation. A mean-field
type equation is derived that describes the dynamics of the distribution of
cultural traits, in terms of the rate of learning, the reproduction rate and
population size. We establish global well-posedness of the initial-boundary
value problem for this equation and give several examples that illustrate the
process of the cultural evolution for some classical games.


  An epidemic carries human and fiscal costs. In the case of imported
pandemics, the first-best solution is to restrict national borders to identify
and isolate infected individuals. However, when that opportunity is not fully
seized and there is no preventative intervention available, second-best options
must be chosen. In this article we develop a system of differential equations
that simulate both the fiscal and human costs associated to different
mitigation measurements. After simulating several scenarios, we conclude that
herd immunity (or unleashing the pandemic) is the worst policy in terms of both
human and fiscal cost. We found that the second-best policy would be a strict
policy (e.g. physical distancing with massive testing) established under the
first 20 days after the pandemic, that lowers the probability of infection by
80%. In the case of the US, this strict policy would save more than 239
thousands lives and almost $170.8 billion to taxpayers when compared to the
herd immunity case.


  This research quantifies the various South African lockdown regimes, for the
SARS-CoV-2 pandemic, in terms of the basic reproduction number, $r_0$. It
further calculates the levels of immunity required for these selfsame lockdown
regimes to begin to work and predicts perceived values, should infections have
been underestimated by a factor of 10. The latter results are compelling. The
first, level-5 lockdown was a valiant attempt to contain the highly infectious,
SARS-CoV-2 virus, based on a limited knowledge. Its basic reproduction number
($r_0 = 1.93$) never came anywhere close to the requirement of being less than
unity. Obviously, it could be anticipated that the same would apply for
subsequent, lower levels of lockdown. The basic reproduction number for the
level-3 lockdown was found to be $2.34$ and that of the level-4 lockdown,
$1.69$.The suggestion is therefore that the level-4 lockdown might have been
marginally `smarter' than the `harder', level-5 lockdown, although its basic
reproduction number may merely reflect an adjustment by the public to the new
normal, or the ever-present error associated with data sets, in general. The
pandemic's basic reproduction number was calculated to be $3.16$, in the
Swedish context. The lockdowns therefore served to ensure that the medical
system was not overwhelmed, bought it valuable time to prepare and provided
useful data. The lockdowns nonetheless failed significantly in meeting any
objective to curtail the pandemic.


  The SIR pandemic model suffers from an unrealistic assumption: The rate of
removal from the infectious class of individuals is assumed to be proportional
to the number of infectious individuals. This means that a change in the rate
of infection is simultaneous with an equal change in the rate of removal. A
more realistic assumption is that an individual is removed at a certain time
interval after having been infected. A simple modified SIR model is proposed
which implements this delay, resulting in a single delay differential equation
which comprises the model. A solution to this equation which is applicable to a
pandemic is of the form A+B L(t) where L(t) is a logistic function, and A and B
are constants. While the classical SIR model is often an oversimplification of
pandemic behavior, it is instructive in that many of the fundamental dynamics
and descriptors of pandemics are clearly and simply defined. The logistic model
is generally used descriptively, dealing as it does with only the susceptible
and infected classes and the rate of transfer between them. The present model
presents a full but modified SIR model with a simpler logistic solution which
is more realistic and equally instructive.


  Why are the epidemic patterns of COVID-19 so different among different cities
or countries which are similar in their populations, medical infrastructures,
and people's behavior? Why are forecasts or predictions made by so-called
experts often grossly wrong, concerning the numbers of people who get infected
or die? The purpose of this study is to better understand the stochastic nature
of an epidemic disease, and answer the above questions. Much of the work on
infectious diseases has been based on "SIR deterministic models," (Kermack and
McKendrick:1927.) We will explore stochastic models that can capture the
essence of the seemingly erratic behavior of an infectious disease. A
stochastic model, in its formulation, takes into account the random nature of
an infectious disease.
  The stochastic model we study here is based on the "birth-and-death process
with immigration" (BDI for short), which was proposed in the study of
population growth or extinction of some biological species. The BDI process
model ,however, has not been investigated by the epidemiology community. The
BDI process is one of a few birth-and-death processes, which we can solve
analytically. Its time-dependent probability distribution function is a
"negative binomial distribution" with its parameter $r$ less than $1$. The
"coefficient of variation" of the process is larger than $\sqrt{1/r} > 1$.
Furthermore, it has a long tail like the zeta distribution. These properties
explain why infection patterns exhibit enormously large variations. The number
of infected predicted by a deterministic model is much greater than the median
of the distribution. This explains why any forecast based on a deterministic
model will fail more often than not.


  This document analyzes the role of data-driven methodologies in Covid-19
pandemic. We provide a SWOT analysis and a roadmap that goes from the access to
data sources to the final decision-making step. We aim to review the available
methodologies while anticipating the difficulties and challenges in the
development of data-driven strategies to combat the Covid-19 pandemic. A
3M-analysis is presented: Monitoring, Modelling and Making decisions. The focus
is on the potential of well-known datadriven schemes to address different
challenges raised by the pandemic: i) monitoring and forecasting the spread of
the epidemic; (ii) assessing the effectiveness of government decisions; (iii)
making timely decisions. Each step of the roadmap is detailed through a review
of consolidated theoretical results and their potential application in the
Covid-19 context. When possible, we provide examples of their applications on
past or present epidemics. We do not provide an exhaustive enumeration of
methodologies, algorithms and applications. We do try to serve as a bridge
between different disciplines required to provide a holistic approach to the
epidemic: data science, epidemiology, controltheory, etc. That is, we highlight
effective data-driven methodologies that have been shown to be successful in
other contexts and that have potential application in the different steps of
the proposed roadmap. To make this document more functional and adapted to the
specifics of each discipline, we encourage researchers and practitioners to
provide feedback. We will update this document regularly.


  We study a resource allocation problem for containing an infectious disease
in a metapopulation subject to resource uncertainty. We propose a two-stage
model where the policy maker seeks to allocate resources in both stages where
the second stage resource is random. Instead of a system of nonlinear
differential equations that governs the epidemic trajectories in the
constraints of the optimization model, we use a data-driven functional form to
model the cumulative number of infected individuals. This flexible data-driven
modeling choice allows us to transform the optimization problem to a tractable
mixed integer linear program. Our flexible approach can handle an online
decision making process, where the decision makers update their decisions for
opening treatment units and allocating beds utilizing the new information about
the epidemic progress. We utilize a detailed simulation model, validated by
real data from the 2014 Ebola epidemic in Sierra Leone. Our results show that
our policies produce about 400 fewer number of infected individuals in Sierra
Leone compared to the policies applied during the actual epidemic. We also
provide a detailed comparison of allocation policies generated by our
optimization framework which sheds light on the optimal resource allocation in
different regions.


  In this paper I show how path integral techniques can be used to put measures
on histories in "susceptible-infectious-recovered" (SIR)-type systems. The
standard SIR solution emerges as the classical saddle point of the action
describing the measure. One can then expand perturbatively around the
background solution, and this paper goes on to work out the covariance of
fluctuations around the background solution. Using a Green's function type
approach, one simply needs to solve additional ordinary differential equations;
an explicit matrix inversion is not required. The computed covariance matrix
should be useful in the construction of fast likelihoods for fitting the
parameters of SIR-type models to data. A comparison of the predictions of the
approach to an ensemble of simulations is presented.


  We present a compartmental SEIRD model aimed at forecasting hospital
occupancy in metropolitan areas during the current COVID-19 outbreak. The model
features asymptomatic and symptomatic infections with detailed hospital
dynamics. We model explicitly branching probabilities and non exponential
residence times in each latent and infected compartments. Using both hospital
admittance confirmed cases and deaths we infer the contact rate and the initial
conditions of the dynamical system, considering break points to model lockdown
interventions. Our Bayesian approach allows us to produce timely probabilistic
forecasts of hospital demand. The model has been used by the federal government
of Mexico to assist public policy, and has been applied for the analysis of
more than 70 metropolitan areas and the 32 states in the country.


  The prototype of a cyclic dominant system is the so-called
rock-scissors-paper game, but similar relation among competing strategies can
be identified in several other models of evolutionary game theory. In this work
we assume that a specific strategy from the available set is reluctant to adopt
alternative states, hence the related learning activity is reduced no matter
which other strategy is considered for adoption. Paradoxically, this
modification of the basic model will primarily elevate the stationary fraction
of another strategy who is the virtual predator of the one with reduced
learning activity. This general reaction of the studied systems is in agreement
with our understanding about Lotka-Volterra type cyclic dominant systems where
lowering the invasion rate between a source and target species promotes the
growth of former population. The observed effect is highly non-linear because
the effective invasion rates between strategies may depend sensitively on the
details of the actual model.


  Nowadays, one of the challenges we face when carrying out modeling of
epidemic spreading is to develop methods to control disease transmission. In
this article we study how the spreading of knowledge of a disease affects the
propagation of that disease in a population of interacting individuals. For
that, we analyze the interaction between two different processes on multiplex
networks: the propagation of an epidemic using the
susceptible-infected-susceptible dynamics and the dissemination of information
about the disease --and its prevention methods-- using the
unaware-aware-unaware dynamics, so that informed individuals are less likely to
be infected. Unlike previous related models where disease and information
spread at the same time scale, we introduce here a parameter that controls the
relative speed between the propagation of the two processes. We study the
behavior of this model using a mean-field approach that gives results in good
agreement with Monte Carlo simulations on homogeneous complex networks. We find
that increasing the rate of information dissemination reduces the disease
prevalence, as one may expect. However, increasing the speed of the information
process as compared to that of the epidemic process has the counter intuitive
effect of increasing the disease prevalence. This result opens an interesting
discussion about the effects of information spreading on disease propagation.


  Background: Network propagation has been widely used for nearly 20 years to
predict gene functions and phenotypes. Despite the popularity of this approach,
little attention has been paid to the question of provenance tracing in this
context, e.g., determining how much any experimental observation in the input
contributes to the score of every prediction.
  Results: We design a network propagation framework with two novel components
and apply it to predict human proteins that directly or indirectly interact
with SARS-CoV-2 proteins. First, we trace the provenance of each prediction to
its experimentally validated sources, which in our case are human proteins
experimentally determined to interact with viral proteins. Second, we design a
technique that helps to reduce the manual adjustment of parameters by users. We
find that for every top-ranking prediction, the highest contribution to its
score arises from a direct neighbor in a human protein-protein interaction
network. We further analyze these results to develop functional insights on
SARS-CoV-2 that expand on known biology such as the connection between
endoplasmic reticulum stress, HSPA5, and anti-clotting agents.
  Conclusions: We examine how our provenance tracing method can be generalized
to a broad class of network-based algorithms. We provide a useful resource for
the SARS-CoV-2 community that implicates many previously undocumented proteins
with putative functional relationships to viral infection. This resource
includes potential drugs that can be opportunistically repositioned to target
these proteins. We also discuss how our overall framework can be extended to
other, newly-emerging viruses.


  We present modeling of the COVID-19 epidemic in Illinois, USA, capturing the
implementation of a Stay-at-Home order and scenarios for its eventual release.
We use a non-Markovian age-of-infection model that is capable of handling long
and variable time delays without changing its model topology. Bayesian
estimation of model parameters is carried out using Markov Chain Monte Carlo
(MCMC) methods. This framework allows us to treat all available input
information, including both the previously published parameters of the epidemic
and available local data, in a uniform manner. To accurately model deaths as
well as demand on the healthcare system, we calibrate our predictions to total
and in-hospital deaths as well as hospital and ICU bed occupancy by COVID-19
patients. We apply this model not only to the state as a whole but also its
sub-regions in order to account for the wide disparities in population size and
density. Without prior information on non-pharmaceutical interventions (NPIs),
the model independently reproduces a mitigation trend closely matching mobility
data reported by Google and Unacast. Forward predictions of the model provide
robust estimates of the peak position and severity and also enable forecasting
the regional-dependent results of releasing Stay-at-Home orders. The resulting
highly constrained narrative of the epidemic is able to provide estimates of
its unseen progression and inform scenarios for sustainable monitoring and
control of the epidemic.


  Accurate forecasts for COVID-19 are necessary for better preparedness and
resource management. Specifically, deciding the response over months or several
months requires accurate long-term forecasts which is particularly challenging
as the model errors accumulate with time. A critical factor that can hinder
accurate long-term forecasts, is the number of unreported/asymptomatic cases.
While there have been early serology tests to estimate this number, more tests
need to be conducted for more reliable results. To identify the number of
unreported/asymptomatic cases, we take an epidemiology data-driven approach. We
show that we can identify lower bounds on this ratio or upper bound on actual
cases as a factor of reported cases. To do so, we propose an extension of our
prior heterogeneous infection rate model, incorporating unreported/asymptomatic
cases. We prove that the number of unreported cases can be reliably estimated
only from a certain time period of the epidemic data. In doing so, we construct
an algorithm called Fixed Infection Rate method, which identifies a reliable
bound on the learned ratio. We also propose two heuristics to learn this ratio
and show their effectiveness on simulated data. We use our approaches to
identify the upper bounds on the ratio of actual to reported cases for New York
City and several US states. Our results demonstrate with high confidence that
the actual number of cases cannot be more than 35 times in New York, 40 times
in Illinois, 38 times in Massachusetts and 29 times in New Jersey, than the
reported cases.


  We analyze the effect of using a screening CT-scan for evaluation of
potential COVID-19 infections in order to isolate and perform contact tracing
based upon a viral pneumonia diagnosis. RT-PCR is then used for continued
isolation based upon a COVID diagnosis. Both the low false negative rates and
rapid results of CT-scans lead to dramatically reduced transmission. The
reduction in cases after 60 days with widespread use of CT-scan screening
compared to PCR by itself is as high as $50\times$, and the reduction of
effective reproduction rate $R(t)$ is $0.20$. Our results imply that much more
rapid extinction of COVID is possible by combining social distancing with
CT-scans and contact tracing.


  This study analyzed the spread and decay durations of the COVID-19 pandemic
in different prefectures of Japan. During the pandemic, affordable healthcare
was widely available in Japan and the medical system did not suffer a collapse,
making accurate comparisons between prefectures possible. For the 16
prefectures included in this study that had daily maximum confirmed cases
exceeding ten, the number of daily confirmed cases follow bell-shape or
log-normal distribution in most prefectures. A good correlation was observed
between the spread and decay durations. However, some exceptions were observed
in areas where travelers returned from foreign countries, which were defined as
the origins of infection clusters. Excluding these prefectures, the population
density was shown to be a major factor affecting the spread and decay patterns,
with R2=0.39 (p<0.05) and 0.42 (p<0.05), respectively, approximately
corresponding to social distancing. The maximum absolute humidity was found to
affect the decay duration normalized by the population density (R2>0.36, p
<0.05). Our findings indicate that the estimated pandemic spread duration,
based on the multivariate analysis of maximum absolute humidity, ambient
temperature, and population density (adjusted R2=0.53, p-value<0.05), could
prove useful for intervention planning during potential future pandemics,
including a second COVID-19 outbreak.


  A key characteristic of the spread of infectious diseases is their ability to
use efficient transmission paths within contact graphs. This enables the
pathogen to maximise infection rates and spread within a target population. In
this work, we devise techniques to localise infections and decrease infection
rates based on a principled analysis of disease transmission paths within
human-contact networks (proximity graphs). Experimental results of disease
spreading shows that that at low visibility rates contact tracing slows disease
spreading. However to stop disease spreading, contact tracing requires both
significant visibility (at least 60%) into the proximity graph and the ability
to place half of the population under isolation. We find that pro-actively
isolating super-links -- key proximity encounters -- has significant benefits:
targeted isolation of a fourth of the population based on 35% visibility into
the proximity graph prevents an epidemic outbreak. It turns out that isolating
super-spreaders is more effective than contact tracing and testing but less
effective than targeting super-links. We highlight the important role of
topology in epidemic outbreaks. We argue that proactive innoculation of a
population by disabling super-links and super-spreaders may have an important
complimentary role alongside contact tracing and testing as part of a
sophisticated public-health response to epidemic outbreaks.


  Genome-scale orthology assignments are usually based on reciprocal best
matches. In the absence of horizontal gene transfer (HGT), every pair of
orthologs forms a reciprocal best match. Incorrect orthology assignments
therefore are always false positives in the reciprocal best match graph. We
consider duplication/loss scenarios and characterize unambiguous false-positive
(u-fp) orthology assignments, that is, edges in the best match graphs (BMGs)
that cannot correspond to orthologs for any gene tree that explains the BMG.
Moreover, we provide a polynomial-time algorithm to identify all u-fp orthology
assignments in a BMG. Simulations show that at least $75\%$ of all incorrect
orthology assignments can be detected in this manner. All results rely only on
the structure of the BMGs and not on any a priori knowledge about underlying
gene or species trees.


  Chemical reaction networks describe interactions between biochemical species.
Once an underlying reaction network is given for a biochemical system, the
system dynamics can be modelled with various mathematical frameworks such as
continuous time Markov processes. In this manuscript, the identifiability of
the underlying network structure with a given stochastic system dynamics is
studied. It is shown that some data types related to the associated stochastic
dynamics can uniquely identify the underlying network structure as well as the
system parameters. The accuracy of the presented network inference is
investigated when given dynamical data is obtained via stochastic simulations.


  The probability of two loci, separated by a certain genome length, being in
contact can be inferred using the Chromosome Conformation Capture (3C) method
and related Hi-C experiments. How to go from the contact map, a matrix listing
the mean contact probabilities between a large number of pairs of loci, to an
ensemble of three-dimensional structures is an open problem. A solution to this
problem, without assuming an assumed energy function, would be the first step
in understanding the way nature has solved the packaging of chromosomes in
tight cellular spaces. We created a theory, based on polymer physics
characteristics of chromosomes and the maximum entropy principles, referred to
as HIPPS (Hi-C-Polymer-Physics-Structures) method, that allows us to calculate
the 3D structures solely from Hi-C contact maps. We created an ensemble of 3D
structures for the 23 chromosomes from lymphoblastoid cells using the measured
contact maps as inputs. The HIPPS method shows that conformations of
chromosomes are heterogeneous even in a single cell type. The differences in
the conformational heterogeneity of the same chromosome in different cell types
(normal as well as cancerous cells) can also be quantitatively discerned using
our theory. We validate the method by showing that the calculated volumes of
the 23 chromosomes from the predicted 3D structures are in good agreement with
experimental estimates. Because the method is general, the 3D structures for
any species may be calculated directly from the contact map without the need to
assume a specific polymer model, as is customarily done.


  We study the target control problem of asynchronous Boolean networks, to
identify a set of nodes, the perturbation of which can drive the dynamics of
the network from any initial state to the desired steady state (or attractor).
We are particularly interested in temporary perturbations, which are applied
for sufficient time and then released to retrieve the original dynamics.
Temporary perturbations have the apparent advantage of averting unforeseen
consequences, which might be induced by permanent perturbations. Despite the
infamous state-space explosion problem, in this work, we develop an efficient
method to compute the temporary target control for a given target attractor of
a Boolean network. We apply our method to a number of real-life biological
networks and compare its performance with the stable motif-based control method
to demonstrate its efficacy and efficiency.


  We present a simple technique to compare the development of the Covid-19
epidemic in different regions, based only on the time series of confirmed
cases. Weekly new infections, taken for every day, are interpreted as infection
potential of Covid-19. We derive a robust time-varying reproduction rate for
the infection potential, including asymptomatic cases, which does not depend on
death rate or testing intensity. It requires few assumptions and shows a more
plausible time course than official reproduction rates in several countries.


  It has been known that the uptake mechanisms of cell-penetrating
peptides(CPPs) depend on the experimental conditions such as concentration of
peptides, lipid composition, temperature, etc. In this study we investigate the
temperature dependence of the penetration of Arg9s into a DOPC/DOPG(4:1)
membrane using molecular dynamics(MD) simulations at two different
temperatures, T = 310 K and T = 288 K. Although it is difficult to identify the
temperature dependence because of having only one single simulation at each
temperature and no evidence of translocation of Arg9s across the membrane at
both temperatures, our simulations suggest that followings are strongly
correlated with the penetration of Arg9s: a number of water molecules
coordinated by Arg9s, electrostatic energy between Arg9s and the lipids
molecules. We also present how Arg9s change a bending rigidity of the membrane
and how a collective behavior between Arg9s enhances the penetration and the
membrane bending. Our analyses can be applicable to any cell-penetrating
peptides(CPPs) to investigate their interactions with various membranes using
MD simulations.


  Researchers across the globe are seeking to rapidly repurpose existing drugs
or discover new drugs to counter the the novel coronavirus disease (COVID-19)
caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). One
promising approach is to train machine learning (ML) and artificial
intelligence (AI) tools to screen large numbers of small molecules. As a
contribution to that effort, we are aggregating numerous small molecules from a
variety of sources, using high-performance computing (HPC) to computer diverse
properties of those molecules, using the computed properties to train ML/AI
models, and then using the resulting models for screening. In this first data
release, we make available 23 datasets collected from community sources
representing over 4.2 B molecules enriched with pre-computed: 1) molecular
fingerprints to aid similarity searches, 2) 2D images of molecules to enable
exploration and application of image-based deep learning methods, and 3) 2D and
3D molecular descriptors to speed development of machine learning models. This
data release encompasses structural information on the 4.2 B molecules and 60
TB of pre-computed data. Future releases will expand the data to include more
detailed molecular simulations, computed models, and other products.


  This study develops a framework for quantification of the impact of changes
in population mobility due to social distancing on the COVID-19 infection
growth rate. Using the Susceptible-Infected-Recovered (SIR) epidemiological
model we establish that under some mild assumptions the growth rate of COVID-19
deaths is a time-delayed approximation of the growth rate of COVID-19
infections. We then hypothesize that the growth rate of COVID-19 infections is
a function of population mobility, which leads to a statistical model that
predicts the growth rate of COVID-19 deaths as a delayed function of population
mobility. The parameters of the statistical model directly reveal the growth
rate of infections, the mobility-dependent transmission rate, the
mobility-independent recovery rate, and the critical mobility, below which
COVID-19 growth rate becomes negative. We fitted the proposed statistical model
on publicly available data from 14 countries where daily death counts exceeded
100 for more than 3 days as of May 6th, 2020. The publicly available Google
Mobility Index (GMI) was used as a measure of population mobility at the
country level. Our results show that the growth rate of COVID-19 deaths can be
accurately estimated 20 days ahead as a quadratic function of the transit
category of GMI (adjusted R-squared = 0.784). The estimated 95% confidence
interval for the critical mobility is in the range between 36.1% and 47.6% of
the pre-COVID-19 mobility. This result indicates that a significant reduction
in population mobility is needed to reverse the growth of COVID-19 epidemic.
Moreover, the quantitative relationship established herein suggests that a
readily available, population-level metric such as GMI can be a useful
indicator of the course of COVID-19 epidemic.


  The present article studies the extension of two deterministic models for
describing the novel coronavirus pandemic crisis, the SIR model and the SEIR
model. The models were studied and compared to real data in order to support
the validity of each description and extract important information regarding
the pandemic, such as the basic reproductive number R0, which might provide
useful information concerning the rate of increase of the pandemic predicted by
each model. We next proceed to making predictions and comparing more complex
models derived from the SEIR model with the SIRD model, in order to find the
most suitable one for describing and predicting the pandemic crisis. Aiming to
answer the question if the simple SIRD model is able to make reliable
predictions and deliver suitable information compared to more complex models.


  Once an epidemic outbreak has been effectively contained through
non-pharmaceutical interventions, a safe protocol is required for the
subsequent release of social distancing restrictions to prevent a disastrous
resurgence of the infection. We report individual-based numerical simulations
of stochastic susceptible-infectious-recovered model variants on four distinct
spatially organized lattice and network architectures wherein contact and
mobility constraints are implemented. We robustly find that the intensity and
spatial spread of the epidemic recurrence wave can be limited to a manageable
extent provided release of these restrictions is delayed sufficiently (for a
duration of at least thrice the time until the peak of the unmitigated
outbreak) and long-distance connections are maintained on a low level (limited
to less than five percent of the overall connectivity).


  We analyze the dynamics of the COVID-19 epidemic taking into account the role
of the unreported cases. After a first section in which we deal with a
framework of very slow test capacity, we turn to the model recently
introduced/implemented by Liu, Magal, Seydi and Webb. First, we prove some
basic structural results for the corresponding ODE, as for instance the
convergence of S(t) to a positive limit. These are similar to those of the
classical SIR model, although the maxima of the corresponding curves are not
necessarily unique. Finally, we implement the model -- but with a variable
transmission rate -- in the Chilean context. A key parameter adjustment
(namely, the fraction of unreported cases) is done via an argument using
mortality rates. We conclude with several conclusions and lines of future
research.


  After their re-emergence in the last decades, dengue fever and other
vector-borne diseases are a potential threat to the lives of millions of
people. Based on a data set of dengue cases in the Brazilian city of Fortaleza,
collected from 2011 to 2016, we study the spatio-temporal characteristics of
dengue outbreaks to characterize epidemic and non-epidemic years. First, we
identify regions that show a high prevalence of dengue cases and mosquito
larvae in different years and also analyze their corresponding correlations.
Our results show that the characteristic correlation length of the epidemic is
of the order of the system size, suggesting that factors such as citizen
mobility may play a major role as a drive for spatial spreading of vector-borne
diseases. Inspired by this observation, we perform a mean-field estimation of
the basic reproduction number and find that our estimated values agree well
with the values reported for other regions, pointing towards similar underlying
spreading mechanisms. These findings provide insights into the spreading
characteristics of dengue in densely populated areas and should be of relevance
for the design of improved disease containment strategies.


  Computer-aided design for synthetic biology promises to accelerate the
rational and robust engineering of biological systems; it requires both
detailed and quantitative mathematical and experimental models of the processes
to (re)design, and software and tools for genetic engineering and DNA assembly.
Ultimately, the increased precision in the design phase will have a dramatic
impact on the production of designer cells and organisms with bespoke functions
and increased modularity. Computer-aided design strategies require quantitative
representations of cells, able to capture multiscale processes and link
genotypes to phenotypes. Here, we present a perspective on how whole-cell,
multiscale models could transform design-build-test-learn cycles in synthetic
biology. We show how these models could significantly aid in the design and
learn phases while reducing experimental testing by presenting case studies
spanning from genome minimization to cell-free systems, and we discuss several
challenges for the realization of our vision. The possibility to describe and
build in silico whole-cells offers an opportunity to develop increasingly
automatized, precise and accessible computer-aided design tools and strategies
throughout novel interdisciplinary collaborations.


  SARS-CoV-2 has disrupted the life of billions of people around the world
since the first outbreak was officially declared in China at the beginning of
2020. Yet, important questions such as how deadly it is or its degree of spread
within different countries remain unanswered. In this work, we exploit the
`universal' growth of the mortality rate with age observed in different
countries since the beginning of their respective outbreaks, combined with the
results of the antibody prevalence tests in the population of Spain, to unveil
both unknowns. We validate these results with an analogous antibody rate survey
in the canton of Geneva, Switzerland. We also argue that the official number of
deaths over 70 years old is importantly underestimated in most of the
countries, and we use the comparison between the official records with the
number of deaths mentioning COVID-19 in the death certificates to quantify by
how much. Using this information, we estimate the fatality infection ratio
(IFR) for the different age segments and the fraction of the population
infected in different countries assuming a uniform exposure to the virus in all
age segments. We also give estimations for the non-uniform IFR using the
sero-epidemiological results of Spain, showing a very similar growth of the
fatality ratio with age. Only for Spain, we estimate the probability (if
infected) of being identified as a case, being hospitalized or admitted in the
intensive care units as function of age. In general, we observe a nearly
exponential growth of the fatality ratio with age, which anticipates large
differences in total IFR in countries with different demographic distributions,
with numbers that range from 1.82\% in Italy, to 0.62\% in China or even 0.14\%
in middle Africa.


  SARS-CoV-2 coronavirus infection is mediated by the binding of its spike
protein to the angiotensin-converting enzyme 2 (ACE2), which plays a pivotal
role in the renin-angiotensin system (RAS). The study of RAS dysregulation due
to SARS-CoV-2 infection is fundamentally important for a better understanding
of the pathogenic mechanisms and risk factors associated with COVID-19
coronavirus disease, and to design effective therapeutic strategies. In this
context, we developed a mathematical model of RAS based on data regarding
protein and peptide concentrations; the model was tested on clinical data from
healthy normotensive and hypertensive individuals. We then used our model to
analyze the impact of SARS-CoV-2 infection on RAS, which we modeled through a
down-regulation of ACE2 as a function of viral load. We also used it to predict
the effect of RAS-targeting drugs, such as RAS-blockers, human recombinant
ACE2, and angiotensin 1-7 peptide, on COVID-19 patients; the model predicted an
improvement of the clinical outcome for some drugs and a worsening for others.


  In December 2019, a newly discovered SARS-CoV-2 virus was emerged from China
and propagated worldwide as a pandemic. In the absence of preventive medicine
or a ready to use vaccine, mathematical models can provide useful scientific
insights about transmission patterns and targets for drug development. In this
study, we propose a within-host mathematical model of SARS-CoV-2 infection
considering innate and adaptive immune responses. We analyze the equilibrium
points of the proposed model and obtain an expression of the basic reproduction
number. We then numerically show the existence of a transcritical bifurcation.
The proposed model is calibrated to real viral load data of two COVID-19
patients. Using the estimated parameters, we perform global sensitivity
analysis with respect to the peak of viral load. Finally, we study the efficacy
of antiviral drugs and vaccination on the dynamics of SARS-CoV-2 infection. Our
results suggest that blocking the production of the virus by infected cells
decreases the viral load more than reducing the infection rate of healthy
cells. Vaccination is also found useful but during the vaccine development
phase, blocking virus production from infected cells can be targeted for
antiviral drug development.


  The interpretation of sampling data plays a crucial role in policy response
to the spread of a disease during an epidemic, such as the COVID-19 epidemic of
2020. However, this is a non-trivial endeavor due to the complexity of real
world conditions and limits to the availability of diagnostic tests, which
necessitate a bias in testing favoring symptomatic individuals. A thorough
understanding of sampling confidence and bias is necessary in order make
accurate conclusions. In this manuscript, we provide a stochastic model of
sampling for assessing confidence in disease metrics such as trend detection,
peak detection, and disease spread estimation. Our model simulates testing for
a disease in an epidemic with known dynamics, allowing us to use Monte-Carlo
sampling to assess metric confidence. This model can provide realistic
simulated data which can be used in the design and calibration of data analysis
and prediction methods. As an example, we use this method to show that trends
in the disease may be identified using under $10000$ biased samples each day,
and an estimate of disease spread can be made with additional $1000-2000$
unbiased samples each day. We also demonstrate that the model can be used to
assess more advanced metrics by finding the precision and recall of a strategy
for finding peaks in the dynamics.


  To slow down the spread of Covid-19, administrative regions within Pakistan
imposed complete and partial lockdown restrictions on socio-economic
activities, religious congregations, and human movement. Here we examine the
impact of regional lockdown strategies on Covid-19 outcomes. After conducting
econometric analyses (Regression Discontinuity and Negative Binomial
Regressions) on official data from the National Institute of Health (NIH)
Pakistan, we find that the strategies did not lead to a similar level of
Covid-19 caseload (positive cases and deaths) in all regions. In terms of
reduction in the overall caseload (positive cases and deaths), compared to no
lockdown, complete and partial lockdown appeared to be effective in four
regions: Balochistan, Gilgit Baltistan (GT), Islamabad Capital Territory (ICT),
and Azad Jammu and Kashmir (AJK). Contrarily, complete and partial lockdowns
did not appear to be effective in containing the virus in the three largest
provinces of Punjab, Sindh, and Khyber Pakhtunkhwa (KPK). The observed regional
heterogeneity in the effectiveness of lockdowns advocates for a careful use of
lockdown strategies based on the demographic, social, and economic factors.


  According to the current perception, symptomatic, presymptomatic, and
asymptomatic infectious persons can infect the healthy population susceptible
to the SARS-Cov-2. More importantly, various reports indicate that the number
of asymptomatic cases can be several-fold higher than the reported symptomatic
cases. In this article, we take the reported cases in India and various states
within the country till September 1, as the specimen to understand the
progression of the COVID-19. Employing a modified SEIRD model, we predict the
spread of COVID-19 by the symptomatic as well as asymptomatic infectious
population. Considering reported infection primarily due to symptomatic we
compare the model predicted results with the available data to estimate the
dynamics of the asymptomatically infected population. Our data indicate that in
the absence of the asymptomatic infectious population, the number of
symptomatic cases would have been much less. Therefore, the current progress of
the symptomatic infection can be reduced by quarantining the asymptomatically
infectious population via extensive or random testing. This study is motivated
strictly towards academic pursuit; this theoretical investigation is not meant
for influencing policy decisions or public health practices.


  Single-molecule F\"{o}rster resonance energy transfer (smFRET) has become a
mainstream technique for probing biomolecular structural dynamics. The rapid
and wide adoption of the technique by an ever-increasing number of groups has
generated many improvements and variations in the technique itself, in methods
for sample preparation and characterization, in analysis of the data from such
experiments, and in analysis codes and algorithms. Recently, several labs that
employ smFRET have joined forces to try to bring the smFRET community together
in adopting a consensus on how to perform experiments and analyze results for
achieving quantitative structural information. These recent efforts include
multi-lab blind-tests to assess the accuracy and precision of smFRET between
different labs using different procedures, the formal assembly of the FRET
community and development of smFRET procedures to be considered for entries in
the wwPDB. Here we delve into the different approaches and viewpoints in the
field. This position paper describes the current "state-of-the field", points
to unresolved methodological issues for quantitative structural studies,
provides a set of 'soft recommendations' about which an emerging consensus
exists, and a list of resources that are openly available. To make further
progress, we strongly encourage 'open science' practices. We hope that this
position paper will provide a roadmap for newcomers to the field, as well as a
reference for seasoned practitioners.


  The novelty of new human coronavirus COVID-19/SARS-CoV-2 and the lack of
effective drugs and vaccines gave rise to a wide variety of strategies employed
to fight this worldwide pandemic. Many of these strategies rely on the
repositioning of existing drugs that could shorten the time and reduce the cost
compared to de novo drug discovery. In this study, we presented a new
network-based algorithm for drug repositioning, called SAveRUNNER (Searching
off-lAbel dRUg aNd NEtwoRk), which predicts drug-disease associations by
quantifying the interplay between the drug targets and the disease-specific
proteins in the human interactome via a novel network-based similarity measure
that prioritizes associations between drugs and diseases locating in the same
network neighborhoods. Specifically, we applied SAveRUNNER on a panel of 14
selected diseases with a consolidated knowledge about their disease-causing
genes and that have been found to be related to COVID-19 for genetic
similarity, comorbidity, or for their association to drugs tentatively
repurposed to treat COVID-19. Focusing specifically on SARS subnetwork, we
identified 282 repurposable drugs, including some the most rumored off-label
drugs for COVID-19 treatments, as well as a new combination therapy of 5 drugs,
actually used in clinical practice. Furthermore, to maximize the efficiency of
putative downstream validation experiments, we prioritized 24 potential
anti-SARS-CoV repurposable drugs based on their network-based similarity
values. These top-ranked drugs include ACE-inhibitors, monoclonal antibodies,
and thrombin inhibitors. Finally, our findings were in-silico validated by
performing a gene set enrichment analysis, which confirmed that most of the
network-predicted repurposable drugs may have a potential treatment effect
against human coronavirus infections.


  We formulate a generalized susceptible exposed infectious recovered (SEIR)
model on a graph, describing the population dynamics of an open crowded place
with an arbitrary topology. As a sample calculation, we discuss three simple
cases, both analytically, and numerically, by means of a cellular automata
simulation of the individual dynamics in the system. As a result, we provide
the infection ratio in the system as a function of controllable parameters,
which allows for quantifying how acting on the human behavior may effectively
lower the disease spread throughout the system.


  Residential colleges and universities face unique challenges in providing
in-person instruction during the COVID-19 pandemic. Administrators are
currently faced with decisions about whether to open during the pandemic and
what modifications of their normal operations might be necessary to protect
students, faculty and staff. There is little information, however, on what
measures are likely to be most effective and whether existing interventions
could contain the spread of an outbreak on campus. We develop a full-scale
stochastic agent-based model to determine whether in-person instruction could
safely continue during the pandemic and evaluate the necessity of various
interventions. Simulation results indicate that large scale randomized testing,
contact-tracing, and quarantining are important components of a successful
strategy for containing campus outbreaks. High test specificity is critical for
keeping the size of the quarantine population manageable. Moving the largest
classes online is also crucial for controlling both the size of outbreaks and
the number of students in quarantine. Increased residential exposure can
significantly impact the size of an outbreak, but it is likely more important
to control non-residential social exposure among students. Finally, necessarily
high quarantine rates even in controlled outbreaks imply significant
absenteeism, indicating a need to plan for remote instruction of quarantined
students.


  Infection can spread easily on networks with heterogeneous degree
distribution. Here, we considered targeted immunization on such networks,
wherein a fraction of individuals with the highest connectivity are immunized.
To quantify the effect of this targeted immunization approach on population
immunity, we proposed a method using the type reproduction number.
Consequently, we derived a precise and simple formula that can yield the
immunization threshold, which, to the best of our knowledge, is the first such
result presented in the literature.


  Metabolism plays a central role in cell physiology because it provides the
molecular machinery for growth. At the genome-scale, metabolism is made up of
thousands of reactions interacting with one another. Untangling this complexity
is key to understand how cells respond to genetic, environmental, or
therapeutic perturbations. Here we discuss the roles of two complementary
strategies for the analysis of genome-scale metabolic models: Flux Balance
Analysis (FBA) and network science. While FBA estimates metabolic flux on the
basis of an optimisation principle, network approaches reveal emergent
properties of the global metabolic connectivity. We highlight how the
integration of both approaches promises to deliver insights on the structure
and function of metabolic systems with wide-ranging implications in discovery
science, precision medicine and industrial biotechnology.


  The nation-wide lockdown starting 25 March 2020, aimed at suppressing the
spread of the COVID-19 disease, was extended until 31 May 2020 in three
subsequent orders by the Government of India. The extended lockdown has had
significant social and economic consequences and `lockdown fatigue' has likely
set in. Phased reopening began from 01 June 2020 onwards. Mumbai, one of the
most crowded cities in the world, has witnessed both the largest number of
cases and deaths among all the cities in India (41986 positive cases and 1368
deaths as of 02 June 2020). Many tough decisions are going to be made on
re-opening in the next few days. In an earlier IISc-TIFR Report, we presented
an agent-based city-scale simulator(ABCS) to model the progression and spread
of the infection in large metropolises like Mumbai and Bengaluru. As discussed
in IISc-TIFR Report 1, ABCS is a useful tool to model interactions of city
residents at an individual level and to capture the impact of
non-pharmaceutical interventions on the infection spread. In this report we
focus on Mumbai. Using our simulator, we consider some plausible scenarios for
phased emergence of Mumbai from the lockdown, 01 June 2020 onwards. These
include phased and gradual opening of the industry, partial opening of public
transportation (modelling of infection spread in suburban trains), impact of
containment zones on controlling infections, and the role of compliance with
respect to various intervention measures including use of masks, case
isolation, home quarantine, etc. The main takeaway of our simulation results is
that a phased opening of workplaces, say at a conservative attendance level of
20 to 33\%, is a good way to restart economic activity while ensuring that the
city's medical care capacity remains adequate to handle the possible rise in
the number of COVID-19 patients in June and July.


  The effect of an oscillating electric field generated from music on yeast
vacuolar proton-ATPase (V-ATPase) activity in its native environment is
reported. An oscillating electric field is generated by electrodes that are
immersed into a dispersion of yeast vacuolar membrane vesicles natively hosting
a high concentration of active V-ATPase. The substantial difference in the ATP
hydrolysing activity of V-ATPase under the most stimulating and inhibiting
music is unprecedented. Since the topic, i.e. an effect of music on
biomolecules, is very attractive for non-scientific, esoteric mystification, we
provide a rational explanation for the observed new phenomenon. Good
correlation is found between changes in the specific activity of the enzyme and
the combined intensity of certain frequency bands of the Fourier spectra of the
music clips. Most prominent identified frequencies are harmonically related to
each other and to the estimated rotation rate of the enzyme. These results lead
to the conclusion that the oscillating electric field interferes with periodic
trans-membrane charge motions in the working enzyme.


  Solar UV$-$C photons do not reach Earth's surface, but are known to be
endowed with germicidal properties that are also effective on viruses. The
effect of softer UV$-$B and UV$-$A photons, which copiously reach the Earth's
surface, on viruses are instead little studied, particularly on
single$-$stranded RNA viruses. Here we combine our measurements of the action
spectrum of Covid$-$19 in response to UV light, Solar irradiation measurements
on Earth during the SARS$-$CoV$-$2 pandemics, worldwide recorded Covid$-$19
mortality data and our 'Solar$-$Pump' diffusive model of epidemics to show that
(a) UV$-$B$/$A photons have a powerful virucidal effect on the
single$-$stranded RNA virus Covid$-$19 and that (b) the Solar radiation that
reaches temperate regions of the Earth at noon during summers, is sufficient to
inactivate 63\perc of virions in open$-$space concentrations (1.5 x 103
TCID50$/$mL, higher than typical aerosol) in less than 2 min. We conclude that
the characteristic seasonality imprint displayed world$-$wide by the
SARS$-$Cov$-$2 mortality time$-$series throughout the diffusion of the outbreak
(with temperate regions showing clear seasonal trends and equatorial regions
suffering, on average, a systematically lower mortality), might have been
efficiently set by the different intensity of UV$-$B$/$A Solar radiation
hitting different Earth's locations at different times of the year. Our results
suggest that Solar UV$-$B$/$A play an important role in planning strategies of
confinement of the epidemics, which should be worked out and set up during
spring$/$summer months and fully implemented during low$-$solar$-$irradiation
periods.


  More than one billion students are out of school because of Covid-19, forced
to a remote learning that has several drawbacks and has been hurriedly
arranged; in addition, most countries are currently uncertain on how to plan
school activities for the 2020-2021 school year; all of this makes learning and
education some of the biggest world issues of the current pandemic.
Unfortunately, due to the length of the incubation period of Covid-19, full
opening of schools seems to be impractical till a vaccine is available. In
order to support the possibility of some in-person learning, we study a
mathematical model of the diffusion of the epidemic due to school opening, and
evaluate plans aimed at containing the extra Covid-19 cases due to school
activities while ensuring an adequate number of in-class learning periods. We
consider a SEAIR model with an external source of infection and a suitable loss
function; after a realistic parameter selection, we numerically determine
optimal school opening strategies by simulated annealing. It turns out that
blended models, with almost periodic alternations of in-class and remote
teaching days or weeks, are generally optimal. Besides containing Covid-19
diffusion, these solutions could be pedagogically acceptable, and could also
become a driving model for the society at large. In a prototypical example, the
optimal strategy results in the school opening 90 days out of 200 with the
number of Covid-19 cases among the individuals related to the school increasing
by about 66%, instead of the about 250% increase that would have been a
consequence of full opening.


  Given the severity of the SARS-CoV-2 pandemic, a major challenge is to
rapidly repurpose existing approved drugs for clinical interventions. While a
number of data-driven and experimental approaches have been suggested in the
context of drug repurposing, a platform that systematically integrates
available transcriptomic, proteomic and structural data is missing. More
importantly, given that SARS-CoV-2 pathogenicity is highly age-dependent, it is
critical to integrate aging signatures into drug discovery platforms. We here
take advantage of large-scale transcriptional drug screens combined with
RNA-seq data of the lung epithelium with SARS-CoV-2 infection as well as the
aging lung. To identify robust druggable protein targets, we propose a
principled causal framework that makes use of multiple data modalities. Our
analysis highlights the importance of serine/threonine and tyrosine kinases as
potential targets that intersect the SARS-CoV-2 and aging pathways. By
integrating transcriptomic, proteomic and structural data that is available for
many diseases, our drug discovery platform is broadly applicable. Rigorous in
vitro experiments as well as clinical trials are needed to validate the
identified candidate drugs.


  We propose a new method to estimate the time-varying effective (or
instantaneous) reproduction number of the novel coronavirus disease (COVID-19).
The method is based on a discrete-time stochastic augmented compartmental model
that describes the virus transmission. A two-stage estimation method, which
combines the Extended Kalman Filter (EKF) to estimate the reported state
variables (active and removed cases) and a low pass filter based on a rational
transfer function to remove short term fluctuations of the reported cases, is
used with case uncertainties that are assumed to follow a Gaussian
distribution. Our method does not require information regarding serial
intervals, which makes the estimation procedure simpler without reducing the
quality of the estimate. We show that the proposed method is comparable to
common approaches, e.g., age-structured and new cases based sequential Bayesian
models. We also apply it to COVID-19 cases in the Scandinavian countries:
Denmark, Sweden, and Norway, where the positive rates were below 5\%
recommended by WHO.


  Many real systems are strongly characterized by collective cooperative
phenomena whose existence and properties still need a satisfactory explanation.
Coherently with their collective nature, they call for new and more accurate
descriptions going beyond pairwise models, such as graphs, in which all the
interactions are considered as involving only two individuals at a time.
Hypergraphs respond to this need, providing a mathematical representation of a
system allowing from pairs to larger groups. In this work, through the use of
different hypergraphs, we study how group interactions influence the evolution
of cooperation in a structured population, by analyzing the evolutionary
dynamics of the public goods game. Here we show that, likewise network
reciprocity, group interactions also promote cooperation. More importantly, by
means of an invasion analysis in which the conditions for a strategy to survive
are studied, we show how, in heterogeneously-structured populations,
reciprocity among players is expected to grow with the increasing of the order
of the interactions. This is due to the heterogeneity of connections and,
particularly, to the presence of individuals standing out as hubs in the
population. Our analysis represents a first step towards the study of
evolutionary dynamics through higher-order interactions, and gives insights
into why cooperation in heterogeneous higher-order structures is enhanced.
Lastly, it also gives clues about the co-existence of cooperative and
non-cooperative behaviors related to the structural properties of the
interaction patterns.


  This paper is devoted to the multidisciplinary modelling of a pandemic
initiated by an aggressive virus, specifically the so-called
\textit{SARS--CoV--2 Severe Acute Respiratory Syndrome, corona virus n.2}. The
study is developed within a multiscale framework accounting for the interaction
of different spatial scales, from the small scale of the virus itself and
cells, to the large scale of individuals and further up to the collective
behaviour of populations. An interdisciplinary vision is developed thanks to
the contributions of epidemiologists, immunologists and economists as well as
those of mathematical modellers. The first part of the contents is devoted to
understanding the complex features of the system and to the design of a
modelling rationale. The modelling approach is treated in the second part of
the paper by showing both how the virus propagates into infected individuals,
successfully and not successfully recovered, and also the spatial patterns,
which are subsequently studied by kinetic and lattice models. The third part
reports the contribution of research in the fields of virology, epidemiology,
immune competition, and economy focused also on social behaviours. Finally, a
critical analysis is proposed looking ahead to research perspectives.


  This contribution aims to shed light on mathematical epidemic dynamics
modelling from the viewpoint of analytical mechanics. To set the stage, it
recasts the basic SIR model of mathematical epidemic dynamics in an analytical
mechanics setting. Thereby, it considers two possible re-parameterizations of
the basic SIR model. On the one hand, it is proposed to re-scale time, while on
the other hand, to transform the coordinates, i.e.\ the independent variables.
In both cases, Hamilton's equations in terms of a suited Hamiltonian as well as
Hamilton's principle in terms of a suited Lagrangian are considered in minimal
and extended phase and state space coordinates, respectively. The corresponding
Legendre transformations relating the various options for the Hamiltonians and
Lagrangians are detailed. Ultimately, this contribution expands on a multitude
of novel vistas on mathematical epidemic dynamics modelling that emerge from
the analytical mechanics viewpoint. As result, it is believed that interesting
and relevant new research avenues open up when exploiting in depth the
analogies between analytical mechanics and mathematical epidemic dynamics
modelling.


  Non-pharmaceutical interventions (NPIs) such as quarantine, self-isolation,
social distancing, and virus-contact tracing can greatly reduce the spread of
the virus during a pandemic. In the wave of the COVID-19 pandemic, many
countries have implemented various NPIs for infection control and mitigation.
However, the stringency of the NPIs and the resulting impact among different
countries remain unclear due to the lack of quantitative factors. In this study
we took a further step to incorporate the effect of the NPIs into the pandemic
dynamics model using the concept of policy intensity factor (PIF). This idea
enables us to characterize the transition rates as time varying quantities
instead of constant values, and thus capturing the dynamical behavior of the
basic reproduction number variation in the pandemic. By leveraging a great
amount of data reported by the governments and the World Health Organization,
we projected the dynamics of the pandemic for the major economies in the world,
including the numbers of infected, susceptible, and recovered cases, as well as
the pandemic durations. It is observed that the proposed variable-rate
susceptible-exposed-infected-recovered (VR-SEIR) model fits and projects the
pandemic dynamics very well. We further showed that the resulting PIFs
correlate with the stringency of NPIs, which allows us to project the final
affected numbers of people in those countries when their current NPIs have been
imposed for 90, 180, 360 days. It provides a quantitative insight into the
effectiveness of the implemented NPIs, and sheds a new light on minimizing both
affected people from COVID-19 and the economic impact.