  The appearance of a new dangerous and contagious disease requires the
development of a drug therapy faster than what is foreseen by usual mechanisms.
Many drug therapy developments consist in investigating through different
clinical trials the effects of different specific drug combinations by
delivering it into a test group of ill patients, meanwhile a placebo treatment
is delivered to the remaining ill patients, known as the control group. We
compare the above technique to a new technique in which all patients receive a
different and reasonable combination of drugs and use this outcome to feed a
Neural Network. By averaging out fluctuations and recognizing different patient
features, the Neural Network learns the pattern that connects the patients
initial state to the outcome of the treatments and therefore can predict the
best drug therapy better than the above method. In contrast to many available
works, we do not study any detail of drugs composition nor interaction, but
instead pose and solve the problem from a phenomenological point of view, which
allows us to compare both methods. Although the conclusion is reached through
mathematical modeling and is stable upon any reasonable model, this is a
proof-of-concept that should be studied within other expertises before
confronting a real scenario. All calculations, tools and scripts have been made
open source for the community to test, modify or expand it. Finally it should
be mentioned that, although the results presented here are in the context of a
new disease in medical sciences, these are useful for any field that requires a
experimental technique with a control group.


  Macroscopic growth laws, solutions of mean field equations, describe in an
effective way an underlying complex dynamics. They are applied to study the
spreading of infections, as in the case of CoviD-19, where the counting of the
cumulated number $N(t)$ of detected infected individuals is a generally
accepted, coarse-grain, variable to understand the epidemic phase. However
$N(t)$ does not take into account the unknown number of asymptomatic, not
detected, cases $A(t)$. Therefore, the question arises if the observed time
series of data of $N(t)$ is a reliable tool for monitoring the evolution of the
infectious disease. We study a system of coupled differential equations which
includes the dynamics of the spreading among symptomatic and asymptomatic
individuals and the strong containment effects due to the social isolation. The
solution is therefore compared with a macroscopic law for the population $N(t)$
coming from a single, non-linear, differential equation with no explicit
reference to $A(t)$, showing the equivalence of the two methods. Indeed, $N(t)$
takes into account a more complex and detailed population dynamics which
permits the evaluation of the number of asymptomatic individuals also. The
model is then applied to Covid-19 spreading in Italy where a transition from an
exponential behavior to a Gompertz growth for $N(t)$ has been observed in more
recent data. Then the information contained in the data analysis of $N(t)$ is
reliable to understand the epidemic phase, although it does not describe the
total infected population. The asymptomatic population is larger than the
symptomatic one in the fast growth phase of the spreading.


  In a previous note we made an analysis of the spreading of the COVID disease
in Italy. We used a model based on the logistic and Hubbert functions, the
analysis we exploited has shown limited usefulness in terms of predictions and
failed in fixing fundamental indications like the point of inflection of the
disease growth. In this note we elaborate on the previous model, using
multi-logistic models and attempt a more realistic analysis.


  The COVID-19 pandemic demands assimilation of all available biomedical
knowledge to decode its mechanisms of pathogenicity and transmission. Despite
the recent renaissance in unsupervised neural networks for decoding
unstructured natural languages, a platform for the real-time synthesis of the
exponentially growing biomedical literature and its comprehensive triangulation
with deep omic insights is not available. Here, we present the nferX platform
for dynamic inference from over 45 quadrillion possible conceptual associations
extracted from unstructured biomedical text, and their triangulation with
Single Cell RNA-sequencing based insights from over 25 tissues. Using this
platform, we identify intersections between the pathologic manifestations of
COVID-19 and the comprehensive expression profile of the SARS-CoV-2 receptor
ACE2. We find that tongue keratinocytes and olfactory epithelial cells are
likely under-appreciated targets of SARS-CoV-2 infection, correlating with
reported loss of sense of taste and smell as early indicators of COVID-19
infection, including in otherwise asymptomatic patients. Airway club cells,
ciliated cells and type II pneumocytes in the lung, and enterocytes of the gut
also express ACE2. This study demonstrates how a holistic data science platform
can leverage unprecedented quantities of structured and unstructured publicly
available data to accelerate the generation of impactful biological insights
and hypotheses.


  A phenomenological model to describe the Corona Virus(covid-19) Pandemic
spread in a given population is developed. It enables the identification of the
key quantities required to form adequate policies for control and mitigation in
terms of observable parameters using the Landau-Stuart equation. It is intended
to be complementary to detailed simulations and methods published recently by
Ferguson and collaborators, March 16, (2020). The results suggest that the
initial growth/spreading rate gamma-c of the disease, and the fraction of
infected persons in the population p-i can be used to define a
`retardation/inhibition coefficient' k-star , which is a measure of the
effectiveness of the control policies adopted.
  The results are obtained analytically and numerically using a simple Python
code. The solutions provide both qualitative and quantitative information. They
substantiate and justify two basic control policies enunciated by WHO and
adopted in many countries: a) Systematic and early intensive testing
individuals for covid-19 and b) Sequestration policies such as `social/physical
distancing' and population density reduction by strict quarantining are
essential for making k-star greater than 1, necessary for suppressing the
pandemic. The model indicates that relaxing such measures when the infection
rate starts to decrease as a result of earlier policies could simply restart
the infection rate in the non-infected population. Presently available
available statistical data in WHO and other reports can be readily used to
determine the the key parameters of the model. Possible extensions to the basic
model to make it more realistic are indicated.


  We analyse the strategy equilibrium of dilemma games considering a payoff
matrix affected by small and random perturbations on the off-diagonal.
  Notably, a recent work [1] reported that, while cooperation is sustained by
perturbations acting on the main diagonal, a less clear scenario emerges when
perturbations act on the off-diagonal. Thus, the second case represents the
core of this investigation, aimed at completing the description of the effects
that payoff perturbations have on the dynamics of evolutionary games. Our
results, achieved by analysing the proposed model under a variety of
configurations, as different update rules, suggest that off-diagonal
perturbations actually constitute a non-trivial form of noise. In particular,
the most interesting effects are detected near the phase transition, as
perturbations tend to move the strategy distribution towards non-ordered states
of equilibrium, supporting cooperation when defection is pervading the
population, and supporting defection in the opposite case. To conclude, we
identified a form of noise that, under controlled conditions, could be used to
enhance cooperation, and greatly delay its extinction.


  We explore a Markov model used in the analysis of gene expression, involving
the bursty production of pre-mRNA, its conversion to mature mRNA, and its
consequent degradation. We demonstrate that the integration used to compute the
solution of the stochastic system can be approximated by the evaluation of
special functions. Furthermore, the form of the special function solution
generalizes to a broader class of burst distributions. In light of the broader
goal of biophysical parameter inference from transcriptomics data, we apply the
method to simulated data, demonstrating effective control of precision and
runtime. Finally, we suggest a non-Bayesian approach to reducing the
computational complexity of parameter inference to linear order in state space
size and number of candidate parameters.


  Omics technologies are powerful tools for analyzing patterns in gene
expression data for thousands of genes. Due to a number of systematic
variations in experiments, the raw gene expression data is often obfuscated by
undesirable technical noises. Various normalization techniques were designed in
an attempt to remove these non-biological errors prior to any statistical
analysis. One of the reasons for normalizing data is the need for recovering
the covariance matrix used in gene network analysis. In this paper, we
introduce a novel normalization technique, called the covariance shift
(C-SHIFT) method. This normalization algorithm uses optimization techniques
together with the blessing of dimensionality philosophy and energy minimization
hypothesis for covariance matrix recovery under additive noise (in biology,
known as the bias). Thus, it is perfectly suited for the analysis of
logarithmic gene expression data. Numerical experiments on synthetic data
demonstrate the method's advantage over the classical normalization techniques.
Namely, the comparison is made with Rank, Quantile, cyclic LOESS (locally
estimated scatterplot smoothing), and MAD (median absolute deviation)
normalization methods. We also evaluate the performance of C-SHIFT algorithm on
real biological data.


  The generalized Lotka-Volterra (gLV) equations model the microbiome as a
collection of interacting ecological species. Here we use a particular
experimentally-derived gLV model of C. difficile infection (CDI) as a case
study to generate methods that are applicable to generic gLV models. We examine
how to transition gLV systems between multiple steady states through the
application of direct control protocols, which alter the state of the system
via the instantaneous addition or subtraction of microbial species. Then, the
geometry of the basins of attraction of point attractors is compressed into an
attractor network, which decomposes a multistable high-dimensional landscape
into web of bistable subsystems. This attractor network is used to identify
efficient (total intervention volume minimizing) protocols that drive the
system from one basin to another. In some cases, the most efficient control
protocol is circuitous and will take the system through intermediate steady
states with sequential interventions. Clinically, the efficient control of the
microbiome has pertinent applications for bacteriotherapies, which seek to
remedy microbiome-affiliated diseases by directly altering the composition of
the gut microbiome.


  We propose a new space of phylogenetic trees which we call wald space. The
motivation is to develop a space suitable for statistical analysis of
phylogenies, but with a geometry based on more biologically principled
assumptions than existing spaces: in wald space, trees are close if they induce
similar distributions on genetic sequence data. As a point set, wald space
contains the previously developed Billera-Holmes-Vogtmann (BHV) tree space; it
also contains disconnected forests, like the edge-product (EP) space but
without certain singularities of the EP space. We investigate two related
geometries on wald space. The first is the geometry of the Fisher information
metric of character distributions induced by the two-state symmetric Markov
substitution process on each tree. Infinitesimally, the metric is proportional
to the Kullback-Leibler divergence, or equivalently, as we show, any to f
-divergence. The second geometry is obtained analogously but using a related
continuous-valued Gaussian process on each tree, and it can be viewed as the
trace metric of the affine-invariant metric for covariance matrices. We derive
a gradient descent algorithm to project from the ambient space of covariance
matrices to wald space. For both geometries we derive computational methods to
compute geodesics in polynomial time and show numerically that the two
information geometries (discrete and continuous) are very similar. In
particular geodesics are approximated extrinsically. Comparison with the BHV
geometry shows that our canonical and biologically motivated space is
substantially different.


  Data assimilation is used to optimally fit a classical epidemiology model to
