  We develop a minimalist compartmental model to study the impact of mobility
restrictions in Italy during the Covid-19 outbreak. We show that an early
lockdown shifts the epidemic in time, while that beyond a critical value of the
lockdown strength, the epidemic tend to restart after lifting the restrictions.
As a consequence, specific mitigation strategies must be introduced. We
characterize the relative importance of different broad strategies by
accounting for two fundamental sources of heterogeneity, i.e. geography and
demography. First, we consider Italian regions as separate administrative
entities, in which social interactions between age classs occur. Due to the
sparsity of the inter-regional mobility matrix, once started the epidemics tend
to develop independently across areas, justifying the adoption of solutions
specific to individual regions or to clusters of regions. Second, we show that
social contacts between age classes play a fundamental role and that measures
which take into account the age structure of the population can provide a
significant contribution to mitigate the rebound effects. Our model is general,
and while it does not analyze specific mitigation strategies, it highlights the
relevance of some key parameters on non-pharmaceutical mitigation mechanisms
for the epidemics.


  We argue that frequent sampling of the fraction of infected people (either by
random testing or by analysis of sewage water), is central to managing the
COVID-19 pandemic because it both measures in real time the key variable
controlled by restrictive measures, and anticipates the load on the healthcare
system due to progression of the disease. Knowledge of random testing outcomes
will (i) significantly improve the predictability of the pandemic, (ii) allow
informed and optimized decisions on how to modify restrictive measures, with
much shorter delay times than the present ones, and (iii) enable the real-time
assessment of the efficiency of new means to reduce transmission rates.
  Here we suggest, irrespective of the size of a suitably homogeneous
population, a conservative estimate of 15000 for the number of randomly tested
people per day which will suffice to obtain reliable data about the current
fraction of infections and its evolution in time, thus enabling close to
real-time assessment of the quantitative effect of restrictive measures. Still
higher testing capacity permits detection of geographical differences in
spreading rates. Furthermore and most importantly, with daily sampling in
place, a reboot could be attempted while the fraction of infected people is
still an order of magnitude higher than the level required for a relaxation of
restrictions with testing focused on symptomatic individuals. This is
demonstrated by considering a feedback and control model of mitigation where
the feed-back is derived from noisy sampling data.


  The abrupt outbreak and transmission of biological diseases has always been a
long-time concern of humankind. For long, mathematical modeling has served as a
simple and yet efficient tool to investigate, predict, and control spread of
communicable diseases through individuals. A myriad of works on epidemic models
and their variants have been reported in the literature. For better prediction
of the dynamics of a particular disease, it is important to adopt the most
suitable model. In this paper, we study some of the widely-appreciated
deterministic epidemic models in which the population is divided into
compartments based on the health status of each individual. In particular, we
provide a demographic classification of such models and study each of them in
terms of mathematical formulation, near equilibrium point stability properties,
and disease outbreak threshold conditions (basic reproduction ratio).
Furthermore, we discuss the various influential factors that need to be
considered during epidemic modeling. The main objective of this article is to
provide a basic understanding of the mathematical complexity incurred in
deterministic epidemic models with the aid of graphical illustrations obtained
through implementation.


  We present a very simple and transparent method to interpret time series of
confirmed Covid-19 cases. Roughly speaking, the analysis of weekly new
infections for each day is a tool for the definition and early detection of the
turning point of the epidemic. In Italy, the growth of Covid-19 activity was
overcome one week ago, Austria and Switzerland followed suit. This note
emphasizes the crucial delay of two weeks between infection and listing of
infection in a central database. Our estimates of the time course of Covid-19
activity and reproduction rates can reduce the information gap by several days.
We show that Spain and Germany are already beyond the turning point. In general
the estimated reproduction rates become undercritical a short time after the
main lockdown measures, and the Covid-19 activity starts to decrease
afterwards. This note gives a very optimistic outlook for all regions which
have taken strict containment measures.


  This paper provides a formal evaluation of the predictive performance of a
model (and its various updates) developed by the Institute for Health Metrics
and Evaluation (IHME) for predicting daily deaths attributed to COVID19 for
each state in the United States. The IHME models have received extensive
attention in social and mass media, and have influenced policy makers at the
highest levels of the United States government. For effective policy making the
accurate assessment of uncertainty, as well as accurate point predictions, are
necessary because the risks inherent in a decision must be taken into account,
especially in the present setting of a novel disease affecting millions of
lives. To assess the accuracy of the IHME models, we examine both forecast
accuracy as well as the predictive performance of the 95% prediction intervals
provided by the IHME models. We find that the initial IHME model underestimates
the uncertainty surrounding the number of daily deaths substantially.
Specifically, the true number of next day deaths fell outside the IHME
prediction intervals as much as 70% of the time, in comparison to the expected
value of 5%. In addition, we note that the performance of the initial model
does not improve with shorter forecast horizons. Regarding the updated models,
our analyses indicate that the later models do not show any improvement in the
accuracy of the point estimate predictions. In fact, there is some evidence
that this accuracy has actually decreased over the initial models. Moreover,
when considering the updated models, while we observe a larger percentage of
states having actual values lying inside the 95% prediction intervals (PI), our
analysis suggests that this observation may be attributed to the widening of
the PIs. The width of these intervals calls into question the usefulness of the
predictions to drive policy making and resource allocation.


  Background: The analysis of the Sars-CoV-2 epidemic is of paramount
importance to understand the dynamics of the coronavirus spread. This can help
health and government authorities take the appropriate measures and implement
suitable politics aimed at fighting and preventing it. Methods: A
time-dependent dynamic SEIR model inspired in a model previously used during
the MERS outbreak in South Korea was used to analyse the time trajectories of
active and hospitalized cases in Portugal. Results: The time evolution of the
virus spread in the country was adequately modelled. The model has changeable
parameters every five days since the onset of mitigation measures. A peak of
about 22,000 active cases is estimated, although the official value for
recovered cases is out of date. Hospitalized cases could reach a peak of about
1,250 cases, of which 200/300 in ICU units. Conclusion: With appropriate
measures, the number of active cases in Portugal can be controlled at about
22,000 people, of which about 1,250 hospitalized and 200/300 in ICU units. This
seems manageable by the country national health service with an estimated 1,140
ventilators.


  We present three data driven model-types for COVID-19 with a minimal number
of parameters to provide insights into the spread of the disease that may be
used for developing policy responses. The first is exponential growth, widely
studied in analysis of early-time data. The second is a self-exciting branching
process model which includes a delay in transmission and recovery. It allows
for meaningful fit to early time stochastic data. The third is the well-known
Susceptible-Infected-Resistant (SIR) model and its cousin, SEIR, with an
"Exposed" component. All three models are related quantitatively, and the SIR
model is used to illustrate the potential effects of short-term distancing
measures in the United States.


  Naloxone, an opioid antagonist, has been widely used to save lives from
opioid overdose, a leading cause for death in the opioid epidemic. However,
naloxone has short brain retention ability, which limits its therapeutic
efficacy. Developing better opioid antagonists is critical in combating the
opioid epidemic.Instead of exhaustively searching in a huge chemical space for
better opioid antagonists, we adopt reinforcement learning which allows
efficient gradient-based search towards molecules with desired physicochemical
and/or biological properties. Specifically, we implement a deep reinforcement
learning framework to discover potential lead compounds as better opioid
antagonists with enhanced brain retention ability. A customized multi-objective
reward function is designed to bias the generation towards molecules with both
sufficient opioid antagonistic effect and enhanced brain retention ability.
Thorough evaluation demonstrates that with this framework, we are able to
identify valid, novel and feasible molecules with multiple desired properties,
which has high potential in drug discovery.


  Some ideas are presented about the physical motivation of the apparent
capacity of generalized logistic equations to describe the outbreak of the
COVID-19 infection, and in general of quite many other epidemics. The main
focuses here are: the complex, possibly fractal, structure of the locus
describing the "contagion event set"; what can be learnt from the models of
trophic webs with "herd behaviour".


  The transcriptomics of cancer tumors are characterized with tens of thousands
of gene expression features. Patient prognosis or tumor stage can be assessed
by machine learning techniques like supervised classification tasks given a
gene expression profile. Feature selection is a useful approach to select the
key genes which helps to classify tumors. In this work we propose a feature
selection method based on Multiple Kernel Learning that results in a reduced
subset of genes and a custom kernel that improves the classification
performance when used in support vector classification. During the feature
selection process this method performs a novel latent regularisation by
relaxing the supervised target problem by introducing unsupervised structure
obtained from the latent space learned by a non linear dimensionality reduction
model. An improvement of the generalization capacity is obtained and assessed
by the tumor classification performance on new unseen test samples when the
classifier is trained with the features selected by the proposed method in
comparison with other supervised feature selection approaches.


  The outbreak of COVID-19 caused by SARS-CoV-2 has rapidly spread worldwide
and has caused over 1,400,000 infections and 80,000 deaths. There are currently
no drugs or vaccines with proven efficacy for its prevention and little
knowledge was known about the pathogenicity mechanism of SARS-CoV-2 infection.
Previous studies showed both virus and host-derived MicroRNAs (miRNAs) played
crucial roles in the pathology of virus infection. In this study, we use
computational approaches to scan the SARS-CoV-2 genome for putative miRNAs and
predict the virus miRNA targets on virus and human genome as well as the host
miRNAs targets on virus genome. Furthermore, we explore miRNAs involved
dysregulation caused by the virus infection. Our results implicated that the
immune response and cytoskeleton organization are two of the most notable
biological processes regulated by the infection-modulated miRNAs. Impressively,
we found hsa-miR-4661-3p was predicted to target the S gene of SARS-CoV-2, and
a virus-encoded miRNA MR147-3p could enhance the expression of TMPRSS2 with the
function of strengthening SARS-CoV-2 infection in the gut. The study may
provide important clues for the mechisms of pathogenesis of SARS-CoV-2.


  To predict if developing human embryos are permissive to coronaviruses, we
