  The cytoskeleton in eukaryotic cells plays several crucial roles. In terms of
intracellular transport, motor proteins use the cytoskeletal filaments as a
backbone along which they can actively transport biological cargos such as
vesicles carrying biochemical reactants. Crossings between such filaments
constitute a key element, as they may serve to alter the destination of such
payload. Although motor proteins are known to display a rich behaviour at such
crossings, the latter have so far only been modelled as simple branching
points. Here we explore a model for a crossing between two microtubules which
retains the individual tracks consisting of protofilaments, and we construct a
schematic representation of the transport paths. We study collective transport
exemplified by the Totally Asymmetric Simple Exclusion Process (TASEP), and
provide a full analysis of the transport features and the associated phase
diagram, by a generic mean-field approach which we confirm through
particle-based stochastic simulations. In particular we show that transport
through such a compound crossing cannot be approximated from a coarse-grained
structure with a simple branching point. Instead, it gives rise to entirely new
and counterintuitive features: the fundamental current-density relation for
traffic flow is no longer a single-valued function, and it furthermore differs
according to whether it is observed upstream or downstream from the crossing.
We argue that these novel features may be directly relevant for interpreting
experimental measurements.


  The novel coronavirus disease (COVID-19) pandemic has resulted in over
200,000 cases in India. Thus far, India has implemented lockdown measures to
curb disease transmission. However, commercial sex work in red-light areas
(RLAs) has potential to lead to COVID-19 resurgence after lockdown. We
developed a model of COVID-19 transmission in RLAs, evaluating the impact of
extended RLA closure compared with RLA reopening on cases, hospitalizations,
and mortality rates within the RLAs of five major Indian cities, within the
cities, and across India. Closure lowered transmission at all scales. More than
90% of cumulative cases and deaths among RLA residents of Kolkata, Pune, and
Nagpur could be averted by the time the epidemic would peak under a re-opening
scenario. Across India, extended closure of RLAs would benefit the population
at large, delaying the peak of COVID-19 cases by 8 to 23 days, and avert 32% to
60.2% of cumulative cases and 43% to 67.6% of cumulative deaths at the peak of
the epidemic. Extended closure of RLAs until better prevention and treatment
strategies are developed would benefit public health in India.


  We introduce an epidemic spreading model on a network using concepts from
percolation theory. The model is motivated by discussing the standard SIR
model, with extensions to describe effects of lockdowns within a population.
The underlying ideas and behavior of the lattice model, implemented using the
same lockdown scheme as for the SIR scheme, are discussed in detail and
illustrated with extensive simulations. A comparison between both models is
presented for the case of COVID-19 data from the USA. Both fits to the
empirical data are very good, but some differences emerge between the two
approaches which indicate the usefulness of having an alternative approach to
the widespread SIR model.


  Introduction: Endpoint choice for randomized controlled trials of treatments
for COVID-19 is complex. A new disease brings many uncertainties, but trials
must start rapidly. COVID-19 is heterogeneous, ranging from mild disease that
improves within days to critical disease that can last weeks and can end in
death. While improvement in mortality would provide unquestionable evidence
about clinical significance of a treatment, sample sizes for a study evaluating
mortality are large and may be impractical. Furthermore, patient states in
between "cure" and "death" represent meaningful distinctions. Clinical severity
scores have been proposed as an alternative. However, the appropriate summary
measure for severity scores has been the subject of debate, particularly in
relating to the uncertainty about the time-course of COVID-19. Outcomes
measured at fixed time-points may risk missing the time of clinical benefit. An
endpoint such as time-to-improvement (or recovery), avoids the timing problem.
However, some have argued that power losses will result from reducing the
ordinal scale to a binary state of "recovered" vs "not recovered."
  Methods: We evaluate statistical power for possible trial endpoints for
COVID-19 treatment trials using simulation models and data from two recent
COVID-19 treatment trials.
  Results: Power for fixed-time point methods depends heavily on the time
selected for evaluation. Time-to-improvement (or recovery) analyses do not
specify a time-point. Time-to-event approaches have reasonable statistical
power, even when compared to a fixed time-point method evaluated at the optimal
time.
  Discussion: Time-to-event analyses methods have advantages in the COVID-19
setting, unless the optimal time for evaluating treatment effect is known in
advance. Even when the optimal time is known, a time-to-event approach may
increase power for interim analyses.


  Under the global health emergency caused by coronavirus disease 2019
(COVID-19), efficient and specific therapies are urgently needed. Compared with
traditional small-molecular drugs, antibody therapies are relatively easy to
develop and as specific as vaccines in targeting severe acute respiratory
syndrome coronavirus 2 (SARS-CoV-2), and thus attract much attention in the
past few months. This work reviews seven existing antibodies for SARS-CoV-2
spike (S) protein with three-dimensional (3D) structures deposited in the
Protein Data Bank. Five antibody structures associated with SARS-CoV are
evaluated for their potential in neutralizing SARS-CoV-2. The interactions of
these antibodies with the S protein receptor-binding domain (RBD) are compared
with those of angiotensin-converting enzyme 2 (ACE2) and RBD complexes. Due to
the orders of magnitude in the discrepancies of experimental binding
affinities, we introduce topological data analysis (TDA), a variety of network
models, and deep learning to analyze the binding strength and therapeutic
potential of the aforementioned fourteen antibody-antigen complexes. The
current COVID-19 antibody clinical trials, which are not limited to the S
protein target, are also reviewed.


  In this paper, we study the effectiveness of the modelling approach on the
pandemic due to the spreading of the novel COVID-19 disease and develop a
susceptible-infected-removed (SIR) model that provides a theoretical framework
to investigate its spread within a community. Here, the model is based upon the
well-known susceptible-infected-removed (SIR) model with the difference that a
total population is not defined or kept constant per se and the number of
susceptible individuals does not decline monotonically. To the contrary, as we
show herein, it can be increased in surge periods! In particular, we
investigate the time evolution of different populations and monitor diverse
significant parameters for the spread of the disease in various communities,
represented by countries and the state of Texas in the USA. The SIR model can
provide us with insights and predictions of the spread of the virus in
communities that the recorded data alone cannot. Our work shows the importance
of modelling the spread of COVID-19 by the SIR model that we propose here, as
it can help to assess the impact of the disease by offering valuable
predictions. Our analysis takes into account data from January to June, 2020,
the period that contains the data before and during the implementation of
strict and control measures. We propose predictions on various parameters
related to the spread of COVID-19 and on the number of susceptible, infected
and removed populations until September 2020. By comparing the recorded data
with the data from our modelling approaches, we deduce that the spread of
COVID-19 can be under control in all communities considered, if proper
restrictions and strong policies are implemented to control the infection rates
early from the spread of the disease.


  The dissociation of ligands from proteins and other biomacromolecules occurs
over a wide range of timescales. For most pharmaceutically relevant inhibitors,
these timescales are far beyond those that are accessible by conventional
molecular dynamics (MD) simulation. Consequently, to explore ligand egress
mechanisms and compute dissociation rates, it is necessary to enhance the
sampling of ligand unbinding. Random Acceleration MD (RAMD) is a simple method
to enhance ligand egress from a macromolecular binding site, which enables the
exploration of ligand egress routes without prior knowledge of the reaction
coordinates. Furthermore, the tauRAMD procedure can be used to compute the
relative residence times of ligands. When combined with a machine-learning
analysis of protein-ligand interaction fingerprints (IFPs), molecular features
that affect ligand unbinding kinetics can be identified. Here, we describe the
implementation of RAMD in GROMACS 2020, which provides significantly improved
computational performance, with scaling to large molecular systems. For the
automated analysis of RAMD results, we developed MD-IFP, a set of tools for the
generation of IFPs along unbinding trajectories and for their use in the
exploration of ligand dynamics. We demonstrate that the analysis of ligand
dissociation trajectories by mapping them onto the IFP space enables the
characterization of ligand dissociation routes and metastable states. The
combined implementation of RAMD and MD-IFP provides a computationally efficient
and freely available workflow that can be applied to hundreds of compounds in a
reasonable computational time and will facilitate the use of tauRAMD in drug
design.


  Coronavirus disease (COVID-19) which is caused by SARS-COV2 has become a
pandemic. This disease is highly infectious and potentially fatal, causing a
global public health concern. To contain the spread of COVID-19, governments
are adopting nationwide interventions, like lockdown, containment and
quarantine, restrictions on travel, cancelling social events and extensive
testing. To understand the effects of these measures on the control of the
epidemic in a data-driven manner, we propose a probabilistic cellular automata
(PCA) based modified SEIQR model. The transitions associated with the model is
driven by data available on chronology, symptoms, pathogenesis and
transmissivity of the virus. By arguing that the lattice-based model captures
the features of the dynamics along with the existing fluctuations, we perform
rigorous computational analyses of the model to take into account of the
spatial dynamics of social distancing measures imposed on the people.
Considering the probabilistic behavioural aspects associated with mitigation
strategies, we study the model considering factors like population density and
testing efficiency. Using the model, we focus on the variability of epidemic
dynamics data for different countries and point out the reasons behind these
contrasting observations. To the best of our knowledge, this is the first
attempt to model COVID-19 spread using PCA that gives us both spatial and
temporal variations of the infection spread with the insight about the
contributions of different infection parameters.


  The COVID-19 pandemic started in Wuhan, China, and caused the worldwide