the Johns Hopkins data of the Covid-19 pandemic. The optimisation is based on
the confirmed cases and confirmed deaths. This is the only data available with
reasonable accuracy. Infection and recovery rates can be infered from the model
as well as the model parameters. The parameters can be linked with government
actions or events like the end of the holiday season. Based on this numbers
predictions for the future can be made and control targets specified. With
other words: We look for a solution to a given model which fits the given data
in an optimal sense. Having that solution, we have all parameters.


  Most causal discovery algorithms find causal structure among a set of
observed variables. Learning the causal structure among latent variables
remains an important open problem, particularly when using high-dimensional
data. In this paper, we address a problem for which it is known that inputs
cause outputs, and these causal relationships are encoded by a causal network
among a set of an unknown number of latent variables. We developed a deep
learning model, which we call a redundant input neural network (RINN), with a
modified architecture and a regularized objective function to find causal
relationships between input, hidden, and output variables. More specifically,
our model allows input variables to directly interact with all latent variables
in a neural network to influence what information the latent variables should
encode in order to generate the output variables accurately. In this setting,
the direct connections between input and latent variables makes the latent
variables partially interpretable; furthermore, the connectivity among the
latent variables in the neural network serves to model their potential causal
relationships to each other and to the output variables. A series of simulation
experiments provide support that the RINN method can successfully recover
latent causal structure between input and output variables.


  Recent advances in modelling VDJ recombination and subsequent selection of T
and B cell receptors provide useful tools to analyze and compare immune
repertoires across time, individuals, and tissues. A suite of tools--IGoR [1],
OLGA [2] and SONIA [3]--have been publicly released to the community that allow
for the inference of generative and selection models from high-throughput
sequencing data. However using these tools requires some scripting or
command-line skills and familiarity with complex datasets. As a result the
application of the above models has not been available to a broad audience. In
this application note we fill this gap by presenting Simple OLGA & SONIA (SOS),
a web-based interface where users with no coding skills can compute the
generation and post-selection probabilities of their sequences, as well as
generate batches of synthetic sequences. The application also functions on
mobile phones.


  We study the epidemic spreading on spatial networks where the probability
that two nodes are connected decays with their distance as a power law. As the
exponent of the distance dependence grows, model networks smoothly transition
from the random network limit to the regular lattice limit. We show that
despite keeping the average number of contacts constant, the increasing
exponent hampers the epidemic spreading by making long-distance connections
less frequent. The spreading dynamics is influenced by the distance-dependence
exponent as well and changes from exponential growth to power-law growth. The
observed power-law growth is compatible with recent analyses of empirical data
on the spreading of COVID-19 in numerous countries.


  In this work we demonstrate how to automate parts of the infectious
disease-control policy-making process via performing inference in existing
epidemiological models. The kind of inference tasks undertaken include
computing the posterior distribution over controllable, via direct
policy-making choices, simulation model parameters that give rise to acceptable
disease progression outcomes. Among other things, we illustrate the use of a
probabilistic programming language that automates inference in existing
simulators. Neither the full capabilities of this tool for automating inference
nor its utility for planning is widely disseminated at the current time. Timely
gains in understanding about how such simulation-based models and inference
automation tools applied in support of policymaking could lead to less
economically damaging policy prescriptions, particularly during the current
COVID-19 pandemic.


  A new coronavirus disease, called COVID-19, appeared in the Chinese region of
Wuhan at the end of last year; since then the virus spread to other countries,
including most of Europe. We propose a differential equation governing the
evolution of the COVID-19. This dynamic equation also describes the evolution
of the number of infected people for 13 common respiratory viruses (including
the SARS-CoV-2). We validate our theoretical predictions with experimental data
for Italy, Belgium and Luxembourg, and compare them with the predictions of the
logistic model. We find that our predictions are in good agreement with the
real world since the beginning of the appearance of the COVID-19; this is not
the case for the logistic model that only applies to the first days. The second
part of the work is devoted to modelling the descending phase, i.e. the
decrease of the number of people tested positive for COVID-19. Also in this
case, we propose a new set of dynamic differential equations that we solved
numerically. We use our differential equations parametrised with experimental
data to make several predictions, such as the date when Italy, Belgium, and
Luxembourg will reach a peak number of SARS-CoV-2 infected people. The
descending curves provide valuable information such as the duration of the
COVID-19 epidemic in a given Country and therefore when it will be possible to
return to normal life. The study of the the dynamics of COVID-19 when the
population have been subject to less restrictive measures is beyond the scope
of this work and it will be matter of future works.


  Objective: To conduct a meta-analysis of current studies that examined sex
differences in severity and mortality in patients with COVID-19, and identify
potential mechanisms underpinning these differences. Methods: We performed a
systematic review to collate data from observational studies examining
associations of sex differences with clinical outcomes of COVID-19. PubMed, Web
of Science and four preprint servers were searched for relevant studies. Data
were extracted and analyzed using meta-analysis where possible, with summary
data presented otherwise. Publicly available bulk RNA sequencing (RNA-seq),
single-cell RNA sequencing (scRNA-seq), and chromatin immunoprecipitation
sequencing (ChIP-seq) data were analyzed to explore the potential mechanisms
underlying the observed association. Results: 39 studies met inclusion
criteria, representing 77932 patients, of which 41510 (53.3%) were males. Men
were at a markedly increased risk of developing severe cases compared with
women. Furthermore, the pooled odds ratio (OR) of mortality for male group
compared with the female group indicated significant higher mortality rate for
male. Data from scRNA-seq suggest that men have a higher amount of
ACE2-expressing pulmonary alveolar type II cells than women. Sex-based
immunological differences exist. The expression of androgen receptor (AR) is
positively correlated with ACE2, and there is evidence that AR may directly
regulate the expression of ACE2. Conclusions: This meta-analysis detected an
increased severity and mortality rate in the male populations with COVID-19,
which might be attributable to the sex-based differences in cellular
compositions and immunological microenvironments of the lung. The host cell
receptor ACE2 is likely regulated by AR signaling pathway, which is identified
as a potential target for prevention and treatment of SARS-Cov-2 infections in
men.


  We study a simple realistic model for describing the diffusion of an
infectious disease on a population of individuals. The dynamics is governed by
a single functional delay differential equation, which, in the case of a large
population, can be solved exactly, even in the presence of a time-dependent
infection rate. This delay model has a higher degree of accuracy than that of
the so-called SIR model, commonly used in epidemiology, which, instead, is
formulated in terms of ordinary differential equations. We apply this model to
describe the outbreak of the new infectious disease, Covid-19, in Italy, taking
into account the containment measures implemented by the government in order to
mitigate the spreading of the virus and the social costs for the population.


  Specific elements of viral genomes regulate interactions within host cells.
Here, we calculated the secondary structure content of >2500 coronaviruses and
computed >100000 human protein interactions with severe acute respiratory
syndrome coronavirus 2 (SARS-CoV-2). We found that the 3 and 5 prime ends are
the most structured elements in the viral genome and the 5 prime end has the
strongest propensity to associate with human proteins. The domain encompassing
nucleotides 23000-24000 is highly conserved both at the sequence and structural
level, while the region upstream varies significantly. These two sequences code
for a domain of the viral protein Spike S that interacts with the human
receptor angiotensin-converting enzyme 2 (ACE2) and has the potential to bind
sialic acids. Our predictions indicate that the first 1000 nucleotides in the 5
prime end can interact with proteins involved in viral RNA processing such as
double-stranded RNA specific editases and ATP-dependent RNA-helicases, in
addition to other high-confidence candidate partners. These interactions,
previously reported to be also implicated in HIV, reveal important information
on host-virus interactions. The list of transcriptional and
post-transcriptional elements recruited by SARS-CoV-2 genome provides clues on
the biological pathways associated with gene expression changes in human cells.


  Human genes required for SARS-CoV-2 entry into human cells, ACE2 and FURIN,
were employed as baits to build genomics-guided maps of up-stream regulatory
elements, their expression and functions in human body, including
pathophysiologically-relevant cell types. Genes acting as repressors and
activators of the ACE2 and FURIN genes were identified based on the analyses of
gene silencing and overexpression experiments as well as relevant transgenic
mouse models. Panels of repressors (VDR; GATA5; SFTPC; HIF1a) and activators
(HMGA2; INSIG1) were then employed to identify existing drugs that could be
repurposed to mitigate the coronavirus infection. Present analyses identify
Vitamin D and Quercetin as promising pandemic mitigation agents. Gene
expression profiles of Vitamin D and Quercetin activities and their established
safety records as over-the-counter medicinal substances suggest that they may
represent viable candidates for further assessment and considerations of their
potential as coronavirus pandemic mitigation agents. Notably, gene set
enrichment analyses and expression profiling experiments identify multiple
drugs, most notably testosterone, dexamethasone, and doxorubicin, smoking, and
many disease conditions that appear to act as putative coronavirus
infection-promoting agents. Discordant patterns of Testosterone versus
Estradiol impacts on SCARS-CoV-2 targets suggest a plausible molecular
explanation of the apparently higher male mortality during coronavirus
pandemic.


  [Abridged] Naturally produced droplets from humans (such as those produced by
breathing, talking, sneezing, and coughing) include several types of cells
(e.g., epithelial cells and cells of the immune system), physiological
electrolytes contained in mucous and saliva (e.g. Na+, K+, Cl-), as well as,
potentially, several infectious agents (e.g. bacteria, fungi, and viruses). In
response to the novel coronavirus SARS-CoV-2 epidemic, which has become a major
public health issue worldwide, we provide a concise overview of airborne germ
transmission as seen from a physics perspective. We study whether coronavirus
aerosols can travel far from the immediate neighborhood and get airborne with
the convective currents developed within confined spaces. We also provide a
recommendation that could help to slow down the spread of the virus.


  We present a simple analytical model to describe the fast increase of deaths