analyzed publicly available single cell RNA-seq datasets of zygotes, 4-cell,
8-cell, morula, inner cell mass, epiblast, primitive endoderm and trophectoderm
for the coronavirus receptors (ACE2, BSG, DPP4 and ANPEP), the Spike protein
cleavage enzymes (TMPRSS2, CTSL). We also analyzed the presence of host genes
involved in viral replication, the endosomal sorting complexes required for
transport (ESCRT) and SARS-Cov-2 interactions. The results reveal that ACE2,
BSG, DPP4 and ANPEP are expressed in the cells of the zygote, to blastocyst
including the trophectodermal lineage. ACE2, TMPRSS, BSG and CTSL are
co-transcribed in a proportion of epiblast cells and most cells of the
trophectoderm. The embryonic and trophectodermal cells also express genes for
proteins ESCRT, viral replication and those that interact with SARS-CoV-2. We
identified 1985 genes in epiblast and 1452 genes in the trophectoderm that are
enriched in the ACE2 and TMPRSS2 co-expressing cells; 216 genes of these are
common in both the cell types. These genes are associated with lipid
metabolism, lysosome, peroxisome and oxidative phosphorylation pathways.
Together our results suggest that developing human embryos could be permissive
to coronavirus entry by both canonical and non-canonical mechanisms and they
also express the genes for proteins involved in viral endocytosis and
replication. This knowledge will be useful for evidence-based patient
management for IVF during the COVID-19 pandemic.


  The outbreak of corona virus disease 2019 (COVID-19), caused by the virus
severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has already
created emergency situations in almost every country of the world. The disease
spreads all over the world within a very short period of time after its first
identification in Wuhan, China in December, 2019. In India, the outbreaks
starts on $2^{nd}$ March, 2020 and after that the cases are increasing
exponentially. Very high population density, the unavailability of specific
medicines or vaccines, insufficient evidences regarding the transmission
mechanism of the disease also make it difficult to fight against the disease
properly in India. Mathematical models have been used to predict the disease
dynamics and also to assess the efficiency of the intervention strategies in
reducing the disease burden. In this work, we propose a mathematical model to
describe the disease transmission mechanism between the individuals. We
consider the initial phase of the outbreak situation in India and our proposed
model is fitted to the daily cumulative new reported cases during the period
$2^{nd}$ March, 2020 to $24^{th}$ March, 2020. We estimate the basic
reproduction number $(R_0)$, effective reproduction number (R(t)) and epidemic
doubling time from the incidence data for the above-mentioned period. We
further assess the effect of preventive measures such as spread of awareness,
lock-down, proper hand sanitization, etc. in reducing the new cases. Two
intervention scenarios are considered depending on the variability of the
intervention strength over the period of implementation. Our study suggests
that higher intervention effort is required to control the disease outbreak
within a shorter period of time in India. Moreover, our analysis reveals that
the strength of the intervention should be strengthened over the time to
eradicate the disease effectively.


  We used Bayesian model inversion to estimate epidemic parameters from the
reported case and death rates from seven countries using data from late January
2020 to April 5th 2020. Two distinct generative model types were employed:
first a continuous time dynamical-systems implementation of a
Susceptible-Exposed-Infectious-Recovered (SEIR) model and second: a partially
observable Markov Decision Process (MDP) or hidden Markov model (HMM)
implementation of an SEIR model. Both models parameterise the size of the
initial susceptible population (S0), as well as epidemic parameters. Parameter
estimation (data fitting) was performed using a standard Bayesian scheme
(variational Laplace) designed to allow for latent unobservable states and
uncertainty in model parameters.
  Both models recapitulated the dynamics of transmissions and disease as given
by case and death rates. The peaks of the current waves were predicted to be in
the past for four countries (Italy, Spain, Germany and Switzerland) and to
emerge in 0.5-2 weeks in Ireland and 1-3 weeks in the UK. For France one model
estimated the peak within the past week and the other in the future in two
weeks. Crucially, Maximum a posteriori (MAP) estimates of S0 for each country
indicated effective population sizes of below 20% (of total population size),
under both the continuous time and HMM models. With a Bayesian weighted average
across all seven countries and both models, we estimated that 6.4% of the total
population would be immune. From the two models the maximum percentage of the
effective population was estimated at 19.6% of the total population for the UK,
16.7% for Ireland, 11.4% for Italy, 12.8% for Spain, 18.8% for France, 4.7% for
Germany and 12.9% for Switzerland.
  Our results indicate that after the current wave, a large proportion of the
total population will remain without immunity.


  In this article, we present QuASeR, a reference-free DNA sequence
reconstruction implementation via de novo assembly on both gate-based and
quantum annealing platforms. Each one of the four steps of the implementation
(TSP, QUBO, Hamiltonians and QAOA) is explained with simple proof-of-concept
examples to target both the genomics research community and quantum application
developers in a self-contained manner. The details of the implementation are
discussed for the various layers of the quantum full-stack accelerator design.
We also highlight the limitations of current classical simulation and available
quantum hardware systems. The implementation is open-source and can be found on
https://github.com/prince-ph0en1x/QuASeR.


  We investigate the problem of the predominance and survival of "weak" species
in the context of the simplest generalization of the spatial stochastic
rock-paper-scissors model to four species by considering models in which one,
two, or three species have a reduced predation probability. We show, using
lattice based spatial stochastic simulations with random initial conditions,
that if only one of the four species has its probability reduced then the most
abundant species is the prey of the "weakest" (assuming that the simulations
are large enough for coexistence to prevail). Also, among the remaining cases,
we present examples in which "weak" and "strong" species have similar average
abundances and others in which either of them dominates -- the most abundant
species being always a prey of a "weak" species with which it maintains a
unidirectional predator-prey interaction. However, in contrast to the
three-species model, we find no systematic difference in the global performance
of "weak" and "strong" species, and we conjecture that the same result will
hold if the number of species is further increased. We also determine the
probability of single species survival and coexistence as a function of the
lattice size, discussing its dependence on initial conditions and on the change
to the dynamics of the model which results from the extinction of one of the
species.


  We study the dynamics of a model of membrane vesicle transport into dendritic
spines, which are bulbous intracellular compartments in neurons driven by
molecular motors. We reduce the lubrication model proposed in [Fai et al.,
Active elastohydrodynamics of vesicles in narrow, blind constrictions. Phys.
Rev. Fluids, 2 (2017), 113601] to a fast-slow system, yielding an analytically
and numerically tractable equation equivalent to the original model in the
overdamped limit. The model's key parameters are the ratio of motors that
prefer to push toward the head of the dendritic spine to the ratio of motors
that prefer to push in the opposite direction. We perform a numerical
bifurcation analysis in these parameters and find that steady-state vesicle
velocities appear and disappear through several saddle-node bifurcations. This
process allows us to identify the region of parameter space in which multiple
stable velocities exist. We show by direct calculations that there can only be
unidirectional motion for sufficiently close vesicle-to-spine diameter ratios.
Our analysis predicts the critical vesicle-to-spine diameter ratio, at which
there is a transition from unidirectional to bidirectional motion, consistent
with experimental observations of vesicle trajectories in the literature.


  The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a
major worldwide public health emergency that has infected over $1.5$ million
people. The partially open state of S1 subunit in spike glycoprotein is
considered vital for its infection with host cell and is represented as a key
target for neutralizing antibodies. However, the mechanism elucidating the
transition from the closed state to the partially open state still remains
unclear. Here, we applied a combination of Markov state model, transition path
theory and random forest to analyze the S1 motion. Our results explored a
promising complete conformational movement of receptor-binding domain, from
buried, partially open, to detached states. We also numerically confirmed the
transition probability between those states. Based on the asymmetry in both the
dynamics behavior and backbone C$\alpha$ importance, we further suggested a
relation between chains in the trimer spike protein, which may help in the
vaccine design and antibody neutralization.


  The correct evaluation of the reproductive number $R$ for COVID-19 -- which
characterizes the average number of secondary cases generated by each typical
primary case -- is central in the quantification of the potential scope of the
pandemic and the selection of an appropriate course of action. In most models,
$R$ is modeled as a universal constant for the virus across outbreak clusters
and individuals -- effectively averaging out the inherent variability of the
transmission process due to varying individual contact rates, population
densities, demographics, or temporal factors amongst many. Yet, due to the
exponential nature of epidemic growth, the error due to this simplification can
be rapidly amplified and lead to inaccurate predictions and/or risk evaluation.
From the statistical modeling perspective, the magnitude of the impact of this
averaging remains an open question: how can this intrinsic variability be
percolated into epidemic models, and how can its impact on uncertainty
quantification and predictive scenarios be better quantified? In this paper, we
propose to study this question through a Bayesian perspective, creating a
bridge between the agent-based and compartmental approaches commonly used in
the literature. After deriving a Bayesian model that captures at scale the
heterogeneity of a population and environmental conditions, we simulate the
spread of the epidemic as well as the impact of different social distancing
strategies, and highlight the strong impact of this added variability on the
reported results. We base our discussion on both synthetic experiments --
thereby quantifying of the reliability and the magnitude of the effects -- and
real COVID-19 data.


  In a recent paper (arXiv:2003.12055), Singh and Adhikari present results of
an analysis of a mathematical model of epidemiology based on which they argue
that a $49$ day lockdown is required in India for containing the pandemic in
India. We assert that as a model study with the stated assumptions, the
analysis presented in the paper is perfectly valid, however, any serious
comparison with real data and attempts at prediction from the model are highly
problematic. The main point of the present paper is to convincingly establish
this assertion while providing a warning that the results and analysis of such
mathematical models should not be taken at face value and need to be used with
great caution by policy makers.


  The ability of random environmental variation to stabilize competitor