spread of the RNA virus SARS-CoV-2, the causative agent of COVID-19. Because of
its mutational rate, wide geographical distribution, and host response variance
this coronavirus is currently evolving into an array of strains with increasing
genetic diversity. Most variants apparently have neutral effects for disease
spread and symptoms severity. However, in the viral Spike protein, which is
responsible for host cell attachment and invasion, an emergent variant,
containing the amino acid substitution D to G in position 614 (D614G), was
suggested to increase viral infection capability. To test whether this variant
has epidemiological impact, the temporal distributions of the SARS-CoV-2
samples bearing D or G at position 614 were compared in the USA, Asia and
Europe. The epidemiological curves were compared at early and late epidemic
stages. At early stages, where containment measures were still not fully
implemented, the viral variants are supposed to be unconstrained and its growth
curves might approximate the free viral dynamics. Our analysis shows that the
D614G prevalence and the growth rates of COVID-19 epidemic curves are
correlated in the USA, Asia and Europe. Our results suggest a selective sweep
that can be explained, at least in part, by a propagation advantage of this
variant, in other words, that the molecular level effects of D614G have
sufficient impact on population transmission dynamics as to be detected by
differences in rate coefficients of epidemic growth curves.


  Mexico, like the rest of the world, is currently facing the The COVID-19
pandemic. Given the size of its territory, the efforts to contain the disease
have involved both national and regional measures. For this work, the curves of
daily new cases of each municipality reported by the federal government were
compared. We found that 114 municipalities form a large network of
statistically dependent epidemic phenomena. Based on the network's modular
structure, these 114 municipalities can be split into four distinct communities
of coordinated epidemic phenomena. These clusters are not limited by
geographical proximity. These findings can be helpful for public health
officials for the evaluation of past strategies and the development of new
directed interventions.


  The spatial variation in the physico-chemical variables and macrobenthic
invertebrates assemblage of a stretch of Ethiope River, Nigeria were evaluated
for 6 months at five study stations. The results showed that temperature,
dissolved oxygen and bicarbonates were similar across the study stations. The
highest values of conductivity and total dissolved solids were recorded at
station 1. The physico-chemical variables showed no significant different (P <
0.05) in the study stations. The abundance of the macrobenthic invertebrates
comprises 61 taxa in 34 families, and the highest and least values occurred at
stations 3 and 1 respectively. The Diptera which accounted for 52.77% of the
overall abundance was the dominant group encountered, while the Decapoda
(10.97%) which was prevalent in station 1 and the Ephemeroptera (8.05%) which
was common in stations 2, 3 and 4 constituted the subdominant group. Over 68%
of the oligochaetes recorded were restricted to station 3 except Dero limnosa
and Nais sp that occurred in all the study stations. The overall abundance of
the macrobenthic invertebrates was not significantly different (p < 0.05) among
the study stations. Based on the canonical correspondence analyses, total
dissolved solids, biochemical oxygen demand, bicarbonates and pH influenced the
macrobenthic invertebrate assemblage. The study stretch of Ethiope River was
reasonably free from gross pollution as indicated by the presence of
Ephemeroptera, Tricoptera, Odonata and Decapoda at the study sites. However,
there were indications of organic contamination as shown by the preponderance
of naidid oligochaetes and certain dipteran species including the chironomids.
This study revealed a high diversity of macrobenthic invertebrates, and this
was attributed to the relative spatial homogeneity of the physico-chemical
variables recorded at the study stations.


  Dry eye disease (DED) affects approximately half of the United States
population. DED is characterized by dryness on the corena surface due to a
variety of causes. This study fills the spatiotemporal gaps in DED epidemiology
by using Google Trends as a novel epidemiological tool for geographically
mapping DED in relation to environmental risk factors. We utilized Google
Trends to extract DED-related queries estimating user intent from 2004-2019 in
the United States. We incorporated national climate data to generate heat maps
comparing geographic, temporal, and environmental relationships of DED.
Multi-variable regression models were constructed to generate quadratic
forecasts predicting DED and control searches. Our results illustrated the
upward trend, seasonal pattern, environmental influence, and spatial
relationship of DED search volume across US geography. Localized patches of DED
interest were visualized along the coastline. There was no significant
difference in DED queries across US census regions. Regression model 1
predicted DED searches over time (R^2=0.97) with significant predictors being
control queries (p=0.0024), time (p=0.001), and seasonality (Winter p=0.0028;
Spring p<0.001; Summer p=0.018). Regression model 2 predicted DED queries per
state (R^2=0.49) with significant predictors being temperature (p=0.0003) and
coastal zone (p=0.025). Importantly, temperature, coastal status, and
seasonality were stronger risk factors of DED searches than humidity, sunshine,
pollution, or region as clinical literature may suggest. Our work paves the way
for future exploration of geographic information systems for locating DED and
other diseases via online search query metrics.


  In order to analyze the effectiveness of three successive nationwide lockdown
enforced in India, we present a data-driven analysis of four key parameters,
reducing the transmission rate, restraining the growth rate, flattening the
epidemic curve and improving the health care system. These were quantified by
the consideration of four different metrics, namely, reproduction rate, growth
rate, doubling time and death to recovery ratio. The incidence data of the
COVID-19 (during the period of 2nd March 2020 to 31st May 2020) outbreak in
India was analyzed for the best fit to the epidemic curve, making use of the
exponential growth, the maximum likelihood estimation, sequential Bayesian
method and estimation of time-dependent reproduction. The best fit (based on
the data considered) was for the time-dependent approach. Accordingly, this
approach was used to assess the impact on the effective reproduction rate. The
period of pre-lockdown to the end of lockdown 3, saw a $45\%$ reduction in the
rate of effective reproduction rate. During the same period the growth rate
reduced from $393\%$ during the pre-lockdown to $33\%$ after lockdown 3,
accompanied by the average doubling time increasing form $4$-$6$ days to
$12$-$14$ days. Finally, the death-to-recovery ratio dropped from $0.28$
(pre-lockdown) to $0.08$ after lockdown 3. In conclusion, all the four metrics
considered to assess the effectiveness of the lockdown, exhibited significant
favourable changes, from the pre-lockdown period to the end of lockdown 3.
Analysis of the data in the post-lockdown period with these metrics will
provide greater clarity with regards to the extent of the success of the
lockdown.


  Genetic sequences are known to possess non-trivial composition together with
symmetries in the frequencies of their components. Recently, it has been shown
that symmetry and structure are hierarchically intertwined in DNA, suggesting a
common origin for both features. However, the mechanism leading to this
relationship is unknown. Here we investigate a biologically motivated dynamics
for the evolution of genetic sequences. We show that a metastable (long-lived)
regime emerges in which sequences have symmetry and structure interlaced in a
way that matches that of extant genomes.


  This paper gives an introduction to rule-based modelling applied to topics in
infectious diseases. Rule-based models generalise reaction-based models with
reagents that have internal state and may be bound together to form complexes,
as in chemistry. Rule-based modelling is directly transferable from molecular
biology to epidemiology and allows us to express a broad class of models for
processes of interest in epidemiology that would not otherwise be feasible in
compartmental models. This includes dynamics commonly found in compartmental
models such as the spread of a virus from an infectious to a susceptible
population, and more complex dynamics outside the typical scope of such models
such as social behaviours and decision-making, testing capacity constraints,
and tracing of people exposed to a virus but not yet symptomatic. We propose
that such dynamics are well-captured with rule-based models, and that doing so
combines intuitiveness and transparency of representation with scalability and
compositionality. We demonstrate this feasibility of our approach using a suite
of seven models to describe a spread of infectious diseases under different
scenarios: wearing masks, infection via fomites and prevention by hand-washing,
the concept of vector-borne diseases, testing and contact tracing
interventions, disease propagation within motif-structured populations with
shared environments such as schools, and superspreading events. The
machine-readable description of these models corresponds closely to the
mathematical description and also functions as a human-readable format so that
one knows readily "what is in the model".


  Epidemic models play a key role in understanding and responding to the
emerging COVID-19 pandemic. Widely used compartmental models are static and are
of limited use to evaluate intervention strategies with the emerging pandemic.
Applying the technology of data assimilation, we propose a Bayesian updating
approach for estimating epidemiological parameters using observable information
for the purpose of assessing the impacts of different intervention strategies.
We adopt a concise renewal model and propose new parameters by disentangling
the reduction of instantaneous reproduction number Rt into mitigation and
suppression factors for quantifying intervention impacts at a finer
granularity. Then we developed a data assimilation framework for estimating
these parameters including constructing an observation function and developing
a Bayesian updating scheme. A statistical analysis framework is then built to
quantify the impact of intervention strategies by monitoring the evolution of
these estimated parameters. By Investigating the impacts of intervention
measures of European countries, the United States and Wuhan with the framework,
we reveal the effects of interventions in these countries and the resurgence
risk in the USA.


  In this paper we afford a quantitative analysis of the sustainability of
current world population growth in relation to the parallel deforestation
process adopting a statistical point of view. We consider a simplified model
based on a stochastic growth process driven by a continuous time random walk,
which depicts the technological evolution of human kind, in conjunction with a
deterministic generalised logistic model for humans-forest interaction and we
evaluate the probability of avoiding the self-destruction of our civilisation.
Based on the current resource consumption rates and best estimate of
technological rate growth our study shows that we have very low probability,
less than 10% in most optimistic estimate, to survive without facing a
catastrophic collapse.


  Recent analysis of early COVID-19 data from China showed that the number of