produced by the corona virus (COVID-19) infections. The 'D' (deaths) model
comes from a simplified version of the SIR (susceptible-infected-recovered)
model known as SI model. It assumes that there is no recovery. In that case the
dynamical equations can be solved analytically and the result is extended to
describe the D-function that depends on three parameters that we can fit to the
data. Results for the data from Spain, Italy and China are presented. The model
is validated by comparing with the data of deaths in China, which are well
described. This allows to make predictions for the development of the disease
in Spain and Italy.


  Previous self-assembly experiments on a model icosahedral plant virus have
shown that, under physiological conditions, capsid proteins initially bind to
the genome through an en masse mechanism and form nucleoprotein complexes in a
disordered state, which raises the questions as to how virions are assembled
into a highly ordered structure in the host cell. Using small-angle X-ray
scattering, we find out that a disorder-order transition occurs under
physiological conditions upon an increase in capsid protein concentrations. Our
cryo-transmission electron microscopy reveals closed spherical shells
containing in vitro transcribed viral RNA even at pH 7.5, in marked contrast
with the previous observations. We use Monte Carlo simulations to explain this
disorder-order transition and find that, as the shell grows, the structures of
disordered intermediates in which the distribution of pentamers does not belong
to the icosahedral subgroups become energetically so unfavorable that the caps
can easily dissociate and reassemble overcoming the energy barriers for the
formation of perfect icosahedral shells. In addition, we monitor the growth of
capsids under the condition that the nucleation and growth is the dominant
pathway and show that the key for the disorder-order transition in both en
masse and nucleation and growth pathways lies in the strength of elastic energy
compared to the other forces in the system including protein-protein
interactions and the chemical potential of free subunits. Our findings explain,
at least in part, why perfect virions with icosahedral order form under
different conditions including physiological ones.


  COVID-19 continues to spread across the country and around the world. Current
strategies for managing the spread of COVID-19 include social distancing. We
present VERA, an interactive AI tool, that first enables users to specify
conceptual models of the impact of social distancing on the spread of COVID-19.
Then, VERA automatically spawns agent-based simulations from the conceptual
models, and, given a data set, automatically fills in the values of the
simulation parameters from the data. Next, the user can view the simulation
results, and, if needed, revise the simulation parameters and run another
experimental trial, or build an alternative conceptual model. We describe the
use VERA to develop a SIR model for the spread of COVID-19 and its relationship
with healthcare capacity.


  The coronavirus disease 2019 (COVID-19) pandemic radically impacts our lives,
while the transmission/infection and recovery dynamics of COVID-19 remain
obscure. A time-dependent Susceptible, Exposed, Infectious, and Recovered
(SEIR) model was proposed and applied to fit and then predict the time series
of COVID-19 evolution observed in the last three months (till 3/22/2020) in
various provinces and metropolises in China. The model results revealed the
space dependent transmission/infection rate and the significant spatiotemporal
variation in the recovery rate, likely due to the continuous improvement of
screening techniques and public hospital systems, as well as full city
lockdowns in China. The validated SEIR model was then applied to predict
COVID-19 evolution in United States, Italy, Japan, and South Korea which have
responded differently to monitoring and mitigating COVID-19 so far, although
these predictions contain high uncertainty due to the intrinsic change of the
maximum infected population and the infection/recovery rates within the
different countries. In addition, a stochastic model based on the random walk
particle tracking scheme, analogous to a mixing-limited bimolecular reaction
model, was developed to evaluate non-pharmaceutical strategies to mitigate
COVID-19 spread. Preliminary tests using the stochastic model showed that
self-quarantine may not be as efficient as strict social distancing in slowing
COVID-19 spread, if not all of the infected people can be promptly diagnosed
and quarantined.


  In the present article we introduce an epidemiological model for the
investigation of the spread of epidemics caused by viruses. The model is
applied specifically to COVID-19, the disease caused by the SARS-Cov-2 virus
(aka "novel coronavirus"). The SIR (Susceptible - Infectious - Recovered) model
is used as a basis for studying the evolution of the epidemic. Nevertheless, we
have modified some of the model hypotheses in order to obtain an estimate of
the contamination free of overestimated predictions. This extended model is
then applied to the case of the recent advance of the epidemic in Brazil. In
this regard, it is possible to obtain the evolution for the number of
infectious significantly close to that provided by current data. Accordingly,
we evaluate possible future scenarios for the disease spread. Regarding the
population susceptibility, we consider different social behaviors in response
to quarantine measures and precautions to avoid contagion. We conclude that the
future scenario of the epidemic depends significantly on the social behavior
adopted to date, as well as on the contagion control measures. The extent of
such measures would be likely to cause thousands, millions or tens of millions
of contaminations in the next few months.


  For preventing the spread of epidemics such as the coronavirus disease
COVID-19, social distancing and the isolation of infected persons are crucial.
However, existing reaction-diffusion equations for epidemic spreading are
incapable of describing these effects. We present an extended model for disease
spread based on combining an SIR model with a dynamical density functional
theory where social distancing and isolation of infected persons are explicitly
taken into account. The model shows interesting nonequilibrium phase separation
associated with a reduction of the number of infections, and allows for new
insights into the control of pandemics.


  The emergence of coronavirus disease 2019 (COVID-19) in the United States has
forced federal and local governments to implement containment measures.
Moreover, the severity of the situation has sparked engagement by both the
research and clinical community with the goal of developing effective
treatments for the disease. This article proposes a time dynamic prediction
model with competing risks for the infected individual and develops a simple
tool for policy makers to compare different strategies in terms of when to
implement the strictest containment measures and how different treatments can
increase or suppress infected cases. Two types of containment strategies are
compared: (1) a constant containment strategy that could satisfy the needs of
citizens for a long period; and (2) an adaptive containment strategy whose
strict level changes across time. We consider how an effective treatment of the
disease can affect the dynamics in a pandemic scenario. For illustration we
consider a region with population 2.8 million and 200 initial infectious cases
assuming a 4% mortality rate compared with a 2% mortality rate if a new drug is
available. Our results show compared with a constant containment strategy,
adaptive containment strategies shorten the outbreak length and reduce maximum
daily number of cases. This, along with an effective treatment plan for the
disease can minimize death rate.


  In the context of the ongoing Covid-19 pandemic, several reports and studies
have attempted to model and predict the spread of the disease. There is also
intense debate about policies for limiting the damage, both to health and to
the economy. On the one hand, the health and safety of the population is the
principal consideration for most countries. On the other hand, we cannot ignore
the potential for long-term economic damage caused by strict nation-wide
lockdowns. In this working paper, we present a quantitative way to compute
lockdown decisions for individual cities or regions, while balancing health and
economic considerations. Furthermore, these policies are learnt automatically
by the proposed algorithm, as a function of disease parameters (infectiousness,
gestation period, duration of symptoms, probability of death) and population
characteristics (density, movement propensity). We account for realistic
considerations such as imperfect lockdowns, and show that the policy obtained
using reinforcement learning is a viable quantitative approach towards
lockdowns.


  We discuss how the presence of a large set of asymptomatic infectives changes
our estimate of the COVID-19 basic reproduction number, also known as $R_0$.


  Different countries -- and sometimes different regions within the same
countries -- have adopted different strategies in trying to contain the ongoing
COVID-19 epidemic; these mix in variable parts social confinement, early
detection and contact tracing. In this paper we discuss the different effects
of these ingredients on the epidemic dynamics; the discussion is conducted with
the help of two simple models, i.e. the classical SIR model and the recently
introduced variant (A-SIR, arXiv:2003.08720) which takes into account the
presence of a large set of asymptomatic infectives.


  The Coronavirus, also known as the COVID-19 virus, has emerged in Wuhan China
since late November 2019. Since that time, it has been spreading at large-scale
until today all around the world. It is currently recognized as the world's
most viral and severe epidemic spread in the last twenty years, as compared to
Ebola 2014, MERS 2012, and SARS 2003. Despite being still in the middle of the
outbreak, there is an urgent need to understand the impact of COVID-19. The
objective is to clarify how it was spread so fast in a short time worldwide in
unprecedented fashion. This paper represents a first initiative to achieve this
goal, and it provides a comprehensive analytical study about the Coronavirus.
The contribution of this paper consists in providing descriptive and predictive
models that give insights into COVID-19 impact through the analysis of
extensive data updated daily for the outbreak in all countries. We aim at
answering several open questions: How does COVID-19 spread around the world?
What is its impact in terms of confirmed and death cases at the continent,
region, and country levels? How does its severity compare with other epidemic
outbreaks, including Ebola 2014, MERS 2012, and SARS 2003? Is there a
correlation between the number of confirmed cases and death cases? We present a
comprehensive analytics visualization to address the questions mentioned above.
To the best of our knowledge, this is the first systematic analytical papers
that pave the way towards a better understanding of COVID-19. The analytical
dashboards and collected data of this study are available online [1].


  Epidemiological data on seasonal influenza show that the growth rate of the