coexistence was pointed out long ago and, in recent years, has received
considerable attention. Analyses have focused on variations in the
log-abundances of species, with mean logarithmic growth rates when rare,
$\mathbb{E}[r]$, used as metrics for persistence. However, invasion
probabilities and the times to extinction are not single-valued functions of
$\mathbb{E}[r]$ and, in some cases, decrease as $\mathbb{E}[r]$ increases.
Here, we present a synthesis of stochasticity-induced stabilization (SIS)
phenomena based on the ratio between the expected arithmetic growth $\mu$ and
its variance $g$. When the diffusion approximation holds, explicit formulas for
invasion probabilities and persistence times are single valued, monotonic
functions of $\mu/g$. The storage effect in the lottery model, together with
other well-known examples drawn from population genetics, microbiology and
ecology (including discrete and continuous dynamics, with overlapping and
non-overlapping generations), are placed together, reviewed, and explained
within this new, transparent theoretical framework. We also clarify the
relationships between life-history strategies and SIS, and study the dynamics
of extinction when SIS fails.


  COVID-19 has spread from China across Europe and the United States and has
become a global pandemic. In countries of the Global South, due to often weaker
socioeconomic options and health care systems, effective local countermeasures
remain debated. We combine large-scale socioeconomic and traffic survey data
with detailed agent-based simulations of local transportation to analyze
COVID-19 spreading in a regional model for the Nelson Mandela Bay Municipality
in South Africa under a range of countermeasure scenarios. The simulations
indicate that any realistic containment strategy, including those similar to
the one ongoing in South Africa, may yield a manifold overload of available
intensive care units (ICUs). Only immediate and the most severe
countermeasures, up to a complete lock-down that essentially inhibits all joint
human activities, can contain the epidemic effectively.
  As South Africa exhibits rather favorable conditions compared to many other
countries of the Global South, our findings constitute rough conservative
estimates and may support identifying strategies towards containing COVID-19 as
well as any major future pandemics in these countries.


  We consider a Spatial Markov Chain model for the spread of viruses. The model
is based on the principle to represent a graph connecting nodes, which
represent humans. The vertices between the nodes represent relations between
humans. In this way, a graph is connected in which the likelihood of infectious
spread from person to person is determined by the intensity of interpersonal
contact. Infectious transfer is determined by chance. The model is extended to
incorporate various lockdown scenarios.


  This article reports a first update on the assesment of the model previously
presented in arXiv:2003.13890v1. New available data have been used to feed the
model and the comparison with real data still shows good agreement. The main
novelty of the model is that it keeps track of the date of infection of a
single individual and uses stochastic distributions to aggregate individuals
who share the same date of infection. In addition, it uses two types of
infections, mild and serious, with a different recovery time. These features
are implemented in a set of differential equations which determine the number
of Carriers, Infections, Recoveries, Hospitalized and Deaths.


  We developed a mathematical model to describe the transmission of new
coronavirus in the S\~ao Paulo State, Brazil. The model divided a community in
subpopulations comprised by young and elder persons, in order to take into
account higher risk of fatality among elder persons with severe CoViD-19. From
data collected in the S\~ao Paulo State, we estimated the transmission and
additional mortality rates, from which we calculated the basic reproduction
number R0. From estimated parameters, estimation of the deaths due to CoViD-19
was three times lower than those found in literature. Considering isolation as
a control mechanism, we varied isolation rates of young and elder persons in
order to assess their epidemiological impacts. The epidemiological scenarios
focused mainly on evaluating the number of severe CoViD-19 cases and deaths due
to this disease when isolation is introduced in a population.


  Infectious diseases are still a major global burden for modern society
causing 13 million deaths annually. One way to reduce the morbidity and
mortality rates from infectious diseases is through preventative or targeted
vaccinations. Current vaccination strategies, however, rely on the highly
specific individual contact information that is difficult and costly to obtain,
in order to identify influential spreading individuals. Current approaches also
focus only on direct contacts between individuals for spreading, and disregard
indirect transmission where a pathogen can spread between one infected
individual and one susceptible individual that visit the same location within a
short time-frame without meeting. This paper presents a novel vaccination
strategy that relies on coarse-grained contact information, both direct and
indirect, that can be easily and efficiently collected. Rather than tracking
exact contact degrees of individuals, our strategy uses the types of places
people visit to estimate a range of contact degrees for individuals,
considering both direct and indirect contacts. We conduct extensive simulations
to evaluate the performance of our strategy in comparison to the state of the
art's vaccination strategies. Results show that our strategy achieves
comparable performance to the oracle approach and outperforms all existing
strategies when considering indirect links.


  Some populations, such as red blood cells (RBCs), exhibit a pattern of
population decline that is closer to linear rather than exponential, which has
proven to be unexpectedly challenging to describe with a single simple
mathematical model. Here we show that a sub-exponential model of population
extinction can approximate very well the experimental curves of RBC extinction,
and that one possible mechanism underlying sub-exponential decay is population
heterogeneity with respect to death rates of individuals. We further show that
a sub-exponential model of population decline can be derived within the
frameworks of frequency-dependent model of population extinction if there
exists heterogeneity with respect to mortality rates such that their initial
distribution is the Gamma distribution. Notably, specific biological mechanisms
that may result in linear pattern of population decay may be different
depending on the specific biological system; however, in the end they must
converge to individual death rates being different within the population, since
uniform death rates would result in exponential population decline. As such,
the proposed model is not intended to describe the complex dynamics of RBC
biology but instead can provide a way to phenomenologically describe linear
decay of population size. We briefly discuss the potential utility of this
model for describing effects of drugs that may cause RBC depletion and conclude
with a suggestion that this tool can provide a lens for discovering linear
extinction patterns in other populations, which may have previously been
overlooked.


  The coronavirus pandemic has rapidly evolved into an unprecedented crisis.
The susceptible-infectious-removed (SIR) model and its variants have been used
for modeling the pandemic. However, time-independent parameters in the
classical models may not capture the dynamic transmission and removal
processes, governed by virus containment strategies taken at various phases of
the epidemic. Moreover, very few models account for possible inaccuracies of
the reported cases. We propose a Poisson model with time-dependent transmission
and removal rates to account for possible random errors in reporting and
estimate a time-dependent disease reproduction number, which may be used to
assess the effectiveness of virus control strategies. We apply our method to
study the pandemic in several severely impacted countries, and analyze and
forecast the evolving spread of the coronavirus. We have developed an
interactive web application to facilitate readers' use of our method.


  Coronavirus disease 2019 (COVID-19) is an ongoing pandemic infectious disease
outbreak that has significantly harmed and threatened the health and lives of
millions or even billions of people. COVID-19 has also negatively impacted the
social and economic activities of many countries significantly. With no
approved vaccine available at this moment, extensive testing of COVID-19
viruses in people are essential for disease diagnosis, virus spread
confinement, contact tracing, and determining right conditions for people to
return to normal economic activities. Identifying people who have antibodies
for COVID-19 can also help select persons who are suitable for undertaking
certain essential activities or returning to workforce. However, the
throughputs of current testing technologies for COVID-19 viruses and antibodies
are often quite limited, which are not sufficient for dealing with COVID-19
viruses' anticipated fast oscillating waves of spread affecting a significant
portion of the earth's population.
  In this paper, we propose to use compressed sensing (group testing can be
seen as a special case of compressed sensing when it is applied to COVID-19
detection) to achieve high-throughput rapid testing of COVID-19 viruses and
antibodies, which can potentially provide tens or even more folds of speedup
compared with current testing technologies. The proposed compressed sensing
system for high-throughput testing can utilize expander graph based compressed
sensing matrices developed by us \cite{Weiyuexpander2007}.


  Based on the classical SIR model, we derive a simple modification for the
dynamics of epidemics with a known incubation period of infection. The model is
described by a system of integro-differential equations. Parameters of our
model directly related to epidemiological data. We derive some analytical
results, as well as perform numerical simulations. We use the proposed model to
analyze COVID-19 epidemic data in Armenia. We propose a strategy: organize a
quarantine, and then conduct extensive testing of risk groups during the
quarantine, evaluating the percentage of the population among risk groups and
people with symptoms.


  Starting with the power law for the total number of detected infections, we
propose differential equations describing the effect of momentum epidemic
management. Our 2-phase formula matches very well the curves of the total
numbers of the Covid-19 infections in many countries; the first phase is
described by Bessel functions. It provides projections for the saturation,
assuming that the management is steady. We discuss Austria, Brazil, Germany,
Japan, India, Israel, Italy, the Netherlands, Sweden, Switzerland, UK, and the
USA, including some analysis of the second waves.


  Testing of symptomatic persons for infection with SARS-CoV-2 is occurring
worldwide. We propose two types of case-control studies that can be carried out
jointly in test-settings for symptomatic persons. The first, the test-negative
case-control design (TND) is the easiest to implement; it only demands
collecting information about potential risk factors for COVID-19 from the
tested symptomatic persons. The second, standard case-control studies with
population controls, requires the collection of data on one or more population
controls for each person who is tested in the test facilities, so that
test-positives and test-negatives can each be compared with population
controls. The TND will detect differences in risk factors between symptomatic
persons who have COVID-19 (test-positives) and those who have other respiratory
infections (test-negatives). However, risk factors with effect sizes of equal
magnitude for both COVID-19 and other respiratory infections will not be
identified by the TND. Therefore, we discuss how to add population controls to
compare with the test-positives and the test-negatives, yielding two additional
case-control studies. We describe two options for population control groups:
one composed of accompanying persons to the test facilities, the other drawn
from existing country-wide health care databases. We also describe other
possibilities for population controls. Combining the TND with population
controls yields a triangulation approach that distinguishes between exposures
that are risk factors for both COVID-19 and other respiratory infections, and
exposures that are risk factors for just COVID-19. This combined design can be
applied to future epidemics, but also to study causes of non-epidemic disease.


  Can public social media data be harnessed to predict COVID-19 case counts? We
