  Patagonian fjords' area is one of the largest estuarine regions in the world.
Every one of its water bodies displays a unique hydrodynamic behavior with
enormous effects on the biogeochemical characteristics of the ecosystems. In
this context, algal blooms are ecological phenomena of major relevance.
Numerical simulation has proved to be a promising tool to understand their
impacts. It has not been used for studying algal blooms in this zone. This
article focuses on proposing a novel numerical model for simulating brief algal
blooms occurring in water bodies of western Patagonia. The proposed model
presents a trade-off between complexity and applicability since field-data
sparsity in the zone discourages using more sophisticated approaches. The model
is based on a two-layer description of the water column. The first layer
represents the euphotic zone where an embedded biogeochemical model of
NPZD-type is used to model a mass-conserving trophic web. High intensity wind
drives the water column mixing, introducing an upward flux of nutrients that
boosts high rates of primary production. A time-dependent Gaussian pulse is
used to describe this process. Mass losses due to detritus sinking are also
included. Then, the ecosystem's dynamics is represented by means of an
externally forced, non-autonomous system of ordinary differential equations
which is characterized by strictly positive trajectories but that it is not
longer mass-conserving. A structure-preserving time integrator based on a
splitting-composition technique is designed for solving the system's equations.
It is cast as a three-steps algorithm and provides an exact estimations of
biomass fluxes. Additionally, a genetic algorithm-based tool is used to
calibrate the model's parameters in realistic scenarios. The proposed model is
applied im a detailed study of a winter bloom in an austral fjord.


  Contact tracing is an effective method to control emerging diseases. Since
the 1980's, mathematical modelers are developing a consistent theory for
contact tracing, with the aim to find effective and efficient implementations
of contact tracing, and to assess the effects of contact tracing on the spread
of an infectious disease. Despite the progress made in the area, there remain
important open questions. In addition, technological developments, especially
in the field of molecular biology (genetic sequencing of pathogens) and modern
communication (digital contact tracing), have posed new challenges for the
modeling community. In the present paper, we discuss modeling approaches for
contact tracing and identify some of the current challenges for the field.


  Somatic hypermutations of immunoglobulin (Ig) genes occuring during affinity
maturation drive B-cell receptors' ability to evolve strong binding to their
antigenic targets. The landscape of these mutations is highly heterogeneous,
with certain regions of the Ig gene being preferentially targeted. However, a
rigorous quantification of this bias has been difficult because of phylogenetic
correlations between sequences and the interference of selective forces. Here,
we present an approach that corrects for these issues, and use it to learn a
model of hypermutation preferences from a recently published large IgH
repertoire dataset. The obtained model predicts mutation profiles accurately
and in a reproducible way, including in the previously uncharacterized
Complementarity Determining Region 3, revealing that both the sequence context
of the mutation and its absolute position along the gene are important. In
addition, we show that hypermutations occurring concomittantly along B-cell
lineages tend to co-localize, suggesting a possible mechanism for accelerating
affinity maturation.


  In this article, we discuss an age-structured SIR model in which disease not
only spread through direct person to person contacts for e.g. infection due to
surface contamination but it can also spread through indirect contacts. It is
evident that age also plays a crucial role in SARS virus infection including
COVID-19 infection. We formulate our model as an abstract semilinear Cauchy
problem in an appropriate Banach space to show the existence of solution and
also show the existence of steady states. It is assumed in this work that the
population is in a demographic stationary state and show that there is no
disease-free equilibrium point as long as there is a transmission of infection
due to the indirect contacts in the environment.


  This paper presents the assessment of time-dependent national-level
restrictions and control actions and their effects in fighting the COVID-19
pandemic. By analysing the transmission dynamics during the first wave of
COVID-19 in the country, the effectiveness of the various levels of control
actions taken to flatten the curve can be better quantified and understood.
This in turn can help the relevant authorities to better plan for and control
the subsequent waves of the pandemic. To achieve this, a deterministic
population model for the pandemic is firstly developed to take into
consideration the time-dependent characteristics of the model parameters,
especially on the ever-evolving value of the reproduction number, which is one
of the critical measures used to describe the transmission dynamics of this
pandemic. The reproduction number alongside other key parameters of the model
can then be estimated by fitting the model to real-world data using numerical
optimisation techniques or by inducing ad-hoc control actions as recorded in
the news platforms. In this paper, the model is verified using a case study
based on the data from the first wave of COVID-19 in the Republic of
Kazakhstan. The model is fitted to provide estimates for two settings in
simulations; time-invariant and time-varying (with bounded constraints)
parameters. Finally, some forecasts are made using four scenarios with
time-dependent control measures so as to determine which would reflect on the
actual situations better.


  Understanding E3 ligase and target substrate interactions are important for