confirmed cases followed a subexponential power-law increase, with a growth
exponent of around 2.2 [B.\,F.~Maier, D.~Brockmann, {\it Science} {\bf 368},
742 (2020)]. The power-law behavior was attributed to a combination of
effective containment and mitigation measures employed as well as behavioral
changes by the population. In this work, we report a random walk Monte Carlo
simulation study of proximity-based infection spread. Control interventions
such as lockdown measures and mobility restrictions are incorporated in the
simulations through a single parameter, the size of each step in the random
walk process. The step size $l$ is taken to be a multiple of $\langle r
\rangle$, which is the average separation between individuals. Three temporal
growth regimes (quadratic, intermediate power-law and exponential) are shown to
emerge naturally from our simulations. For $l = \langle r \rangle$, we get
intermediate power-law growth exponents that are in general agreement with
available data from China. On the other hand, we obtain a quadratic growth for
smaller step sizes $l \lesssim \langle r \rangle/2 $, while for large $l$ the
growth is found to be exponential. %Together with available data, these results
suggest that the early containment of the disease within China was close to
optimal. We further performed a comparative case study of early fatality data
(under varying levels of lockdown conditions) from three other countries,
India, Brazil and South Africa. We show that reasonable agreement with these
data can be obtained by incorporating small-world-like connections in our
simulations.


  In this paper, we estimate the impact of national lockdown on COVID-19
related total and daily deaths, per million people, in select European
countries. In particular, we compare countries that imposed a nationwide
lockdown (Treatment group); Belgium, Denmark, France, Germany, Italy, Norway,
Spain, United Kingdom (UK), and the US, to Sweden (Control group) that did not
impose national lockdown using a changes-in-changes (CIC) estimation model. The
key advantage of the CIC model as compared to the standard
difference-in-difference model is that CIC allows for mean and variance of the
outcomes to change over time in the absence of any policy intervention, and CIC
accounts for endogeneity in the choice of policy intervention. Our results
indicate that in contrast to Sweden, which did not impose a national lockdown,
Germany, and to some extent, the US were the two countries where nationwide
lockdown had a significant impact on the reduction in COVID-19 related total
and daily deaths per million people. In Norway and Denmark, there was no
significant impact on total and daily deaths per million people relative to
Sweden. Whereas in other countries; Belgium, France, Italy, Spain, and the UK,
the effect of the lockdown was in the opposite direction, that is, they
experienced significantly higher COVID-19 related total and daily deaths per
million people, post the lockdown as compared to Sweden. Our results suggest
that the impact of nationwide lockdown on COVID-19 related total and daily
deaths per million people varied from one country to another.


  Binding of anthracycline antibiotic carminomycin (CM) to synthetic
polyribonucleotides poly(A) and poly(U) was studied in solution of low ionic
strength in a wide phosphate-to-dye (P/D) range using absorption and polarized
fluorescence spectroscopy. Two different modes of CM binding to the
ss-polynucleotides have been identified. The first of them dominating at low
phosphate-to-dye (P/D) ratios is self-assembly of the heterocyclic dye on the
polymer surface driven by to cooperative electrostatic binding of amino group
of CM sugar moiety to negatively charged polynucleotide phosphate groups with
the chromophores self-stacking. At high P/D values, the stacking-associates
disintegrate and monomeric binding of ligand to nucleic bases become prevalent.
Thermodynamic parameters of binding were estimated for the both cases.


  We consider the SIR epidemiological model applied to the evolution of
COVID-19 with two approaches. In the first place we fit a global SIR model,
with time delay, and constant parameters throughout the outbreak, including the
contagion rate. The contention measures are reflected on an effective reduced
susceptible population $N_{eff}$. In the second approach we consider a
time-dependent contagion rate that reflects the contention measures either
through a step by step fitting process or by following an exponential decay. In
this last model the population is considered the one of the country $N$. In the
linear region of the differential equations, when the total population $N$ is
large the predictions are independent of $N$. We apply these methodologies to
study the spread of the pandemic in Argentina, Brazil, Colombia, Mexico, and
South Africa for which the infection peaks are yet to be reached. In all of
these cases we provide estimates for the reproduction and recovery rates. The
scenario for a time varying contagion rate is optimistic, considering that
reasonable measures are taken such that the reproduction factor $R_0$ decreases
exponentially. The measured values for the recovery rate $\gamma$ are reported
finding a universality of this parameter over various countries. We discuss the
correspondence between the global SIR with effective population $N_{eff}$ and
the evolution of the time local SIR.


  A deterministic multi-stage malaria model with a non-therapeutic control
measure, the use of mosquito bednet is formulated and analyzed. The model basic
reproduction number is derived, and analytical results show that the models
equilibria are locally and globally asymptotically stable when certain
threshold conditions are satisfied. Pontryagin's Maximum Principle with respect
to a time dependent constant is used to derive the necessary conditions for the
optimal usage of the Long-Lasting Insecticide-treated bednets(LLINs) to
mitigate the malaria transmission dynamics. This is accomplished by introducing
biologically admissible control and e-approximate sub-optimal control. The
results from this study could help public health planners and policy
decision-makers to design reachable and more practical malaria prevention
programs "close" to the optimal strategy.


  Accurately forecasting county level COVID-19 confirmed cases is crucial to
optimizing medical resources. Forecasting emerging outbreaks pose a particular
challenge because many existing forecasting techniques learn from historical
seasons trends. Recurrent neural networks (RNNs) with LSTM-based cells are a
logical choice of model due to their ability to learn temporal dynamics. In
this paper, we adapt the state and county level influenza model, TDEFSI-LONLY,
proposed in Wang et a. [l2020] to national and county level COVID-19 data. We
show that this model poorly forecasts the current pandemic. We analyze the two
week ahead forecasting capabilities of the TDEFSI-LONLY model with combinations
of regularization techniques. Effective training of the TDEFSI-LONLY model
requires data augmentation, to overcome this challenge we utilize an SEIR model
and present an inter-county mixing extension to this model to simulate
sufficient training data. Further, we propose an alternate forecast model, {\it
County Level Epidemiological Inference Recurrent Network} (\alg{}) that trains
an LSTM backbone on national confirmed cases to learn a low dimensional time
pattern and utilizes a time distributed dense layer to learn individual county
confirmed case changes each day for a two weeks forecast. We show that the
best, worst, and median state forecasts made using CLEIR-Net model are
respectively New York, South Carolina, and Montana.


  The paper formulates and solves the problem of identification of unknown
parameters of mathematical models of the spread of COVID-19 coronavirus
infection, based on SEIR type models, based on additional information about the
number of detected cases, mortality, self-isolation coefficient and tests
performed for the Moscow city and the Novosibirsk Region from 03.23.2020.
Within the framework of the models used, the population is divided into seven
(SEIR-HCD) and five (SEIR-D) groups with similar characteristics with
transition probabilities between groups depending on a specific region.
Identifiability analysis of the SEIR-HCD mathematical model was carried out,
which revealed the least sensitive unknown parameters to additional
measurements. The tasks of refining the parameters are reduced to minimizing
the corresponding target functionals, which were solved using stochastic
methods (simulating annealing, differential evolution, genetic algorithm,
etc.). For a different amount of tested data, a prognostic scenario for the
development of the disease in the city of Moscow and the Novosibirsk region was
developed, the peak is predicted the development of the epidemic in Moscow with
an error of 2 days and 174 detected cases, and an analysis of the applicability
of the developed models was carried out.


  Polyampholyte field theory and explicit-chain molecular dynamics models of
sequence-specific phase separation of a system with two intrinsically
disordered protein (IDP) species indicate consistently that a substantial
polymer excluded volume and a significant mismatch of the IDP sequence charge
patterns can act in concert, but not in isolation, to demix the two IDP species
upon condensation. This finding reveals an energetic-geometric interplay in a
stochastic, "fuzzy" molecular recognition mechanism that may facilitate
subcompartmentalization of membraneless organelles.


  Analysing and understanding the transmission and evolution of the COVID-19
pandemic is mandatory to be able to design the best social and medical
policies, foresee their outcomes and deal with all the subsequent
socio-economic effects. We address this important problem from a computational
and machine learning perspective. More specifically, we want to statistically
estimate all the relevant parameters for the new coronavirus COVID-19, such as
the reproduction number, fatality rate or length of infectiousness period,
based on Romanian patients, as well as be able to predict future outcomes. This
endeavor is important, since it is well known that these factors vary across
the globe, and might be dependent on many causes, including social, medical,
age and genetic factors. We use a recently published improved version of SEIR,
which is the classic, established model for infectious diseases. We want to
infer all the parameters of the model, which govern the evolution of the
pandemic in Romania, based on the only reliable, true measurement, which is the
number of deaths. Once the model parameters are estimated, we are able to
predict all the other relevant measures, such as the number of exposed and
infectious people. To this end, we propose a self-supervised approach to train
a deep convolutional network to guess the correct set of Modified-SEIR model
parameters, given the observed number of daily fatalities. Then, we refine the
solution with a stochastic coordinate descent approach. We compare our deep
learning optimization scheme with the classic grid search approach and show
great improvement in both computational time and prediction accuracy. We find
an optimistic result in the case fatality rate for Romania which may be around
0.3% and we also demonstrate that our model is able to correctly predict the
number of daily fatalities for up to three weeks in the future.


  A molecular and cellular understanding of how SARS-CoV-2 variably infects and
causes severe COVID-19 remains a bottleneck in developing interventions to end
the pandemic. We sought to use deep learning to study the biology of SARS-CoV-2
infection and COVID-19 severity by identifying transcriptomic patterns and cell
types associated with SARS-CoV-2 infection and COVID-19 severity. To do this,
we developed a new approach to generating self-supervised edge features. We
propose a model that builds on Graph Attention Networks (GAT), creates edge
features using self-supervised learning, and ingests these edge features via a
Set Transformer. This model achieves significant improvements in predicting the
disease state of individual cells, given their transcriptome. We apply our
model to single-cell RNA sequencing datasets of SARS-CoV-2 infected lung
organoids and bronchoalveolar lavage fluid samples of patients with COVID-19,
achieving state-of-the-art performance on both datasets with our model. We then
borrow from the field of explainable AI (XAI) to identify the features (genes)
and cell types that discriminate bystander vs. infected cells across time and
moderate vs. severe COVID-19 disease. To the best of our knowledge, this
represents the first application of deep learning to identifying the molecular
and cellular determinants of SARS-CoV-2 infection and COVID-19 severity using
single-cell omics data.


  Combinations of intense non-pharmaceutical interventions ('lockdowns') were