analyzed approximately 15 million COVID-19 related posts on Weibo, a popular
Twitter-like social media platform in China, from November 1, 2019 to March 31,
2020. We developed a machine learning classifier to identify "sick posts,"
which are reports of one's own and other people's symptoms and diagnosis
related to COVID-19. We then modeled the predictive power of sick posts and
other COVID-19 posts on daily case counts. We found that reports of symptoms
and diagnosis of COVID-19 significantly predicted daily case counts, up to 14
days ahead of official statistics. But other COVID-19 posts did not have
similar predictive power. For a subset of geotagged posts (3.10% of all
retrieved posts), we found that the predictive pattern held true for both Hubei
province and the rest of mainland China, regardless of unequal distribution of
healthcare resources and outbreak timeline. Researchers and disease control
agencies should pay close attention to the social media infosphere regarding
COVID-19. On top of monitoring overall search and posting activities, it is
crucial to sift through the contents and efficiently identify true signals from
noise.


  Mathematical models of epidemiological systems enable investigation of and
predictions about potential disease outbreaks. However, commonly used models
are often highly simplified representations of incredibly complex systems.
Because of these simplifications, the model output, of say new cases of a
disease over time, or when an epidemic will occur, may be inconsistent with
available data. In this case, we must improve the model, especially if we plan
to make decisions based on it that could affect human health and safety, but
direct improvements are often beyond our reach. In this work, we explore this
problem through a case study of the Zika outbreak in Brazil in 2016. We propose
an embedded discrepancy operator---a modification to the model equations that
requires modest information about the system and is calibrated by all relevant
data. We show that the new enriched model demonstrates greatly increased
consistency with real data. Moreover, the method is general enough to easily
apply to many other mathematical models in epidemiology.


  The novel coronavirus SARS-CoV-2 and resulting COVID-19 disease have had an
unprecedented spread and continue to cause an increasing number of fatalities
worldwide. While vaccines are still under development, social distancing,
extensive testing, and quarantining of confirmed infected subjects remain the
most effective measures to contain the pandemic. These measures carry a
significant socioeconomic cost. In this work, we introduce a novel
optimization-based decision-making framework for managing the COVID-19 outbreak
in the US. This includes modeling the dynamics of affected populations,
estimating the model parameters and hidden states from data, and an optimal
control strategy for sequencing social distancing and testing events such that
the number of infections is minimized. The analysis of our extensive
computational efforts reveals that social distancing and quarantining are most
effective when implemented early, with quarantining of confirmed infected
subjects having a much higher impact. Further, we find that "on-off" policies
alternating between strict social distancing and relaxing such restrictions can
be effective at "flattening" the curve while likely minimizing social and
economic cost.


  We measure the effect of different public health regulations to the spread of
COVID-19, based on a SEIRA model -- a SEIR model including asymptomatic
transmissions. The cumulative confirmed cases and death show nonlinear positive
relationship with the value of asymptomatic rate. Based on this model, we
analyze the inhibit effects to COVID-19 of three types of public health
policies, i.e. isolation of laboratory confirmed cases, general personal
protection and quarantine (lock-down). The simulations conclude that the
isolation display limited effects to the asymptomatic viral carriers. The
general personal protection and quarantine perform similar effects when the
their percentages of participants are same. When the total proportion of
asymptomatic, mild symptomatic and neglected patients is 40%, only depends on
isolation policy may lead to an additional 75% infections, compared with
general personal protection or quarantine with an efficiency 80%. At end, we
provide seven recommendations of public health intervention before and during
an aerial transmitted epidemic (COVID-19).


  Cographs are exactly the hereditarily well-colored graphs, i.e., the graphs
for which a greedy vertex coloring of every induced subgraph uses only the
minimally necessary number of colors $\chi(G)$. We show that greedy colorings
are a special case of the more general hierarchical vertex colorings, which
recently were introduced in phylogenetic combinatorics. Replacing cotrees by
modular decomposition trees generalizes the concept of hierarchical colorings
to arbitrary graphs. We show that every graph has a modularly-minimal coloring
$\sigma$ satisfying $|\sigma(M)|=\chi(M)$ for every strong module $M$ of $G$.
This, in particular, shows that modularly-minimal colorings provide a useful
device to design efficient coloring algorithms for certain hereditary graph
classes. For cographs, the hierarchical colorings coincide with the
modularly-minimal coloring. As a by-product, we obtain a simple linear-time
algorithm to compute a modularly-minimal coloring of $P_4$-sparse graphs.


  Presently, India bears amongst the highest burden of non-communicable
diseases such as diabetes mellitus (DM), hypertension (HT), and cardio vascular
disease (CVD) and thus represents a vulnerable target to the
SARS-CoV-2/COVID-19 pandemic. Involvement of the angiotensin converting enzyme
2 (ACE2) in susceptibility to infection and pathogenesis by SARS-CoV-2 is
currently an actively pursued research area. An increased susceptibility to
infection in individuals with DM, HT and CVD together with higher levels of
circulating ACE2 in these settings presents a scenario where interaction with
soluble ACE2 may result in disseminated virus-receptor complexes that could
enhance virus acquisition and pathogenesis. Thus, understanding the SARS-CoV-2
receptor binding domain-ACE2 interaction, both membrane bound and in the cell
free context may contribute to elucidating the role of co-morbidities in
increased susceptibility to infection and pathogenesis. Both Azithromycin and
Hydroxychloroquine (HCQ) have shown efficacy in mitigating viral carriage in
infected individuals. Furthermore, each of these compounds generate active
metabolites which in turn may also modulate virus-receptor interaction and thus
influence clinical outcomes. In this study, we model the structural interaction
of S1 with both full-length and soluble ACE2. Additionally, therapeutic drugs
and their active metabolites were docked with soluble ACE2 protein. Our results
show that S1 from either of the reported Indian sequences can bind both
full-length and soluble ACE2, albeit with varying affinity that can be
attributed to a reported substitution in the RBD. Furthermore, both
Azythromycin and HCQ together with their active metabolites can allosterically
affect, to a range of extents, binding of S1 to ACE2.


  The emergence of COVID-19 has severely compromised the arsenal of antiviral
and antibiotic drugs. Drug discovery is a multistep process with a high failure
rate, high cost and it takes approximately 10-12 years for the development of
new molecules into the clinical candidate. On the other side, drug repurposing
also called old drugs for new uses, is an attractive alternative approach for a
new application of marketed FDA approved or investigational drugs. In the
current pandemic situation raised due to COVID-19, repurposing of existing FDA
approved drugs are emerging as the first line of the treatment. The causative
viral agent of this highly contagious disease and acute respiratory syndrome
coronavirus (SARS-CoV) shares high nucleotide similarity. Therefore, many
existing viral targets are structurally expected to be similar to SARS-CoV and
likely to be inhibited by the same compounds. Here, we selected three viral key
proteins based on their vital role in viral life cycle: ACE2 (helps in entry
into the human host), viral nonstructural proteins RNA-dependent RNA polymerase
(RdRp) NSP12, and NSP16 which helps in replication, and viral latency (invasion
from immunity). The FDA approved drugs chloroquine (CQ), hydroxychloroquine
(HCQ), remdesivir (RDV) and arbidol (ABD) are emerging as promising agents to
combat COVID-19. Our hypothesis behind the docking studies is to determine the
binding affinities of these drugs and identify the key amino acid residues
playing a key role in their mechanism of action. The docking studies were
carried out through Autodock and online COVID-19 docking server. Further
studies on a broad range of FDA approved drugs including few more protein
targets, molecular dynamics studies, in-vitro and in-vivo biological evaluation
are required to identify the combination therapy targeting various stages of
the viral life cycle.


  Recent statistics of malaria shows that over 200 million cases and estimated
deaths of nearly half a million occur globally. Africa alone accounts for
almost 90% of the cases. Several studies have been conducted to understand the
disease transmission dynamics. In particular, mathematical methods have been
frequently used to model and understand the disease dynamics and outbreak
patterns. Although, mathematical methods have provided good results for
homogeneous populations, these methods impose significant limitations for
studying malaria dynamics in heterogeneous populations, a result of various
factors, e.g. spatial and temporal fluctuations, social networks, human
movements pattern etc. This paper proposes an agent-based modelling approach
that permits modelling and analysing malaria dynamics for heterogenous
populations. Our approach is illustrated using the climate and demographic data
for the Tripura, Limpopo and Benin cities. Our agent-based simulation has been
validated against the reported cases of malaria collected in the cities
mentioned. Furthermore, the efficiency of the proposed model has been compared
with the mathematical model used as benchmark. A statistical test confirms the
proposed model is robust and has potential for predicting the peak seasons of
malaria. This potentially makes our methods a useful tool as an intervention
mechanism, which will have impact on hospitals, healthcare providers, health
organisations.


  Studies in malaria control cover many areas such as medicine, sociology,
biology, mathematic, physic, computer science and so forth. Researches in the
realm of mathematic are conducted to predict the occurrence of the disease and
to support the eradication process. Basically, the modeling methodology is
predominantly deterministic and differential equation based while selecting
clinical and biological features that seem to be important. Yet, if the
individual characteristics matter when modeling the disease, the overall
estimation of the malaria is not done based on the health status of each
individual but in a non-specified percentage of the global population. The goal
of this paper is to propose a model that relies on a daily evolution of the
individual's state, which depends on their mobility and the characteristics of
the area they visit. Thus, the mobility data of a single person moving from one
area to another, gathered thanks to mobile networks, is the essential building
block to predict the outcome of the disease. We implement our solution and
demonstrate its effectiveness through empirical experiments. The results show
how promising the model is in providing possible insights into the failure of
the disease control in the Kedougou region.


  For the description of a pandemic mathematical models could be interesting.