number of infected individuals can increase passing from one exponential growth
rate to another one with a larger exponent. Such behavior is not described by
conventional epidemiological models. In this work an immuno-epidemiological
model is proposed in order to describe this two-stage growth. It takes into
account that the growth in the number of infected individuals increases the
initial viral load and provides a passage from the first stage of epidemic
where only people with weak immune response are infected to the second stage
where people with strong immune response are also infected. This scenario may
be viewed as an increase of the effective number of susceptible increasing the
effective growth rate of infected.


  Since the first outbreak of the COVID-19 epidemic at the end of 2019, data
has been made available on the number of infections, deaths and recoveries for
all countries of the World, and that data can be used for statistical analysis.
The primary interest of this paper is how well the logistic equation can
predict the outcome of COVID-19 epidemic in any regions of the World assuming
that the methodology of the testing process, namely the data collection method
and social behavior is not changing over the course of time. Besides the social
relevance, this study has two scientific purposes: we investigate if a simple
saturation model can describe the trend of the COVID-19 epidemic and if so, we
would like to determine, from which point during the epidemic the fitting
parameters provide reliable predictions. We also give estimations for the
outcome of this epidemic in several countries based on the logistic model and
the data available on 27 March, 2020. Based on the saturated cases in China, we
have managed to find some criteria to judge the reliability of the predictions.


  Proteins are key building blocks of virtually all life, providing the
material foundation of spider silk, cells, and hair, but also offering other
functions from enzymes to drugs, and pathogens like viruses. Based on a
nanomechanical analysis of the structure and motions of atoms and molecules at
multiple scales, we report sonified versions of the coronavirus spike protein
of the pathogen of COVID-19, 2019-nCoV. The audio signal, created using a novel
nanomechanical sonification method, features an overlay of the vibrational
signatures of the protein's primary, secondary and higher-order structures.
Presenting musical encoding in two versions - one in the amino-acid scale and
one based on equal temperament tuning - the method allows for expressing
protein structures in audible space, offering novel avenues to represent,
analyze and design architectural features across length- and time-scales. We
further report a hierarchical frequency spectrum analysis of five distinct
protein structures, which offer insights into how genetic mutations, and the
binding of the virus spike protein to the human ACE2 cell receptor directly
influence the audio. Applications of the approach may include the development
of de novo antibodies by designing protein sequences that match, through
melodic counterpoints, the binding sites in the spike protein. Other
applications of audible coding of matter include material design by
manipulating sound, detecting mutations, and offering a way to reach out to
broader communities to explain the physics of proteins. It also forms a
physics-based compositional technique to create new art, referred to as
materiomusic, which is akin to finding a new palette of colors for a painter.
Here, the nanomechanical structure of matter, reflected in an oscillatory
framework, presents a new palette for sound generation, and can complement or
support human creativity.


  The data on number of registered cases of COVID-19 disease in twenty European
countries is analyzed by the least-squares fitting procedure with generic
analytic functions. Three regimes of the expansion of the disease are
identified and quantified -- early exponential expansion, damped exponential,
and linear expansion. Differences among countries in the early expansion period
are quantified. The velocity of the expansion in the exponential regime lies
within one standard deviation from the average value for 11 countries. The
number of infected individuals at the initial time is excessively high for
Italy, 7 standard deviations from the average value.
  Method for predicting the expansion based on extrapolation in the parametric
space is presented. One-week predictions based on extrapolations have average
precision of 18% and 29% during the later period of the damped exponential
expansion for the case of Italy and Czechia, respectively. The method based on
extrapolations in the parametric space may provide an elementary method to
quantify the impact of restrictive measures on the spreading of the disease.


  Current available data of the worldwide impact of the COVID-19 pandemic has
been analyzed using dimensional analysis and self-similarity hypotheses. We
show that the time series of infected population and deaths of the most
impacted and unprepared countries exhibits an asymptotic power law behavior,
compatible with the propagation of a signal in a fractal network. We propose a
model which predicts an asymptotically self-similar expansion of deaths in time
before containment, and the final death toll under total containment measures,
as a function of the delay in taking those measures after the expansion is
observed. The physics of the model resembles the expansion of a flame in a
homogeneous domain with a fractal dimension 3.75. After containment measures
are taken, the natural fractal structure of the network is drastically altered
and a secondary evolution is observed. This evolution, akin to the homogeneous
combustion in a static isolated enclosure with a final quenching, has a
characteristic time of 20.1 days, according to available data of the pandemic
behavior in China. The proposed model is remarkably consistent with available
data, which supports the simplifying hypotheses made in the model. A universal
formulation for a quarantine as a function of that delay is also proposed.


  We model and forecast the early evolution of the COVID-19 pandemic in Brazil
using Brazilian recent data from February 25, 2020 to March 30, 2020. This
early period accounts for unawareness of the epidemiological characteristics of
the disease in a new territory, sub-notification of the real numbers of
infected people and the timely introduction of social distancing policies to
flatten the spread of the disease. We use two variations of the SIR model and
we include a parameter that comprises the effects of social distancing
measures. Short and long term forecasts show that the social distancing policy
imposed by the government is able to flatten the pattern of infection of the
COVID-19. However, our results also show that if this policy does not last
enough time, it is only able to shift the peak of infection into the future
keeping the value of the peak in almost the same value. Furthermore, our long
term simulations forecast the optimal date to end the policy. Finally, we show
that the proportion of asymptomatic individuals affects the amplitude of the
peak of symptomatic infected, suggesting that it is important to test the
population.


  The COVID-19 pandemic has affected all countries of the world producing a
substantial number of fatalities accompanied by a major disruption in their
social, financial, and educational organization. The strict disciplinary
measures implemented by China were very effective and thus were subsequently
adopted by most world countries to various degrees. The infection duration and
number of infected persons are of critical importance for the battle against
the pandemic. We use the quantitative landscape of the disease spreading in
China as a benchmark and utilize infection data from eight countries to
estimate the complete evolution of the infection in each of these countries.
The analysis predicts successfully both the expected number of daily infections
per country and, perhaps more importantly, the duration of the epidemic in each
country. Our quantitative approach is based on a Gaussian spreading hypothesis
that is shown to arise as a result of imposed measures in a simple dynamical
infection model. This may have consequences and shed light in the efficiency of
policies once the phenomenon is over.


  When testing for infections, the standard method is to test each subject
individually. If testing methodology is such that samples from multiple
subjects can be efficiently combined and tested at once, yielding a positive
results if any one subject in the subgroup is positive, then one can often
identify the infected sub-population with a considerably lower number of tests
compared to the number of test subjects. We present two such methods that allow
an increase in testing efficiency (in terms of total number of test performed)
by a factor of $\approx$ 10 if population infection rate is $10^{-2}$ and a
factor of $\approx$50 when it is $10^{-3}$. Such methods could be useful when
testing large fractions of the total population, as will be perhaps required
during the current coronavirus pandemic.


  Uncompetitive antagonists of the N-methyl D-aspartate receptor (NMDAR) have
demonstrated therapeutic benefit in the treatment of neurological diseases such
as Parkinson's and Alzheimer's, but some also cause dissociative effects that
have led to the synthesis of illicit drugs. The ability to generate NMDAR
antagonists in silico is therefore desirable both for new medication
development and for preempting and identifying new designer drugs. Recently,
generative deep learning models have been applied to de novo drug design as a
means to expand the amount of chemical space that can be explored for potential
drug-like compounds. In this study, we assess the application of a generative
model to the NMDAR to achieve two primary objectives: (i) the creation and
release of a comprehensive library of experimentally validated NMDAR
phencyclidine (PCP) site antagonists to assist the drug discovery community and
(ii) an analysis of both the advantages conferred by applying such generative
artificial intelligence models to drug design and the current limitations of
the approach. We apply, and provide source code for, a variety of ligand- and
structure-based assessment techniques used in standard drug discovery analyses
to the deep learning-generated compounds. We present twelve candidate
antagonists that are not available in existing chemical databases to provide an
example of what this type of workflow can achieve, though synthesis and
experimental validation of these compounds is still required.


  We study the perception of COVID-2019 epidemic in Polish society using
quantitative analysis of its digital footprints on the Internet (on Twitter,
Google, YouTube, Wikipedia and electronic media represented by Event Registry)
from January 2020 to 12.03.2020 (before and after official introduction to
Poland on 04.03.2020). To this end we utilize data mining, social network
analysis, natural language processing techniques. Each examined internet
platform was analyzed for representativeness and composition of the target
group. We identified three temporal major cluster of the interest before
disease introduction on the topic COVID-2019: China- and Italy-related peaks on
all platforms, as well as a peak on social media related to the recent special
law on combating COVID-2019. Besides, there was a peak in interest on the day
of officially confirmed introduction as well as an exponential increase of
interest when the Polish government declared war against disease with a massive
mitigation program. From sociolingistic perspective, we found that concepts and
issues of threat, fear and prevention prevailed before introduction. After
introduction, practical concepts about disease and epidemic dominate. We have
found out that Twitter reflected the structural division of the Polish
political sphere. We were able to identify clear communities of governing
party, mainstream oppostition and protestant group and potential sources of
misinformation. We have also detected bluring boundaries between comminities
after disease introduction.


  In this work we analyse the growth of the cumulative number of confirmed