introduced in countries worldwide to reduce SARS-CoV-2 transmission. Many
governments have begun to implement lockdown exit strategies that allow
restrictions to be relaxed while attempting to control the risk of a surge in
cases. Mathematical modelling has played a central role in guiding
interventions, but the challenge of designing optimal exit strategies in the
face of ongoing transmission is unprecedented. Here, we report discussions from
the Isaac Newton Institute 'Models for an exit strategy' workshop (11-15 May
2020). A diverse community of modellers who are providing evidence to
governments worldwide were asked to identify the main questions that, if
answered, will allow for more accurate predictions of the effects of different
exit strategies. Based on these questions, we propose a roadmap to facilitate
the development of reliable models to guide exit strategies. The roadmap
requires a global collaborative effort from the scientific community and
policy-makers, and is made up of three parts: i) improve estimation of key
epidemiological parameters; ii) understand sources of heterogeneity in
populations; iii) focus on requirements for data collection, particularly in
Low-to-Middle-Income countries. This will provide important information for
planning exit strategies that balance socio-economic benefits with public
health.


  Metastasis, the spread of cancer cells from a primary tumor to secondary
location(s) in the human organism, is the ultimate cause of death for the
majority of cancer patients. That is why, it is crucial to understand
metastases evolution in order to successfully combat the disease. We consider a
metastasized cancer cell population after medical treatment (e.g.
chemotherapy). Arriving in a different environment the cancer cells may change
their lifespan and reproduction, thus they may proliferate into different
types. If the treatment is effective, in the context of branching processes it
means, the reproduction of cancer cells is such that the mean offspring of each
cell is less than one. However, it is possible mutations to occur during cell
division cycle. These mutations can produce a new cancer cell type, which is
resistant to the treatment. Cancer cells from this new type may lead to the
rise of a non-extinction branching process. The above scenario leads us to the
choice of a reducible multi-type age-dependent branching process as a relevant
framework for studying the asymptotic behavior of such complex structures. Our
previous theoretical results are related to the asymptotic behavior of the
waiting time until the first occurrence of a mutant starting a non-extinction
process and the modified hazard function as a measure of immediate recurrence
of cancer disease. In the present paper these asymptotic results are used for
developing numerical schemes and algorithms implemented in Python via the NumPy
package for approximate calculation of the corresponding quantities. In
conclusion, our conjecture is that this methodology can be advantageous in
revealing the role of the lifespan distribution of the cancer cells in the
context of cancer disease evolution and other complex cell population systems,
in general.


  We show that logic computational circuits in gene regulatory networks arise
from a fibration symmetry breaking in the network structure. From this idea we
implement a constructive procedure that reveals a hierarchy of genetic
circuits, ubiquitous across species, that are surprising analogues to the
emblematic circuits of solid-state electronics: starting from the transistor
and progressing to ring oscillators, current-mirror circuits to toggle switches
and flip-flops. These canonical variants serve fundamental operations of
synchronization and clocks (in their symmetric states) and memory storage (in
their broken symmetry states). These conclusions introduce a theoretically
principled strategy to search for computational building blocks in biological
networks, and present a systematic route to design synthetic biological
circuits.


  In this paper we consider the epidemic competition between two generic
diffusion processes, where each competing side is represented by a different
state of a stochastic process. For this setting, we present the Generalized
Largest Reduction in Infectious Edges (gLRIE) dynamic resource allocation
strategy to advantage the preferred state against the other. Motivated by
social epidemics, we apply this method to a generic continuous-time SIS-like
diffusion model where we allow for: i) arbitrary node transition rate functions
that describe the dynamics of propagation depending on the network state, and
ii) competition between the healthy (positive) and infected (negative) states,
which are both diffusive at the same time, yet mutually exclusive on each node.
Finally we use simulations to compare empirically the proposed gLRIE against
competitive approaches from literature.


  Immune repertoires provide a unique fingerprint reflecting the immune history
of individuals, with potential applications in precision medicine. However, the
question of how personal that information is and how it can be used to identify
individuals has not been explored. Here, we show that individuals can be
uniquely identified from repertoires of just a few thousands lymphocytes. We
present "Immprint," a classifier using an information-theoretic measure of
repertoire similarity to distinguish pairs of repertoire samples coming from
the same versus different individuals. Using published T-cell receptor
repertoires and statistical modeling, we tested its ability to identify
individuals with great accuracy, including identical twins, by computing false
positive and false negative rates $< 10^{-6}$ from samples composed of 10,000
T-cells. We verified through longitudinal datasets and simulations that the
method is robust to acute infections and the passage of time. These results
emphasize the private and personal nature of repertoire data.


  The Covid-19 pandemic has made clear the need to improve modern multivariate
time-series forecasting models. Current state of the art predictions of future
daily deaths and, especially, hospital resource usage have confidence intervals
that are unacceptably wide. Policy makers and hospitals require accurate
forecasts to make informed decisions on passing legislation and allocating
resources. We used US county-level data on daily deaths and population
statistics to forecast future deaths. We extended the SIR epidemiological model
to a novel model we call the PECAIQR model. It adds several new variables and
parameters to the naive SIR model by taking into account the ramifications of
the partial quarantining implemented in the US. We fitted data to the model
parameters with numerical integration. Because of the fit degeneracy in
parameter space and non-constant nature of the parameters, we developed several
methods to optimize our fit, such as training on the data tail and training on
specific policy regimes. We use cross-validation to tune our hyper parameters
at the county level and generate a CDF for future daily deaths. For predictions
made from training data up to May 25th, we consistently obtained an averaged
pinball loss score of 0.096 on a 14 day forecast. We finally present examples
of possible avenues for utility from our model. We generate longer-time horizon
predictions over various 1-month windows in the past, forecast how many medical
resources such as ventilators and ICU beds will be needed in counties, and
evaluate the efficacy of our model in other countries.


  We discuss the generation of various reproduction ratios or numbers that are
very useful to monitor an ongoing epidemic like Covid-19 and examine the
effects of intervention measures. A detailed SEIR algorithm is described for
their computation, with applications given to the current Covid-19 outbreaks in
a number of countries (Argentina, Brazil, France, Italy, Mexico, Spain, UK and
USA). The corresponding matlab script, complete and ready to use, is provided
for free downloading.


  The COVID-19 pandemic is the profoundest health crisis of the 21rst century.
The SARS-CoV-2 virus arrived in Brazil around March, 2020 and its social and
economical backlashes are catastrophic. In this paper, it is investigated how
Model Predictive Control (MPC) could be used to plan appropriate social
distancing policies to mitigate the pandemic effects in Bahia and Santa
Catarina, two states of different regions, culture, and population demography
in Brazil. In addition, the parameters of
Susceptible-Infected-Recovered-Deceased (SIRD) models for these two states are
identified using an optimization procedure. The control input to the process is
a social isolation guideline passed to the population. Two MPC strategies are
designed: a) a centralized MPC, which coordinates a single control policy for
both states; and b) a decentralized strategy, for which one optimization is
solved for each state. Simulation results are shown to illustrate and compare
both control strategies. The framework serves as guidelines to deals with such
pandemic phenomena.


  We investigate the effect of school closure and subsequent reopening on the
transmission of COVID-19, by considering Denmark, Norway, Sweden, and German
states as case studies. By comparing the growth rates in daily hospitalisations
or confirmed cases under different interventions, we provide evidence that the
effect of school closure is visible as a reduction in the growth rate
approximately 9 days after implementation. Limited school attendance, such as
older students sitting exams or the partial return of younger year groups, does
not appear to significantly affect community transmission. A large-scale
reopening of schools while controlling or suppressing the epidemic appears
feasible in countries such as Denmark or Norway, where community transmission
is generally low. However, school reopening can contribute to significant
increases in the growth rate in countries like Germany, where community
transmission is relatively high. Our findings underscore the need for a
cautious evaluation of reopening strategies that ensure low classroom occupancy
and a solid infrastructure to quickly identify and isolate new infections.


  Many quantities we are interested in predicting are geometric tensors; we
refer to this class of problems as geometric prediction. Attempts to perform
geometric prediction in real-world scenarios have been limited to approximating
them through scalar predictions, leading to losses in data efficiency. In this
work, we demonstrate that equivariant networks have the capability to predict
real-world geometric tensors without the need for such approximations. We show
the applicability of this method to the prediction of force fields and then
propose a novel formulation of an important task, biomolecular structure
refinement, as a geometric prediction problem, improving state-of-the-art
structural candidates. In both settings, we find that our equivariant network
is able to generalize to unseen systems, despite having been trained on small
sets of examples. This novel and data-efficient ability to predict real-world
geometric tensors opens the door to addressing many problems through the lens
of geometric prediction, in areas such as 3D vision, robotics, and molecular
and structural biology.


  Currently, many countries are considering the introduction of tracing
software on mobile smartphones with the main purpose to inform and alarm the
mobile app user. Here, we demonstrate that, in addition to alarming and
informing, mobile tracing can detect nearly all users that are infected by
SARS-CoV-2. Our algorithm BETIS (Bayesian Estimation for Tracing Infection
States) makes use of self-reports of the user's health status. Then, BETIS
guarantees that almost all SARS-CoV-2 infections of the group of users can be
detected. Furthermore, BETIS estimates the virus prevalence in the whole
population, consisting of users and non-users. BETIS is based on a hidden
Markov epidemic model and recursive Bayesian filtering. The potential that
mobile tracing apps, in addition to medical testing and quarantining, can
eradicate COVID-19 may persuade citizens to trade-off privacy against public
health.


  The recent coronavirus disease (COVID-19) outbreak has dramatically increased