Both for physicians and politicians as a base for decisions to treat the
disease. The responsible estimation of parameters is a main issue of
mathematical pandemic models. Especially a good choice of $\beta$ as the number
of others that one infected person encounters per unit time (per day)
influences the adequateness of the results of the model. For the actual
COVID-19 pandemic some aspects of the parameter choice will be discussed.
Because of the incompatibility of the data of the Johns-Hopkins-University to
the data of the German Robert-Koch-Institut we use the COVID-19 data of the
European Centre for Disease Prevention and Control (ECDC) as a base for the
parameter estimation. Two different mathematical methods for the data analysis
will be discussed in this paper and possible sources of trouble will be shown.
As example of the parameter choice serve the data of the USA and the UK. The
resulting parameters will be used estimated and used in W.\,O. Kermack and
A.\,G. McKendrick's SIR model. Strategies for the commencing and ending of
social and economic shutdown measures are discussed. The numerical solution of
the ordinary differential equation system of the modified SIR model is being
done with a Runge-Kutta integration method of fourth order. At the end the
applicability of the SIR model could be shown essentially. Suggestions about
appropriate points in time at which to commence with lockdown measures based on
the acceleration rate of infections conclude the paper.


  Within the context that Greece promises a success story in the fight against
the disease, this paper proposes a novel method to study the evolution of the
Greek COVID-19 infection-curve in relation to the anti-COVID-19 policies
applied to control the pandemic. Based on the ongoing spreading of COVID-19 and
the insufficient data for applying classic time-series approaches, the analysis
builds on the visibility graph algorithm to study the Greek COVID-19
infection-curve as a complex network. By using the modularity optimization
algorithm, the generated visibility graph is divided into communities defining
periods of different connectivity in the time-series body. These periods reveal
a sequence of different typologies in the evolution of the disease, starting
with a power pattern, where a second order polynomial (U-shaped) pattern
intermediates, being followed by a couple of exponential patterns, and ending
up with a current logarithmic pattern revealing that the evolution of the Greek
COVID-19 infection-curve tends into saturation. The network analysis also
illustrates stability of hubs and instability of medium and low-degree nodes,
implying a low probability to meet maximum (infection) values at the future and
high uncertainty in the variability of other values below the average. The
overall approach contributes to the scientific research by proposing a novel
method for the structural decomposition of a time-series into periods, which
allows removing from the series the disconnected past-data facilitating better
forecasting, and provides insights of good policy and decision-making practices
and management that may help other countries improve their performance in the
war against COVID-19.


  We study the problem of optimal control of the stochastic SIR model. Models
of this type are used in mathematical epidemiology to capture the time
evolution of highly infectious diseases such as COVID-19. Our approach relies
on reformulating the Hamilton-Jacobi-Bellman equation as a stochastic minimum
principle. This results in a system of forward backward stochastic differential
equations, which is amenable to numerical solution via Monte Carlo simulations.
We present a number of numerical solutions of the system under a variety of
scenarios.


  Relapse, the recurrence of a disorder following a symptomatic remission, is a
frequent outcome in substance abuse disorders. Some of our prior results
suggested that relapse, in the context of abusive drinking, is likely an
"unbeatable" force as long as recovered individuals continue to interact in the
environments that lead to and/or reinforce the persistence of abusive drinking
behaviors. Our earlier results were obtained via a deterministic model that
ignored differences between individuals, that is, in a rather simple "social"
setting. In this paper, we address the role of relapse on drinking dynamics but
use models that incorporate the role of "chance", or a high degree of "social"
heterogeneity, or both. Our focus is primarily on situations where relapse
rates are high. We first use a Markov chain model to simulate the effect of
relapse on drinking dynamics. These simulations reinforce the conclusions
obtained before, with the usual caveats that arise when the outcomes of
deterministic and stochastic models are compared. However, the simulation
results generated from stochastic realizations of an "equivalent" drinking
process in populations "living" in small world networks, parameterized via a
disorder parameter $p$, show that there is no social structure within this
family capable of reducing the impact of high relapse rates on drinking
prevalence, even if we drastically limit the interactions between individuals
($p\approx 0$).


  The observed sequence variation at a locus informs about the evolutionary
history of the sample and past population size dynamics. The Kingman coalescent
is used in a generative model of molecular sequence variation to infer
evolutionary parameters. However, it is well understood that inference under
this model does not scale well with sample size. Here, we build on recent work
based on a lower resolution coalescent process, the Tajima coalescent, to model
longitudinal samples. While the Kingman coalescent models the ancestry of
labeled individuals, the heterochronous Tajima coalescent models the ancestry
of individuals labeled by their sampling time. We propose a new inference
scheme for the reconstruction of effective population size trajectories based
on this model with the potential to improve computational efficiency. Modeling
of longitudinal samples is necessary for applications (e.g. ancient DNA and RNA
from rapidly evolving pathogens like viruses) and statistically desirable
(variance reduction and parameter identifiability). We propose an efficient
algorithm to calculate the likelihood and employ a Bayesian nonparametric
procedure to infer the population size trajectory. We provide a new MCMC
sampler to explore the space of heterochronous Tajima's genealogies and model
parameters. We compare our procedure with state-of-the-art methodologies in
simulations and applications.


  We consider a single outbreak susceptible-infected-recovered (SIR) model and
corresponding estimation procedures for the effective reproductive number
$\mathcal{R}(t)$. We discuss the estimation of the underlying SIR parameters
with a generalized least squares (GLS) estimation technique. We do this in the
context of appropriate statistical models for the measurement process. We use
asymptotic statistical theories to derive the mean and variance of the limiting
(Gaussian) sampling distribution and to perform post statistical analysis of
the inverse problems. We illustrate the ideas and pitfalls (e.g., large
condition numbers on the corresponding Fisher information matrix) with both
synthetic and influenza incidence data sets.


  In this paper, we propose a new real-time differential virus transmission
model, which can give more accurate and robust short-term predictions of
COVID-19 transmitted infectious disease with benefits of near-term trend
projection. Different from the existing Susceptible-Exposed-Infected-Removed
(SEIR) based virus transmission models, which fits well for pandemic modeling
with sufficient historical data, the new model, which is also SEIR based, uses
short history data to find the trend of the changing disease dynamics for the
infected, the dead and the recovered so that it can naturally accommodate the
adaptive real-time changes of disease mitigation, business activity and social
behavior of populations. As the parameters of the improved SEIR models are
trained by short history window data for accurate trend prediction, our
differential epidemic model, essentially are window-based time-varying SEIR
model. Since SEIR model still is a physics-based disease transmission model,
its near-term (like one month) projection can still be very instrumental for
policy makers to guide their decision for disease mitigation and business
activity policy change in a real-time. This is especially useful if the
pandemic lasts more than one year with different phases across the world like
1918 flu pandemic. Numerical results on the recent COVID-19 data from China,
Italy and US, California and New York states have been analyzed.


  This paper repurposes the classic insight from network theory that
long-distance connections drive disease propagation into a strategy for
controlling a second wave of Covid-19. We simulate a scenario in which a
lockdown is first imposed on a population and then partly lifted while
long-range transmission is kept at a minimum. Simulated spreading patterns
resemble contemporary distributions of Covid-19 across nations, regions, and
provinces, providing some model validation. Results suggest that the proposed
strategy may significantly flatten a second wave. We also find that
post-lockdown flare-ups remain local longer, aiding geographical containment.
Public policy may target long ties by heavily focusing medical testing and
mobility tracking efforts on traffic and transport. This policy can be
communicated to the general public as a simple and reasonable principle: Stay
nearby or get checked.


  The current COVID-19 pandemic is affecting different countries in different
ways. The assortment of reporting techniques alongside other issues, such as
underreporting and budgetary constraints, makes predicting the spread and
lethality of the virus a challenging task. This work attempts to gain a better
understanding of how COVID-19 will affect one of the least studied countries,
namely Brazil. Currently, several Brazilian states are in a state of lock-down.
However, there is political pressure for this type of measures to be lifted.
This work considers the impact that such a termination would have on how the
virus evolves locally. This was done by extending the SEIR model with an on /
off strategy. Given the simplicity of SEIR we also attempted to gain more
insight by developing a neural regressor. We chose to employ features that
current clinical studies have pinpointed has having a connection to the
lethality of COVID-19. We discuss how this data can be processed in order to
obtain a robust assessment.


  Machine learning encompasses a set of tools and algorithms which are now
becoming popular in almost all scientific and technological fields. This is
true for molecular dynamics as well, where machine learning offers promises of
extracting valuable information from the enormous amounts of data generated by
simulation of complex systems. We provide here a review of our current
understanding of goals, benefits, and limitations of machine learning
techniques for computational studies on atomistic systems, focusing on the
construction of empirical force fields from ab-initio databases and the
determination of reaction coordinates for free energy computation and enhanced
sampling.


  Social distancing and isolation have been introduced widely to counter the
COVID-19 pandemic. However, more moderate contact reduction policies become
desirable owing to adverse social, psychological, and economic consequences of
a complete or near-complete lockdown. Adopting a social network approach, we
evaluate the effectiveness of three targeted distancing strategies designed to
'keep the curve flat' and aid compliance in a post-lockdown world. These are
limiting interaction to a few repeated contacts, seeking similarity across
contacts, and strengthening communities via triadic strategies. We simulate
stochastic infection curves that incorporate core elements from infection
models, ideal-type social network models, and statistical relational event
models. We demonstrate that strategic reduction of contact can strongly
increase the efficiency of social distancing measures, introducing the
possibility of allowing some social contact while keeping risks low. This
approach provides nuanced insights to policy makers for effective social
distancing that can mitigate negative consequences of social isolation.


  The moral decisions we make during this period, such as deciding whether to