cell biology and therapeutic development. However, experimental identification
of E3 target relationships is not an easy task due to the labor-intensive
nature of the experiments. In this article, a sequence-based E3-target
prediction model is proposed for the first time. The proposed framework
utilizes composition of k-spaced amino acid pairs (CKSAAP) to learn the
relationship between E3 ligases and their target protein. A class separable
latent space encoding scheme is also devised that provides a compressed
representation of feature space. A thorough ablation study is performed to
identify an optimal gap size for CKSAAP and the number of latent variables that
can represent the E3-target relationship successfully. The proposed scheme is
evaluated on an independent dataset for a variety of standard quantitative
measures. In particular, it achieves an average accuracy of $70.63\%$ on an
independent dataset. The source code and datasets used in the study are
available at the author's GitHub page
(https://github.com/psychemistz/E3targetPred).


  Microorganisms live in environments that inevitably fluctuate between mild
and harsh conditions. As harsh conditions may cause extinctions, the rate at
which fluctuations occur can shape microbial communities and their diversity,
but we still lack an intuition on how. Here, we build a mathematical model
describing two microbial species living in an environment where substrate
supplies randomly switch between abundant and scarce. We then vary the rate of
switching as well as different properties of the interacting species, and
measure the probability of the weaker species driving the stronger one extinct.
We find that this probability increases with the strength of demographic noise
under harsh conditions and peaks at either low, high, or intermediate switching
rates depending on both species' ability to withstand the harsh environment.
This complex relationship shows why finding patterns between environmental
fluctuations and diversity has historically been difficult. In parameter ranges
where the fittest species was most likely to be excluded, however, the beta
diversity in larger communities also peaked. In sum, how environmental
fluctuations affect interactions between a few species pairs predicts their
effect on the beta diversity of the whole community.


  The new Coronavirus (COVID-19) is an emerging disease responsible for
infecting millions of people since the first notification until nowadays.
Developing efficient short-term forecasting models allow knowing the number of
future cases. In this context, it is possible to develop strategic planning in
the public health system to avoid deaths. In this paper, autoregressive
integrated moving average (ARIMA), cubist (CUBIST), random forest (RF), ridge
regression (RIDGE), support vector regression (SVR), and stacking-ensemble
learning are evaluated in the task of time series forecasting with one, three,
and six-days ahead the COVID-19 cumulative confirmed cases in ten Brazilian
states with a high daily incidence. In the stacking learning approach, the
cubist, RF, RIDGE, and SVR models are adopted as base-learners and Gaussian
process (GP) as meta-learner. The models' effectiveness is evaluated based on
the improvement index, mean absolute error, and symmetric mean absolute
percentage error criteria. In most of the cases, the SVR and stacking ensemble
learning reach a better performance regarding adopted criteria than compared
models. In general, the developed models can generate accurate forecasting,
achieving errors in a range of 0.87% - 3.51%, 1.02% - 5.63%, and 0.95% - 6.90%
in one, three, and six-days-ahead, respectively. The ranking of models in all
scenarios is SVR, stacking ensemble learning, ARIMA, CUBIST, RIDGE, and RF
models. The use of evaluated models is recommended to forecasting and monitor
the ongoing growth of COVID-19 cases, once these models can assist the managers
in the decision-making support systems.


  To increase situational awareness and support evidence-based policy-making,
we formulated two types of mathematical models for COVID-19 transmission within
a regional population. One is a fitting function that can be calibrated to
reproduce an epidemic curve with two timescales (e.g., fast growth and slow
decay). The other is a compartmental model that accounts for quarantine,
self-isolation, social distancing, a non-exponentially distributed incubation
period, asymptomatic individuals, and mild and severe forms of symptomatic
disease. Using Bayesian inference, we have been calibrating our models daily
for consistency with new reports of confirmed cases from the 15 most populous
metropolitan statistical areas in the United States and quantifying uncertainty
in parameter estimates and predictions of future case reports. This online
learning approach allows for early identification of new trends despite
considerable variability in case reporting. We infer new significant upward
trends for five of the metropolitan areas starting between 19-April-2020 and
12-June-2020.


  In this paper we propose an eco-friendly optimization of banana or plantain
yield by the control of the pest burrowing nematode \textit{Radopholus
similis}. This control relies on fallow deployment, with greater respect for
the environment than chemical methods. The optimization is based on a
multi-seasonal model in which fallow periods follow cropping seasons. The aim
is to find the best way, in terms of profit, to allocate the durations of
fallow periods between the cropping seasons, over a fixed time horizon spanning
several seasons. The existence of an optimal allocation is proven and an
adaptive random search algorithm is proposed to solve the optimization problem.
For a relatively long time horizon, deploying one season less than the maximum
possible number of cropping seasons allows to increase the fallow period
durations and results in a better multi-seasonal profit. For regular fallow
durations, the profit is lower than the optimal solution, but the final soil
infestation is also lower.


  Background: The impact of individual non-pharmaceutical interventions (NPI)
such as state-wide stay-at-home orders, school closures and gathering size
limitations, on the COVID-19 epidemic is unknown. Understanding the impact that
above listed NPI have on disease transmission is critical for policy makers,
particularly as case counts increase again in some areas.
  Methods: Using a Bayesian framework, we reconstructed the incidence and
time-varying reproductive number (Rt) curves to investigate the relationship
between Rt, individual mobility as measured by Google Community Mobility
Reports, and NPI.
  Results: We found a strong relationship between reproductive number and
mobility, with each 10% drop in mobility being associated with an expected
10.2% reduction in Rt compared to baseline. The effects of limitations on the
size of gatherings, school and business closures, and stay-at-home orders were
dominated by the trend over time, which was associated with a 48% decrease in
the reproductive number, adjusting for the NPI.
  Conclusions: We found that the decrease in mobility associated with time may
be due to individuals changing their behavior in response to perceived risk or
external factors.


  In both criminal cases and civil cases there is an increasing demand for the
analysis of DNA mixtures involving relationships. The goal might be, for
example, to identify the contributors to a DNA mixture where the donors may be
related, or to infer the relationship between individuals based on a DNA
mixture. This paper applies a recent approach to modelling and computation for
DNA mixtures involving contributors with arbitrarily complex relationships to
two real cases from the Spanish Forensic Police.


  The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been
mutating since it was first sequenced in early January 2020. The genetic
variants have developed into a few distinct clusters with different properties.
Since the United States (US) has the highest number of viral infected patients
globally, it is essential to understand the US SARS-CoV-2. Using genotyping,
sequence-alignment, time-evolution, $k$-means clustering, protein-folding
stability, algebraic topology, and network theory, we reveal that the US
SARS-CoV-2 has four substrains and five top US SARS-CoV-2 mutations were first
detected in China (2 cases), Singapore (2 cases), and the United Kingdom (1
case). The next three top US SARS-CoV-2 mutations were first detected in the
US. These eight top mutations belong to two disconnected groups. The first
group consisting of 5 concurrent mutations is prevailing, while the other group
with three concurrent mutations gradually fades out. Our analysis suggests that
female immune systems are more active than those of males in responding to
SARS-CoV-2 infections. We identify that one of the top mutations,
27964C$>$T-(S24L) on ORF8, has an unusually strong gender dependence. Based on
the analysis of all mutations on the spike protein, we further uncover that
three of four US SASR-CoV-2 substrains become more infectious. Our study calls
for effective viral control and containing strategies in the US.


  Species extinction is a core process that affects the diversity of life on
Earth. Competition between species in a population is considered by ecological
niche-based theories as a key factor leading to different severity of species
extinctions. There are population dynamics models that describe a simple and
easily understandable mechanism for resource competition. However, these models
can not efficiently characterize and quantify new emergent extinctions in a
large population appearing due to environmental forcing. To address this issue
we develop a stochastic physics-inspired approach to analyze how environmental
forcing influences the severity of species extinctions in such models. This
approach is based on the large deviations theory of stochastic processes (the
Freidlin-Wentzell theory). We show that there are three possible fundamentally
different scenarios of extinctions, which we call catastrophic extinctions,
asymmetric ones, and extinctions with exponentially small probabilities. The
realization of those scenarios depends on environmental noise properties and
the boundaries of niches, which define the domain, where species survive.
Furthermore, we describe a hysteresis effect in species extinction showing that
fluctuations can lead to dramatic consequences even if an averaged resource
supply is sufficient to support population survival. Our stochastic
physics-inspired approach generalizes niche theory by accounting environmental
forcing and will be useful to find, by available data, which environmental
perturbations may induce extinctions.


  Protein function prediction may be framed as predicting subgraphs (with
certain closure properties) of a directed acyclic graph describing the
hierarchy of protein functions. Graph neural networks (GNNs), with their
built-in inductive bias for relational data, are hence naturally suited for
this task. However, in contrast with most GNN applications, the graph is not
related to the input, but to the label space. Accordingly, we propose
Tail-GNNs, neural networks which naturally compose with the output space of any
neural network for multi-task prediction, to provide relationally-reinforced
labels. For protein function prediction, we combine a Tail-GNN with a dilated
convolutional network which learns representations of the protein sequence,
making significant improvement in F_1 score and demonstrating the ability of
Tail-GNNs to learn useful representations of labels and exploit them in
real-world problem solving.


  The SIQR model is reformulated where compartments for infected and
quarantined are redefined so as to be appropriate to COVID-19, and exact
properties of the model are presented. It is shown that the maximum number of
infected at large depends strongly on the quarantine rate and that the
quarantine measure is more effective than the lockdown measure in controlling
the pandemic. The peak of the number of quarantined patients is shown to appear
some time later than the time that the number of infected becomes maximum. On
the basis of the expected utility theory, a theoretical framework to find out
an optimum strategy in the space of lockdown measure and quarantine measure is
proposed for minimizing the maximum number of infected and for controlling the
outbreak of pandemic at its early stage.


  We use analytical methods to investigate a continuous vaccination strategy
effects on the infectious disease dynamics in a closed population and a
demographically open population. The methodology and key assumptions are based
on Breda et al (2012). We show that the cumulative force of infection for the
closed population and the endemic force of infection in the demographically
open population can be reduced significantly by combining two factors: the
vaccine effectiveness and the vaccination rate. The impact of these factors on
the force of infection can transform an endemic steady state into a
disease-free state. Keywords: Force of infection, Cumulative force of
infection, Scalar-renewal equation, Per capita death rate, Lambert function,
adaptive vaccination strategy


  A second wave pandemic constitutes an imminent threat to society, with a
potentially immense toll in terms of human lives and a devastating economic
impact. We employ the epidemic renormalisation group approach to pandemics,
together with the first wave data for COVID-19, to efficiently simulate the
dynamics of disease transmission and spreading across different European
countries. The framework allows us to model, not only inter and extra European
border control effects, but also the impact of social distancing for each
country. We perform statistical analyses averaging on different level of human
interaction across Europe and with the rest of the world. Our results are
neatly summarised as an animation reporting the time evolution of the first and
second waves of the European COVID-19 pandemic. Our temporal playbook of the
second wave pandemic can be used by governments, financial markets, the
industries and individual citizens, to efficiently time, prepare and implement
local and global measures.


  Nascent adhesions are submicron transient structures promoting the early
adhesion of cells to the extracellular matrix. Nascent adhesions typically
consist of several tens of integrins, and serve as platforms for the
recruitment and activation of proteins to build mature focal adhesions. They
are also associated with early stage signalling and the mechanoresponse.
Despite their crucial role in sampling the local extracellular matrix, very
little is known about the mechanism of their formation. Consequently, there is
a strong scientific activity focused on elucidating the physical and
biochemical foundation of their development and function. Precisely the results
of this effort will be summarized in this article.


  To what extent are effectiveness estimates of nonpharmaceutical interventions
(NPIs) against COVID-19 influenced by the assumptions our models make? To
answer this question, we investigate 2 state-of-the-art NPI effectiveness
models and propose 6 variants that make different structural assumptions. In
particular, we investigate how well NPI effectiveness estimates generalise to
unseen countries, and their sensitivity to unobserved factors. Models that
account for noise in disease transmission compare favourably. We further
evaluate how robust estimates are to different choices of epidemiological
parameters and data. Focusing on models that assume transmission noise, we find
that previously published results are remarkably robust across these variables.
Finally, we mathematically ground the interpretation of NPI effectiveness
estimates when certain common assumptions do not hold.


  Recent events leading to the worldwide pandemic of COVID-19 have demonstrated
the effective use of genomic sequencing technologies to establish the genetic
sequence of this virus. In contrast, the COVID-19 pandemic has demonstrated the
absence of computational approaches to understand the molecular basis of this
infection rapidly. Here we present an integrated approach to the study of the
nsp1 protein in SARS-CoV-1, which plays an essential role in maintaining the
expression of viral proteins and further disabling the host protein expression,
also known as the host shutoff mechanism. We present three independent methods
of evaluating two potential binding sites speculated to participate in host
shutoff by nsp1. We have combined results from computed models of nsp1, with
deep mining of all existing protein structures (using PDBMine), and binding
site recognition (using msTALI) to examine the two sites consisting of residues
55-59 and 73-80. Based on our preliminary results, we conclude that the
residues 73-80 appear as the regions that facilitate the critical initial steps
in the function of nsp1. Given the 90% sequence identity between nsp1 from
SARS-CoV-1 and SARS-CoV-2, we conjecture the same critical initiation step in
the function of COVID-19 nsp1.


  In this paper we propose a three stages analysis of the evolution of Covid19
in Romania.
  There are two main issues when it comes to pandemic prediction. The first one
is the fact that the numbers reported of infected and recovered are unreliable,
however the number of deaths is more accurate. The second issue is that there
were many factors which affected the evolution of the pandemic.
  In this paper we propose an analysis in three stages. The first stage is
based on the classical SIR model which we do using a neural network. This
provides a first set of daily parameters.
  In the second stage we propose a refinement of the SIR model in which we
separate the deceased into a distinct category. By using the first estimate and
a grid search, we give a daily estimation of the parameters.
  The third stage is used to define a notion of turning points (local extremes)
for the parameters. We call a regime the time between these points.
  We outline a general way based on time varying parameters of SIRD to make
predictions.


  A central mechanism in machine learning is to identify, store, and recognize
patterns. How to learn, access, and retrieve such patterns is crucial in
Hopfield networks and the more recent transformer architectures. We show that
the attention mechanism of transformer architectures is actually the update
rule of modern Hopfield networks that can store exponentially many patterns. We
exploit this high storage capacity of modern Hopfield networks to solve a
challenging multiple instance learning (MIL) problem in computational biology:
immune repertoire classification. Accurate and interpretable machine learning
methods solving this problem could pave the way towards new vaccines and
therapies, which is currently a very relevant research topic intensified by the
COVID-19 crisis. Immune repertoire classification based on the vast number of
immunosequences of an individual is a MIL problem with an unprecedentedly
massive number of instances, two orders of magnitude larger than currently
considered problems, and with an extremely low witness rate. In this work, we
present our novel method DeepRC that integrates transformer-like attention, or
equivalently modern Hopfield networks, into deep learning architectures for
massive MIL such as immune repertoire classification. We demonstrate that
DeepRC outperforms all other methods with respect to predictive performance on
large-scale experiments, including simulated and real-world virus infection
data, and enables the extraction of sequence motifs that are connected to a
given disease class. Source code and datasets: https://github.com/ml-jku/DeepRC


  Hepatitis C is a viral infection that appears as a result of the Hepatitis C
Virus (HCV), and it has been recognized as the main reason for liver diseases.
HCV incidence is growing as an important issue in the epidemiology of
infectious diseases. In the present study, a mathematical model is employed for
simulating the dynamics of HCV outbreak in a population. The total population
is divided into five compartments, including unaware and aware susceptible,
acutely and chronically infected, and treated classes. Then, a Lyapunov-based
nonlinear adaptive method is proposed for the first time to control the HCV
epidemic considering modelling uncertainties. A positive definite Lyapunov
candidate function is suggested, and adaptation and control laws are attained
based on that. The main goal of the proposed control strategy is to decrease
the population of unaware susceptible and chronically infected compartments by
pursuing appropriate treatment scenarios. As a consequence of this decrease in
the mentioned compartments, the population of aware susceptible individuals
increases and the population of acutely infected and treated humans decreases.
The Lyapunov stability theorem and Barbalat's lemma are employed in order to
prove the tracking convergence to desired population reduction scenarios. Based
on the acquired numerical results, the proposed nonlinear adaptive controller
can achieve the above-mentioned objective by adjusting the inputs (rates of
informing the susceptible people and treatment of chronically infected ones)
and estimating uncertain parameter values based on the designed control and
adaptation laws, respectively. Moreover, the proposed strategy is designed to
be robust in the presence of different levels of parametric uncertainties.


  The study carries out predictive modeling based on publicly available
COVID-19 data for the duration 01 April to 20 June 2020 pertaining to India and
five of its most infected states: Maharashtra, Tamil Nadu, Delhi, Gujarat, and
Rajasthan using susceptible, infected, recovered, and dead (SIRD) model. The
basic reproduction number R0 is derived by exponential growth method using
RStudio package R0. The differential equations reflecting SIRD model have been
solved using Python 3.7.4 on Jupyter Notebook platform. For visualization,
Python Matplotlib 3.2.1 package is used. The study offers insights on
peak-date, peak number of COVID-19 infections, and end-date pertaining to India
and five of its states. The results could be leveraged by political leadership,
health authorities, and industry doyens for policy planning and execution.


  We have developed a globally applicable diagnostic Covid-19 model by
augmenting the classical SIR epidemiological model with a neural network
module. Our model does not rely upon previous epidemics like SARS/MERS and all
parameters are optimized via machine learning algorithms employed on publicly
available Covid-19 data. The model decomposes the contributions to the
infection timeseries to analyze and compare the role of quarantine control
policies employed in highly affected regions of Europe, North America, South
America and Asia in controlling the spread of the virus. For all continents
considered, our results show a generally strong correlation between
strengthening of the quarantine controls as learnt by the model and actions
taken by the regions' respective governments. Finally, we have hosted our
quarantine diagnosis results for the top 70 affected countries worldwide, on a
public platform, which can be used for informed decision making by public
health officials and researchers alike.


  Identifying the relative importance of the different transmission routes of
the SARS-CoV-2 virus is an urgent research priority. To that end, the different
transmission routes, and their role in determining the evolution of the
Covid-19 pandemic are analyzed in this work. Probability of infection caused by
inhaling virus-laden droplets (initial, ejection diameters between $0.5-750\mu
m$) and the corresponding desiccated nuclei that mostly encapsulate the virions
post droplet evaporation, are individually calculated. At typical,
air-conditioned yet quiescent indoor space, for average viral loading, cough
droplets of initial diameter between $10-50 \mu m$ have the highest infection
probability. However, by the time they are inhaled, the diameters reduce to
about $1/6^{th}$ of their initial diameters. While the initially near unity
infection probability due to droplets rapidly decays within the first $25s$,
the small yet persistent infection probability of desiccated nuclei decays
appreciably only by $\mathcal{O} (1000s)$, assuming the virus sustains equally
well within the dried droplet nuclei as in the droplets. Combined with
molecular collision theory adapted to calculate frequency of contact between
the susceptible population and the droplet/nuclei cloud, infection rate
constants are derived ab-initio, leading to a SEIR model applicable for any
respiratory event - vector combination. Viral load, minimum infectious dose,
sensitivity of the virus half-life to the phase of its vector and dilution of
the respiratory jet/puff by the entraining air are shown to mechanistically
determine specific physical modes of transmission and variation in the basic
reproduction number $\mathcal{R}_0$, from first principle calculations.


  We study the problem of finding a temporal hybridization network for a set of
phylogenetic trees that minimizes the number of reticulations. First, we
introduce an FPT algorithm for this problem on an arbitrary set of $m$ binary
trees with $n$ leaves each with a running time of $O(5^k\cdot n\cdot m)$, where
$k$ is the minimum temporal hybridization number. We also present the concept
of temporal distance, which is a measure for how close a tree-child network is
to being temporal. Then we introduce an algorithm for computing a tree-child
network with temporal distance at most $d$ and at most $k$ reticulations in
$O((8k)^d5^ k\cdot n\cdot m)$ time. Lastly, we introduce a $O(6^kk!\cdot k\cdot
n^2)$ time algorithm for computing a minimum temporal hybridization network for
a set of two nonbinary trees. We also provide an implementation of all
algorithms and an experimental analysis on their performance.


  The novel coronavirus SARS-CoV-2 has rapidly emerged as a significant threat
to global public health, in particular because -- as is not uncommon with novel
pathogens -- there is no effective pharmaceutical treatment or prophylaxis to
the viral syndrome it causes. In the absence of such specific treatment
modalities, the mainstay of public health response rests on non-pharmaceutical
interventions (NPIs), such as social distancing. This paper contributes to the
understanding of social distancing against SARS-CoV-2 by quantitatively
analysing the statistical dynamics of disease propagation as a differential
game, and estimating the relative costs of distancing versus not distancing,
identifying marginal utility of distancing based on known population
epidemiological data about SARS-CoV-2 and concluding that unless the costs of
distancing vastly exceed the cost of illness per unit time, social distancing
remains a dominant strategy. These findings can assist in solidly anchoring
public health responses based on social distancing within a quantitative
framework attesting to their effectiveness.


  The COVID-19 pandemic has posed a policy making crisis where efforts to slow
down or end the pandemic conflict with economic priorities. This paper provides
mathematical analysis of optimal disease control policies with idealized
compartmental models for disease propagation and simplistic models of social
and economic costs. Two locally optimal control strategies are found and
categorized as `suppression' and `mitigation' strategies. We analyze how these
strategies change when we incorporate vaccination into the model and find a new
optimal `delay-mitigation' strategy.


  Evolutionary Game Theory (EGT) simulations are used to model populations
undergoing biological and cultural evolution in a range of fields, from biology
to economics to linguistics. In this paper we present DyPy, an open source
Python package that can perform evolutionary simulations for any matrix form
game for three common evolutionary dynamics: Moran, Wright-Fisher and
Replicator. We discuss the basic components of this package and illustrate how
it can be used to run a variety of simulations. Our package allows a user to
run such simulations fairly easily without much prior Python knowledge. We hope
that this will be a great asset to researchers in a number of different fields.


  COVID-19 is now the largest pandemic crisis of this century, with over 16
million registered cases worldwide. African countries have now begun
registering an increasing number of cases, yet, not many models developed focus
in specific African countries. In our study we use a simple SEIR model to
evaluate and predict future scenarios regarding the pandemic crisis in
Mozambique. We compare the effect of different policies on the infection curve
and estimate epidemiological parameters such as the current infection
reproduction number Rt and the growth rate g. We have found a low value for Rt,
ranging from 1.11 to 1.48 and a positive growth rate, between g = 0.22 to 0.27.
Our simulations also suggest that a lockdown shows potential for reducing the
infection peak height in 28%, on average, ranging from 20 to 36%.


  This study analyzed the morbidity and mortality rates of the COVID-19
pandemic in different prefectures of Japan. Under the constraint that daily
maximum confirmed deaths and daily maximum cases should exceed 4 and 10,
respectively, 14 prefectures were included, and cofactors affecting the
morbidity and mortality rates were evaluated. In particular, the number of
confirmed deaths was assessed excluding the cases of nosocomial infections and
nursing home patients. A mild correlation was observed between morbidity rate
and population density (R2=0.394). In addition, the percentage of the elderly
per population was also found to be non-negligible. Among weather parameters,
the maximum temperature and absolute humidity averaged over the duration were
found to be in modest correlation with the morbidity and mortality rates,
excluding the cases of nosocomial infections. The lower morbidity and mortality
are observed for higher temperature and absolute humidity. Multivariate
analysis considering these factors showed that determination coefficients for
the spread, decay, and combined stages were 0.708, 0.785, and 0.615,
respectively. These findings could be useful for intervention planning during
future pandemics, including a potential second COVID-19 outbreak.


  Lockdown procedures have been proven successful in mitigating the spread of
the viruses in this COVID-19 pandemic, but they also have devastating impact on
the economy. We use a modified Susceptible-Infectious-Recovered-Deceased model
with time dependent infection rate to simulate how the infection is spread
under lockdown. The economic cost due to the loss of workforce and incurred
medical expenses is evaluated with a simple model. We find the best strategy,
meaning the smallest economic cost for the entire course of the pandemic, is to
keep the strict lockdown as long as possible.


  We propose a robust parameter estimation method for dynamical systems based
on Statistical Learning techniques which aims to estimate a set of parameters
that well fit the dynamics in order to obtain robust evidences about the
qualitative behaviour of its trajectory. The method is quite general and
flexible, since it does not rely on any specific property of the dynamical
system, and represents a reinterpretation of Approximate Bayesian Computation
methods through the lens of Statistical Learning. The method is specially
useful for estimating parameters in epidemiological compartmental models in
order to obtain qualitative properties of a disease evolution. We apply it to
simulated and real data about COVID-19 spread in the US in order to evaluate
qualitatively its evolution over time, showing how one may assess the
effectiveness of measures implemented to slow the spread and some qualitative
features of the disease current and future evolution.


  Background: Northern Italy was one of the most impacted areas by COVID. It is
now widely assumed that the virus was silently spreading for at least 2 weeks
before the first patient was identified. During this silent phase, and in the
following weeks when the hospital system was overburdened, data collection was
not performed in an accurate enough way to estimate an epidemic curve. With the
aim of assessing both the dynamics of the introduction of the virus and the
effectiveness of containment measures introduced, we try to reconstruct the
epidemic curve using all cause mortality data. Methods: we collected all cause
mortality data stratified by age from the national institute of statistics,
together with COVID-related deaths data released by other government
structures. Using a SEIR model together with estimates of the exposure to death
time distribution, we fitted the reproduction number in different phases of the
spread at regional level. Results: We estimate a reproduction number of
2.6+/-0.1 before case 1 was identified. School closures in Lombardy lowered it
to 1.3. Soft lockdown measures resulted in R<0.8 and no further reductions were
observed when a hard lockdown was introduced (e.g. Emilia-Romagna soft lockdown
0.67 +/-0.07, hard lockdown 0.69+/-0.071). Reproduction number for the >75 age
range during hard lockdown are consistently higher than for the rest of the
population (e.g. 0.98 vs 0.71 in Milan province), suggesting outbreaks in
retirement facilities. Reproduction numbers in Bergamo and Brescia provinces
starting from March 7th are markedly lower than in other areas with the same
strict lockdown measures (Nearby provinces: 0.73, Brescia: 0.52, Bergamo 0.43)
supporting the hypothesis that in those provinces a large percentage of the
population had already been infected by the beginning of March.


  Coronaviruses are enveloped, non-segmented positive-sense RNA viruses that
have the largest genome among RNA viruses. The genome contains a large
replicase ORF encodes nonstructural proteins (NSPs), structural and accessory
genes. NSP15 is a nidoviral RNA uridylate-specific endoribonuclease (NendoU)
has C-terminal catalytic domain. The endoribonuclease activity of NSP15
interferes with the innate immune response of the host. Here, we screened
Selleckchem Natural product database of compounds against the NSP15,
Thymopentin and Oleuropein showed highest binding energies. The binding of
these molecules was further validated by Molecular dynamic simulation and found
very stable complexes. These drugs might serve as effective counter molecules
in the reduction of virulence of this virus. Future validation of both these
inhibitors are worth consideration for patients being treated for COVID -19.


  We present a compartmental meta-population model for the spread of Covid-19
in India. Our model simulates populations at a district or state level using an
epidemiological model that is appropriate to Covid-19. Different districts are
connected by a transportation matrix developed using available census data. We
introduce uncertainties in the testing rates into the model that takes into
account the disparate responses of the different states to the epidemic and
also factors in the state of the public healthcare system. Our model allows us
to generate qualitative projections of Covid-19 spread in India, and further
allows us to investigate the effects of different proposed interventions. By
building in heterogeneity at geographical and infrastructural levels and in
local responses, our model aims to capture some of the complexity of
epidemiological modeling appropriate to a diverse country such as India.


  Whether aneuploid cells with diverse karyotypes have any properties in common
has a been a subject of intense interest. A recent study by Terhorst et al. (1)
reinvestigated the common aneuploidy gene expression (CAGE), disputing the
conclusion of our recent work (2). In this short article, which has been
submitted to PNAS as a Letter to the Editor, we explain our major concerns
about Terhorst et al. and why we believe that our previous conclusion stands
valid.


  We introduce a theoretical framework that highlights the impact of physical
distancing variables such as human mobility and physical proximity on the
evolution of epidemics and, crucially, on the reproduction number. In
particular, in response to the coronavirus disease (CoViD-19) pandemic,
countries have introduced various levels of 'lockdown' to reduce the number of
new infections. Specifically we use a collisional approach to an infection-age
structured model described by a renewal equation for the time homogeneous
evolution of epidemics. As a result, we show how various contributions of the
lockdown policies, namely physical proximity and human mobility, reduce the
impact of SARS-CoV-2 and mitigate the risk of disease resurgence. We check our
theoretical framework using real-world data on physical distancing with two
different data repositories, obtaining consistent results. Finally, we propose
an equation for the effective reproduction number which takes into account
types of interactions among people, which may help policy makers to improve
remote-working organizational structure.


  Molecular transport of biomolecules plays a pivotal role in the machinery of
life. Yet, this role is poorly understood due the lack of quantitative
information. Here, the role and properties of the C-terminal region of
Escherichia coli Hfq is reported, involved in controlling the flow of a DNA
solution. A combination of experimental methodologies has been used to probe
the interaction of Hfq with DNA and to measure the rheological properties of
the complex. A physical gel with a temperature reversible elasticity modulus is
formed due to formation of non-covalent crosslinks. The mechanical response of
the complexes shows that they are inhomogeneous soft solids. Our experiments
indicate that Hfq C-terminal region could contribute to genome mechanical
response. The reported viscoelasticity of the DNA-protein complex might have
implications for cellular pro-cesses involving molecular transport of DNA or
segments thereof.


  The Colless index for bifurcating phylogenetic trees, introduced by Colless
(1982), is defined as the sum, over all internal nodes $v$ of the tree, of the
absolute value of the difference of the sizes of the clades defined by the
children of $v$. It is one of the most popular phylogenetic balance indices,
because, in addition to measuring the balance of a tree in a very simple and
intuitive way, it turns out to be one of the most powerful and discriminating
phylogenetic shape indices. But it has some drawbacks. On the one hand,
although its minimum value is reached at the so-called maximally balanced
trees, it is almost always reached also at trees that are not maximally
balanced. On the other hand, its definition as a sum of absolute values of
differences makes it difficult to study analytically its distribution under
probabilistic models of bifurcating phylogenetic trees. In this paper we show
that if we replace in its definition the absolute values of the differences of
clade sizes by the squares of these differences, all these drawbacks are
overcome and the resulting index is still more powerful and discriminating than
the original Colless index.


  Energy-based modelling brings engineering insight to the understanding of
biomolecular systems. It is shown how well-established control engineering
concepts, such as loop-gain, arise from energy feedback loops and are therefore
amenable to control engineering insight. In particular, a novel method is
introduced to allow the transfer function based approach of classical linear
control to be utilised in the analysis of feedback systems modelled by network
thermodynamics and thus amalgamate energy-based modelling with control systems
analysis. The approach is illustrated using a class of metabolic cycles with
activation and inhibition leading the concept of Cyclic Flow Modulation.


  Disease and information spread over social and information networks.
Understanding the spread phenomena in networks requires paying attention not
only to the degree distribution but also to the degree correlation. However, it
is considered difficult to analytically deal with the effect of degree
correlation on spread phenomena. Here, we introduce degree correlation using a
simple method and present the theoretical formulas of the outbreak threshold
and basic reproduction number. We theoretically clarify the effect of the
degree correlation.


  We consider a global (location independent) model of pandemic growth which
generalizes the SIR model to accommodate important features of the COVID-19
pandemic, notably the implementation of pandemic reduction measures. This
"SHIR" model is applied to COVID-19 data, and shows promise as a simple,
tractable formalism with few parameters that can be used to model pandemic case
numbers. As an example we show that the average time dependence of new COVID-19
cases per day from 15 Central and Western European countries is in good
agreement with the analytic, parameter-free prediction of the model


  Viral transmission pathways have profound implications for public safety; it
is thus imperative to establish a complete understanding of viable infectious
avenues. Mounting evidence suggests SARS-CoV-2 can be transmitted via the air;
however, this has not yet been demonstrated. Here we quantitatively analyze
virion accumulation by accounting for aerosolized virion emission and
destabilization. Reported superspreading events analyzed within this framework
point towards aerosol mediated transmission of SARS-CoV-2. Virion exposure
calculated for these events is found to trace out a single value, suggesting a
universal minimum infective dose (MID) via aerosol that is comparable to the
MIDs measured for other respiratory viruses; thus, the consistent infectious
exposure levels and their commensurability to known aerosol-MIDs establishes
the plausibility of aerosol transmission of SARS-CoV-2. Using filtration at a
rate exceeding the destabilization rate of aerosolized SARS-CoV-2 can reduce
exposure below this infective dose.


  The nucleotide sequence representation of DNA can be inadequate for resolving
protein-DNA binding sites and regulatory substrates, such as those involved in
gene expression and horizontal gene transfer. Considering that sequence-like
representations are algorithmically very useful, here we fused over 60
currently available DNA physicochemical and conformational variables into
compact structural representations that can encode single DNA binding sites to
whole regulatory regions. We find that the main structural components reflect
key properties of protein-DNA interactions and can be condensed to the amount
of information found in a single nucleotide position. The most accurate
structural representations compress functional DNA sequence variants by 30% to
50%, as each instance encodes from tens to thousands of sequences. We show that
a structural distance function discriminates among groups of DNA substrates
more accurately than nucleotide sequence-based metrics. As this opens up a
variety of implementation possibilities, we develop and test a distance-based
alignment algorithm, demonstrating the potential of using the structural
representations to enhance sequence-based algorithms. Due to the bias of most
current bioinformatic methods to nucleotide sequence representations, it is
possible that considerable performance increases might still be achievable with
such solutions.


  Interactions among individuals in natural populations often occur in a
dynamically changing environment. Understanding the role of environmental
variation in population dynamics has long been a central topic in theoretical
ecology and population biology. However, the key question of how individuals,
in the middle of challenging social dilemmas (e.g., the "tragedy of the
commons"), modulate their behaviors to adapt to the fluctuation of the
environment has not yet been addressed satisfactorily. Utilizing evolutionary
game theory and stochastic games, we develop a game-theoretical framework that
incorporates the adaptive mechanism of reinforcement learning to investigate
whether cooperative behaviors can evolve in the ever-changing group interaction
environment. When the action choices of players are just slightly influenced by
past reinforcements, we construct an analytical condition to determine whether
cooperation can be favored over defection. Intuitively, this condition reveals
why and how the environment can mediate cooperative dilemmas. Under our model
architecture, we also compare this learning mechanism with two non-learning
decision rules, and we find that learning significantly improves the propensity
for cooperation in weak social dilemmas, and, in sharp contrast, hinders
cooperation in strong social dilemmas. Our results suggest that in complex
social-ecological dilemmas, learning enables the adaptation of individuals to
varying environments.


  To shorten the time required to find effective new drugs, like antivirals, a
key parameter to consider is membrane permeability, as a compound intended for
an intracellular target with poor permeability will have low efficacy. Here, we
present a computational model that considers both drug characteristics and
membrane properties for the rapid assessment of drugs permeability through the
coronavirus envelope and various cellular membranes. We analyze 79 drugs that
are considered as potential candidates for the treatment of SARS-CoV-2 and
determine their time of permeation in different organelle membranes grouped by
viral baits and mammalian processes. The computational results are correlated
with experimental data, present in the literature, on bioavailability of the
drugs, showing a negative correlation between fast permeation and most
promising drugs. This model represents an important tool capable of evaluating
how permeability affects the ability of compounds to reach both intended and
unintended intracellular targets in an accurate and rapid way. The method is
general and flexible and can be employed for a variety of molecules, from small
drugs to nanoparticles, as well to a variety of biological membranes.


  Containment measures implemented by some countries to suppress the spread of
COVID-19 have resulted in a slowdown of the epidemic characterized by time
series of daily infections plateauing over extended periods of time. We prove
that such a dynamical pattern is compatible with critical
Susceptible-Infected-Removed (SIR) dynamics. In traditional analyses of the
critical SIR model, the critical dynamical regime is started from a single
infected node. The application of containment measures to an ongoing epidemic,
however, has the effect to make the system enter in its critical regime with a
number of infected individuals potentially large. We describe how such
non-trivial starting conditions affect the critical behavior of the SIR model.
We perform a theoretical and large-scale numerical investigation of the model.
We show that the expected outbreak size is an increasing function of the
initial number of infected individuals, while the expected duration of the
outbreak is a non-monotonic function of the initial number of infected
individuals. Also, we precisely characterize the magnitude of the fluctuations
associated with the size and duration of the outbreak in critical SIR dynamics
with non-trivial initial conditions. Far from heard immunity, fluctuations are
much larger than average values, thus indicating that predictions of plateauing
time series may be particularly challenging.


  An activated process consists of energy activation and barrier crossing; the
former is a prerequisite for the latter. Barrier crossing has been studied
extensively, but energy activation has been overlooked due to a lack of means
to gauge its progress. We define reaction stability as the probability that
reactive trajectories pass a vicinity in phase space; it enabled us to analyze
energy activation of a biomolecular isomerization. This process follows a
mechanism fundamentally different from presumed mechanisms in standard reaction
rate theories: it features accumulation of high kinetic energy in reaction
coordinates, achieved by precise synergy between them coordinated by momentum
space.


  Ever since the outbreak of the COVID-19 epidemic, various public health
control strategies have been proposed and tested against the coronavirus
SARS-CoV-2. We study three specific COVID-19 epidemic control models: the
susceptible, exposed, infectious, recovered (SEIR) model with vaccination
control; the SEIR model with shield immunity control; and the susceptible,
un-quarantined infected, quarantined infected, confirmed infected (SUQC) model
with quarantine control. We express the control requirement in metric temporal
logic (MTL) formulas (a type of formal specification languages) which can
specify the expected control outcomes such as "the deaths from the infection
should never exceed one thousand per day within the next three months" or "the
population immune from the disease should eventually exceed 200 thousand within
the next 100 to 120 days". We then develop methods for synthesizing control
strategies with MTL specifications. To the best of our knowledge, this is the
first paper to systematically synthesize control strategies based on the
COVID-19 epidemic models with formal specifications. We provide simulation
results in three different case studies: vaccination control for the COVID-19
epidemic with model parameters estimated from data in Lombardy, Italy; shield
immunity control for the COVID-19 epidemic with model parameters estimated from
data in Lombardy, Italy; and quarantine control for the COVID-19 epidemic with
model parameters estimated from data in Wuhan, China. The results show that the
proposed synthesis approach can generate control inputs such that the
time-varying numbers of individuals in each category (e.g., infectious, immune)
satisfy the MTL specifications. The results also show that early intervention
is essential in mitigating the spread of COVID-19, and more control effort is
needed for more stringent MTL specifications.


  The ongoing SARS-CoV-2 pandemic has been holding the world hostage for more
than a year now. Mobility is key to viral spreading and its restriction is the
main non-pharmaceutical interventions to fight the virus expansion. Previous
works have shown a connection between the structural organization of cities and
the movement patterns of their residents. This puts urban centers in the focus
of epidemic surveillance and interventions. Here we show that the organization
of urban flows has a tremendous impact on disease spreading and on the
amenability of different mitigation strategies. By studying anonymous and
aggregated intra-urban flows in a variety of cities in the United States and
other countries, and a combination of empirical analysis and analytical
methods, we demonstrate that the response of cities to epidemic spreading can
be roughly classified in two major types according to the overall organization
of those flows. Hierarchical cities, where flows are concentrated primarily
between mobility hotspots, are particularly vulnerable to the rapid spread of
epidemics. Nevertheless, mobility restrictions in such types of cities are very
effective in mitigating the spread of a virus. Conversely, in sprawled cities
which present many centers of activity, the spread of an epidemic is much
slower, but the response to mobility restrictions is much weaker and less
effective. Investing resources on early monitoring and prompt ad-hoc
interventions in more vulnerable cities may prove helpful in containing and
reducing the impact of future pandemics.


  This paper is based on the observation that, during Covid-19 epidemic, the
choice of which individuals should be tested has an important impact on the
effectiveness of selective confinement measures. This decision problem is
closely related to the problem of optimal sensor selection, which is a very
active research subject in control engineering. The goal of this paper is to
propose a policy to smartly select the individuals to be tested. The main idea
is to model the epidemics as a stochastic dynamic system and to select the
individual to be tested accordingly to some optimality criteria, e.g. to
minimize the probability of undetected asymptomatic cases. Every day, the
probability of infection of the different individuals is updated making use of
the stochastic model of the phenomenon and of the information collected in the
previous days. Simulations for a closed community of 10000 individuals show
that the proposed technique, coupled with a selective confinement policy, can
reduce the spread of the disease while limiting the number of individuals
confined if compared to the simple contact tracing of positive and to an
off-line test selection strategy based on the number of contacts.


  Chemical Reaction Networks (CRNs) provide a useful abstraction of molecular
interaction networks in which molecular structures as well as mass conservation
principles are abstracted away to focus on the main dynamical properties of the
network structure. In their interpretation by ordinary differential equations,
we say that a CRN with distinguished input and output species computes a
positive real function $f : R+ $\rightarrow$ R+$, if for any initial
concentration x of the input species, the concentration of the output molecular
species stabilizes at concentration f (x). The Turing-completeness of that
notion of chemical analog computation has been established by proving that any
computable real function can be computed by a CRN over a finite set of
molecular species. Rate-independent CRNs form a restricted class of CRNs of
high practical value since they enjoy a form of absolute robustness in the
sense that the result is completely independent of the reaction rates and
depends solely on the input concentrations. The functions computed by
rate-independent CRNs have been characterized mathematically as the set of
piecewise linear functions from input species. However, this does not provide a
mean to decide whether a given CRN is rate-independent. In this paper, we
provide graphical conditions on the Petri Net structure of a CRN which entail
the rate-independence property either for all species or for some output
species. We show that in the curated part of the Biomodels repository, among
the 590 reaction models tested, 2 reaction graphs were found to satisfy our
rate-independence conditions for all species, 94 for some output species, among
which 29 for some non-trivial output species. Our graphical conditions are
based on a non-standard use of the Petri net notions of place-invariants and
siphons which are computed by constraint programming techniques for efficiency
reasons.


  A number of epidemics, including the SARS-CoV-1 epidemic of 2002-2004, have
been known to exhibit superspreading, in which a small fraction of infected
individuals is responsible for the majority of new infections. The existence of
superspreading implies a fat-tailed distribution of infectiousness (new
secondary infections caused per day) among different individuals. Here, we
present a simple method to estimate the variation in infectiousness by
examining the variation in early-time growth rates of new cases among different
subpopulations. We use this method to estimate the mean and variance in the
infectiousness, $\beta$, for SARS-CoV-2 transmission during the early stages of
the pandemic within the United States. We find that $\sigma_\beta/\mu_\beta
\gtrsim 3.2$, where $\mu_\beta$ is the mean infectiousness and $\sigma_\beta$
its standard deviation, which implies pervasive superspreading. This result
allows us to estimate that in the early stages of the pandemic in the USA, over
81% of new cases were a result of the top 10% of most infectious individuals.


  The main goal of this work is to estimate the actual number of cases of
Covid-19 in Spain in the period 01-31-2020 / 06-01-2020 by Autonomous
Communities. Based on these estimates, this work allows us to accurately
re-estimate the lethality of the disease in Spain, taking into account
unreported cases. A hierarchical Bayesian model recently proposed in the
literature has been adapted to model the actual number of Covid-19 cases in
Spain. The results of this work show that the real load of Covid-19 in Spain in
the period considered is well above the data registered by the public health
system. Specifically, the model estimates show that, cumulatively until June
1st, 2020, there were 2,425,930 cases of Covid-19 in Spain with characteristics
similar to those reported (95\% credibility interval: 2,148,261 - 2,813,864),
from which were actually registered only 518,664. Considering the results
obtained from the second wave of the Spanish seroprevalence study, which
estimates 2,350,324 cases of Covid-19 produced in Spain, in the period of time
considered, it can be seen that the estimates provided by the model are quite
good. This work clearly shows the key importance of having good quality data to
optimize decision-making in the critical context of dealing with a pandemic.


  Single-molecule experimental techniques track the real-time dynamics of
molecules by recording a small number of experimental observables. Following
these observables provides a coarse-grained, low-dimensional representation of
the conformational dynamics but does not furnish an atomistic representation of
the instantaneous molecular structure. Takens' Delay Embedding Theorem asserts
that, under quite general conditions, these low-dimensional time series can
contain sufficient information to reconstruct the full molecular configuration
of the system up to an a priori unknown transformation. By combining Takens'
Theorem with tools from statistical thermodynamics, manifold learning,
artificial neural networks, and rigid graph theory, we establish an approach
Single-molecule TAkens Reconstruction (STAR) to learn this transformation and
reconstruct molecular configurations from time series in
experimentally-measurable observables such as intramolecular distances
accessible to single molecule F\"orster resonance energy transfer. We
demonstrate the approach in applications to molecular dynamics simulations of a
C24H50 polymer chain and the artificial mini-protein Chignolin. The trained
models reconstruct molecular configurations from synthetic time series data in
the head-to-tail molecular distances with atomistic root mean squared deviation
accuracies better than 0.2 nm. This work demonstrates that it is possible to
accurately reconstruct protein structures from time series in
experimentally-measurable observables and establishes the theoretical and
algorithmic foundations to do so in applications to real experimental data.


  Whole-cell modelling is constrained by the laws of nature in general and the
laws of thermodynamics in particular. This paper shows how one prolific source
of information, stoichiometric models of biomolecular systems, can be
integrated with thermodynamic principles using the bond graph approach to
network thermodynamics.


  This paper is concerned with the robustness of the sustained oscillations
predicted by an epidemic ODE model defined on contact networks. The model
incorporates the spread of awareness among individuals and, moreover, a small
inflow of imported cases. These cases prevent stochastic extinctions when we
simulate the epidemics and, hence, they allow to check whether the average
dynamics for the fraction of infected individuals are accurately predicted by
the ODE model. Stochastic simulations confirm the existence of sustained
oscillations for different types of random networks, with a sharp transition
from a non-oscillatory asymptotic regime to a periodic one as the alerting rate
of susceptible individuals increases from very small values. This abrupt
transition to periodic epidemics of high amplitude is quite accurately
predicted by the Hopf-bifurcation curve computed from the ODE model using the
alerting rate and the infection transmission rate for aware individuals as
tuning parameters.


  The global SARS-CoV-2 pandemic suggests a novel type of disease spread
dynamics. WHO states that there is currently no evidence that people who have
recovered from COVID-19 and have antibodies are immune from a second infection
[WHO]. Conventional mathematical models consider cases for which a recovered
individual either becomes susceptible again or develops an immunity. Here, we
study the case where infected agents recover and only develop immunity if they
are continuously infected for some time. Otherwise, they become susceptible
again. We show that field theory bounds the peak of the infectious rate.
Consequently, the theory's phases characterise the disease dynamics: (i) a
pandemic phase and (ii) a response regime. The model excellently describes the
epidemic spread of the SARS-CoV-2 outbreak in the city of Wuhan, China. We find
that only 30% of the recovered agents have developed an immunity. We anticipate
our paper to influence the decision making upon balancing the economic impact
and the pandemic impact on society. As long as disease controlling measures
keep the disease dynamics in the "response regime", a pandemic escalation
('second wave') is ruled out.


  Modern computational organic chemistry is becoming increasingly data-driven.
There remain a large number of important unsolved problems in this area such as
product prediction given reactants, drug discovery, and metric-optimized
molecule synthesis, but efforts to solve these problems using machine learning
have also increased in recent years. In this work, we propose the use of
attention to study functional groups and other property-impacting molecular
substructures from a data-driven perspective, using a transformer-based model
(BERT) on datasets of string representations of molecules and analyzing the
behavior of its attention heads. We then apply the representations of
functional groups and atoms learned by the model to tackle problems of
toxicity, solubility, drug-likeness, and synthesis accessibility on smaller
datasets using the learned representations as features for graph convolution
and attention models on the graph structure of molecules, as well as
fine-tuning of BERT. Finally, we propose the use of attention visualization as
a helpful tool for chemistry practitioners and students to quickly identify
important substructures in various chemical properties.


  The exact analytical solution in closed form of a modified SIR system where
recovered individuals are removed from the population is presented. In this
dynamical system the populations $S(t)$ and $R(t)$ of susceptible and recovered
individuals are found to be generalized logistic functions, while infective
ones $I(t)$ are given by a generalized logistic function times an exponential,
all of them with the same characteristic time. The dynamics of this modified
SIR system is analyzed and the exact computation of some epidemiologically
relevant quantities is performed. The main differences between this modified
SIR model and original SIR one are presented and explained in terms of the
zeroes of their respective conserved quantities. Moreover, it is shown that the
modified SIR model with time-dependent transmission rate can be also solved in
closed form for certain realistic transmission rate functions.


  Without vaccines and treatments, societies must rely on non-pharmaceutical
intervention strategies to control the spread of emerging diseases such as
COVID-19. Though complete lockdown is epidemiologically effective, because it
eliminates infectious contacts, it comes with significant costs. Several recent
studies have suggested that a plausible compromise strategy for minimizing
epidemic risk is periodic closure, in which populations oscillate between
wide-spread social restrictions and relaxation. However, no underlying theory
has been proposed to predict and explain optimal closure periods as a function
of epidemiological and social parameters. In this work we develop such an
analytical theory for SEIR-like model diseases, showing how characteristic
closure periods emerge that minimize the total outbreak, and increase
predictably with the reproductive number and incubation periods of a disease,
as long as both are within predictable limits. Using our approach we
demonstrate a sweet-spot effect in which optimal periodic closure is maximally
effective for diseases with similar incubation and recovery periods. Our
results compare well to numerical simulations, including in COVID-19 models
where infectivity and recovery show significant variability.


  Computational protein structure determination involves optimization in a
problem space much too large to exhaustively search. Existing approaches
include optimization algorithms such as gradient descent and simulated
annealing, but these typically only find local minima. One novel approach
implemented in REDcRAFT is to instead of folding a protein all at the same
time, fold it residue by residue. This simulates a protein folding as each
residue exits from the generating ribosome. While REDcRAFT exponentially
reduces the problem space so it can be explored in polynomial time, it is still
extremely computationally demanding. This algorithm does have the advantage
that most of the execution time is spent in inherently parallelizable code.
However, preliminary results from parallel execution indicate that
approximately two-thirds of execution time is dedicated to system overhead.
Additionally, by carefully analyzing and timing the structure of the program
the major bottlenecks can be identified. After addressing these issues,
REDcRAFT becomes a scalable parallel application with nearly two orders of
magnitude improvement.


  A novel coronavirus disease (COVID-19) is sweeping the world and has taken
away thousands of lives. As the current epicenter, the United States has the
largest number of confirmed and death cases of COVID-19. Meteorological factors
have been found associated with many respiratory diseases in the past studies.
In order to understand that how and during which period of time do the
meteorological factors have the strongest association with the transmission and
fatality of COVID-19, we analyze the correlation between each meteorological
factor during different time periods within the incubation window and the
confirmation and fatality rate, and develop statistic models to quantify the
effects at county level. Results show that meteorological variables except
maximum wind speed during the day 13 - 0 before current day shows the most
significant correlation (P < 0.05) with the daily confirmed rate, while
temperature during the day 13 - 8 before are most significantly correlated (P <
0.05) with the daily fatality rate. Temperature is the only meteorological
factor showing dramatic positive association nationally, particularly in the
southeast US where the current outbreak most intensive. The influence of
temperature is remarkable on the confirmed rate with an increase of over 5 pmp
in many counties, but not as much on the fatality rate (mostly within 0.01%).
Findings in this study will help understanding the role of meteorological
factors in the spreading of COVID-19 and provide insights for public and
individual in fighting against this global epidemic.


  We provide a probabilistic SIRD model for the COVID-19 pandemic in Italy,
where we allow the infection, recovery and death rates to be random. In
particular, the underlying random factor is driven by a fractional Brownian
motion. Our model is simple and needs only some few parameters to be
calibrated.


  We study a SEIR model considered by Gomes et al. \cite{Gomes2020} and Aguas
et al. \cite{Aguas2020} where different individuals are assumed to have
different levels of susceptibility or exposure to infection. Under this
heterogeneity assumption, epidemic growth is effectively suppressed when the
percentage of population having acquired immunity surpasses a critical level -
the herd immunity threshold - that is lower than in homogeneous populations. We
find explicit formulas to calculate herd immunity thresholds and stable
configuration, and explore extensions of the model.


  (1) Background: RNA viruses and especially coronaviruses could act inside
host cells not only by building their own proteins, but also by perturbing the
cell metabolism. We show the possibility of miRNA-like inhibitions by the
SARS-CoV-2 concerning for example the hemoglobin and type I interferons
syntheses, hence highly perturbing oxygen distribution in vital organs and
immune response as described by clinicians; (2) Methods: We compare RNA
subsequences of SARS-CoV-2 protein S and RNA-dependent RNA polymerase genes to
mRNA sequences of beta-globin and type I interferons; (3) Results: RNA
subsequences longer than eight nucleotides from SARS-CoV-2 genome could
hybridize subsequences of the mRNA of beta-globin and of type I interferons;
(4) Conclusions: Beyond viral protein production, Covid-19 might affect vital
processes like host oxygen transport and immune response.


  As organisms are faced with intense rapidly changing selective pressures, new
genetic material is required to facilitate adaptation. Among sources of genetic
novelty, gene duplications and transposable elements (TEs) offer new genes or
new regulatory patterns that can facilitate evolutionary change. With advances
in genome sequencing it is possible to gain a broader view of how gene family
proliferation and TE content evolve in non-model species when populations
become threatened. Freshwater bivalves (Unionidae) currently face severe
anthropogenic challenges. Over 70% of species in the United States are
threatened, endangered or extinct due to pollution, damming of waterways, and
overfishing. We have created a reference genome for M. nervosa to determine how
genome content has evolved in the face of these widespread environmental
challenges. We observe a burst of recent transposable element proliferation
causing a 382 Mb expansion in genome content. Gene family expansion is common,
with a duplication rate of 1.16 x 10^-8 per gene per generation. Cytochrome
P450, ABC transporters, Hsp70 genes, von Willebrand proteins, chitin metabolism
genes, mitochondria eating proteins, and opsin gene families have experienced
significantly greater amplification and show signatures of selection. We use
evolutionary theory to assess the relative contribution of SNPs and
duplications in evolutionary change. Estimates suggest that gene family
evolution may offer an exceptional substrate of genetic variation in M.
nervosa, with Psgv=0.185 compared with Psgv=0.067 for single nucleotide
changes. Hence, we suggest that gene family evolution is a source of "hopeful
monsters" within the genome that facilitate adaptation.


  Repeated interaction promotes cooperation among rational individuals under
the shadow of future, but it is hard to maintain cooperation when a large
number of error-prone individuals are involved. One way to construct a
cooperative Nash equilibrium is to find a `friendly-rivalry' strategy, which
aims at full cooperation but never allows the co-players to be better off.
Recently it has been shown that for the iterated Prisoner's Dilemma in the
presence of error, a friendly rival can be designed with the following five
rules: Cooperate if everyone did, accept punishment for your own mistake,
punish defection, recover cooperation if you find a chance, and defect in all
the other circumstances. In this work, we construct such a friendly-rivalry
strategy for the iterated $n$-person public-goods game by generalizing those
five rules. The resulting strategy makes a decision with referring to the
previous $m=2n-1$ rounds. A friendly-rivalry strategy for $n=2$ inherently has
evolutionary robustness in the sense that no mutant strategy has higher
fixation probability in this population than that of a neutral mutant. Our
evolutionary simulation indeed shows excellent performance of the proposed
strategy in a broad range of environmental conditions when $n= 2$ and $3$.


  Motivation: One of the branches of Systems Biology is focused on a deep
understanding of underlying regulatory networks through the analysis of the
biomolecules oscillations and their interplay. Synthetic Biology exploits gene
or/and protein regulatory networks towards the design of oscillatory networks
for producing useful compounds. Therefore, at different levels of application
and for different purposes, the study of biomolecular oscillations can lead to
different clues about the mechanisms underlying living cells. It is known that
network-level interactions involve more than one type of biomolecule as well as
biological processes operating at multiple omic levels. Combining
network/pathway-level information with genetic information it is possible to
describe well-understood or unknown bacterial mechanisms and organism-specific
dynamics. Results: Network multi-omic integration has led to the discovery of
interesting oscillatory signals. Following the methodologies used in signal
processing and communication engineering, a new methodology is introduced to
identify and quantify the extent of the multi-omic oscillations of the signal.
New signal metrics are designed to allow further biotechnological explanations
and provide important clues about the oscillatory nature of the pathways and
their regulatory circuits. Our algorithms designed for the analysis of
multi-omic signals are tested and validated on 11 different bacteria for
thousands of multi-omic signals perturbed at the network level by different
experimental conditions. Information on the order of genes, codon usage, gene
expression, and protein molecular weight is integrated at three different
functional levels. Oscillations show interesting evidence that network-level
multi-omic signals present a synchronized response to perturbations and
evolutionary relations along with taxa.


  We study and predict the evolution of Covid-19 in six US states from the
period May 1 through August 31 using a discrete compartment-based model and
prescribe active intervention policies, like lockdowns, on the basis of
minimizing a loss function, within the broad framework of partially observed
Markov decision processes. For each state, Covid-19 data for 40 days (starting
from May 1 for two northern states and June 1 for four southern states) are
analyzed to estimate the transition probabilities between compartments and
other parameters associated with the evolution of the epidemic. These
quantities are then used to predict the course of the epidemic in the given
state for the next 50 days (test period) under various policy allocations,
leading to different values of the loss function over the training horizon. The
optimal policy allocation is the one corresponding to the smallest loss. Our
analysis shows that none of the six states need lockdowns over the test period,
though the no lockdown prescription is to be interpreted with caution:
responsible mask use and social distancing of course need to be continued. The
caveats involved in modeling epidemic propagation of this sort are discussed at
length. A sketch of a non-Markovian formulation of Covid-19 propagation (and
more general epidemic propagation) is presented as an attractive avenue for
future research in this area.


  Empirical studies show that epidemiological models based on an epidemic's
initial spread rate often fail to predict the true scale of that epidemic. Most
epidemics with a rapid early rise die out before affecting a significant
fraction of the population, whereas the early pace of some pandemics is rather
modest. Recent models suggest that this could be due to the heterogeneity of
the target population's susceptibility. We study a computer malware ecosystem
exhibiting spread mechanisms resembling those of biological systems while
offering details unavailable for human epidemics. Rather than comparing models,
we directly estimate reach from a new and vastly more complete data from a
parallel domain, that offers superior details and insight as concerns
biological outbreaks. We find a highly heterogeneous distribution of computer
susceptibilities, with nearly all outbreaks initially over-affecting the tail
of the distribution, then collapsing quickly once this tail is depleted. This
mechanism restricts the correlation between an epidemic's initial growth rate
and its total reach, thus preventing the majority of epidemics, including
initially fast-growing outbreaks, from reaching a macroscopic fraction of the
population. The few pervasive malwares distinguish themselves early on via the
following key trait: they avoid infecting the tail, while preferentially
targeting computers unaffected by typical malware.


  We generalize the Susceptible-Infected-Removed model for epidemics to take
into account generic effects of heterogeneity in the degree of susceptibility
to infection in the population. We introduce a single new parameter
corresponding to a power-law exponent of the susceptibility distribution that
characterizes the population heterogeneity. We show that our generalized model
is as simple as the original model which is contained as a limiting case.
Because of this simplicity, numerical solutions can be generated easily and key
properties of the epidemic wave can still be obtained exactly. In particular,
we present exact expressions for the herd immunity level, the final size of the
epidemic, as well as for the shape of the wave and for observables that can be
quantified during an epidemic. We find that in strongly heterogeneous
populations the epidemic reaches only a small fraction of the population. This
implies that the herd immunity level can be much lower than in commonly used
models with homogeneous populations. Using our model to analyze data for the
SARS-CoV-2 epidemic in Germany shows that the reported time course is
consistent with several scenarios characterized by different levels of
immunity. These scenarios differ in population heterogeneity and in the time
course of the infection rate, for example due to mitigation efforts or
seasonality. Our analysis reveals that quantifying the effects of mitigation
requires knowledge on the degree of heterogeneity in the population. Our work
shows that key effects of population heterogeneity can be captured without
increasing the complexity of the model. We show that information about
population heterogeneity will be key to understand how far an epidemic has
progressed and what can be expected for its future course.


  The gradual accumulation of damage and dysregulation during the aging of
living organisms can be quantified. Even so, the aging process is complex and
has multiple interacting physiological scales -- from the molecular to cellular
to whole tissues. In the face of this complexity, we can significantly advance
our understanding of aging with the use of computational models that simulate
realistic individual trajectories of health as well as mortality. To do so,
they must be systems-level models that incorporate interactions between
measurable aspects of age-associated changes. To incorporate individual
variability in the aging process, models must be stochastic. To be useful they
should also be predictive, and so must be fit or parameterized by data from
large populations of aging individuals. In this perspective, we outline where
we have been, where we are, and where we hope to go with such computational
models of aging. Our focus is on data-driven systems-level models, and on their
great potential in aging research.


  From patients and physicians perspectives, the clinical definition of a tumor
malignant phenotype could be restricted to the early diagnosis of sub-types of
malignancies with the increased risk of existing therapy failure and high
likelihood of death from cancer. It is the viewpoint from which the
understanding of malignant regulatory signatures is considered in this
contribution. Analyses from this perspective of experimental and clinical
observations revealed the pivotal role of human stem cell-associated retroviral
sequences (SCARS) in the origin and pathophysiology of clinically-lethal
malignancies. SCARS represent evolutionary- and biologically-related family of
genomic regulatory sequences, the principal physiological function of which is
to create and maintain the stemness phenotype during human preimplantation
embryogenesis. SCARS expression must be silenced during cellular
differentiation and SCARS activity remains silent in most
terminally-differentiated human cells performing specialized functions in the
human body. De-repression and sustained activation of SCARS result in
differentiation-defective phenotypes, tissue- and organ-specific clinical
manifestations of which are diagnosed as pathological conditions defined by a
consensus of pathomorphological, molecular, and genetic examinations as the
malignant growth. Contemporary evidence are presented that high-fidelity
molecular signals of continuing activities of SCARS in association with genomic
regulatory networks of thousands functionally-active enhancers triggering
engagements of down-stream genetic loci may serve as both reliable diagnostic
tools and druggable molecular targets readily amenable for diagnosis and
efficient therapeutic management of clinically-lethal malignancies.


  We present a compartmental mathematical model with demography for the spread
of the COVID-19 disease, considering also asymptomatic infectious individuals.
We compute the basic reproductive ratio of the model and study the local and
global stability for it. We solve the model numerically based on the case of
Italy. We propose a vaccination model and we derive threshold conditions for
preventing infection spread in the case of imperfect vaccines.


  This manuscript proposes a model of $n$ distinct populations interaction
structured SEIR to describe the spread of COVID-19 pandemic diseases. The
proposed model has the flexibility to include geographically separated
communities as well as taking into account aging population groups and their
interactions. We show that distinct assumptions on the dynamics of the proposed
model lead to a plateau-like curve of the infected population, reflecting
collected data from large countries like Brazil. Such observations aim to the
following conjecture: " Covid-19 diseased diffusion from the capitals to the
Brazil interior, as reflected by the proposed model, is responsible for
plateau-like reported cases in the country". We present numerical simulations
of some scenarios and comparison with reported data from Brazil that
corroborate with the aforementioned conclusions.


  Genome analysis fundamentally starts with a process known as read mapping,
where sequenced fragments of an organism's genome are compared against a
reference genome. Read mapping is currently a major bottleneck in the entire
genome analysis pipeline, because state-of-the-art genome sequencing
technologies are able to sequence a genome much faster than the computational
techniques employed to analyze the genome. We describe the ongoing journey in
significantly improving the performance of read mapping. We explain
state-of-the-art algorithmic methods and hardware-based acceleration
approaches. Algorithmic approaches exploit the structure of the genome as well
as the structure of the underlying hardware. Hardware-based acceleration
approaches exploit specialized microarchitectures or various execution
paradigms (e.g., processing inside or near memory). We conclude with the
challenges of adopting these hardware-accelerated read mappers.


  A mathematical model for estimating the risk of airborne transmission of a
respiratory infection such as COVID-19, is presented. The model employs basic
concepts from fluid dynamics and incorporates the known scope of factors
involved in the airborne transmission of such diseases. Simplicity in the
mathematical form of the model is by design, so that it can serve not only as a
common basis for scientific inquiry across disciplinary boundaries, but also be
understandable by a broad audience outside science and academia. The caveats
and limitations of the model are discussed in detail. The model is used to
assess the protection from transmission afforded by face coverings made from a
variety of fabrics. The reduction in transmission risk associated with
increased physical distance between the host and susceptible is also quantified
by coupling the model with available data on scalar dispersion in canonical
flows. Finally, the effect of the level of physical activity (or exercise
intensity) of the host and the susceptible in enhancing transmission risk, is
also assessed.


  Proton nuclear magnetic resonance (1H-NMR) is a widely used tool for chemical
structural analysis. However, 1H-NMR spectra suffer from natural aberrations
that render computer-assisted automated identification of these spectra
difficult, and at times impossible. Previous efforts have successfully
implemented instrument dependent or conditional identification of these
spectra. In this paper, we report the first instrument independent
computer-assisted automated identification system for a group of complex
carbohydrates known as the xyloglucan oligosaccharides. The developed system is
also implemented on the world wide web (http://www.ccrc.uga.edu) as part of an
identification package called the CCRC-Net and is intended to recognize any
submitted 1H-NMR spectrum of these structures with reasonable signal-to-noise
ratio, recorded on any 500 MHz NMR instrument. The system uses Artificial
Neural Networks (ANNs) technology and is insensitive to the instrument and
environment-dependent variations in 1H-NMR spectroscopy. In this paper,
comparative results of the ANN engine versus a multidimensional Bayes'
classifier is also presented.


  To capture the death rates and strong weekly, biweekly and probably monthly
patterns in the Canada COVID-19, we utilize the generalized additive models in
the absence of direct statistically based measurement of infection rates. By
examining the death rates of Canada in general and Quebec, Ontario and Alberta
in particular, one can easily figured out that there are substantial
overdispersion relative to the Poisson so that the negative binomial
distribution is an appropriate choice for the analysis. Generalized additive
models (GAMs) are one of the main modeling tools for data analysis. GAMs can
efficiently combine different types of fixed, random and smooth terms in the
linear predictor of a regression model to account for different types of
effects. GAMs are a semi-parametric extension of the generalized linear models
(GLMs), used often for the case when there is no a priori reason for choosing a
particular response function such as linear, quadratic, etc. and need the data
to 'speak for themselves'. GAMs do this via the smoothing functions and take
each predictor variable in the model and separate it into sections delimited by
'knots', and then fit polynomial functions to each section separately, with the
constraint that there are no links at the knots - second derivatives of the
separate functions are equal at the knots.


  Mediterranean ecosystems such as those found in California, Central Chile,
Southern Europe, and Southwest Australia host numerous, diverse, fire-adapted
micro-ecosystems. These micro-ecosystems are as diverse as mountainous conifer
to desert-like chaparral communities. Over the last few centuries, human
intervention, invasive species, and climate warming have drastically affected
the composition and health of Mediterranean ecosystems on almost every
continent. Increased fuel load from fire suppression policies and the continued
range expansion of non-native insects and plants, some driven by long-term
drought, produced the deadliest wildfire season on record in 2018. As a
consequence of these fires, a large number of structures are destroyed,
releasing household chemicals into the environment as uncontrolled toxins. The
mobilization of these materials can lead to health risks and disruption in both
human and natural systems. This article identifies drivers that led to a
structural weakening of the mosaic of fire-adapted ecosystems in California,
and subsequently increased the risk of destructive and explosive wildfires
throughout the state. Under a new climate regime, managing the impacts on
systems moving out-of-phase with natural processes may protect lives and ensure
the stability of ecosystem services.


  Gene regulation is an important fundamental biological process. The
regulation of gene expression is managed through a variety of methods including
epigenetic processes (e.g., DNA methylation). Understanding the role of
epigenetic changes in gene expression is a fundamental question of molecular
biology. Predictions of gene expression values from epigenetic data have
tremendous research and clinical potential. Despite active research, studies to
date have focused on using statistical models to predict gene expression from
methylation data. In contrast, dynamical systems can be used to generate a
model to predict gene expression using epigenetic data and a gene regulatory
network (GRN) which can also serve as a mechanistic hypothesis. Here we present
a novel stochastic dynamical systems model that predicts gene expression levels
from methylation data of genes in a given GRN. We provide an evaluation of the
model using real patient data and a GRN created from robust reference sources.
Software for dataset preparation, model parameter fitting and prediction
generation, and reporting are available at
\verb|https://github.com/kordk/stoch_epi_lib|.


  Background: One of the most challenging keys to understand COVID-19 evolution
is to have a measure on those mild cases which are never tested because their
few symptoms are soft and/or fade away soon. The problem is not only that they
are difficult to identify and test, but also that it is believed that they may
constitute the bulk of the cases and could be crucial in the pandemic equation.
  Methods: We present a novel algorithm to extract the number of these mild
cases by correlating a COVID-line calls to reported cases in given districts.
The key assumption is to realize that, being a highly contagious disease, the
number of calls by mild cases should be proportional to the number of reported
cases. Whereas a background of calls not related to infected people should be
proportional to the district population.
  Results: We find that for Buenos Aires Province, in addition to the
background, there are in signal 6.6 +/- 0.4 calls per each reported COVID-19
case. Using this we estimate in Buenos Aires Province 20 +/- 2 COVID-19
symptomatic cases for each one reported.
  Conclusions: A very simple algorithm that models the COVID-line calls as sum
of signal plus background allows to estimate the crucial number of the rate of
symptomatic to reported COVID-19 cases in a given district. The result from
this method is an early and inexpensive estimate and should be contrasted to
other methods such as serology and/or massive testing.


  In this article, we present a modification to the network-based stochastic
SEIR epidemic model which allows for modifications to the underlying contact
network to account for the effects of quarantine. We also discuss the changes
needed to the model to incorporate situations where some proportion of the
individuals who are infected remain asymptomatic throughout the course of the
disease. Using a generic network model where every potential contact exists
with the same common probability, we conduct a simulation study in which we
vary four key model parameters (transmission rate, probability of remaining
asymptomatic, and the mean lengths of time spent in the Exposed and Infectious
disease states) and examine the resulting impacts on various metrics of
epidemic severity, including the effective reproduction number. We find that
the mean length of time spent in the Infectious state and the transmission rate
are the most important model parameters, while the mean length of time spent in
the Exposed state and the probability of remaining asymptomatic are less
important.