infected cases by the COVID-19 until March 27th, 2020, from countries of Asia,
Europe, North and South America. Our results show (i) that power-law growth is
observed for all countries; (ii) by using the distance correlation, that the
power-law curves between countries are statistically highly correlated,
suggesting the universality of such curves around the World; and (iii) that
soft quarantine strategies are inefficient to flatten the growth curves.
Furthermore, we present a model and strategies which allow the government to
reach the flattening of the power-law curves. We found that, besides the social
distance of individuals, of well known relevance, the strategy of identifying
and isolating infected individuals in a large daily rate can help to flatten
the power-laws. These are essentially the strategies used in the Republic of
Korea. The high correlation between the power-law curves of different countries
strongly indicate that the government containment measures can be applied with
success around the whole World. These measures must be scathing and applied as
soon as possible.


  Motivated by the rapid spread of COVID-19 all across the globe, we have
performed simulations of a system dynamic epidemic spread model in different
possible situations. The simulation, not only captures the model dynamic of the
spread of the virus, but also, takes care of population and mobility data. The
model is calibrated based on epidemic data and events specifically of Pakistan,
which can easily be generalized. The simulation results are quite disturbing,
indicating that, even with stringent social distancing and testing strategies
and for a quite long time (even beyond one year), the spread would be
significant (in tens of thousands). The real alarm is when some of these
measures got leaked for a short time within this duration, which may result in
catastrophic situation when millions of people would be infected.


  The concept of $k$-compatibility measures how many phylogenetic trees it
would take to display all splits in a given set. A set of trees that display
every single possible split is termed a \textit{universal tree set}. In this
note, we find $A(n)$, the minimal size of a universal tree set for $n$ taxa. By
normalising the $k$-compatibility using $A(n)$, one can then compare
incompatibility of split systems across different taxa sizes. We demonstrate
this application by comparing two SplitsTree networks of different sizes
derived from archaeal genomes.


  The epidemic curve and the final extent of the COVID-19 pandemic are usually
predicted from the rate of early exponential raising using the SIR model. These
predictions implicitly assume a full social mixing, which is not plausible
generally. Here I am showing a counterexample to the these predictions, based
on random propagation of an epidemic in Barab\'asi--Albert scale-free network
models. The start of the epidemic suggests $R_0=2.6$, but unlike $\Omega\approx
70\%{}$ predicted by the SIR model, they reach a final extent of only
$\Omega\approx 4\%{}$ without external mitigation and $\Omega\approx
0.5$--$1.5\%{}$ with mitigation. Daily infection rate at the top is also 1--1.5
orders of magnitude less than in SIR models. Quarantining only the 1.5\%{} most
active superspreaders has similar effect on extent and top infection rate as
blind quarantining a random 50\%{} of the full community.


  We present in this work some results from analysing the spread of Covid-19 in
different countries and regions around the world and the potential relations
with climate, geographical location, and GDP. While the situation remains
dynamic, we believe this analysis has the potential to uncover certain
underlying trends. We primarily intend the results to drive further, more
detailed analysis of the relevant data by other researchers that would help us
gather a better understanding of the situation, aiding our preparedness. Our
observations show that countries in high latitudes, with temperate and/or
continental climate, and with reasonably "open" economies are the most
vulnerable to this outbreak.


  Early assessments of the spreading rate of COVID-19 were subject to
significant uncertainty, as expected with limited data and difficulties in case
ascertainment, but more reliable inferences can now be made. Here, we estimate
from European data that COVID-19 cases are expected to double initially every
three days, until social distancing interventions slow this growth, and that
the impact of such measures is typically only seen nine days - i.e. three
doubling times - after their implementation. We argue that such temporal
patterns are more critical than precise estimates of the basic reproduction
number for initiating interventions. This observation has particular
implications for the low- and middle-income countries currently in the early
stages of their local epidemics.


  Stem cells can precisely and robustly undergo cellular differentiation and
lineage commitment, referred to as stemness. However, how the gene network
underlying stemness regulation reliably specifies cell fates is not well
understood. To address this question, we applied a recently developed
computational method, Random Circuit Perturbation (RACIPE), to a nine-component
gene regulatory network (GRN) governing stemness, from which we identified
fifteen robust gene states. Among them, four out of the five most probable gene
states exhibit gene expression patterns observed in single mouse embryonic
cells at 32-cell and 64-cell stages. These gene states can be robustly
predicted by the stemness GRN but not by randomized versions of the stemness
GRN. Strikingly, we found a hierarchical structure of the GRN with the
Oct4/Cdx2 motif functioning as the first decision-making module followed by
Gata6/Nanog. We propose that stem cell populations, instead of being viewed as
all having a specific cellular state, can be regarded as a heterogeneous
mixture including cells in various states. Upon perturbations by external
signals, stem cells lose the capacity to access certain cellular states,
thereby becoming differentiated. The findings demonstrate that the functions of
the stemness GRN is mainly determined by its well-evolved network topology
rather than by detailed kinetic parameters.


  COVID-19, a member of corona virus family is spreading its tentacles across
the world due to lack of drugs at present. However, the main viral proteinase
(Mpro/3CLpro) has recently been regarded as a suitable target for drug design
against SARS infection due to its vital role in polyproteins processing
necessary for coronavirus reproduction. The present in silico study was
designed to evaluate the effect of Jensenone, a essential oil component from
eucalyptus oil, on Mpro by docking study. In the present study, molecular
docking studies were conducted by using 1-click dock and swiss dock tools.
Protein interaction mode was calculated by Protein Interactions Calculator.The
calculated parameters such as binding energy, and binding site similarity
indicated effective binding of Jensenone to COVID-19 proteinase. Active site
prediction further validated the role of active site residues in ligand
binding. PIC results indicated that, Mpro/ Jensenone complexes forms
hydrophobic interactions, hydrogen bond interactions and strong ionic
interactions. Therefore, Jensenone may represent potential treatment potential
to act as COVID-19 Mpro inhibitor. However, further research is necessary to
investigate their potential medicinal use.


  Direct reciprocity is one of the key mechanisms accounting for cooperation in
our social life. According to recent understanding, most of classical
strategies for direct reciprocity fall into one of two classes, `partners' or
`rivals'. A `partner' is a generous strategy achieving mutual cooperation, and
a `rival' never lets the co-player become better off. They have different
working conditions: For example, partners show good performance in a large
population, whereas rivals do in head-to-head matches. By means of exhaustive
enumeration, we demonstrate the existence of strategies that act as both
partners and rivals. Among them, we focus on a human-interpretable strategy,
named `CAPRI' after its five characteristic ingredients, i.e., cooperate,
accept, punish, recover, and defect otherwise. Our evolutionary simulation
shows excellent performance of CAPRI in a broad range of environmental
conditions.


  Feature selection for predictive analytics is the problem of identifying a
minimal-size subset of features that is maximally predictive of an outcome of
interest. To apply to molecular data, feature selection algorithms need to be
scalable to tens of thousands of available features. In this paper, we propose
gOMP, a highly-scalable generalisation of the Orthogonal Matching Pursuit
feature selection algorithm to several directions: (a) different types of
outcomes, such as continuous, binary, nominal, and time-to-event, (b) different
types of predictive models (e.g., linear least squares, logistic regression),
(c) different types of predictive features (continuous, categorical), and (d)
different, statistical-based stopping criteria. We compare the proposed
algorithm against LASSO, a prototypical, widely used algorithm for
high-dimensional data. On dozens of simulated datasets, as well as, real gene
expression datasets, gOMP is on par, or outperforms LASSO for case-control
binary classification, quantified outcomes (regression), and (censored)
survival times (time-to-event) analysis. gOMP has also several theoretical
advantages that are discussed. While gOMP is based on quite simple and basic
statistical ideas, easy to implement and to generalize, we also show in an
extensive evaluation that it is also quite effective in bioinformatics analysis
settings.


  A key challenge to make effective use of evolutionary algorithms is to choose
appropriate settings for their parameters. However, the appropriate parameter
setting generally depends on the structure of the optimisation problem, which
is often unknown to the user. Non-deterministic parameter control mechanisms
adjust parameters using information obtained from the evolutionary process.
Self-adaptation -- where parameter settings are encoded in the chromosomes of
individuals and evolve through mutation and crossover -- is a popular parameter
control mechanism in evolutionary strategies. However, there is little
theoretical evidence that self-adaptation is effective, and self-adaptation has
largely been ignored by the discrete evolutionary computation community.
  Here we show through a theoretical runtime analysis that a non-elitist,
discrete evolutionary algorithm which self-adapts its mutation rate not only
outperforms EAs which use static mutation rates on \leadingones, but also
improves asymptotically on an EA using a state-of-the-art control mechanism.
The structure of this problem depends on a parameter $k$, which is \emph{a
priori} unknown to the algorithm, and which is needed to appropriately set a
fixed mutation rate. The self-adaptive EA achieves the same asymptotic runtime
as if this parameter was known to the algorithm beforehand, which is an
asymptotic speedup for this problem compared to all other EAs previously
studied. An experimental study of how the mutation-rates evolve show that they
respond adequately to a diverse range of problem structures.
  These results suggest that self-adaptation should be adopted more broadly as
a parameter control mechanism in discrete, non-elitist evolutionary algorithms.


  RNA function crucially depends on its structure. Thermodynamic models
currently used for secondary structure prediction rely on computing the
partition function of folding ensembles, and can thus estimate minimum
free-energy structures and ensemble populations. These models sometimes fail in
identifying native structures unless complemented by auxiliary experimental
data. Here, we build a set of models that combine thermodynamic parameters,
chemical probing data (DMS, SHAPE), and co-evolutionary data (Direct Coupling
Analysis, DCA) through a network that outputs perturbations to the ensemble
free energy. Perturbations are trained to increase the ensemble populations of
a representative set of known native RNA structures. In the chemical probing
nodes of the network, a convolutional window combines neighboring reactivities,
enlightening their structural information content and the contribution of local
conformational ensembles. Regularization is used to limit overfitting and
improve transferability. The most transferable model is selected through a
cross-validation strategy that estimates the performance of models on systems
on which they are not trained. With the selected model we obtain increased
ensemble populations for native structures and more accurate predictions in an
independent validation set. The flexibility of the approach allows the model to
be easily retrained and adapted to incorporate arbitrary experimental
information.


  Coronavirus disease (COVID-2019) is a severe ongoing novel pandemic that is