comply with quarantine rules, have unprecedented societal effects. We simulate
the "escape from Milan" that occurred on March 7th-8th 2020, when many
travelers moved from a high-risk zone (Milan) to southern regions of Italy
(Campania and Lazio) immediately after an imminent lockdown was announced. Our
simulations show that fewer than 50 active cases might have caused the sudden
spread of the virus observed afterwards in these regions. The surprising
influence of the actions of few individuals on societal dynamics challenges our
cognitive expectations -- as in normal conditions, collective dynamics are
rather robust to the decisions of few "cheaters". This situation therefore
requires novel educational strategies that increase our awareness and
understanding of the unprecedented effects of our individual moral decisions.


  The pandemic prevalence of COVID-19 has become a very serious global health
issue. Scientists all over the world have been heavily invested in the
discovery of a drug to combat SARS-CoV-2. It has been found that RNA-dependent
RNA Polymerase (RdRp) plays a crucial role in SARS-CoV-2 replication, and thus
could be a potential drug target. Here, comprehensive computational approaches
including drug repurposing and molecular docking were employed to predict an
effective drug candidate targeting RdRp of SARS-CoV-2. This study revealed that
Rifabutin, Rifapentine, Fidaxomicin,
7-methyl-guanosine-5'-triphosphate-5'-guanosine and Ivermectin have a potential
inhibitory interaction with RdRp of SARS-CoV-2, and could be effective drugs
for COVID-19. In addition, virtual screening of the compounds from ZINC
database also allowed the prediction of two compounds (ZINC09128258 and ZINC
09883305) with pharmacophore features that interact effectively with RdRp of
SARS-CoV-2; indicating their potentiality as effective inhibitors of the
enzyme. Furthermore, ADME analysis along with analysis of toxicity was also
investigated to check the pharmacokinetics and drug-likeness properties of the
two compounds. Comparative structural analysis of protein-inhibitor complexes
revealed that positions of the amino acid Y32, K47, Y122, Y129, H133, N138,
D140, T141, S709 and N781 are crucial for drug surface hotspot in the RdRp of
SARS-CoV-2.


  We are in the midst of a global viral pandemic, one with no cure and a high
mortality rate. The Human Leukocyte Antigen (HLA) gene complex plays a critical
role in host immunity. We predicted HLA class I and II alleles from the
transcriptome sequencing data prepared from the bronchoalveolar lavage fluid
samples of five patients at the early stage of the COVID-19 outbreak. We
identified the HLA-I allele A*24:02 in four out of five patients, which is
higher than the expected frequency (17.2%) in the South Han Chinese population.
The difference is statistically significant with a p-value less than $10^{-4}$.
Our analysis results may help provide future insights on disease
susceptibility.


  Much scientific enquiry across disciplines is founded upon a mechanistic
treatment of dynamic systems that ties form to function. A highly visible
instance of this is in molecular biology, where an important goal is to
determine functionally-relevant forms/structures that a protein molecule
employs to interact with molecular partners in the living cell. This goal is
typically pursued under the umbrella of stochastic optimization with algorithms
that optimize a scoring function. Research repeatedly shows that current
scoring function, though steadily improving, correlate weakly with molecular
activity. Inspired by recent momentum in generative deep learning, this paper
proposes and evaluates an alternative approach to generating
functionally-relevant three-dimensional structures of a protein. Though
typically deep generative models struggle with highly-structured data, the work
presented here circumvents this challenge via graph-generative models. A
comprehensive evaluation of several deep architectures shows the promise of
generative models in directly revealing the latent space for sampling novel
tertiary structures, as well as in highlighting axes/factors that carry
structural meaning and open the black box often associated with deep models.
The work presented here is a first step towards interpretative, deep generative
models becoming viable and informative complementary approaches to protein
structure prediction.


  It is of vital importance to understand and track the dynamics of rapidly
unfolding epidemics. The health and economic consequences of the current
COVID-19 pandemic provide a poignant case. Here we point out that since they
are based on differential equations, the most widely used models of epidemic
spread are plagued by an approximation that is not justified in the case of the
current COVID-19 pandemic. Taking the example of data from New York City, we
show that currently used models significantly underestimate the initial basic
reproduction number ($R_0$). The correct description, based on integral
equations, can be implemented in most of the reported models and it much more
accurately accounts for the dynamics of the epidemic after sharp changes in
$R_0$ due to restrictive public congregation measures. It also provides a novel
way to determine the incubation period, and most importantly, as we demonstrate
for several countries, this method allows an accurate monitoring of $R_0$ and
thus a fine-tuning of any restrictive measures. Integral equation based models
do not only provide the conceptually correct description, they also have more
predictive power than differential equation based models, therefore we do not
see any reason for using the latter.


  The fatality rate of Covid-19 escalates with age and is larger in men than
women. I show that these variations correlate strongly with the level of the
viral receptor protein ACE2 in rat lungs, which is consistent with the still
limited and apparently contradictory data on human ACE2. Surprisingly, lower
levels of the receptor correlate with higher fatality. However, a previous
mathematical model predicts that the speed of viral progression in the organism
has a maximum and then declines with the receptor level. Moreover, many
manifestations of severe CoViD-19, such as severe lung injury, exacerbated
inflammatory response and thrombotic problems may derive from increased
Angiotensin II (Ang-II) level that results from degradation of ACE2 by the
virus. I present here a mathematical model based on the influence of ACE2 on
viral propagation and disease severity. The model fits Covid-19 fatality rate
across age and sex with high accuracy ($r^2>0.9$) under the hypothesis that
SARS-CoV-2 infections are in the dynamical regimes in which increased receptor
slows down viral propagation. Moreover, rescaling the model parameters by the
ratio of the binding rates of the spike proteins of SARS-CoV and SARS-CoV-2
allows predicting the fatality rate of SARS-CoV across age and sex, thus
linking the molecular and epidemiological levels. The presented model opposes
the fear that angiotensin receptor blockers (ARB), suggested as a therapy
against the most adverse effects of CoViD-19, may favour viral propagation, and
suggests that Ang-II and ACE2 are candidate prognostic factors for detecting
population that needs stronger protection.


  The current pandemic has highlighted the need for methodologies that can
quickly and reliably prioritize clinically approved compounds for their
potential effectiveness for SARS-CoV-2 infections. In the past decade, network
medicine has developed and validated multiple predictive algorithms for drug
repurposing, exploiting the sub-cellular network-based relationship between a
drug's targets and disease genes. Here, we deployed algorithms relying on
artificial intelligence, network diffusion, and network proximity, tasking each
of them to rank 6,340 drugs for their expected efficacy against SARS-CoV-2. To
test the predictions, we used as ground truth 918 drugs that had been
experimentally screened in VeroE6 cells, and the list of drugs under clinical
trial, that capture the medical community's assessment of drugs with potential
COVID-19 efficacy. We find that while most algorithms offer predictive power
for these ground truth data, no single method offers consistently reliable
outcomes across all datasets and metrics. This prompted us to develop a
multimodal approach that fuses the predictions of all algorithms, showing that
a consensus among the different predictive methods consistently exceeds the
performance of the best individual pipelines. We find that 76 of the 77 drugs
that successfully reduced viral infection do not bind the proteins targeted by
SARS-CoV-2, indicating that these drugs rely on network-based actions that
cannot be identified using docking-based strategies. These advances offer a
methodological pathway to identify repurposable drugs for future pathogens and
neglected diseases underserved by the costs and extended timeline of de novo
drug development.


  A population can be immune to epidemics even if not all of its individual
members are immune to the disease, so long as sufficiently many are immune -
this is the traditional notion of herd immunity. In the smartphone era a
population can be immune to epidemics even if not a single one of its members
is immune to the disease - a notion we call "digital herd immunity", which is
similarly an emergent characteristic of the population. This immunity arises
because contact-tracing protocols based on smartphone capabilities can lead to
highly efficient quarantining of infected population members and thus the
extinguishing of nascent epidemics. When the disease characteristics are
favorable and smartphone usage is high enough, the population is in this immune
phase. As usage decreases there is a novel "contact-tracing phase transition"
to an epidemic phase. We present and study a simple branching-process model for
COVID-19 and show that digital immunity is possible regardless of the
proportion of non-symptomatic transmission.


  This report provides a visual examination of Covid-19 case and death data. In
particular, it shows that country specific differences can too a large extend
be explained by two easily interpreted parameters. Namely, the delay between
reported cases and deaths and the fraction of cases observed. Furthermore, this
allows to lower bound the actual total number of people already infected.


  A novel Coronavirus pandemic emerged in December of 2019, causing devastating
public health impact across the world. In the absence of a safe and effective
vaccine or antiviral, strategies for mitigating the burden of the pandemic are
focused on non-pharmaceutical interventions, such as social-distancing,
contact-tracing, quarantine, isolation and the use of face-masks in public. We
develop a new mathematical model for assessing the population-level impact of
these mitigation strategies. Simulations of the model, using data relevant to
COVID-19 transmission in New York state and the entire US, show that the
pandemic will peak in mid and late April, respectively. The worst-case scenario
projections for cumulative mortality (based on the baseline levels of
anti-COVID non-pharmaceutical interventions considered in the study) in New
York State and the entire US decrease dramatically by 80% and 64%,
respectively, if the strict social-distancing measures implemented are
maintained until the end of May or June, 2020. This study shows that early
termination of strict social-distancing could trigger a devastating second wave
with burden similar to that projected before the onset of strict
social-distance. The use of efficacious face-masks (efficacy greater than 70%)
could lead to the elimination of the pandemic if at least 70% of the residents
of New York state use such masks consistently (nationwide, a compliance of at
least 80% will be required using such masks). The use of low efficacy masks,
such as cloth masks (of efficacy less than 30%), could also lead to significant
reduction of COVID-19 burden (albeit, they are not able to lead to
elimination). Combining low efficacy masks with improved levels of other
anti-COVID-19 intervention measures can lead to elimination of the pandemic.
The mask coverage needed to eliminate COVID-19 decreases if mask-use is
combined with strict social-distancing.


  During decades, mathematical models have been used to predict the behavior of