the public awareness and appreciation of the utility of dynamic models. At the
same time, the dissemination of contradictory model predictions has highlighted
their limitations. If some parameters and/or state variables of a model cannot
be determined from output measurements, its ability to yield correct insights
-- as well as the possibility of controlling the system -- may be compromised.
Epidemic dynamics are commonly analysed using compartmental models, and many
variations of such models have been used for analysing and predicting the
evolution of the COVID-19 pandemic. In this paper we survey the different
models proposed in the literature, assembling a list of 36 model structures and
assessing their ability to provide reliable information. We address the problem
using the control theoretic concepts of structural identifiability and
observability. Since some parameters can vary during the course of an epidemic,
we consider both the constant and time-varying parameter assumptions. We
analyse the structural identifiability and observability of all of the models,
considering all plausible choices of outputs and time-varying parameters, which
leads us to analyse 255 different model versions. We classify the models
according to their structural identifiability and observability under the
different assumptions and discuss the implications of the results. We also
illustrate with an example several alternative ways of remedying the lack of
observability of a model. Our analyses provide guidelines for choosing the most
informative model for each purpose, taking into account the available knowledge
and measurements.


  Management performance of five alternative stock assessment methods was
evaluated by using them to set harvest levels targeting multi-species maximum
yield in a multi-species flatfish fishery, including single-species and
hierarchical multi-species models, and methods that pooled data across species
and spatial strata, with catch outcomes of each method under three data
scenarios compared to catch under an omniscient manager simulation. Operating
models included technical interactions between species intended to produce
choke effects often observed in output controlled multi-species fisheries.
Hierarchical multi-species models outperformed all other methods under
data-poor and data-moderate scenarios, and outperformed single-species models
under the data-rich scenario. Hierarchical models were least sensitive to prior
precision, sometimes improving in performance when prior precision was reduced.
Choke effects were found to both positive and negative effects, sometimes
leading to underfishing of non-choke species, but at other times preventing
overfishing of non-choke species. We highlight the importance of including
technical interactions in multi-species assessment models and management
objectives, how choke species can indicate mismatches between management
objectives and system dynamics, and recommend hierarchical multi-species models
for multi-species fishery management systems.


  We are currently facing a highly critical case of a world-wide pandemic. The
novel coronavirus (SARS-CoV-2, a.k.a. COVID-19) has proved to be extremely
contagious and the original outbreak from Asia has now spread to all
continents. This situation will fruitfully profit from the study in regards of
the spread of the virus, assessing effective countermeasures to weight the
impact of the adopted strategies. The standard Susceptible-Infectious-Recovered
(SIR) model is a very successful and widely used mathematical model for
predicting the spread of an epidemic. We adopt the SIR model on a random
network and extend the model to include control strategies {\em lockdown} and
{\em testing} -- two often employed mitigation strategies. The ability of these
strategies in controlling the pandemic spread is investigated by varying the
effectiveness with which they are implemented. The possibility of a second
outbreak is evaluated in detail after the mitigation strategies are withdrawn.
We notice that, in any case, a sudden interruption of such mitigation
strategies will likely induce a resurgence of a second outbreak, whose peak
will be correlated to the number of susceptible individuals. In fact, we find
that a population will remain vulnerable to the infection until the herd
immunity is achieved. We also test our model with real statistics and
information on the epidemic spread in South Korea, Germany, and New York and
find a remarkable agreement with the simulation data.


  The present report briefly reviews the basic features of the current strategy
for the protection against the health harms caused by solar ultraviolet (UV,
295 - 400 nm). The emphasis has been made upon the erythema as being the
best-studied UV harm and the ability of UV irradiance to damage the
deoxyribonucleic acid (DNA) molecules, which leads to carcinogenesis. The
erythemally weighted UV irradiance that determines the ultraviolet index (UVI),
represents a common measure of the solar UV radiation level at the Earth's
surface and the current protective messages have been made by using UVI as a
basic parameter. However, such an approach seems insufficiently grounded to be
used also in the messages against the skin cancer, bearing in mind the
different nature of the erythema and DNA lesions. In this context, an
enlargement of the strategy basis by including additional biological effects of
UV radiation studied during the past years has been discussed. For instance,
the weight of the spectral UV-A (315 - 400 nm) band that in practice had been
neglected by UVI definition can be enhanced since it was found to play an
important role in DNA damaging. In addition, features of the contemporaneous
life style can be taken into account together with some people habits.


  The ongoing novel coronavirus epidemic has been announced a pandemic by the
World Health Organization on March 11, 2020, and the Govt. of India has
declared a nationwide lockdown from March 25, 2020, to prevent community
transmission of COVID-19. Due to absence of specific antivirals or vaccine,
mathematical modeling play an important role to better understand the disease
dynamics and designing strategies to control rapidly spreading infectious
diseases. In our study, we developed a new compartmental model that explains
the transmission dynamics of COVID-19. We calibrated our proposed model with
daily COVID-19 data for the four Indian provinces, namely Jharkhand, Gujarat,
Andhra Pradesh, and Chandigarh. We study the qualitative properties of the
model including feasible equilibria and their stability with respect to the
basic reproduction number $\mathcal{R}_0$. The disease-free equilibrium becomes
stable and the endemic equilibrium becomes unstable when the recovery rate of
infected individuals increased but if the disease transmission rate remains
higher then the endemic equilibrium always remain stable. For the estimated
model parameters, $\mathcal{R}_0 >1$ for all the four provinces, which suggests
the significant outbreak of COVID-19. Short-time prediction shows the
increasing trend of daily and cumulative cases of COVID-19 for the four
provinces of India.


  The maintenance of cooperation in the presence of spatial restrictions has
been studied extensively. It is well-established that the underlying graph
topology can significantly influence the outcome of games on graphs.
Maintenance of cooperation could be difficult, especially in the absence of
spatial restrictions. The evolution of cooperation would naturally depend on
payoffs. However, payoffs are generally considered to be invariant in a given
game. A natural yet unexplored question is whether the topology of the
underlying structures on which the games are played, possesses no role
whatsoever in the determination of payoffs. Herein, we introduce the notion of
cooperator graphs and defector graphs as well as a new form of game payoff,
which is weakly dependent on the underlying network topology. These concepts
are inspired by the well-known microbial phenomenon of quorum sensing. We
demonstrate that even with such a weak dependence, the fundamental game
dynamics and indeed the very nature of the game may be altered. Such changes in
the nature of a game have been well-reported in theoretical and experimental
studies.


  The reproductive number R_0 (and its value after initial disease emergence R)
has long been used to predict the likelihood of pathogen invasion, to gauge the
potential severity of an epidemic, and to set policy around interventions.
However, often ignored complexities have generated confusion around use of the
metric. This is particularly apparent with the emergent pandemic virus
SARS-CoV-2, the causative agent of COVID-19. We address some of these
misconceptions, namely, how R changes over time, varies over space, and relates
to epidemic size by referencing the mathematical definition of R and examples
from the current pandemic. We hope that a better appreciation of the uses,
nuances, and limitations of R facilitates a better understanding of epidemic
spread, epidemic severity, and the effects of interventions in the context of
SARS-CoV-2.


  COVID-19 is a pandemic disease that began to rapidly spread in the US with
the first case detected on January 19, 2020, in Washington State. March 9,
2020, and then increased rapidly with total cases of 25,739 as of April 20,
2020. The Covid-19 pandemic is so unnerving that it is difficult to understand
how any person is affected by the virus. Although most people with coronavirus
81%, according to the U.S. Centers for Disease Control and Prevention (CDC),
will have little to mild symptoms, others may rely on a ventilator to breathe
or not at all. SEIR models have broad applicability in predicting the outcome
of the population with a variety of diseases. However, many researchers use
these models without validating the necessary hypotheses. Far too many
researchers often "overfit" the data by using too many predictor variables and
small sample sizes to create models. Models thus developed are unlikely to
stand validity check on a separate group of population and regions. The
researcher remains unaware that overfitting has occurred, without attempting
such validation. In the paper, we present a combination algorithm that combines
similar days features selection based on the region using Xgboost, K Means, and
long short-term memory (LSTM) neural networks to construct a prediction model
(i.e., K-Means-LSTM) for short-term COVID-19 cases forecasting in Louisana
state USA. The weighted k-means algorithm based on extreme gradient boosting is
used to evaluate the similarity between the forecasts and past days. The
results show that the method with K-Means-LSTM has a higher accuracy with an
RMSE of 601.20 whereas the SEIR model with an RMSE of 3615.83.


  Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the strain of
coronavirus that causes coronavirus disease 2019 (COVID-19), the respiratory
illness responsible for the COVID-19 pandemic. Currently there is no known
vaccine or specific antiviral treatment for COVID-19 and so, there is an urgent
need for expedite discovery of new therapeutics to combat the disease until a
vaccine will be available worldwide. Drug repurposing is a strategy for
identifying new uses for approved drugs that has the advantage (over
conventional approaches that attempt to develop a drug from scratch) that time
frame of the overall process can be significantly reduced because of the few
number of clinical trial required. In this work, a virtual screening of
FDA-approved drugs was performed for repositioning as potential inhibitors of
the main protease Mpro of SARS-CoV-2. As a result of this study, 12 drugs are
proposed as candidates for inhibitors of the Mpro enzyme. Some of the selected
compounds are antiviral drugs that are already being tested in COVID-19
clinical trials (i.e. ribavirin) or are used to alleviate symptoms of the
disease (i.e. codeine). Surprisingly, the most promising candidate is the
naturally occurring broad spectrum antibiotic oxytetracycline. This compound
has largely outperformed the remaining selected candidates along all filtering
steps of our virtual screening protocol. If the activity of any of these drugs
is experimentally corroborated, they could be used directly in clinical trials
without the need for pre-clinical testing or safety evaluation since they are
already used as drugs for other diseases.


  We analyse the infection-age-dependent SIR model from a numerical point of