spreading quickly across the world. Italy, that is widely considered one of the
main epicenters of the pandemic, has registered the highest COVID-2019 death
rates and death toll in the world, to the present day. In this article I
estimate an autoregressive integrated moving average (ARIMA) model to forecast
the epidemic trend over the period after April 4, 2020, by using the Italian
epidemiological data at national and regional level. The data refer to the
number of daily confirmed cases officially registered by the Italian Ministry
of Health (www.salute.gov.it) for the period February 20 to April 4, 2020. The
main advantage of this model is that it is easy to manage and fit. Moreover, it
may give a first understanding of the basic trends, by suggesting the
hypothetic epidemic's inflection point and final size.


  Key traits of unicellular species, like cell size, often follow scale-free or
self-similar distributions, hinting at the possibility of an underlying
critical process. However, linking such empirical scaling laws to the critical
regime of realistic individual-based model classes is difficult. Here we reveal
new empirical scaling evidence associated with a transition in the population
and chlorophyll dynamics of phytoplankton. We offer a possible explanation for
these observations by deriving scaling laws in the vicinity of the critical
point of a new universality class of non-local cell growth and division models.
This `criticality hypothesis' can be tested through new scaling predictions
derived for our model class, for the response of chlorophyll distributions to
perturbations. The derived scaling laws may also be generalized to other
cellular traits and environmental drivers relevant to phytoplankton ecology.


  The outbreak of COVID-19 caused by SARS-CoV-2 in Wuhan and other cities in
China in 2019 has become a global pandemic as declared by World Health
Organization (WHO) in the first quarter of 2020 . The delay in diagnosis,
limited hospital resources and other treatment resources leads to rapid spread
of COVID-19. In this article, we consider an optimal control COVID-19
transmission model and assess the impact of some control measures that can lead
to the reduction of exposed and infectious individuals in the population. We
investigate three control strategies for this deadly infectious disease using
personal protection, treatment with early diagnosis, treatment with delay
diagnosis and spraying of virus in the environment as time-dependent control
functions in our dynamical model to curb the disease spread.


  The rapid spread of the Coronavirus (COVID-19) confronts policy makers with
the problem of measuring the effectiveness of containment strategies, balancing
public health considerations with the economic costs of social distancing
measures. We introduce a modified epidemic model that we name the
controlled-SIR model, in which the disease reproduction rate evolves
dynamically in response to political and societal reactions. An analytic
solution is presented. The model reproduces official COVID-19 cases counts of a
large number of regions and countries that surpassed the first peak of the
outbreak. A single unbiased feedback parameter is extracted from field data and
used to formulate an index that measures the efficiency of containment
strategies (the CEI index). CEI values for a range of countries are given. For
two variants of the controlled-SIR model, detailed estimates of the total
medical and socio-economic costs are evaluated over the entire course of the
epidemic. Costs comprise medical care cost, the economic cost of social
distancing, as well as the economic value of lives saved. Under plausible
parameters, strict measures fare better than a hands-off policy. Strategies
based on current case numbers lead to substantially higher total costs than
strategies based on the overall history of the epidemic.


  The classification of the long-term behavior of dynamical systems is a
fundamental problem in mathematics. For both deterministic and stochastic
dynamics specific classes of models verify Palis' conjecture: the long-term
behavior is determined by a finite number of stationary distributions. In this
paper we consider the classification problem for stochastic models of
interacting species. For a large class of three-species, stochastic
differential equation models, we prove a variant of Palis' conjecture: the
long-term statistical behavior is determined by a finite number of stationary
distributions and, generically, three general types of behavior are possible:
1) convergence to a unique stationary distribution that supports all species,
2) convergence to one of a finite number of stationary distributions supporting
two or fewer species, 3) convergence to convex combinations of single species,
stationary distributions due to a rock-paper-scissors type of dynamic.
Moreover, we prove that the classification reduces to computing Lyapunov
exponents (external Lyapunov exponents) that correspond to the average
per-capita growth rate of species when rare. Our results stand in contrast to
the deterministic setting where the classification is incomplete even for
three-dimensional, competitive Lotka--Volterra systems. For these SDE models,
our results also provide a rigorous foundation for ecology's modern coexistence
theory (MCT) which assumes the external Lyapunov exponents determine long-term
ecological outcomes.


  Collective behavior, and swarm formation in particular, has been studied from
several perspectives within a large variety of fields, ranging from biology to
physics. In this work, we apply Projective Simulation to model each individual
as an artificial learning agent that interacts with its neighbors and
surroundings in order to make decisions and learn from them. Within a
reinforcement learning framework, we discuss one-dimensional learning scenarios
where agents need to get to food resources to be rewarded. We observe how
different types of collective motion emerge depending on the distance the
agents need to travel to reach the resources. For instance, strongly aligned
swarms emerge when the food source is placed far away from the region where
agents are situated initially. In addition, we study the properties of the
individual trajectories that occur within the different types of emergent
collective dynamics. Agents trained to find distant resources exhibit
individual trajectories with L\'evy-like characteristics as a consequence of
the collective motion, whereas agents trained to reach nearby resources present
Brownian-like trajectories.


  The aim of this paper consists in the application of a recent epidemiological
model, namely SEIR with Social Distancing (SEIR--SD), extended here through the
definition of a social distancing function varying over time, to assess the
situation related to the spreading of the coronavirus Covid--19 in Italy and in
two of its most important regions, i.e., Lombardy and Campania. To profitably
use this model, the most suitable values of its parameters must be found. The
estimation of the SEIR--SD model parameters takes place here through the use of
Differential Evolution, a heuristic optimization technique. In this way, we are
able to evaluate for each of the three above-mentioned scenarios the daily
number of infectious cases from today until the end of virus spreading, the
day(s) in which this number will be at its highest peak, and the day in which
the infected cases will become very close to zero.


  Unfolding of a coarse grained COVN protein from its native configuration
shows a linear response with increasing temperature followed by a nonmonotonic
double peaks in its radius of gyration. The protein conforms to a random coil
of folded segments in native state with increasing tenuous and globular
structures in specific temperature regimes where the effective dimensions of
corresponding structures D is about 1.6 to 2.4. Thermal agitation alone is not
sufficient to fully eradicate its segmental folding as few folds are found to
persist around such residues as 65W, 110Y, 224L, 374P even at high
temperatures.


  The scope of this work is to serve as a guiding tool against subjective
estimations on real pandemic situations (mainly due to the inability to acquire
objective real data over whole populations). The previously introduced model of
closed self-organized criticality (SOC), is adapted in the case of a
virus-induced epidemic. In this version this physical model can distinguish the
virus spread according to the virus aggressiveness. The study presented,
highlights the critical value of virus density over a population. For low
values of the initial virus density (lower than the critical value) it is
proved that the virus-diffusion behavior is safe and quantitatively similar to
usual real epidemical data. However, it is revealed that very close to the
critical point, the critical slowing-down (CSD) phenomenon, introduced by the
theory of critical phenomena, emerges, leading to a tremendous increase of both
the percentage of active carriers and the duration of the epidemic. A behavior
of the epidemic obeying to a second order phase transition, also occurs. For
virus density values higher than the critical value, the epidemic duration
becomes extremely prolonged. Additionally, the effect of the closed system
population size revealed interesting properties. All these results, together
with an investigation of the effectiveness of applying physical contact
restriction measures, document scientifically their worthiness, while they also
demonstrate the limits for which herd immunity holds safely. Finally, the model
has been compared against real epidemic data in the case of Greece, which
imposed restrictive measures consistently and in time.


  This work covers the analysis of the COVID 19 spread in different countries
and dealing the main feature of COVID 19 growth, which is the spread due to the
social-contact structure, which is governed by the parameter \b{eta}. The
dependency of this parameter \b{eta} on the transmission level in society gives
a sense of the effectiveness of the measures taken for social distancing. A
separate algorithm is hardcoded in python using Scipy which learns the
social-contact structure and gives a suitable value for \b{eta}, which has a
major impact on the outcome of the result. Forecasting for the epidemic spread
in India was done, and it was found that the strictness at which social
distancing in India is done, is insufficient for the growth of COVID 19.


  Coronaviruses have recently caused world-wide severe outbreaks: SARS (Severe
Acute Respiratory Syndrome) in 2002 and MERS (Middle-East Respiratory Syndrome)
in 2012. At the end of 2019, a new coronavirus outbreak appeared in Wuhan
(China) seafood market as first focus of infection, becoming a pandemics in
2020, spreading mainly into Europe and Asia. Although the virus family is
well-known, this specific virus presents considerable differences, as higher
transmission rates, being a challenge for diagnostic methods, treatments and
vaccines. Coronavirus(C++).pro is a C++ application which simulates Coronavirus
replication cycle. This software has identified virus type in short times and
provided FASTA files of virus proteins, a list of mRNA sequences and secondary
structures. Furthermore, the software has identified a list of structural,
non-structural and accessory proteins in 2019-nCoV virus genome more similar to
SARS than to MERS, as several fusion proteins characteristics of this virus
type. These results are useful as a first step in order to develop diagnostic
methods, new vaccines or antiviral drugs, which could avoid virus replication
in any stage: fusion inhibitors, RdRp inhibitors and PL2pro/3CLpro protease
inhibitors.


  As the coronavirus disease 2019 (COVID-19) continues to be a global pandemic,
policy makers have enacted and reversed non-pharmaceutical interventions with
various levels of restrictions to limit its spread. Data driven approaches that
analyze temporal characteristics of the pandemic and its dependence on regional
conditions might supply information to support the implementation of mitigation
and suppression strategies. To facilitate research in this direction on the
example of the United States, we present a machine-readable dataset that
aggregates relevant data from governmental, journalistic, and academic sources
on the U.S. county level. In addition to county-level time-series data from the
JHU CSSE COVID-19 Dashboard, our dataset contains more than 300 variables that
summarize population estimates, demographics, ethnicity, housing, education,
employment and income, climate, transit scores, and healthcare system-related
metrics. Furthermore, we present aggregated out-of-home activity information
for various points of interest for each county, including grocery stores and
hospitals, summarizing data from SafeGraph and Google mobility reports. We
compile information from IHME, state and county-level government, and
newspapers for dates of the enactment and reversal of non-pharmaceutical
interventions. By collecting these data, as well as providing tools to read
them, we hope to accelerate research that investigates how the disease spreads
and why spread may be different across regions. Our dataset and associated code
are available at github.com/JieYingWu/COVID-19_US_County-level_Summaries.


  Critical thinking and skepticism are fundamental mechanisms that one may use
to prevent the spreading of rumors, fake-news and misinformation. We consider a
simple model in which agents without previous contact with the rumor, being
skeptically oriented, may convince spreaders to stop their activity or, once
exposed to the rumor, decide not to propagate it as a consequence, for example,
of fact-checking. We extend a previous, mean-field analysis of the combined
effect of these two mechanisms, active and passive skepticism, to include
spatial correlations. This can be done either analytically, through the pair
approximation, or simulating an agent-based version on diverse networks. Our
results show that while in mean-field there is no coexistence between spreaders
and susceptibles (although, depending on the parameters, there may be
bistability depending on the initial conditions), when spatial correlations are
included, because of the protective effect of the isolation provided by removed
agents, coexistence is possible.


  Targeted immunization or attacks of large-scale networks has attracted
significant attention by the scientific community. However, in real-world
scenarios, knowledge and observations of the network may be limited thereby
precluding a full assessment of the optimal nodes to immunize (or remove) in
order to avoid epidemic spreading such as that of current COVID-19 epidemic.
Here, we study a novel immunization strategy where only $n$ nodes are observed
at a time and the most central between these $n$ nodes is immunized (or
attacked). This process is continued repeatedly until $1-p$ fraction of nodes
are immunized (or attacked). We develop an analytical framework for this
approach and determine the critical percolation threshold $p_c$ and the size of
the giant component $P_{\infty}$ for networks with arbitrary degree
distributions $P(k)$. In the limit of $n\to\infty$ we recover prior work on
targeted attack, whereas for $n=1$ we recover the known case of random failure.
Between these two extremes, we observe that as $n$ increases, $p_c$ increases
quickly towards its optimal value under targeted immunization (attack) with
complete information. In particular, we find a new scaling relationship between
$|p_c(\infty)-p_c(n)|$ and $n$ as $|p_c(\infty)-p_c(n)|\sim n^{-1}\exp(-\alpha
n)$. For Scale-free (SF) networks, where $P(k)\sim k^{-\gamma}, 2<\gamma<3$, we
find that $p_c$ has a transition from zero to non-zero when $n$ increases from
$n=1$ to order of $\log N$ ($N$ is the size of network). Thus, for SF networks,
knowledge of order of $\log N$ nodes and immunizing them can reduce
dramatically an epidemics.


  In population genetics, extant samples are usually used for inference of past
population genetic forces. With the Kingman coalescent and the backward
diffusion equation, inference of the marginal likelihood proceeds from an
extant sample backward in time. Conditional on an extant sample, the Moran
model can also be used backward in time with identical results, up to a scaling
of time. In particular, all three approaches -- the coalescent, the backward
diffusion, and the Moran model -- lead to the identical marginal likelihood of
the sample. If probabilities of ancestral states are also inferred, either of
discrete ancestral allele particle configurations, as in the coalescent, or of
ancestral population allele proportions, as in the backward diffusion, the
backward algorithm needs to be combined with the corresponding forward
algorithm to the forward-backward algorithm. Generally orthogonal polynomials,
solving the diffusion equation, are numerically simpler than the other
approaches: they implicitly sum over many intermediate ancestral particle
configurations; furthermore, while the Moran model requires iterative matrix
multiplication with a transition matrix of a dimension of the population size
squared, expansion of the polynomials is only necessary up to the sample size.
For discrete samples, forward-in-time moving pure birth processes similar to
the Polya- or Hoppe-urn models complement the backward-looking coalescent.
Because, the sample size is a random variable forward in time, pure-birth
processes are unsuited to model population demography given extant samples.
With orthogonal polynomials, however, not only ancestral allele proportions but
also probabilities of ancestral particle configurations can be calculated
easily. Assuming only mutation and drift, the use of orthogonal polynomials is
numerically advantageous over alternative strategies.


  Emerging infectious diseases are existential threats to human health and
global stability. The recent outbreaks of the novel coronavirus COVID-19 have
rapidly formed a global pandemic, causing hundreds of thousands of infections
and huge economic loss. The WHO declares that more precise measures to track,
detect and isolate infected people are among the most effective means to
quickly contain the outbreak. Based on trajectory provided by the big data and
the mean field theory, we establish an aggregated risk mean field that contains
information of all risk-spreading particles by proposing a spatio-temporal
model named HiRES risk map. It has dynamic fine spatial resolution and high
computation efficiency enabling fast update. We then propose an objective
individual epidemic risk scoring model named HiRES-p based on HiRES risk maps,
and use it to develop statistical inference and machine learning methods for
detecting suspected epidemic-infected individuals. We conduct numerical
experiments by applying the proposed methods to study the early outbreak of
COVID-19 in China. Results show that the HiRES risk map has strong ability in
capturing global trend and local variability of the epidemic risk, thus can be
applied to monitor epidemic risk at country, province, city and community
levels, as well as at specific high-risk locations such as hospital and
station. HiRES-p score seems to be an effective measurement of personal
epidemic risk. The accuracy of both detecting methods are above 90\% when the
population infection rate is under 20\%, which indicates great application
potential in epidemic risk prevention and control practice.


  Deterministic continuum models formulated in terms of non-local partial
differential equations for the evolutionary dynamics of populations structured
by phenotypic traits have been used recently to address open questions
concerning the adaptation of asexual species to periodically fluctuating
environmental conditions. These deterministic continuum models are usually
defined on the basis of population-scale phenomenological assumptions and
cannot capture adaptive phenomena that are driven by stochastic variability in
the evolutionary paths of single individuals. In this paper, we develop a
stochastic individual-based model for the coevolution between two competing
phenotype-structured cell populations that are exposed to time-varying nutrient
levels and undergo spontaneous, heritable phenotypic variations with different
probabilities. The evolution of every cell is described by a set of rules that
result in a discrete-time branching random walk on the space of phenotypic
states. We formally show that the deterministic continuum counterpart of this
model comprises a system of non-local partial differential equations for the
cell population density functions coupled with an ordinary differential
equation for the nutrient concentration. We compare the individual-based model
and its continuum analogue, focussing on scenarios whereby the predictions of
the two models differ. Our results clarify the conditions under which
significant differences between the two models can emerge due to stochastic
effects associated with small population levels. These differences arise in the
presence of low probabilities of phenotypic variation, and become more apparent
when the two populations are characterised by less fit initial mean phenotypes
and smaller initial levels of phenotypic heterogeneity.


  COVID-19 is an emerging respiratory infectious disease caused by the
coronavirus SARS-CoV-2. It was first reported on in early December 2019 in
Wuhan, China and within three month spread as a pandemic around the whole
globe. Here, we study macro-epidemiological patterns along the time course of
the pandemic. We compute the distribution of confirmed COVID-19 cases and
deaths for countries worldwide and for counties in the US, and show that both
distributions follow a truncated power-law over five orders of magnitude. We
are able to explain the origin of this scaling behavior as a dual-scale
process: the large-scale spread of the virus between countries and the
small-scale accumulation of case numbers within each country. Assuming
exponential growth on both scales, the critical exponent of the power-law is
determined by the ratio of large-scale to small-scale growth rates. We confirm
this theory in numerical simulations in a simple meta-population model,
describing the epidemic spread in a network of interconnected countries. Our
theory gives a mechanistic explanation why most COVID-19 cases occurred within
a few epicenters, at least in the initial phase of the outbreak. Assessing how
well a simple dual-scale model predicts the early spread of epidemics, despite
the huge contrasts between countries, could help identify critical temporal and
spatial scales of response in which to mitigate future epidemic threats.


  COVID-19 pandemic has become a major threat to the country. Till date, well
tested medication or antidote is not available to cure this disease. According
to WHO reports, COVID-19 is a severe acute respiratory syndrome which is
transmitted through respiratory droplets and contact routes. Analysis of this
disease requires major attention by the Government to take necessary steps in
reducing the effect of this global pandemic. In this study, outbreak of this
disease has been analysed for India till 30th March 2020 and predictions have
been made for the number of cases for the next 2 weeks. SEIR model and
Regression model have been used for predictions based on the data collected
from John Hopkins University repository in the time period of 30th January 2020
to 30th March 2020. The performance of the models was evaluated using RMSLE and
achieved 1.52 for SEIR model and 1.75 for the regression model. The RMSLE error
rate between SEIR model and Regression model was found to be 2.01. Also, the
value of R0 which is the spread of the disease was calculated to be 2.02.
Expected cases may rise between 5000-6000 in the next two weeks of time. This
study will help the Government and doctors in preparing their plans for the
next two weeks. Based on the predictions for short-term interval, these models
can be tuned for forecasting in long-term intervals.


  The recent worldwide outbreak of the novel coronavirus (COVID-19) has opened
up new challenges to the research community. Artificial intelligence (AI)
driven methods can be useful to predict the parameters, risks, and effects of
such an epidemic. Such predictions can be helpful to control and prevent the
spread of such diseases. The main challenges of applying AI is the small volume
of data and the uncertain nature. Here, we propose a shallow long short-term
memory (LSTM) based neural network to predict the risk category of a country.
We have used a Bayesian optimization framework to optimize and automatically
design country-specific networks. The results show that the proposed pipeline
outperforms state-of-the-art methods for data of 180 countries and can be a
useful tool for such risk categorization. We have also experimented with the
trend data and weather data combined for the prediction. The outcome shows that
the weather does not have a significant role. The tool can be used to predict
long-duration outbreak of such an epidemic such that we can take preventive
steps earlier


  Due to the current severe acute respiratory syndrome coronavirus 2
(SARS-CoV-2) pandemic, there is an urgent need for novel therapies and drugs.
We conducted a large-scale virtual screening for small molecules that are
potential CoV-2 inhibitors. To this end, we utilized "ChemAI", a deep neural
network trained on more than 220M data points across 3.6M molecules from three
public drug-discovery databases. With ChemAI, we screened and ranked one
billion molecules from the ZINC database for favourable effects against CoV-2.
We then reduced the result to the 30,000 top-ranked compounds, which are
readily accessible and purchasable via the ZINC database. Additionally, we
screened the DrugBank using ChemAI to allow for drug repurposing, which would
be a fast way towards a therapy. We provide these top-ranked compounds of ZINC
and DrugBank as a library for further screening with bioassays at
https://github.com/ml-jku/sars-cov-inhibitors-chemai.


  Biological membranes are shaped by various proteins that either generate
inward or outward membrane curvature. In this article, we investigate the
membrane morphologies induced by mixtures of arc-shaped particles with
coarse-grained modeling and simulations. The particles bind to the membranes
either with their inward, concave side or their outward, convex side and, thus,
generate membrane curvature of opposite sign. We find that small fractions of
convex-binding particles can stabilize three-way junctions of membrane tubules,
as suggested for the protein lunapark in the endoplasmic reticulum of cells.
For comparable fractions of concave-binding and convex-binding particles, we
observe lines of particles of the same type, and diverse membrane morphologies
with grooves and bulges induced by these particle lines. The alignment and
segregation of the particles is driven by indirect, membrane-mediated
interactions.


  With the booming of next generation sequencing technology and its
implementation in clinical practice and life science research, the need for
faster and more efficient data analysis methods becomes pressing in the field
of sequencing. Here we report on the evaluation of an optimized germline
mutation calling pipeline, HummingBird, by assessing its performance against
the widely accepted BWA-GATK pipeline. We found that the HummingBird pipeline
can significantly reduce the running time of the primary data analysis for
whole genome sequencing and whole exome sequencing while without significantly
sacrificing the variant calling accuracy. Thus, we conclude that expansion of
such software usage will help to improve the primary data analysis efficiency
for next generation sequencing.


  Heterologous gene expression draws resources from host cells. These resources
include vital components to sustain growth and replication, and the resulting
cellular burden is a widely recognised bottleneck in the design of robust
circuits. In this tutorial we discuss the use of computational models that
integrate gene circuits and the physiology of host cells. Through various use
cases, we illustrate the power of host-circuit models to predict the impact of
design parameters on both burden and circuit functionality. Our approach relies
on a new generation of computational models for microbial growth that can
flexibly accommodate resource bottlenecks encountered in gene circuit design.
Adoption of this modelling paradigm can facilitate fast and robust design
cycles in synthetic biology.


  This chapter gives a synopsis of recent approaches to model and analyse the
evolution of microbial populations under selection. The first part reviews two
population genetic models of Lenski's long-term evolution experiment with
Escherichia coli, where models aim at explaining the observed curve of the
evolution of the mean fitness. The second part describes a model of a
host-pathogen system where the population of pathogenes experiences balancing
selection, migration, and mutation, as motivated by observations of the genetic
diversity of HCMV (the human cytomegalovirus) across hosts.


  As COVID-19 is rapidly spreading across the globe, short-term modeling
forecasts provide time-critical information for decisions on containment and
mitigation strategies. A main challenge for short-term forecasts is the
assessment of key epidemiological parameters and how they change when first
interventions show an effect. By combining an established epidemiological model
with Bayesian inference, we analyze the time dependence of the effective growth
rate of new infections. Focusing on the COVID-19 spread in Germany, we detect
change points in the effective growth rate that correlate well with the times
of publicly announced interventions. Thereby, we can quantify the effect of
interventions, and we can incorporate the corresponding change points into
forecasts of future scenarios and case numbers. Our code is freely available
and can be readily adapted to any country or region.


  The life of a cell is governed by highly dynamical microscopic processes. Two
notable examples are the diffusion of membrane receptors and the kinetics of
transcription factors governing the rates of gene expression. Different
fluorescence imaging techniques have emerged to study molecular dynamics. Among
them, fluorescence correlation spectroscopy (FCS) and single-particle tracking
(SPT) have proven to be instrumental to our understanding of cell dynamics and
function. The analysis of SPT and FCS is an ongoing effort, and despite decades
of work, much progress remains to be done. In this paper, we give a quick
overview of the existing techniques used to analyze anomalous diffusion in
cells and propose a collaborative challenge to foster the development of
state-of-the-art analysis algorithms. We propose to provide labelled (training)
and unlabelled (evaluation) simulated data to competitors all over the world in
an open and fair challenge. The goal is to offer unified data benchmarks based
on biologically-relevant metrics in order to compare the diffusion analysis
software available for the community.


  Protein folding is a central challenge in computational biology, with
important applications in molecular biology, drug discovery and catalyst
design. As a hard combinatorial optimisation problem, it has been studied as a
potential target problem for quantum annealing. Although several experimental
implementations have been discussed in the literature, the computational
scaling of these approaches has not been elucidated. In this article, we
present a numerical study of quantum annealing applied to a large number of
small peptide folding problems, aiming to infer useful insights for near-term
applications. We present two conclusions: that even naive quantum annealing,
when applied to protein lattice folding, has the potential to outperform
classical approaches, and that careful engineering of the Hamiltonians and
schedules involved can deliver notable relative improvements for this problem.
Overall, our results suggest that quantum algorithms may well offer
improvements for problems in the protein folding and structure prediction
realm.


  We study the spread of information in finite and infinite inhomogeneous
spatial random graphs. We assume that each edge has a transmission cost that is
a product of an i.i.d. random variable L and a penalty factor: edges between
vertices of expected degrees w_1 and w_2 are penalised by a factor of
(w_1w_2)^\mu for all \mu >0. We study this process for scale-free percolation,
for (finite and infinite) Geometric Inhomogeneous Random Graphs, and for
Hyperbolic Random Graphs, all with power law degree distributions with exponent
\tau > 1. For \tau < 3, we find a threshold behaviour, depending on how fast
the cumulative distribution function of L decays at zero. If it decays at most
polynomially with exponent smaller than (3-\tau)/(2\mu) then explosion happens,
i.e., with positive probability we can reach infinitely many vertices with
finite cost (for the infinite models), or reach a linear fraction of all
vertices with bounded costs (for the finite models). On the other hand, if the
cdf of L decays at zero at least polynomially with exponent larger than
(3-\tau)/(2\mu), then no explosion happens. This behaviour is arguably a better
representation of information spreading processes in social networks than the
case without penalising factor, in which explosion always happens unless the
cdf of L is doubly exponentially flat around zero. Finally, we extend the
results to other penalty functions, including arbitrary polynomials in w_1 and
w_2. In some cases the interesting phenomenon occurs that the model changes
behaviour (from explosive to conservative and vice versa) when we reverse the
role of w_1 and w_2. Intuitively, this could corresponds to reversing the flow
of information: gathering information might take much longer than sending it
out.


  In this work, a proposal for the estimation of the populations using logistic
curve fitting is presented. This type of curve is used to study population
growth, in this case population of people infected with the Covid-19 virus; and
it can also be used to approximate the survival curve used in actuarial and
similar studies in Spanish: En este trabajos se presenta una propuesta para la
estimaci\'on de la poblaciones usando ajuste de curvas del tipo log\'istica.
Este tipo de curvas se utilizan para el estudio de crecimiento de poblaciones,
en este casos poblaci\'on de personas infectadas por el virus Covid-19; y
tambi\'en se puede utilizar para aproximar la curva de supervivencia que se
utiliza en estudios actuariales y otras similares


  In this study, we use US county-level COVID-19 case data from January
21-March 25, 2020 to study the exponential behavior of case growth at the
metropolitan scale. In particular, we assume that all localized outbreaks are
in an early stage (either undergoing exponential growth in the number of cases,
or are effectively contained) and compare the explanatory performance of
different simple exponential and linear growth models for different
metropolitan areas. While we find no relationship between city size and
exponential growth rate (directly related to $R0$, which denotes average the
number of cases an infected individual infects), we do find that larger cities
seem to begin exponential spreading earlier and are thus in a more advanced
stage of the pandemic at the time of submission. We also use more recent data
to compute prediction errors given our models, and find that in many cities,
exponential growth models trained on data before March 26 are poor predictors
for case numbers in this more recent period (March 26-30), likely indicating a
reduction in the number of new cases facilitated through social distancing.