physical and biologic systems, and to define strategies aiming the minimization
of the effects regarding different types of diseases. In the present days, the
development of mathematical models to simulate the dynamic behavior of novel
coronavirus disease (COVID-19) is considered an important theme due to the
quantity of infected people worldwide. In this work, the aim is to determine an
optimal control strategy for vaccine administration in COVID-19 pandemic
treatment considering real data from China. For this purpose, an inverse
problem is formulated and solved in order to determine the parameters of the
compartmental SIR (Susceptible-Infectious-Recovered) model. To solve such
inverse problem, the Differential Evolution (DE) algorithm is employed. After
this step, two optimal control problems (mono- and multi-objective) to
determine the optimal strategy for vaccine administration in COVID-19 pandemic
treatment are proposed. The first consists of minimizing the quantity of
infected individuals during the treatment. The second considers minimizing
together the quantity of infected individuals and the prescribed vaccine
concentration during the treatment, i.e., a multi-objective optimal control
problem. The solution of each optimal control problems is obtained using DE and
Multi-Objective Differential Evolution (MODE) algorithms, respectively. The
results regarding the proposed multi-objective optimal control problem provides
a set of evidences from which an optimal strategy for vaccine administration
can be chosen, according to a given criterion.


  We propose a mathematical model for the transmission dynamics of some strains
of the bacterium Vibrio cholerae, responsible for the cholera disease in
humans. We prove that, when the basic reproduction number is equal to one, a
transcritical bifurcation occurs for which the endemic equilibrium emanates
from the disease-free point. A control function is introduced into the model,
representing the distribution of chlorine water tablets for water purification.
An optimal control problem is then proposed and analyzed, where the goal is to
determine the fraction of susceptible individuals who should have access to
chlorine water tablets in order to minimize the total number of new infections
plus the total cost associated with the distribution of chlorine water tablets,
over the considered period of time. Finally, we consider real data of the
cholera outbreak in Yemen, from 27 April 2017 to 15 April 2018, choosing the
values of the parameters of the uncontrolled model that fit the real data.
Using our optimal control results, we show, numerically, that the distribution
of chlorine water tablets could have stopped, in a fast way, the worst cholera
outbreak that ever occurred in human history. Due to the critical situation of
Yemen, we also simulate the case where only a small percentage of susceptible
individuals has access to chlorine water tablets and obtain an optimal control
solution that decreases, substantially, the maximum number of infective
individuals attained at the outbreak.


  Eukaryotic cells transmit extracellular signal information to cellular
interiors through the formation of a ternary complex made up of a ligand (or
agonist), G-protein, and G-protein coupled receptor (GPCR). Previously
formalized theories of ternary complex formation have mainly assumed that
observable states of receptors can only take the form of monomers. Here, we
propose a multiary complex model of GPCR signaling activations via the vector
representation of various unobserved aggregated receptor states. Our results
from model simulations imply that receptor aggregation processes can govern
cooperative effects in a regime inaccessible by previous theories. In
particular, we show how the affinity of ligand-receptor binding can be largely
varied by various oligomer formations in the low concentration range of
G-protein stimulus.


  This work deals with a mesoscopic and deterministic theory for the collective
dynamics of Water Coherence Domains, and it represents a continuation of a
previous work published in Plos ONE
https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0188753 This
approach qualitatively reproduces some features of the experimental
phenomenology, and it could sustain an evolutionary theory based on layered
synchronization processes. Its relevance for the description of the processes
behind the emergence of life is discussed.


  In recent work we reported the vibrational spectrum of more than 100,000
known protein structures, and a self-consistent sonification method to render
the spectrum in the audible range of frequencies (Extreme Mechanics Letters,
2019). Here we present a method to transform these molecular vibrations into
materialized vibrations of thin water films using acoustic actuators, leading
to complex patterns of surface waves, and using the resulting macroscopic
images in further processing using deep convolutional neural networks.
Specifically, the patterns of water surface waves for each protein structure is
used to build training sets for neural networks, aimed to classify and further
process the patterns. Once trained, the neural network model is capable of
discerning different proteins solely by analyzing the macroscopic surface wave
patterns in the water film. Not only can the method distinguish different types
of proteins (e.g. alpha-helix vs hybrids of alpha-helices and beta-sheets), but
it is also capable of determining different folding states of the same protein,
or the binding events of proteins to ligands. Using the DeepDream algorithm,
instances of key features of the deep neural network can be made visible in a
range of images, allowing us to explore the inner workings of protein surface
wave patter neural networks, as well as the creation of new images by finding
and highlighting features of protein molecular spectra in a range of
photographic input. The integration of the water-focused realization of
cymatics, combined with neural networks and especially generative methods,
offer a new direction to realize materiomusical "Inceptionism" as a possible
direction in nano-inspired art. The method could have applications for
detecting different protein structures, the effect of mutations, or uses in
medical imaging and diagnostics, with broad impact in nano-to-macro
transitions.


  Witnessing its first case in late January 2020 India has seen a sharp rise in
the number of positive cases of COVID-19. 34 States/UT (s) of the country have
been found to be affected by the pandemic to date. We in this work, study the
progress of COVID-19 pandemic in India. We aim to create transmission network
visualization (s) of COVID-19 in India and perform analysis upon them. Using
the transmission networks obtained we attempt to find the possible Super
Spreader Individual (s) and Super Spreader Events (SSE) responsible for the
outbreak in their respective regions. We discuss the potentials of network
analysis in mitigating the further spread of the disease. This is one of the
initial studies of the outbreak in India and the first attempt to study the
pandemic in the country from a transmission network perspective.


  Multiple lines of evidence strongly suggest that infection hotspots, where a
single individual infects many others, play a key role in the transmission
dynamics of COVID-19. However, most of the existing epidemiological models fail
to capture this aspect by neither representing the sites visited by individuals
explicitly nor characterizing disease transmission as a function of individual
mobility patterns. In this work, we introduce a temporal point process modeling
framework that specifically represents visits to the sites where individuals
get in contact and infect each other. Under our model, the number of infections
caused by an infectious individual naturally emerges to be overdispersed. Using
an efficient sampling algorithm, we demonstrate how to estimate the
transmission rate of infectious individuals at the sites they visit and in
their households using Bayesian optimization and longitudinal case data.
Simulations using fine-grained and publicly available demographic data and site
locations from Bern, Switzerland showcase the flexibility of our framework. To
facilitate research and analyses of other cities and regions, we release an
open-source implementation of our framework.


  Spectral analysis characterises oscillatory time series behaviours such as
cycles, but accurate estimation requires reasonable numbers of observations.
Current COVID-19 time series for many countries are short: pre- and
post-lockdown series are shorter still. Accurate estimation of potentially
interesting cycles within such series seems beyond reach. We solve the problem
of obtaining accurate estimates from short time series by using recent Bayesian
spectral fusion methods. Here we show that transformed new daily COVID-19 cases
for many countries generally contain three cycles operating at wavelengths of
around 2.7, 4.1 and 6.7 days (weekly). We show that the shorter cycles are
suppressed after lockdown. The pre- and post lockdown differences suggest that
the weekly effect is at least partly due to non-epidemic factors, whereas the
two shorter cycles seem intrinsic to the epidemic. Unconstrained, new cases
grow exponentially, but the internal cyclic structure causes periodic falls in
cases. This suggests that lockdown success might only be indicated by four or
more daily falls in cases. Spectral learning for epidemic time series
contributes to the understanding of the epidemic process, helping evaluate
interventions and assists with forecasting. Spectral fusion is a general
technique that is able to fuse spectra recorded at different sampling rates,
which can be applied to a wide range of time series from many disciplines.


  The COVID-19 crisis called for rapid reaction from all the fields of
biomedical research. Traditional drug development involves time consuming
pipelines that conflict with the urgence of identifying effective therapies
during a health and economic emergency. Drug repositioning, that is the
discovery of new clinical applications for drugs already approved for different
therapeutic contexts, could provide an effective shortcut to bring COVID-19
treatments to the bedside in a timely manner. Moreover, computational
approaches can help accelerate the process even further. Here we present the
application of computational drug repositioning tools based on transcriptomics
data to identify drugs that are potentially able to counteract SARS-CoV-2
infection, and also to provide insights on their mode of action. We believe
that mucolytics and HDAC inhibitors warrant further investigation. In addition,
we found that the DNA Mismatch repair pathway is strongly modulated by drugs
with experimental in vitro activity against SARS-CoV-2 infection. Both full
results and methods are publicly available.


  Calibration of a SIR (Susceptibles-Infected-Recovered) model with official
international data for the COVID-19 pandemics provides a good example of the
difficulties inherent the solution of inverse problems. Inverse modeling is set
up in a framework of discrete inverse problems, which explicitly considers the
role and the relevance of data. Together with a physical vision of the model,
the present work addresses numerically the issue of parameters calibration in
SIR models, it discusses the uncertainties in the data provided by
international authorities, how they influence the reliability of calibrated
model parameters and, ultimately, of model predictions.


  In the absence of any available vaccines or drugs, prevention of the spread