view. First, we present an algorithm for calculating the solution the
infection-age-structured SIR model without demography of the background host.
Second, we examine how and under which conditions, the conventional SIR model
(without infection-age) serves as a practical approximation to the
infection-age SIR model. Special emphasis is given on the effective
reproduction number.


  The water dynamics, as characterized by the local hydrophobicity (LH), is
investigated for tetrameric hemoglobin and dimeric melittin. For the T0 to R0
transition in Hb it is found that LH provides additional molecular-level
insight into the Perutz mechanism, i.e., the breaking and formation of salt
bridges at the alpha1 / beta2 and alpha2 / beta1 interface is accompanied by
changes in LH. For Hb in cubic water boxes with 90 Aengstroem and 120
Aengstroem edge length it is observed that following a decrease in LH as a
consequence of reduced water density or change of water orientation at the
protein/water interface the alpha / beta interfaces are destabilized; this is a
hallmark of the Perutz stereochemical model for the T to R transition in Hb.
The present work thus provides a dynamical view of the classical structural
model relevant to the molecular foundations of Hb function. For dimeric
melittin, earlier results by Cheng and Rossky (Nature, 1998, 392, 696-699) are
confirmed and interpreted on the basis of LH from simulations in which the
protein structure is frozen. For the flexible melittin dimer the changes in the
local hydration can be as much as 30 % than for the rigid dimer, reflecting the
fact that protein and water dynamics are coupled.


  One of the critical measures to control infectious diseases is a lockdown.
Once past the lockdown stage in many parts of the world, the crucial question
now concerns the effects of relaxing the lockdown and finding the best ways to
implement further lockdown(s), if required, to control the spread. With the
relaxation of lockdown, people migrate to different cities and enhance the
spread of the disease. This work presents the population migration model for
n-cities and applies the model for migration between two and three cities. The
reproduction number is calculated, and the effect of the migration rate is
analyzed. A punctuated lockdown is implemented to simulate a protocol of
repeated lockdowns that limits the resurgence of infections. A damped
oscillatory behavior is observed with multiple peaks over a period.


  Originating in Wuhan, the novel coronavirus, severe acute respiratory
syndrome 2 (SARS-CoV-2), has astonished health-care systems across globe due to
its rapid and simultaneous spread to the neighboring and distantly located
countries. To gain the systems level understanding of the role of global
transmission routes in the COVID-19 spread, in this study, we have developed
the first, empirical, global, index-case transmission network of SARS-CoV-2
termed as C19-TraNet. We manually curated the travel history of country wise
index-cases using government press releases, their official social media
handles and online news reports to construct this C19-TraNet that is a
spatio-temporal, sparse, growing network comprising of 187 nodes and 199 edges
and follows a power-law degree distribution. To model the growing C19-TraNet, a
novel stochastic scale free (SSF) algorithm is proposed that accounts for
stochastic addition of both nodes as well as edges at each time step. A
peculiar connectivity pattern in C19-TraNet is observed, characterized by a
fourth degree polynomial growth curve, that significantly diverges from the
average random connectivity pattern obtained from an ensemble of its 1,000 SSF
realizations. Partitioning the C19-TraNet, using edge betweenness, it is found
that most of the large communities are comprised of a heterogeneous mixture of
countries belonging to different world regions suggesting that there are no
spatial constraints on the spread of disease. This work characterizes the
superspreaders that have very quickly transported the virus, through multiple
transmission routes, to long range geographical locations alongwith their local
neighborhoods.


  Transformer architectures have proven to learn useful representations for
protein classification and generation tasks. However, these representations
present challenges in interpretability. In this work, we demonstrate a set of
methods for analyzing protein Transformer models through the lens of attention.
We show that attention: (1) captures the folding structure of proteins,
connecting amino acids that are far apart in the underlying sequence, but
spatially close in the three-dimensional structure, (2) targets binding sites,
a key functional component of proteins, and (3) focuses on progressively more
complex biophysical properties with increasing layer depth. We find this
behavior to be consistent across three Transformer architectures (BERT, ALBERT,
XLNet) and two distinct protein datasets. We also present a three-dimensional
visualization of the interaction between attention and protein structure. Code
for visualization and analysis is available at
https://github.com/salesforce/provis.


  The COVID-19 pandemic (SARS-CoV-2 virus) is the defying global health crisis
of our time. The absence of mass testing and the relevant presence of
asymptomatic individuals causes the available data of the COVID-19 pandemic in
Brazil to be largely under-reported regarding the number of infected
individuals and deaths. We propose an adapted Susceptible-Infected-Recovered
(SIR) model which explicitly incorporates the under-reporting and the response
of the population to public policies (such as confinement measures, widespread
use of masks, etc) to cast short-term and long-term predictions. Large amounts
of uncertainty could provide misleading models and predictions. In this paper,
we discuss the role of uncertainty in these prediction, which is illustrated
regarding three key aspects. First, assuming that the number of infected
individuals is under-reported, we demonstrate an anticipation regarding the
peak of infection. Furthermore, while a model with a single class of infected
individuals yields forecasts with increased peaks, a model that considers both
symptomatic and asymptomatic infected individuals suggests a decrease of the
peak of symptomatic. Second, considering that the actual amount of deaths is
larger than what is being register, then demonstrate the increase of the
mortality rates. Third, when consider generally under-reported data, we
demonstrate how the transmission and recovery rate model parameters change
qualitatively and quantitatively. We also investigate the effect of the
"COVID-19 under-reporting tripod", i.e. the under-reporting in terms of
infected individuals, of deaths and the true mortality rate. If two of these
factors are known, the remainder can be inferred, as long as proportions are
kept constant. The proposed approach allows one to determine the margins of
uncertainty by assessments on the observed and true mortality rates.


  A novel predictive modeling framework for the spread of infectious diseases
using high dimensional partial differential equations is developed and
implemented. A scalar function representing the infected population is defined
on a high-dimensional space and its evolution over all directions is described
by a population balance equation (PBE). New infections are introduced among the
susceptible population from non-quarantined infected population based on their
interaction, adherence to distancing norms, hygiene levels and any other
societal interventions. Moreover, recovery, death, immunity and all
aforementioned parameters are modeled on the high-dimensional space. To
epitomize the capabilities and features of the above framework, prognostic
estimates of Covid-19 spread using a six-dimensional (time, 2D space, infection
severity, duration of infection, and population age) PBE is presented. Further,
scenario analysis for different policy interventions and population behavior is
presented, throwing more insights into the spatio-temporal spread of infections
across disease age, intensity and age of population. These insights could be
used for science-informed policy planning.


  After the breakout of the disease caused by the new virus COVID-19, the
mitigation stage has been reached in most of the countries in the world. During
this stage, a more accurate data analysis of the daily reported cases and other
parameters became possible for the European countries and has been performed in
this work. Based on a proposed parametrization model appropriate for
implementation to an epidemic in a large population, we focused on the disease
spread and we studied the obtained curves, as well as, we investigated probable
correlations between the country's characteristics and the parameters of the
parametrization. We have also developed a methodology for coupling our model to
the SIR-based models determining the basic and the effective reproductive
number referring to the parameter space. The obtained results and conclusions
could be useful in the case of a recurrence of this repulsive disease in the
future.


  Background: Psoriasis is a multifactorial chronic inflammatory disorder of
the skin with significant morbidity, characterized by hyper proliferation of
the epidermis. Even though psoriasis etiology is not fully understood, it is
believed to be multifactorial with numerous key components. Methods: In order
to cast light on the complex molecular interactions in psoriasis vulgaris at
both protein-protein interactions and transcriptomics levels, we analyzed a set
of microarray gene expression analysis consisting of 170 paired lesional and
non-lesional samples. Afterwards, a network analysis was conducted on
protein-protein interaction network of differentially expressed genes based on
micro- and macro-level network metrics at a systemic level standpoint. Results:
We found 17 top communicative genes, all of which experimentally proven to be
pivotal in psoriasis were identified in two modules, namely, cell cycle and
immune system. Intra- and inter-gene interaction subnetworks from the top
communicative genes might provide further insight into the corresponding
characteristic mechanisms. Conclusions: Potential gene combinations for
therapeutic/diagnostics purposes were identified. Moreover, our proposed
pipeline could be of interest to a broader range of biological network analysis
studies.


  We develop an agent-based model to assess the cumulative number of deaths
during hypothetical Covid-19-like epidemics for various non-pharmaceutical
intervention strategies. We consider local and non-local modes of disease
transmission. The first simulates transmission through social contacts in the
vicinity of the place of residence while the second through social contacts in
public places: schools, hospitals, airports, etc., where many people meet, who
live in remote geographic locations. Epidemic spreading is modeled as a
discrete-time stochastic process on random geometric networks. We use the
Monte-Carlo method in the simulations. The~following assumptions are made. The
basic reproduction number is 2.5 and the infectious period lasts approximately
ten days. Infections lead to SARS in about one percent of cases, which are
likely to lead to respiratory default and death, unless the patient receives an
appropriate medical treatment. The~healthcare system capacity is simulated by
the availability of respiratory ventilators or intensive care beds. Some
parameters of the model, like mortality rates or the number of respiratory
ventilators per 100000 inhabitants, are chosen to simulate the real values for
the USA and Poland. In the simulations we compare `do-nothing' strategy with
mitigation strategies based on social distancing and reducing social mixing. We
study epidemics in the pre-vaccine era, where immunity is obtained only by
infection. The model applies only to epidemics for which reinfections are rare
and can be neglected. The results of the simulations show that strategies that
slow the development of an epidemic too much in the early stages do not
significantly reduce the overall number of deaths in the long term, but
increase the duration of the epidemic. In particular, a~hybrid strategy where
lockdown is held for some time and is then completely released, is inefficient.


  The effect of masking the general population on a COVID-19 epidemic is
estimated by computer simulation using two separate state-of-the-art web-based
softwares, one of them calibrated for the SARS-CoV-2 virus. The questions
addressed are these: 1. Can mask use by the general population limit the spread
of SARS-CoV-2 in a country? 2. What types of masks exist, and how elaborate
must a mask be to be effective against COVID-19? 3. Does the mask have to be
applied early in an epidemic? 4. A brief general discussion of masks and some
possible future research questions regarding masks and SARS-CoV-2. Results are
as follows: (1) The results indicate that any type of mask, even simple
home-made ones, may be effective. Masks use seems to have an effect in lowering
new patients even the protective effect of each mask (here dubbed "one-mask
protection") is low. Strict adherence to mask use does not appear to be
critical. However, increasing the one-mask protection to > 50% was found to be
advantageous. Masks seemed able to reduce overflow of capacity, e.g. of
intensive care. As the default parameters of the software included another
intervention, it seems possible to combine mask and other interventions. (2)
Masks do seem to reduce the number of new cases even if introduced at a late
stage in an epidemic. However, early implementation helps reduce the cumulative
and total number of cases. (3) The simulations suggest that it might be
possible to eliminate a COVID-19 outbreak by widespread mask use during a
limited period. The results from these simulations are encouraging, but do not
necessarily represent the real-life situation, so it is suggested that clinical
trials of masks are now carried out while continuously monitoring effects and
side-effects.


  This is a pedagogical paper on estimating the number of people that can be
infected by one infectious person during an epidemic outbreak, known as the
reproduction number. Knowing the number is crucial for developing policy
responses. There are generally two types of such a number, i.e., basic and
effective (or instantaneous). While basic reproduction number is the average
expected number of cases directly generated by one case in a population where
all individuals are susceptible, effective reproduction number is the number of
cases generated in the current state of a population. In this paper, we exploit
the deterministic susceptible-infected-removed (SIR) model to estimate them
through three different numerical approximations. We apply the methods to the
pandemic COVID-19 in Italy to provide insights into the spread of the disease
in the country. We see that the effect of the national lockdown in slowing down
the disease exponential growth appeared about two weeks after the
implementation date. We also discuss available improvements to the simple (and
naive) methods that have been made by researchers in the field.
  Authors of this paper are members of the SimcovID (Simulasi dan Pemodelan
COVID-19 Indonesia) collaboration.


  Without a vaccine, the fight against the spreading of the coronavirus has
focused on maintaining physical distance. To study the impact of such measures
on inter-municipality traffic, we analyze a mobile dataset with the daily flow
of people in Portugal in March and April 2020. We find that the reduction in
inter-municipality traffic depends strongly on its initial outflow. In
municipalities where the mobility is low, the outflow reduced by $10-20\%$ and
this decrease was independent of the population size. Whereas, for
municipalities of high mobility, the reduction was a monotonic increasing
function of the population size and it even exceeded $60\%$ for the largest
municipalities. As a consequence of such heterogeneities, there were
significant structural changes on the most probable paths for the spreading of
the virus, which must be considered when modeling the impact of control
measures.


  The current review focuses on important aspects of pyroptosis, such as
morphological features of inflammasomes, as well as on knowledge about
coronavirus infection mechanisms and the association between SARS-CoV-2
infection and cell pyroptosis. The application of immunomodulation therapies to
Covid-19 patients is the most promising treatment under investigation nowadays.
The possibilities of assistance to help protecting patients from infections are
also addressed in the present study.


  We apply the cavity master equation (CME) approach to epidemics models. We
explore mostly the susceptible-infectious-susceptible (SIS) model, which can be
readily treated with the CME as a two-state. We show that this approach is more
accurate than individual based and pair based mean field methods, and a
previously published dynamic message passing scheme. We explore average case
predictions and extend the cavity master equation to SIR and SIRS models.


  Forecasting models have been influential in shaping decision-making in the
COVID-19 pandemic. However, there is concern that their predictions may have
been misleading. Here, we dissect the predictions made by four models for the
daily COVID-19 death counts between March 25 and June 5 in New York state, as
well as the predictions of ICU bed utilisation made by the influential IHME
model. We evaluated the accuracy of the point estimates and the accuracy of the
uncertainty estimates of the model predictions. First, we compared the "ground
truth" data sources on daily deaths against which these models were trained.
Three different data sources were used by these models, and these had
substantial differences in recorded daily death counts. Two additional data
sources that we examined also provided different death counts per day. For
accuracy of prediction, all models fared very poorly. Only 10.2% of the
predictions fell within 10% of their training ground truth, irrespective of
distance into the future. For accurate assessment of uncertainty, only one
model matched relatively well the nominal 95% coverage, but that model did not
start predictions until April 16, thus had no impact on early, major decisions.
For ICU bed utilisation, the IHME model was highly inaccurate; the point
estimates only started to match ground truth after the pandemic wave had
started to wane. We conclude that trustworthy models require trustworthy input
data to be trained upon. Moreover, models need to be subjected to prespecified
real time performance tests, before their results are provided to policy makers
and public health officials.


  Ecological systems, as is often noted, are complex. Equally notable is the
generalization that complex systems tend to be oscillatory, whether Huygens
simple patterns of pendulum entrainment or the twisted chaotic orbits of Lorenz
convection rolls. The analytics of oscillators may thus provide insight into
the structure of ecological systems. One of the most popular analytical tools
for such study is the Kuramoto model of coupled oscillators. Using a
well-studied system of pests and their enemies in an agroecosystem, we apply
this model as a stylized vision of the dynamics of that real system, to ask
whether its actual natural history is reflected in the dynamics of the
qualitatively instantiated Kuramoto model. Emerging from the model is a series
of synchrony groups generally corresponding to subnetworks of the natural
system, with an overlying chimeric structure, depending on the strength of the
inter-oscillator coupling. We conclude that the Kuramoto model presents a novel
window through which interesting questions about the structure of ecological
systems may emerge.


  A mathematical model for the COVID-19 pandemic spread, which integrates
age-structured Susceptible-Exposed-Infected-Recovered-Deceased dynamics with
real mobile phone data accounting for the population mobility, is presented.
The dynamical model adjustment is performed via Approximate Bayesian
Computation. Optimal lockdown and exit strategies are determined based on
nonlinear model predictive control, constrained to public-health and
socio-economic factors. Through an extensive computational validation of the
methodology, it is shown that it is possible to compute robust exit strategies
with realistic reduced mobility values to inform public policy making, and we
exemplify the applicability of the methodology using datasets from England and
France. Code implementing the described experiments is available at
https://github.com/OptimalLockdown.


  Ecological, environmental and geophysical time series consistently exhibit
the characteristics of coloured (1/f^\b{eta}) noise. Here we briefly survey the
literature on coloured noise, population persistence and related evolutionary
dynamics, before introducing coloured noise as an appropriate model for
environmental variation in artificial evolutionary systems. To illustrate and
explore the effects of different noise colours, a simple evolutionary model
that examines the trade-off between specialism and generalism in fluctuating
environments is applied. The results of the model clearly demonstrate a need
for greater generalism as environmental variability becomes `whiter', whilst
specialisation is favoured as environmental variability becomes `redder'. Pink
noise, sitting midway between white and red noise, is shown to be the point at
which the pressures for generalism and specialism balance, providing some
insight in to why `pinker' noise is increasingly being seen as an appropriate
model of typical environmental variability. We go on to discuss how the results
presented here feed in to a wider discussion on evolutionary responses to
fluctuating environments. Ultimately we argue that Artificial Life as a field
should embrace the use of coloured noise to produce models of environmental
variability.


  In this article, we perform an analysis of COVID-19 on one of the South
Saharan countries (hot zone), the Senegal (West Africa). Many questions remain
unanswered: why the African continent is not very contaminated compared to
other continents. Factors of cross immunity, temperature, population density,
youth, etc. are taken into account for an analysis of the contamination
factors. Numerical simulations are carried out for a prediction over the coming
week.


  To curb the spread of COVID-19, many governments around the world have
implemented tiered lockdowns with varying degrees of stringency. Lockdown
levels are typically increased when the disease spreads and reduced when the
disease abates. A predictive control approach is used to develop optimized
lockdown strategies for curbing the spread of COVID-19. The strategies are then
applied to South African data. The South African case is of interest as the
South African government has defined five distinct levels of lockdown, which
serves as a discrete control input. An epidemiological model for the spread of
COVID-19 in South Africa was previously developed, and is used in conjunction
with a hybrid model predictive controller to optimize lockdown management under
different policy scenarios. Scenarios considered include how to flatten the
curve to a level that the healthcare system can cope with, how to balance lives
and livelihoods, and what impact the compliance of the population to the
lockdown measures has on the spread of COVID-19. The main purpose of this paper
is to show what the optimal lockdown level should be given the policy that is
in place, as determined by the closed-loop feedback controller.


  The current COVID-19 pandemic and subsequent lockdowns have highlighted the