of Coronavirus Disease 2019 (COVID-19) pandemic is being achieved by putting
many mitigation measures in place. It is indispensable to have robust and
reliable ways of evaluating the effectiveness of these measures. In this work,
we assume that, at a very coarse-grained level of description, the overall
effect of all the mitigation measures is that we can still describe the spread
of the pandemic using the most basic Susceptible-Exposed-Infectious-Removed
($SEIR$) model but with an "effective contact rate" ($\beta$) which is
time-dependent. We then use the time series data of the number of infected
individuals in the population to extract the instantaneous effective contact
rate which is the result of various social interventions put in place. This
approach has the potential to be significantly useful while evaluating the
impact of mitigation measures on the spread of COVID-19 in near future.


  With a two-layer contact-dispersion model and data in China, we analyze the
cost-effectiveness of three types of antiepidemic measures for COVID-19:
regular epidemiological control, local social interaction control, and
inter-city travel restriction. We find that: 1) intercity travel restriction
has minimal or even negative effect compared to the other two at the national
level; 2) the time of reaching turning point is independent of the current
number of cases, and only related to the enforcement stringency of
epidemiological control and social interaction control measures; 3) strong
enforcement at the early stage is the only opportunity to maximize both
antiepidemic effectiveness and cost-effectiveness; 4) mediocre stringency of
social interaction measures is the worst choice. Subsequently, we cluster
countries/regions into four groups based on their control measures and provide
situation assessment and policy suggestions for each group.


  Combination therapy has shown to improve therapeutic efficacy while reducing
side effects. Importantly, it has become an indispensable strategy to overcome
resistance in antibiotics, anti-microbials, and anti-cancer drugs. Facing
enormous chemical space and unclear design principles for small-molecule
combinations, the computational drug-combination design has not seen generative
models to meet its potential to accelerate resistance-overcoming drug
combination discovery. We have developed the first deep generative model for
drug combination design, by jointly embedding graph-structured domain knowledge
and iteratively training a reinforcement learning-based chemical graph-set
designer. First, we have developed Hierarchical Variational Graph Auto-Encoders
(HVGAE) trained end-to-end to jointly embed gene-gene, gene-disease, and
disease-disease networks. Novel attentional pooling is introduced here for
learning disease-representations from associated genes' representations.
Second, targeting diseases in learned representations, we have recast the
drug-combination design problem as graph-set generation and developed a deep
learning-based model with novel rewards. Specifically, besides chemical
validity rewards, we have introduced a novel generative adversarial award,
being generalized sliced Wasserstein, for chemically diverse molecules with
distributions similar to known drugs. We have also designed a network
principle-based reward for drug combinations. Numerical results indicate that,
compared to graph embedding methods, HVGAE learns more informative and
generalizable disease representations. Case studies on four diseases show that
network-principled drug combinations tend to have low toxicity. The generated
drug combinations collectively cover the disease module similar to FDA-approved
drug combinations and could potentially suggest novel systems-pharmacology
strategies.


  Roccia San Sebastiano is a tectonic-karstic cave located at the foot of the
southern slope of Mt. Massico, in the territory of Mondragone (Caserta) in
Campania (southern Italy). Systematic excavation has been carried out since
2001, leading to the partial exploration of an important Pleistocene deposit,
extraordinarily rich in lithic and faunal remains. The aim of this paper is to
(1) present the stratigraphic sequence of Roccia San Sebastiano, and (2)
technologically describe the lithic materials of squares F14 t18, t19, t20; E16
t16, t17, t18 recently recognised as Uluzzian. The stratigraphic sequence is
more than 3 metres thick and dates from the Middle to the Upper Palaeolithic.
It contains different techno-complexes: Gravettian, Aurignacian, Uluzzian and
Mousterian. In the Uluzzian lithic assemblage mostly local pebbles of chert
were used in order to produce small-sized objects. The concept of debitage
mainly deals with unidirectional debitage with absent or fairly accurate
management of the convexities and angles; the striking platforms are usually
natural or made by one stroke. It is attested the use of both direct freehand
percussion and bipolar technique on anvil in the same reduction sequence.
Amongst the retouched tools the presence of two lunates is of note. This study
of the Roccia San Sebastiano Uluzzian lithic complexes is significant for
understanding the dynamics of the transition from Middle to Upper Palaeolithic
in the Tyrrhenian margin of southern Italy.


  We discuss methods for {\em a priori} selection of parameters to be estimated
in inverse problem formulations (such as Maximum Likelihood, Ordinary and
Generalized Least Squares) for dynamical systems with numerous state variables
and an even larger number of parameters. We illustrate the ideas with an
in-host model for HIV dynamics which has been successfully validated with
clinical data and used for prediction.


  We use daily data from Lombardy, the Italian region most affected by the
COVID-19 outbreak, to calibrate a SIR model individually on each municipality.
These are all covered by the same health system and, in the post-lockdown phase
we focus on, all subject to the same social distancing regulations. We find
that municipalities with a higher number of cases at the beginning of the
period analyzed have a lower rate of diffusion, which cannot be imputed to herd
immunity. In particular, there is a robust and strongly significant negative
correlation between the estimated basic reproduction number ($R_0$) and the
initial outbreak size, in contrast with the role of $R_0$ as a \emph{predictor}
of outbreak size. We explore different possible explanations for this
phenomenon and conclude that a higher number of cases causes changes of
behavior, such as a more strict adoption of social distancing measures among
the population, that reduce the spread. This result calls for a transparent,
real-time distribution of detailed epidemiological data, as such data affects
the behavior of populations in areas affected by the outbreak.


  An accurate closed-form solution is obtained to the SIR Epidemic Model
through the use of Asymptotic Approximants (Barlow et. al, 2017, Q. Jl Mech.
Appl. Math, 70 (1), 21-48). The solution is created by analytically continuing
the divergent power series solution such that it matches the long-time
asymptotic behavior of the epidemic model. The utility of the analytical form
is demonstrated through its application to the COVID-19 pandemic.


  COVID-19, a novel coronavirus, is currently a major worldwide threat. It has
infected more than a million people globally leading to hundred-thousands of
deaths. In such grave circumstances, it is very important to predict the future
infected cases to support prevention of the disease and aid in the healthcare
service preparation. Following that notion, we have developed a model and then
employed it for forecasting future COVID-19 cases in India. The study indicates
an ascending trend for the cases in the coming days. A time series analysis
also presents an exponential increase in the number of cases. It is supposed
that the present prediction models will assist the government and medical
personnel to be prepared for the upcoming conditions and have more readiness in
healthcare systems.


  We propose four novel mathematical models, describing the microscopic
mechanisms of force generation in the cardiac muscle tissue, which are suitable
for multiscale numerical simulations of cardiac electromechanics. Such models
are based on a biophysically accurate representation of the regulatory and
contractile proteins in the sarcomeres. Our models, unlike most of the
sarcomere dynamics models that are available in the literature and that feature
a comparable richness of detail, do not require the time-consuming Monte Carlo
method for their numerical approximation. Conversely, the models that we
propose only require the solution of a system of PDEs and/or ODEs (the most
reduced of the four only involving 20 ODEs), thus entailing a significant
computational efficiency. By focusing on the two models that feature the best
trade-off between detail of description and identifiability of parameters, we
propose a pipeline to calibrate such parameters starting from experimental
measurements available in literature. Thanks to this pipeline, we calibrate
these models for room-temperature rat and for body-temperature human cells. We
show, by means of numerical simulations, that the proposed models correctly
predict the main features of force generation, including the steady-state
force-calcium and force-length relationships, the length-dependent prolongation
of twitches and increase of peak force, the force-velocity relationship.
Moreover, they correctly reproduce the Frank-Starling effect, when employed in
multiscale 3D numerical simulation of cardiac electromechanics.


  Aggregation of ubiquitinated cargo by oligomers of the protein p62 is an
important preparatory step in cellular autophagy. In this work a mathematical
model for the dynamics of these heterogeneous aggregates in the form of a
system of ordinary differential equations is derived and analyzed. Three
different parameter regimes are identified, where either aggregates are
unstable, or their size saturates at a finite value, or their size grows
indefinitely as long as free particles are abundant. The boundaries of these
regimes as well as the finite size in the second case can be computed
explicitly. The growth in the third case (quadratic in time) can also be made
explicit by formal asymptotic methods. The qualitative results are illustrated
by numerical simulations. A comparison with recent experimental results permits
a partial parametrization of the model.


  We present a modified diffusive epidemic process that has a finite threshold
on scale-free graphs. The diffusive epidemic process describes the epidemic
spreading in a non-sedentary population, and it is a reaction-diffusion
process. In the diffusion stage, the individuals can jump between connected
nodes, according to different diffusive rates for the infected and susceptible
individuals. In the reaction stage, the contagion can happen if there is an
infected individual sharing the same node, and infected individuals can
spontaneously recover. Our main modification is to turn the number of
individuals' interactions independent on the population size by using Gillespie
algorithm with a reaction time $t_\mathrm{max}$, exponentially distributed with
mean inversely proportional to the node concentration. Our simulation results
of the modified model on Barabasi-Albert networks are compatible with a
continuous phase transition with a finite threshold from an absorbing phase to
an active phase when increasing the concentration. The transition obeys the
mean-field critical exponents of the order parameter, its fluctuations and the
spatial correlation length, whose values are $\beta=1$, $\gamma'=0$ and
$\nu_\perp=1/2$, respectively. In addition, the system presents logarithmic
corrections with pseudo-exponents $\widehat{\beta}=\widehat{\gamma}'=-3/2$ on
the order parameter and its fluctuations, respectively. The most evident
implication of our simulation results is if the individuals avoid social
interactions in order to not spread a disease, this leads the system to have a
finite threshold in scale-free graphs, allowing for epidemic control.