close and delicate relationship between a country's public health and economic
health. Macroeconomic models that use preexisting epidemic models to calculate
the impacts of a disease outbreak are therefore extremely useful for
policymakers seeking to evaluate the best course of action in such a crisis. We
develop an SIR model of the COVID-19 pandemic that explicitly considers herd
immunity, behavior-dependent transmission rates, remote workers, and indirect
externalities of lockdown. This model is presented as an exit time control
problem where lockdown ends when the population achieves herd immunity, either
naturally or via a vaccine. A social planner prescribes separate levels of
lockdown for two separate sections of the adult population: low-risk (ages
20-64) and high-risk (ages 65 and over). These levels are determined via
optimization of an objective function which assigns a macroeconomic cost to the
level of lockdown and the number of deaths. We find that, by ending lockdowns
once herd immunity is reached, high-risk individuals are able to leave lockdown
significantly before the arrival of a vaccine without causing large increases
in mortality. Moreover, if we incorporate a behavior-dependent transmission
rate which represents increased personal caution in response to increased
infection levels, both output loss and total mortality are lowered. Lockdown
efficacy is further increased when there is less interaction between low- and
high-risk individuals, and increased remote work decreases output losses.
Overall, our model predicts that a lockdown which ends at the arrival of herd
immunity, combined with individual actions to slow virus transmission, can
reduce total mortality to one-third of the no-lockdown level, while allowing
high-risk individuals to leave lockdown well before vaccine arrival.


  By assuming that tau protein can be in seven kinetic states, we developed a
model of tau protein transport in the axon and in the axon initial segment
(AIS). Two separate sets of kinetic constants were determined, one in the axon
and the other in the AIS. This was done by fitting the model predictions in the
axon with experimental results and by fitting the model predictions in the AIS
with the assumed linear increase of the total tau concentration in the AIS. The
calibrated model was used to make predictions about tau transport in the axon
and in the AIS. To the best of our knowledge, this is the first paper that
presents a mathematical model of tau transport in the AIS. Our modeling results
suggest that binding of free tau to MTs creates a negative gradient of free tau
in the AIS. This leads to diffusion-driven tau transport from the soma into the
AIS. The model further suggests that slow axonal transport and diffusion-driven
transport of tau work together in the AIS, moving tau anterogradely. Our
numerical results predict an interplay between these two mechanisms: as the
distance from the soma increases, the diffusion-driven transport decreases,
while motor-driven transport becomes larger. Thus, the machinery in the AIS
works as a pump, moving tau into the axon.


  Renewal processes are a popular approach used in modelling infectious disease
outbreaks. In a renewal process, previous infections give rise to future
infections. However, while this formulation seems sensible, its application to
infectious disease can be difficult to justify from first principles. It has
been shown from the seminal work of Bellman and Harris that the renewal
equation arises as the expectation of an age-dependent branching process. In
this paper we provide a detailed derivation of the original Bellman Harris
process. We introduce generalisations, that allow for time-varying reproduction
numbers and the accounting of exogenous events, such as importations. We show
how inference on the renewal equation is easy to accomplish within a Bayesian
hierarchical framework. Using off the shelf MCMC packages, we fit to South
Korea COVID-19 case data to estimate reproduction numbers and importations. Our
derivation provides the mathematical fundamentals and assumptions underpinning
the use of the renewal equation for modelling outbreaks.


  The COVID-19 pandemic has created unprecedented challenges worldwide.
Strained healthcare providers make difficult decisions on patient triage,
treatment and care management on a daily basis. Policy makers have imposed
social distancing measures to slow the disease, at a steep economic price. We
design analytical tools to support these decisions and combat the pandemic.
Specifically, we propose a comprehensive data-driven approach to understand the
clinical characteristics of COVID-19, predict its mortality, forecast its
evolution, and ultimately alleviate its impact. By leveraging cohort-level
clinical data, patient-level hospital data, and census-level epidemiological
data, we develop an integrated four-step approach, combining descriptive,
predictive and prescriptive analytics. First, we aggregate hundreds of clinical
studies into the most comprehensive database on COVID-19 to paint a new
macroscopic picture of the disease. Second, we build personalized calculators
to predict the risk of infection and mortality as a function of demographics,
symptoms, comorbidities, and lab values. Third, we develop a novel
epidemiological model to project the pandemic's spread and inform social
distancing policies. Fourth, we propose an optimization model to re-allocate
ventilators and alleviate shortages. Our results have been used at the clinical
level by several hospitals to triage patients, guide care management, plan ICU
capacity, and re-distribute ventilators. At the policy level, they are
currently supporting safe back-to-work policies at a major institution and
equitable vaccine distribution planning at a major pharmaceutical company, and
have been integrated into the US Center for Disease Control's pandemic
forecast.


  We proposed a three-species food chain model with hunting co-operation among
the middle predator. In this model, third species prey on the middle species
and the middle prey on the first species. The hunting cooperation among the
middle predator affects interestingly on the numbers of both the predators and
the prey. We examined the linear stability of the model theoretically and
numerically. We conducted the two-parameter numerical analysis to check the
long-term behavior and the change in the number of species with respect to
hunting co-operation. Our findings supported the postulates from the two
species food chain model with hunting co-operation.


  We consider a population evolving due to mutation, selection and
recombination, where selection includes single-locus terms (additive fitness)
and two-loci terms (pairwise epistatic fitness). We further consider the
problem of inferring fitness in the evolutionary dynamics from one or several
snap-shots of the distribution of genotypes in the population. In the recent
literature this has been done by applying the Quasi-Linkage Equilibrium (QLE)
regime first obtained by Kimura in the limit of high recombination. Here we
show that the approach also works in the interesting regime where the effects
of mutations are comparable to or larger than recombination. This leads to a
modified main epistatic fitness inference formula where the rates of mutation
and recombination occur together. We also derive this formula using by a
previously developed Gaussian closure that formally remains valid when
recombination is absent. The findings are validated through numerical
simulations.


  Transport in crowded, complex environments occurs across many spatial scales.
Geometric restrictions can hinder the motion of individuals and, combined with
crowding between individuals, can have drastic effects on global transport
phenomena. However, in general, the interplay between crowding and geometry in
complex real-life environments is poorly understood. Existing analytical
methodologies are not always readily extendable to heterogeneous environments:
in these situations predictions of crowded transport behaviour within
heterogeneous environments rely on computationally intensive mesh-based
approaches. Here, we take a different approach by employing networked
representations of complex environments to provide an efficient framework
within which the interactions between networked geometry and crowding can be
explored. We demonstrate how the framework can be used to: extract detailed
information at the level of the whole population or an individual within it;
identify the topological features of environments that enable accurate
prediction of transport phenomena; and, provide insights into the design of
optimal environments.


  The ongoing COVID-19 pandemic calls for a multi-faceted public health
response comprising complementary interventions to control the spread of the
disease while vaccines and therapies are developed. Many of these interventions
need to be informed by epidemic risk predictions given available data,
including symptoms, contact patterns, and environmental factors. Here we
propose a novel probabilistic formalism based on Individual-Level Models (ILMs)
that offers rigorous formulas for the probability of infection of individuals,
which can be parameterised via Maximum Likelihood Estimation (MLE) applied on
compartmental models defined at the population level. We describe an approach
where individual data collected in real-time is integrated with overall case
counts to update the a predictor of the susceptibility of infection as a
function of individual risk factors.


  Best match graphs (BMGs) are vertex-colored directed graphs that were
introduced to model the relationships of genes (vertices) from different
species (colors) given an underlying evolutionary tree that is assumed to be
unknown. In real-life applications, BMGs are estimated from sequence similarity
data. Measurement noise and approximation errors usually result in empirically
determined graphs that in general violate characteristic properties of BMGs.
The arc modification problems for BMGs aim at correcting such violations and
thus provide a means to improve the initial estimates of best match data. We
show here that the arc deletion, arc completion and arc editing problems for
BMGs are NP-complete and that they can be formulated and solved as integer
linear programs. To this end, we provide a novel characterization of BMGs in
terms of triples (binary trees on three leaves) and a characterization of BMGs
with two colors in terms of forbidden subgraphs.


  The outbreak of COVID-19 disrupts the life of many people in the world. The
state of Arizona in the U.S. emerges as one of the country's newest COVID-19
hot spots. Accurate forecasting for COVID-19 cases will help governments to
implement necessary measures and convince more people to take personal
precautions to combat the virus. It is difficult to accurately predict the
COVID-19 cases due to many human factors involved. This paper aims to provide a
forecasting model for COVID-19 cases with the help of human activity data from
the Google Community Mobility Reports. To achieve this goal, a specific partial
differential equation (PDE) is developed and validated with the COVID-19 data
from the New York Times at the county level in the state of Arizona in the U.S.
The proposed model describes the combined effects of transboundary spread among
county clusters in Arizona and human actives on the transmission of COVID-19.
The results show that the prediction accuracy of this model is well acceptable
(above 94\%). Furthermore, we study the effectiveness of personal precautions
such as wearing face masks and practicing social distancing on COVID-19 cases
at the local level. The localized analytical results can be used to help to
slow the spread of COVID-19 in Arizona. To the best of our knowledge, this work
is the first attempt to apply PDE models on COVID-19 prediction with the Google
Community Mobility Reports.


  Designing compounds with desired properties is a key element of the drug
discovery process. However, measuring progress in the field has been
challenging due to the lack of realistic retrospective benchmarks, and the
large cost of prospective validation. To close this gap, we propose a benchmark
based on docking, a popular computational method for assessing molecule binding
to a protein. Concretely, the goal is to generate drug-like molecules that are
scored highly by SMINA, a popular docking software. We observe that popular
graph-based generative models fail to generate molecules with a high docking
score when trained using a realistically sized training set. This suggests a
limitation of the current incarnation of models for de novo drug design.
Finally, we propose a simplified version of the benchmark based on a simpler
scoring function, and show that the tested models are able to partially solve
it. We release the benchmark as an easy to use package available at
https://github.com/cieplinski-tobiasz/smina-docking-benchmark. We hope that our
benchmark will serve as a stepping stone towards the goal of automatically
generating promising drug candidates.