  Many pairwise additive force fields are in active use for intrinsically
disordered proteins (IDPs) and regions (IDRs), some of which modify energetic
terms to improve description of IDPs/IDRs, but are largely in disagreement with
solution experiments for the disordered states. We have evaluated
representative pairwise and many-body protein and water force fields against
experimental data on representative IDPs and IDRs, a peptide that undergoes a
disorder-to-order transition, and for seven globular proteins ranging in size
from 130-266 amino acids. We find that force fields with the largest
statistical fluctuations consistent with the radius of gyration and universal
Lindemann values for folded states simultaneously better describe IDPs and IDRs
and disorder to order transitions. Hence the crux of what a force field should
exhibit to well describe IDRs/IDPs is not just the balance between protein and
water energetics, but the balance between energetic effects and configurational
entropy of folded states of globular proteins.


  Proteins perform a large variety of functions in living organisms, thus
playing a key role in biology. As of now, available learning algorithms to
process protein data do not consider several particularities of such data
and/or do not scale well for large protein conformations. To fill this gap, we
propose two new learning operations enabling deep 3D analysis of large-scale
protein data. First, we introduce a novel convolution operator which considers
both, the intrinsic (invariant under protein folding) as well as extrinsic
(invariant under bonding) structure, by using $n$-D convolutions defined on
both the Euclidean distance, as well as multiple geodesic distances between
atoms in a multi-graph. Second, we enable a multi-scale protein analysis by
introducing hierarchical pooling operators, exploiting the fact that proteins
are a recombination of a finite set of amino acids, which can be pooled using
shared pooling matrices. Lastly, we evaluate the accuracy of our algorithms on
several large-scale data sets for common protein analysis tasks, where we
outperform state-of-the-art methods.


  We study the effectiveness of tracking and testing in mitigating or
suppressing epidemic outbreaks, in combination with or as an alternative to
quarantines and global lockdowns. We study these intervention methods on a
network-based SEIR model, augmented with an additional probability to model
symptomatic, asymptomatic and pre-symptomatic cases. Our focus is on the basic
trade-offs between economic costs and human lives lost, and how these
trade-offs change under different lockdown, quarantine, tracking and testing
policies.
  Our main findings are as follows: (i) Tests combined with patient quarantines
reduce both economic costs and mortality, but require a large-scale testing
capacity to achieve a significant improvement; (ii) Tracking significantly
reduces both economic costs and mortality; (iii) Tracking combined with a
limited number of tests can achieve containment without lockdowns; (iv) If
there is a small flow of new incoming infections, dynamic "On-Off" lockdowns
are more efficient than fixed lockdowns.
  Our simulation results underline the extreme effectiveness of tracking and
testing policies in reducing both economic costs and mortality and their
potential to contain epidemic outbreaks without imposing social distancing
restrictions. This highlights the difficult social question of trading-off
these gains with the privacy loss that tracking necessarily entails.


  Among the many aspects that characterize the COVID-19 pandemic, two seem
particularly challenging to understand: (i) the great geographical differences
in the degree of virus contagiousness and lethality which were found in the
different phases of the epidemic progression, and (ii) the potential role of
the infected people's blood type in both the virus infectivity and the
progression of the disease. A recent hypothesis could shed some light on both
aspects. Specifically, it has been proposed that in the subject-to-subject
transfer SARS-CoV-2 conserves on its capsid the erythrocytes' antigens of the
source subject. Thus these conserved antigens can potentially cause an immune
reaction in a receiving subject that has previously acquired specific
antibodies for the source subject antigens. This hypothesis implies a blood
type-dependent infection rate. The strong geographical dependence of the blood
type distribution could be, therefore, one of the factors at the origin of the
observed heterogeneity in the epidemics spread. Here, we present an
epidemiological deterministic model where the infection rules based on blood
types are taken into account and compare our model outcomes with the exiting
worldwide infection progression data. We found an overall good agreement, which
strengthens the hypothesis that blood types do play a role in the COVID-19
infection.


  Regular spatial structures emerge in a wide range of different dynamics
characterized by local and/or nonlocal coupling terms. In several research
fields this has spurred the study of many models, which can explain pattern
formation. The modulations of patterns, occurring on long spatial and temporal
scales, can not be captured by linear approximation analysis. Here, we show
that, starting from a general model with long range couplings displaying
patterns, the spatio-temporal evolution of large scale modulations at the onset
of instability is ruled by the well-known Ginzburg-Landau equation,
independently of the details of the dynamics. Hence, we demonstrate the
validity of such equation in the description of the behavior of a wide class of
systems. We introduce a novel mathematical framework that is also able to
retrieve the analytical expressions of the coefficients appearing in the
Ginzburg-Landau equation as functions of the model parameters. Such framework
can include higher order nonlocal interactions and has much larger
applicability than the model considered here, possibly including pattern
formation in models with very different physical features.


  Bovine tuberculosis, a disease that affects cattle and badgers in Ireland,
was studied via stochastic epidemic modeling using incidence data from the Four
Area Project (Griffin et al., 2005). The Four Area Project was a large scale
field trial conducted in four diverse farming regions of Ireland over a
five-year period (1997-2002) to evaluate the impact of badger culling on bovine
tuberculosis incidence in cattle herds.
  Based on the comparison of several models, the model with no between-herd
transmission and badger-to-herd transmission proportional to the total number
of infected badgers culled was best supported by the data.
  Detailed model validation was conducted via model prediction, identifiability
checks and sensitivity analysis.
  The results suggest that badger-to-cattle transmission is of more importance
than between-herd transmission and that if there was no badger-to-herd
transmission, levels of bovine tuberculosis in cattle herds in Ireland could
decrease considerably.


  The rapid rise of antibiotic resistance is a serious threat to global public
health. Without further incentives, pharmaceutical companies have little
interest in developing antibiotics, since the success probability is low and
development costs are huge. The situation is exacerbated by the "antibiotics
dilemma": Developing narrow-spectrum antibiotics against resistant bacteria is
most beneficial for society, but least attractive for companies since their
usage is more limited than for broad-spectrum drugs and thus sales are low.
Starting from a general mathematical framework for the study of
antibiotic-resistance dynamics with an arbitrary number of antibiotics, we
identify efficient treatment protocols and introduce a market-based refunding
scheme that incentivizes pharmaceutical companies to develop narrow-spectrum
antibiotics: Successful companies can claim a refund from a newly established
antibiotics fund that partially covers their development costs. The proposed
refund involves a fixed and variable part. The latter (i) increases with the
use of the new antibiotic for currently resistant strains in comparison with
other newly developed antibiotics for this purpose---the resistance
premium---and (ii) decreases with the use of this antibiotic for non-resistant
bacteria. We outline how such a refunding scheme can solve the antibiotics
dilemma and cope with various sources of uncertainty inherent in antibiotic
R\&D. Finally, connecting our refunding approach to the recently established
antimicrobial resistance (AMR) action fund, we discuss how the antibiotics fund
can be financed.


  In this short technical report we model, within the Bayesian framework, the
rate of positive tests reported by the the State of Indiana, accounting also
for the substantial variability (and overdispeartion) in the daily count of the
tests performed. The approach we take, results with a simple procedure for
prediction, a posteriori, of this rate of 'positivity' and allows for an easy
and a straightforward adaptation by any agency tracking daily results of
COVID-19 tests. The numerical results provided herein were obtained via an
updatable R Markdown document.


  The infectivity profile of an individual with COVID-19 is attributed to the
paper Temporal dynamics in viral shedding and transmissibility of COVID-19 by
He et al., published in Nature Medicine in April 2020. However, the analysis
within this paper contains a mistake such that the published infectivity
profile is incorrect and the conclusion that infectiousness begins 2.3 days
before symptom onset is no longer supported. In this document we discuss the
error and compute the correct infectivity profile. We also establish confidence
intervals on this profile, quantify the difference between the published and
the corrected profiles, and discuss an issue of normalisation when fitting
serial interval data. This infectivity profile plays a central role in policy
and decision making, thus it is crucial that this issue is corrected with the
utmost urgency to prevent the propagation of this error into further studies
and policies. We hope that this preprint will reach all researchers and policy
makers who are using the incorrect infectivity profile to inform their work.


  Protein tertiary structure prediction has improved dramatically in recent
years. A considerable fraction of various proteomes can be modelled in the
absence of structural templates. We ask whether our DMPfold method can model
all the proteins without templates in the JCVI-syn3.0 minimal genome, which
contains 438 proteins. We find that a useful tertiary structure annotation can
be provided for all but 10 proteins. The models may help annotate function in
cases where it is unknown, and provide coverage for 29 predicted
protein-protein interactions which lacked monomer models. We also show that
DMPfold performs well on proteins with structures released since initial
publication. It is likely that the minimal genome will have complete structural
coverage within a few years.


  We consider species tree estimation under a standard stochastic model of gene
tree evolution that incorporates incomplete lineage sorting (as modeled by a
coalescent process) and gene duplication and loss (as modeled by a branching
process). Through a probabilistic analysis of the model, we derive sample
complexity bounds for widely used quartet-based inference methods that
highlight the effect of the duplication and loss rates in both subcritical and
supercritical regimes.


  Protein conformational changes are activated processes essential for protein
functions. Activation in a protein differs from activation in a small molecule
in that it involves directed and systematic energy flows through preferred
channels encoded in the protein structure. Understanding the nature of these
energy flow channels and how energy flows through them during activation is
critical for understanding protein conformational changes. We recently
developed a rigorous statistical mechanical framework for understanding
potential energy flows. Here we complete this theoretical framework with a
rigorous theory for kinetic energy flows: potential and kinetic energy
inter-convert when impressed forces oppose inertial forces whereas kinetic
energy transfers directly from one coordinate to another when inertial forces
oppose each other. This theory is applied to analyzing a prototypic system for
biomolecular conformational dynamics: the isomerization of an alanine
dipeptide. Among the two essential energy flow channels for this process,
dihedral phi confronts the activation barrier, whereas dihedral theta receives
energy from potential energy flows. Intriguingly, theta helps phi to cross the
activation barrier by transferring to phi via direct kinetic energy flow all
the energy it received: increase in theta caused by potential energy flow
converts into increase in phi. As a compensation, theta receives kinetic energy
from bond angle alpha via direct mechanism and bond angle beta via indirect
mechanism.


  COVID-19 patients show varying severity of the disease ranging from
asymptomatic to requiring intensive care. Although a number of SARS-CoV-2
specific monoclonal antibodies have been identified, we still lack an
understanding of the overall landscape of B-cell receptor (BCR) repertoires in
COVID-19 patients. Here, we used high-throughput sequencing of bulk and plasma
B-cells collected over multiple time points during infection to characterize
signatures of B-cell response to SARS-CoV-2 in 19 patients. Using principled
statistical approaches, we determined differential features of BCRs associated
with different disease severity. We identified 38 significantly expanded clonal
lineages shared among patients as candidates for specific responses to
SARS-CoV-2. Using single-cell sequencing, we verified reactivity of BCRs shared
among individuals to SARS-CoV-2 epitopes. Moreover, we identified natural
emergence of a BCR with cross-reactivity to SARS-CoV-1 and SARS-CoV-2 in a
number of patients. Our results provide important insights for development of
rational therapies and vaccines against COVID-19.


  Single particle, cryogenic electron microscopy (cryo-EM) experiments now
routinely produce high-resolution data for large proteins and their complexes.
Building an atomic model into a cryo-EM density map is challenging,
particularly when no structure for the target protein is known a priori.
Existing protocols for this type of task often rely on significant human
intervention and can take hours to many days to produce an output. Here, we
present a fully automated, template-free model building approach that is based
entirely on neural networks. We use a graph convolutional network (GCN) to
generate an embedding from a set of rotamer-based amino acid identities and
candidate 3-dimensional C$\alpha$ locations. Starting from this embedding, we
use a bidirectional long short-term memory (LSTM) module to order and label the
candidate identities and atomic locations consistent with the input protein
sequence to obtain a structural model. Our approach paves the way for
determining protein structures from cryo-EM densities at a fraction of the time
of existing approaches and without the need for human intervention.


  The contact structure of a population plays an important role in transmission
of infection. Many ``structured models'' capture aspects of the contact
structure through an underlying network or a mixing matrix. An important
observation in such models, is that once a fraction $1-1/\mathcal{R}_0$ has
been infected, the residual susceptible population can no longer sustain an
epidemic. A recent observation of some structured models is that this threshold
can be crossed with a smaller fraction of infected individuals, because the
disease acts like a targeted vaccine, preferentially immunizing higher-risk
individuals who play a greater role in transmission. Therefore, a limited
``first wave'' may leave behind a residual population that cannot support a
second wave once interventions are lifted. In this paper, we systematically
analyse a number of mean-field models for networks and other structured
populations to address issues relevant to the Covid-19 pandemic. In particular,
we consider herd-immunity under several scenarios. We confirm that, in networks
with high degree heterogeneity, the first wave confers herd-immunity with
significantly fewer infections than equivalent models with lower degree
heterogeneity. However, if modelling the intervention as a change in the
contact network, then this effect might become more subtle. Indeed, modifying
the structure can shield highly connected nodes from becoming infected during
the first wave and make the second wave more substantial. We confirm this
finding by using an age-structured compartmental model parameterised with real
data and comparing lockdown periods implemented either as a global scaling of
the mixing matrix or age-specific structural changes. We find that results
regarding herd immunity levels are strongly dependent on the model, the
duration of lockdown and how lockdown is implemented.


  Computer-aided drug discovery is an essential component of modern drug
development. Therein, deep learning has become an important tool for rapid
screening of billions of molecules in silico for potential hits containing
desired chemical features. Despite its importance, substantial challenges
persist in training these models, such as severe class imbalance, high decision
thresholds, and lack of ground truth labels in some datasets. In this work we
argue in favor of directly optimizing the receiver operating characteristic
(ROC) in such cases, due to its robustness to class imbalance, its ability to
compromise over different decision thresholds, certain freedom to influence the
relative weights in this compromise, fidelity to typical benchmarking measures,
and equivalence to positive/unlabeled learning. We also propose new training
schemes (coherent mini-batch arrangement, and usage of out-of-batch samples)
for cost functions based on the ROC, as well as a cost function based on the
logAUC metric that facilitates early enrichment (i.e. improves performance at
high decision thresholds, as often desired when synthesizing predicted hit
compounds). We demonstrate that these approaches outperform standard deep
learning approaches on a series of PubChem high-throughput screening datasets
that represent realistic and diverse drug discovery campaigns on major drug
target families.


  Countries highly exposed to incoming traffic from China were expected to be
at the highest risk of COVID-19 spread. However, COVID-19 case numbers
(infection levels) are negatively correlated with incoming traffic-level.
Moreover, infection levels are positively correlated with population-size,
while the latter should only affect infection-level once herd immunity is
reached. These could be explained if a low-virulence strain (LVS) began
spreading a few months earlier from China, providing immunity from the later
emerging known SARS-CoV-2 high-virulence strain (HVS). We find that the
dynamics of the COVID-19 pandemic depend on the LVS and HVS spread
doubling-times and the delay between their initial onsets. We find that LVS
doubling-time to be $T_L\sim1.59\pm0.17$ times slower than the HVS ($T_H$), but
its earlier onset allowed its global wide-spread to the levels required for
herd-immunity. In countries exposed earlier to the LVS and/or having smaller
population-size, the LVS achieved herd-immunity earlier, allowing less time for
the spread of the HVS, and giving rise to lower HVS-infection levels. Such
model accurately predicts a country's infection-level ({\rm R^{2}=0.74};
p-value of {\rm 5.2\times10^{-13}}), given only its population-size and
incoming-traffic from China. It explains the negative correlation with
incoming-traffic ($c_{exp}$), the positive correlation with the population size
(n_{pop}) and their specific relations (${\rm N}_{{\rm cases}}\propto
n_{pop}^{{\rm T_{L}/{\rm T_{H}}}}\times c_{exp}^{{\rm T_{L}/{\rm T_{H}-1}}}$).
We find that most countries should have already achieved herd-immunity. Further
COVID-19-spread in these countries is limited and is not expected to rise by
more than a factor of 2-3. We suggest tests/predictions to further verify the
model and biologically identify the LVS, and discuss the implications.


  Since the beginning of the epidemic, daily reports of CoViD-19 cases,
hospitalizations, and deaths from around the world have been publicly
available. This paper describes methods to characterize broad features of the
spread of the disease, with relatively long periods of constant transmission
rates, using a new population modeling framework based on discrete-time
difference equations. Comparative parameters are chosen for their weak
dependence on model assumptions. Approaches for their point and interval
estimation, accounting for additional sources of variance in the case data, are
presented. These methods provide a basis to quantitatively assess the impact of
changes to social distancing policies using publicly available data. As
examples, data from Ontario and German states are analyzed using this
framework. German case data show a small increase in transmission rates
following the relaxation of lock-down rules on May 6, 2020. By combining case
and death data from Germany, the mean and standard deviation of the time from
infection to death are estimated.


  Background: Wuhan, China was the epicenter of COVID-19 pandemic. The goal of
current study is to understand the infection transmission dynamics before
intervention measures were taken. Methods: Data and key events were searched
through pubmed and internet. Epidemiological data were calculated using data
extracted from a variety of data sources. Results: We established a timeline
showing by January 1, 2020, Chinese authorities had been presented convincing
evidence of human-to-human transmission; however, it was not until January 20,
2020 that this information was shared with the public. Our study estimated that
there would have been 10989 total infected cases if interventions were taken on
January 2, 2020, versus 239875 cases when lockdown was put in place on January
23, 2020. Conclusions: China's withholding of key information about the 2020
COVID-19 outbreak and its delayed response ultimately led to the largest public
health crisis of this century and could have been avoided with earlier
countermeasures.


  One of the central tools to control the COVID-19 pandemics is the knowledge
of its spreading dynamics. Here we develop a fractal model capable of describe
this dynamics, in term of daily new cases, and provide quantitative criteria
for some predictions. We propose a fractal dynamical model using conformed
derivative and fractal time scale. A Burr-XII shaped solution of the
fractal-like equation is obtained. The model is tested using data from several
countries, showing that a single function is able to describe very different
shapes of the outbreak. The diverse behavior of the outbreak on those countries
is presented and discussed. Moreover, a criterion to determine the existence of
the pandemic peak and a expression to find the time to reach herd immunity are
also obtained.


  Despite the huge effort to contain the infection, the novel SARS-CoV-2
coronavirus has rapidly become pandemics, mainly due to its extremely high
human-to-human transmission capability, and a surprisingly high viral charge of
symptom-less people. While the seek of a vaccine is still ongoing, promising
results have been obtained with antiviral compounds. In particular, lactoferrin
is found to have beneficial effects both in preventing and soothing the
infection. Here, we explore the possible molecular mechanisms with which
lactoferrin interferes with SARS-CoV-2 cell invasion, preventing attachment
and/or entry of the virus. To this aim, we search for possible interactions
lactoferrin may have with virus structural proteins and host receptors.
Representing the molecular iso-electron surface of proteins in terms of
2D-Zernike descriptors, we (i) identified putative regions on the lactoferrin
surface able to bind sialic acid receptors on the host cell membrane,
sheltering the cell from the virus attachment; (ii) showed that no significant
shape complementarity is present between lactoferrin and the ACE2 receptor,
while (iii) two high complementarity regions are found on the N- and C-terminal
domains of the SARS-CoV-2 spike protein, hinting at a possible competition
between lactoferrin and ACE2 for the binding to the spike protein.


  COVID-19 testing has become a standard approach for estimating prevalence
which then assist in public health decision making to contain and mitigate the
spread of the disease. The sampling designs used are often biased in that they
do not reflect the true underlying populations. For instance, individuals with
strong symptoms are more likely to be tested than those with no symptoms. This
results in biased estimates of prevalence (too high). Typical post-sampling
corrections are not always possible. Here we present a simple bias correction
methodology derived and adapted from a correction for publication bias in meta
analysis studies. The methodology is general enough to allow a wide variety of
customization making it more useful in practice. Implementation is easily done
using already collected information. Via a simulation and two real datasets, we
show that the bias corrections can provide dramatic reductions in estimation
error.


  Mathematical disease modelling has long operated under the assumption that
any one infectious disease is caused by one transmissible pathogen spreading
among a population. This paradigm has been useful in simplifying the biological
reality of epidemics and has allowed the modelling community to focus on the
complexity of other factors such as population structure and interventions.
However, there is an increasing amount of evidence that the strain diversity of
pathogens, and their interplay with the host immune system, can play a large
role in shaping the dynamics of epidemics. Here, we introduce a disease model
with an underlying genotype network to account for two important mechanisms.
One, the disease can mutate along network pathways as it spreads in a host
population. Two, the genotype network allows us to define a genetic distance
across strains and therefore to model the transcendence of immunity often
observed in real world pathogens. We study the emergence of epidemics in this
model, through its epidemic phase transitions, and highlight the role of the
genotype network in driving cyclicity of diseases, large scale fluctuations,
sequential epidemic transitions, as well as localization around specific
strains of the associated pathogen. More generally, our model illustrates the
richness of behaviours that are possible even in well-mixed host populations
once we consider strain diversity and go beyond the "one disease equals one
pathogen" paradigm.


  We develop an agent-based model of disease transmission in remote communities
in Western Australia. Despite extreme isolation, we show that the movement of
people amongst a large number of small but isolated communities has the effect
of causing transmission to spread quickly. Significant movement between remote
communities, and regional and urban centres allows for infection to quickly
spread to and then among these remote communities. Our conclusions are based on
two characteristic features of remote communities in Western Australia: (1)
high mobility of people amongst these communities, and (2) relatively high
proportion of travellers from very small communities to major population
centres. In models of infection initiated in the state capital, Perth, these
remote communities are collectively and uniquely vulnerable. Our model and
analysis does not account for possibly heightened impact due to preexisting
conditions, such additional assumptions would only make the projections of this
model more dire. We advocate stringent monitoring and control of movement to
prevent significant impact on the indigenous population of Western Australia.


  The coronavirus disease 2019 (COVID-19) had caused more that 8 million
infections as of middle June 2020. Recently, Brazil has become a new epicentre
of COVID-19, while India and African region are potential epicentres. This
study aims to predict the inflection point and outbreak size of these
new/potential epicentres at the early phase of the epidemics by borrowing
information from more `mature' curves from other countries. We modeled the
cumulative cases to the well-known sigmoid growth curves to describe the
epidemic trends under the mixed-effect models and using the four-parameter
logistic model after power transformations. African region is predicted to have
the largest total outbreak size of 3.9 million cases (2.2 to 6 million), and
the inflection will come around September 13, 2020. Brazil and India are
predicted to have a similar final outbreak size of around 2.5 million cases
(1.1 to 4.3 million), with the inflection points arriving June 23 and July 26,
respectively. We conclude in Brazil, India, and African the epidemics of COVI19
have not yet passed the inflection points; these regions potentially can take
over USA in terms of outbreak size


  Membrane tubulation is a ubiquitous process that occurs both at the plasma
membrane and on the membranes of intracellular organelles. These tubulation
events are known to be mediated by forces applied on the membrane either due to
motor proteins, by polymerization of the cytoskeleton, or due to the
interactions between membrane proteins binding onto the membrane. The numerous
experimental observations of tube formation have been amply supported by
mathematical modeling of the associated membrane mechanics and have provided
insights into the force-displacement relationships of membrane tubes. Recent
advances in quantitative biophysical measurements of membrane-protein
interactions and tubule formation have necessitated the need for advances in
modeling that will account for the interplay of multiple aspects of physics
that occur simultaneously. Here, we present a comprehensive review of
experimental observations of tubule formation and provide context from the
framework of continuum modeling. Finally, we explore the scope for future
research in this area with an emphasis on iterative modeling and experimental
measurements that will enable us to expand our mechanistic understanding of
tubulation processes in cells.


  Here, we consider an SIS epidemic model where the individuals are distributed
on several distinct patches. We construct a stochastic model and then prove
that it converges to a deterministic model as the total population size tends
to infinity. Furthermore, we show the existence and the global stability of a
unique endemic equilibrium provided that the migration rates of susceptible and
infectious individuals are equal. Finally, we compare the equilibra with those
of the homogeneous model, and with those of isolated patches.


  3CL-Pro (or M-Pro) is the SARS-CoV-2 main protease, acting as a homodimer, is
responsible for the cleavage of the large polyprotein 1ab transcript in
proteins acting on viral growth and replication. 3CL-Pro has been one of the
most studied SARS-CoV-2 proteins and the subject of therapeutic interventions,
targeting its catalytic domain. A number of drug candidates have been reported,
including some natural products. Here, we investigated in silico, through
binding and molecular dynamics simulations, the natural product space for the
identification of candidates of 3CL-Pro dimerization inhibitors. We report that
fortunellin (acacetin 7-O-neohesperidoside), a natural flavonoid O-glycoside,
is a potent inhibitor of 3CL-Pro dimerization. A search of the ZINC natural
products database identified another 16 related molecules, including apilin and
rhoifolin, with interesting pharmacological properties. We propose that
fortunellin and its structural analogs might be the basis of novel
pharmaceuticals and dietary supplements against SARS-CoV-2 induced COVID-19
disease.


  An incoherent feed-forward loop (IFFL) is a network motif known for its
ability to accelerate responses and generate pulses. Though functions of IFFLs
are well studied, most previous computational analysis of IFFLs used ordinary
differential equation (ODE) models where retroactivity, the effect downstream
binding sites exert on the dynamics of an upstream transcription factor (TF),
was not considered. It remains an open question to understand the behavior of
IFFLs in contexts with high levels of retroactivity, e.g., in cells
transformed/transfected with high-copy plasmids, or in eukaryotic cells where a
TF binds to numerous high-affinity binding sites in addition to one or more
functional target sites. Here we study the behavior of IFFLs by simulating and
comparing ODE models with different levels of retroactivity. We find that
increasing retroactivity in an IFFL can increase, decrease, or keep the
network's response time and pulse amplitude constant. This suggests that
increasing retroactivity, traditionally considered as an impediment to
designing robust synthetic systems, could be exploited to improve the
performance of IFFLs. We compare the behaviors of IFFLs to negative
autoregulatory loops, another sign-sensitive response-accelerating network
motif, and find that increasing retroactivity in a negative autoregulated
circuit can only slow the response. The inability of a negative autoregulatory
loop to flexibly handle retroactivity may have contributed to its lower
abundance in eukaryotic relative to bacterial regulatory networks, a sharp
contrast to the significant abundance of IFFLs in both cell types.


  A new epidemic model for Covid-19 has been constructed and simulated for
eight US states. The coefficients for this model, based on seven coupled
differential equations, are carefully evaluated against recorded data on cases
and deaths. These projections reveal that Covid-19 will become endemic,
spreading for more than two years. If stay-at-home orders are relaxed, most
states may experience a secondary peak in 2021. The number of Covid-19 deaths
could have been significantly lower in most states that opened up, if lockdowns
had been maintained. Additionally, our model predicts that decreasing contact
rate by 10%, or increasing testing by approximately 15%, or doubling lockdown
compliance (from the current $\sim$ 15%) will eradicate infections in Texas
within a year. Applied to the entire US, the predictions based on the current
situation indicate about 11 million total infections (including undetected), 8
million cumulative confirmed cases, and 630,000 cumulative deaths by November
1, 2020.


  COVID-19--a viral infectious disease--has quickly emerged as a global
pandemic infecting millions of people with a significant number of deaths
across the globe. The symptoms of this disease vary widely. Depending on the
symptoms an infected person is broadly classified into two categories namely,
asymptomatic and symptomatic. Asymptomatic individuals display mild or no
symptoms but continue to transmit the infection to otherwise healthy
individuals. This particular aspect of asymptomatic infection poses a major
obstacle in managing and controlling the transmission of the infectious
disease. In this paper, we attempt to mathematically model the spread of
COVID-19 in India under various intervention strategies. We consider SEIR type
epidemiological models, incorporated with India specific social contact matrix
representing contact structures among different age groups of the population.
Impact of various factors such as presence of asymptotic individuals, lockdown
strategies, social distancing practices, quarantine, and hospitalization on the
disease transmission is extensively studied. Numerical simulation of our model
is matched with the real COVID-19 data of India till May 15, 2020 for the
purpose of estimating the model parameters. Our model with zone-wise lockdown
is seen to give a decent prediction for July 20, 2020.


  Raw data on the cumulative number of deaths at a country level generally
indicate a spatially variable distribution of the incidence of COVID-19
disease. An important issue is to determine whether this spatial pattern is a
consequence of environmental heterogeneities, such as the climatic conditions,
during the course of the outbreak. Another fundamental issue is to understand
the spatial spreading of COVID-19. To address these questions, we consider four
candidate epidemiological models with varying complexity in terms of initial
conditions, contact rates and non-local transmissions, and we fit them to
French mortality data with a mixed probabilistic-ODE approach. Using standard
statistical criteria, we select the model with non-local transmission
corresponding to a diffusion on the graph of counties that depends on the
geographic proximity, with time-dependent contact rate and spatially constant
parameters. This original spatially parsimonious model suggests that in a
geographically middle size centralized country such as France, once the
epidemic is established, the effect of global processes such as restriction
policies, sanitary measures and social distancing overwhelms the effect of
local factors. Additionally, this modeling approach reveals the latent
epidemiological dynamics including the local level of immunity, and allows us
to evaluate the role of non-local interactions on the future spread of the
disease. In view of its theoretical and numerical simplicity and its ability to
accurately track the COVID-19 epidemic curves, the framework we develop here,
in particular the non-local model and the associated estimation procedure, is
of general interest in studying spatial dynamics of epidemics.


  We propose a Bayesian Heterogeneity Learning approach for
Susceptible-Infected-Removal-Susceptible (SIRS) model that allows underlying
clustering patterns for transmission rate, recovery rate, and loss of immunity
rate for the latest coronavirus (COVID-19) among different regions. Our
proposed method provides simultaneously inference on parameter estimation and
clustering information which contains both number of clusters and cluster
configurations. Specifically, our key idea is to formulates the SIRS model into
a hierarchical form and assign the Mixture of Finite mixtures priors for
heterogeneity learning. The properties of the proposed models are examined and
a Markov chain Monte Carlo sampling algorithm is used to sample from the
posterior distribution. Extensive simulation studies are carried out to examine
empirical performance of the proposed methods. We further apply the proposed
methodology to analyze the state level COVID-19 data in U.S.


  Genetic circuits need a cellular environment to operate in, which naturally
couples the circuit function with the overall functionality of gene regulatory
network. To execute their functions all gene circuits draw resources in the
form of RNA polymerases, ribosomes, and tRNAs. Recent experiments pointed out
that the role of resource competition on synthetic circuit outputs could be
immense. However, the effect of complexity of the circuit architecture on
resource sharing dynamics is yet unexplored. In this paper, we employ
mathematical modelling and in-silico experiments to identify the sources of
resource trade-off and to quantify its impact on the function of a genetic
circuit, keeping our focus on regulation of immediate downstream proteins. We
take the example of the fluorescent reporters, which are often used as protein
read-outs. We show that estimating gene expression dynamics from readings of
downstream protein data might be unreliable when the resource is limited and
ribosome affinities are asymmetric. We focus on the impact of mRNA copy number
and RBS strength on the nonlinear isocline that emerges with two regimes,
prominently separated by a tipping point, and study how correlation and
competition dominate each other depending on various circuit parameters.
Focusing further on genetic toggle circuit, we have identified major effects of
resource competition in this model motif, and quantified the observations. The
observations are testable in wet-lab experiments, as all the parameters chosen
are experimentally relevant.


  Deep learning is catalyzing a scientific revolution fueled by big data,
accessible toolkits, and powerful computational resources, impacting many
fields including protein structural modeling. Protein structural modeling, such
as predicting structure from amino acid sequence and evolutionary information,
designing proteins toward desirable functionality, or predicting properties or
behavior of a protein, is critical to understand and engineer biological
systems at the molecular level. In this review, we summarize the recent
advances in applying deep learning techniques to tackle problems in protein
structural modeling and design. We dissect the emerging approaches using deep
learning techniques for protein structural modeling, and discuss advances and
challenges that must be addressed. We argue for the central importance of
structure, following the "sequence -> structure -> function" paradigm. This
review is directed to help both computational biologists to gain familiarity
with the deep learning methods applied in protein modeling, and computer
scientists to gain perspective on the biologically meaningful problems that may
benefit from deep learning techniques.


  Phylogenetics uses alignments of molecular sequence data to learn about
evolutionary trees. Substitutions in sequences are modelled through a
continuous-time Markov process, characterised by an instantaneous rate matrix,
which standard models assume is time-reversible and stationary. These
assumptions are biologically questionable and induce a likelihood function
which is invariant to a tree's root position. This hampers inference because a
tree's biological interpretation depends critically on where it is rooted.
Relaxing both assumptions, we introduce a model whose likelihood can
distinguish between rooted trees. The model is non-stationary, with step
changes in the instantaneous rate matrix at each speciation event. Exploiting
recent theoretical work, each rate matrix belongs to a non-reversible family of
Lie Markov models. These models are closed under matrix multiplication, so our
extension offers the conceptually appealing property that a tree and all its
sub-trees could have arisen from the same family of non-stationary models.
  We adopt a Bayesian approach, describe an MCMC algorithm for posterior
inference and provide software. The biological insight that our model can
provide is illustrated through an analysis in which non-reversible but
stationary, and non-stationary but reversible models cannot identify a
plausible root.


  We interpret attitudes towards science and pseudosciences as cultural traits
that diffuse in society through communication efforts exerted by agents. We
present a tractable model that allows us to study the interaction among the
diffusion of an epidemic, vaccination choices, and the dynamics of cultural
traits. We apply it to study the impact of homophily between pro-vaxxers and
anti-vaxxers on the total number of cases (the cumulative infection). We show
that, during the outbreak of a disease, homophily has the direct effect of
decreasing the speed of recovery. Hence, it may increase the number of cases
and make the disease endemic. The dynamics of the shares of the two cultural
traits in the population is crucial in determining the sign of the total effect
on the cumulative infection: more homophily is beneficial if agents are not too
flexible in changing their cultural trait, is detrimental otherwise.


  Phylogenetic networks can represent evolutionary events that cannot be
described by phylogenetic trees. These networks are able to incorporate
reticulate evolutionary events such as hybridization, introgression, and
lateral gene transfer. Recently, network-based Markov models of DNA sequence
evolution have been introduced along with model-based methods for
reconstructing phylogenetic networks. For these methods to be consistent, the
network parameter needs to be identifiable from data generated under the model.
Here, we show that the semi-directed network parameter of a triangle-free,
level-1 network model with any fixed number of reticulation vertices is
generically identifiable under the Jukes-Cantor, Kimura 2-parameter, or Kimura
3-parameter constraints.


  Approaches to the calculation of the full state vector of a larger
epidemiological model for the spread of COVID-19 in Sweden at the initial time
instant from available data and with a simplified dynamical model are proposed
and evaluated. The larger epidemiological model is based on a continuous Markov
chain and captures the demographic composition of and the transport flows
between the counties of Sweden. Its intended use is to predict the outbreak
development in temporal and spatial coordinates as well as across the
demographic groups. It can also support evaluating and comparing of prospective
intervention strategies in terms of e.g. lockdown in certain areas or isolation
of specific age groups. The simplified model is a discrete time-invariant
linear system that has cumulative infectious incidence, infected population,
asymptomatic population, exposed population, and infectious pressure as the
state variables. Since the system matrix of the model depends on a number
transition rates, structural properties of the model are investigated for
suitable parameter ranges. It is concluded that the model becomes unobservable
for some parameter values. Two contrasting approaches to the initial state
estimation are considered. One is a version of Rauch-Tung-Striebel smoother and
another is based on solving a batch nonlinear optimization problem. The
benefits and shortcomings of the considered estimation techniques are analyzed
and compared on synthetic data for several Swedish counties.


  We analyze a general theory for coexistence and extinction of ecological
communities that are influenced by stochastic temporal environmental
fluctuations. The results apply to discrete time (stochastic difference
equations), continuous time (stochastic differential equations), compact and
non-compact state spaces and degenerate or non-degenerate noise. In addition,
we can also include in the dynamics auxiliary variables that model
environmental fluctuations, population structure, eco-environmental feedbacks
or other internal or external factors.
  We are able to significantly generalize the recent discrete time results by
Benaim and Schreiber (Journal of Mathematical Biology '19) to non-compact state
spaces, and we provide stronger persistence and extinction results. The
continuous time results by Hening and Nguyen (Annals of Applied Probability
'18) are strengthened to include degenerate noise and auxiliary variables.
  Using the general theory, we work out several examples. In discrete time, we
classify the dynamics when there are one or two species, and look at the Ricker
model, Log-normally distributed offspring models, lottery models, discrete
Lotka-Volterra models as well as models of perennial and annual organisms. For
the continuous time setting we explore models with a resource variable,
stochastic replicator models, and three dimensional Lotka-Volterra models.


  Sterols have been ascribed a major role in the organization of biological
membranes, in particular for the formation of liquid ordered domains in complex
lipid mixtures. Here, we employed molecular dynamics simulations to compare the
effects of cholesterol and ergosterol as the major sterol of mammalian and
fungal cells, respectively, on binary mixtures with
1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) as a proxy for saturated
lipids. In agreement with previous work, we observe that the addition of sterol
molecules modifies the order of DPPC both in the gel phase and in the liquid
phase. When disentangling the overall tilt angle and the structure of the tail
imposed by trans/gauche configurations of torsion angles in the tail,
respectively, a more detailed picture of the impact of sterols can be
formulated, revealing, e.g., an approximate temperature-concentration
superposition ranging from the liquid to the gel phase. Furthermore, a new
quantitative measure to identify the presence of collective sterol effects is
discussed. Moreover, when comparing both types of sterols, addition of
cholesterol has a noticeably stronger impact on phospholipid properties than of
ergosterol. The observed differences can be attributed to higher planarity of
the cholesterol ring system. This planarity combined with an inherent asymmetry
in its molecular interactions leads to better alignment and hence stronger
interaction with saturated acyl chains. Our results suggest that the high order
demonstrated for ergosterol in fungal plasma membranes must therefore be
generated via additional mechanisms.


  Non-genetic heterogeneity is key to cellular decisions, as even genetically
identical cells respond in very different ways to the same external stimulus,
e.g., during cell differentiation or therapeutic treatment of disease. Strong
heterogeneity is typically already observed at the level of signaling pathways
that are the first sensors of external inputs and transmit information to the
nucleus where decisions are made. Since heterogeneity arises from random
fluctuations of cellular components, mathematical models are required to fully
describe the phenomenon and to understand the dynamics of heterogeneous cell
populations. Here, we review the experimental and theoretical literature on
cellular signaling heterogeneity, with special focus on the TGFb/SMAD signaling
pathway.


  Protein interactions are important in a broad range of biological processes.
Traditionally, computational methods have been developed to automatically
predict protein interface from hand-crafted features. Recent approaches employ
deep neural networks and predict the interaction of each amino acid pair
independently. However, these methods do not incorporate the important
sequential information from amino acid chains and the high-order pairwise
interactions. Intuitively, the prediction of an amino acid pair should depend
on both their features and the information of other amino acid pairs. In this
work, we propose to formulate the protein interface prediction as a 2D dense
prediction problem. In addition, we propose a novel deep model to incorporate
the sequential information and high-order pairwise interactions to perform
interface predictions. We represent proteins as graphs and employ graph neural
networks to learn node features. Then we propose the sequential modeling method
to incorporate the sequential information and reorder the feature matrix. Next,
we incorporate high-order pairwise interactions to generate a 3D tensor
containing different pairwise interactions. Finally, we employ convolutional
neural networks to perform 2D dense predictions. Experimental results on
multiple benchmarks demonstrate that our proposed method can consistently
improve the protein interface prediction performance.


  Spectral clustering is one of the fundamental unsupervised learning methods
widely used in data analysis. Sparse spectral clustering (SSC) imposes sparsity
to the spectral clustering and it improves the interpretability of the model.
This paper considers a widely adopted model for SSC, which can be formulated as
an optimization problem over the Stiefel manifold with nonsmooth and nonconvex
objective. Such an optimization problem is very challenging to solve. Existing
methods usually solve its convex relaxation or need to smooth its nonsmooth
part using certain smoothing techniques. In this paper, we propose a manifold
proximal linear method (ManPL) that solves the original SSC formulation. We
also extend the algorithm to solve the multiple-kernel SSC problems, for which
an alternating ManPL algorithm is proposed. Convergence and iteration
complexity results of the proposed methods are established. We demonstrate the
advantage of our proposed methods over existing methods via the single-cell RNA
sequencing data analysis.


  The COVID-19 pandemic, caused by the novel coronavirus SARS-CoV-2, is
currently spreading rapidly throughout the world, causing many deaths and
severe economic damage. It is believed that hot and humid conditions do not
favor the novel coronavirus, yet this is still under debate due to many
uncertainties associated with the COVID-19 data. Here we propose surrogate data
tests to examine the preference of this virus to spread under different climate
conditions. We find that the disease is significantly (above the 95% confidence
level) more common when the temperature is $\sim$10$^\circ$C, the specific
humidity is $\sim$5 g/kg, and the ultraviolet (UV) radiation is $\sim$80
kJ/m$^2$. The significance of relative humidity is below the 95% confidence
level and does not show a preferred value. The results are supported using
global and regional data, spanning the time period from January to July 2020.
The COVID-19 data includes the daily reported new cases and daily death cases;
for both, the population size is either taken into account or ignored.


  We present a unifying, tractable approach for studying the spread of viruses
causing complex diseases requiring to be modeled using a large number of types
(e.g., infective stage, clinical state, risk factor class). We show that
recording each infected individual's infection age, i.e., the time elapsed
since infection, has three benefits.
  First, regardless of the number of types, the age distribution of the
population can be described by means of a first-order, one-dimensional partial
differential equation (PDE) known as the McKendrick-von Foerster equation. The
frequency of type $i$ is simply obtained by integrating the probability of
being in state $i$ at a given age against the age distribution. This
representation induces a simple methodology based on the additional assumption
of Poisson sampling to infer and forecast the epidemic. We illustrate this
technique using French data from the COVID-19 epidemic.
  Second, our approach generalizes and simplifies standard compartmental models
using high-dimensional systems of ordinary differential equations (ODEs) to
account for disease complexity. We show that such models can always be
rewritten in our framework, thus, providing a low-dimensional yet equivalent
representation of these complex models.
  Third, beyond the simplicity of the approach, we show that our population
model naturally appears as a universal scaling limit of a large class of fully
stochastic individual-based epidemic models, where the initial condition of the
PDE emerges as the limiting age structure of an exponentially growing
population starting from a single individual.


  We consider an optimization problem in which a (single) bat aims to exploit
the nectar in a set of $n$ cacti with the objective of maximizing the expected
total amount of nectar it drinks. Each cactus $i \in [n]$ is characterized by a
parameter $r_{i} > 0$ that determines the rate in which nectar accumulates in
$i$. In every round, the bat can visit one cactus and drink all the nectar
accumulated there since its previous visit. Furthermore, competition with other
bats, that may also visit some cacti and drink their nectar, is modeled by
means of a stochastic process in which cactus $i$ is emptied in each round
(independently) with probability $0 < s_i < 1$. Our attention is restricted to
purely-stochastic strategies that are characterized by a probability vector
$(p_1, \ldots, p_n)$ determining the probability $p_i$ that the bat visits
cactus $i$ in each round. We prove that for every $\epsilon > 0$, there exists
a purely-stochastic strategy that approximates the optimal purely-stochastic
strategy to within a multiplicative factor of $1 + \epsilon$, while exploiting
only a small core of cacti. Specifically, we show that it suffices to include
at most $\frac{2 (1 - \sigma)}{\epsilon \cdot \sigma}$ cacti in the core, where
$\sigma = \min_{i \in [n]} s_{i}$. We also show that this upper bound on core
size is asymptotically optimal as a core of a significantly smaller size cannot
provide a $(1 + \epsilon)$-approximation of the optimal purely-stochastic
strategy. This means that when the competition is more intense (i.e., $\sigma$
is larger), a strategy based on exploiting smaller cores will be favorable.


  The identification of essential genes/proteins is a critical step towards a
better understanding of human biology and pathology. Computational approaches
helped to mitigate experimental constraints by exploring machine learning (ML)
methods and the correlation of essentiality with biological information,
especially protein-protein interaction (PPI) networks, to predict essential
genes. Nonetheless, their performance is still limited, as network-based
centralities are not exclusive proxies of essentiality, and traditional ML
methods are unable to learn from non-Euclidean domains such as graphs. Given
these limitations, we proposed EPGAT, an approach for essentiality prediction
based on Graph Attention Networks (GATs), which are attention-based Graph
Neural Networks (GNNs) that operate on graph-structured data. Our model
directly learns patterns of gene essentiality from PPI networks, integrating
additional evidence from multiomics data encoded as node attributes. We
benchmarked EPGAT for four organisms, including humans, accurately predicting
gene essentiality with AUC score ranging from 0.78 to 0.97. Our model
significantly outperformed network-based and shallow ML-based methods and
achieved a very competitive performance against the state-of-the-art node2vec
embedding method. Notably, EPGAT was the most robust approach in scenarios with
limited and imbalanced training data. Thus, the proposed approach offers a
powerful and effective way to identify essential genes and proteins.


  In this paper, we formulate a Nonlinear Model Predictive Control (NMPC) to
plan appropriate social distancing measures (and relaxations) in order to
mitigate the COVID-19 pandemic effects, considering the contagion development
in Brazil. The NMPC strategy is designed upon an adapted data-driven
Susceptible-Infected-Recovered-Deceased (SIRD) contagion model, which takes
into account the effects of social distancing. Furthermore, the adapted SIRD
model includes time-varying auto-regressive contagion parameters, which
dynamically converge according to the stage of the pandemic. This new model is
identified through a three-layered procedures, with analytical regressions,
Least-Squares optimization runs and auto-regressive model fits. The data-driven
model is validated and shown to adequately describe the contagion curves over
large forecast horizons. In this model, control input is defined as finitely
parametrized values for social distancing guidelines, which directly affect the
transmission and infection rates of the SARS-CoV-2 virus. The NMPC strategy
generates piece-wise constant quarantine guidelines which can be
relaxed/strengthen as each week passes. The implementation of the method is
pursued through a search mechanism, since the control is finitely parametrized
and, thus, there exist a finite number of possible control sequences.
Simulation essays are shown to illustrate the results obtained with the
proposed closed-loop NMPC strategy, which is able to mitigate the number of
infections and progressively loosen social distancing measures. With respect to
an "open-loop"/no control condition, the number of deaths still could be
reduced in up to 30 %. The forecast preview an infection peak to September 2nd,
2020, which could lead to over 1.5 million deaths if no coordinate health
policy is enacted. The framework serves as guidelines for possible public
health policies in Brazil.


  Although extensive behavioral changes often exist between closely related
animal species, our understanding of the genetic basis underlying the evolution
of behavior has remained limited. Here, we propose a new framework to study
behavioral evolution by computational estimation of ancestral behavioral
repertoires. We measured the behaviors of individuals from six species of fruit
flies using unsupervised techniques and identified suites of stereotyped
movements exhibited by each species. We then fit a Generalized Linear Mixed
Model to estimate the suites of behaviors exhibited by ancestral species, as
well as the intra- and inter-species behavioral covariances. We found that much
of intraspecific behavioral variation is explained by differences between
individuals in the status of their behavioral hidden states, what might be
called their "mood." Lastly, we propose a method to identify groups of
behaviors that appear to have evolved together, illustrating how sets of
behaviors, rather than individual behaviors, likely evolved. Our approach
provides a new framework for identifying co-evolving behaviors and may provide
new opportunities to study the genetic basis of behavioral evolution.


  Tree-child networks are a recently-described class of directed acyclic graphs
that have risen to prominence in phylogenetics (the study of evolutionary trees
and networks). Although these networks have a number of attractive mathematical
properties, many combinatorial questions concerning them remain intractable. In
this paper, we show that endowing these networks with a biologically relevant
ranking structure yields mathematically tractable objects, which we term ranked
tree-child networks (RTCNs). We explain how to derive exact and explicit
combinatorial results concerning the enumeration and generation of these
networks. We also explore probabilistic questions concerning the properties of
RTCNs when they are sampled uniformly at random. These questions include the
lengths of random walks between the root and leaves (both from the root to the
leaves and from a leaf to the root); the distribution of the number of cherries
in the network; and sampling RTCNs conditional on displaying a given tree. We
also formulate a conjecture regarding the scaling limit of the process that
counts the number of lineages in the ancestry of a leaf. The main idea in this
paper, namely using ranking as a way to achieve combinatorial tractability, may
also extend to other classes of networks.


  A multilayer network approach combines different network layers, which are
connected by interlayer edges, to create a single mathematical object. These
networks can contain a variety of information types and represent different
aspects of a system. However, the process for selecting which information to
include is not always straightforward. Using data on two agonistic behaviors in
a captive population of monk parakeets (Myiopsitta monachus), we developed a
framework for investigating how pooling or splitting behaviors at the scale of
dyadic relationships (between two individuals) affects individual- and
group-level social properties. We designed two reference models to test whether
randomizing the number of interactions across behavior types results in similar
structural patterns as the observed data. Although the behaviors were
correlated, the first reference model suggests that the two behaviors convey
different information about some social properties and should therefore not be
pooled. However, once we controlled for data sparsity, we found that the
observed measures corresponded with those from the second reference model.
Hence, our initial result may have been due to the unequal frequencies of each
behavior. Overall, our findings support pooling the two behaviors. Awareness of
how selected measurements can be affected by data properties is warranted, but
nonetheless our framework disentangles these efforts and as a result can be
used for myriad types of behaviors and questions. This framework will help
researchers make informed and data-driven decisions about which behaviors to
pool or separate, prior to using the data in subsequent multilayer network
analyses.


  The novel coronavirus disease 2019 COVID-19 has been leading the world into a
prominent crisis. As of May 19, 2020, the virus had spread to 215 countries
with more than 4,622,001 confirmed cases and 311,916 reported deaths worldwide,
including Algeria with 7201 cases and 555 deaths. Analyze and forecast COVID-19
cases and deaths growth could be useful in many ways, governments could
estimate medical equipment and take appropriate policy responses, and experts
could approximate the peak and the end of the disease. In this work, we first
train a time series Prophet model to analyze and forecast the number of
COVID-19 cases and deaths in Algeria based on the previously reported numbers.
Then, to better understand the spread and the properties of the COVID-19, we
include external factors that may contribute to accelerate/slow the spread of
the virus, construct a dataset from reliable sources, and conduct a large-scale
data analysis considering 82 countries worldwide. The evaluation results show
that the time series Prophet model accurately predicts the number of cases and
deaths in Algeria with low RMSE scores of 218.87 and 4.79 respectively, while
the forecast suggests that the total number of cases and deaths are expected to
increase in the coming weeks. Moreover, the worldwide data-driven analysis
reveals several correlations between the increase/decrease in the number of
cases and deaths and external factors that may contribute to accelerate/slow
the spread of the virus such as geographic, climatic, health, economic, and
demographic factors.


  COVID-19, a viral respiratory pandemic, has rapidly spread throughout the
globe. Large scale and rapid testing of the population is required to contain
the disease, but such testing is prohibitive in terms of resources, cost and
time. Recently RT-PCR based pooled testing has emerged as a promising way to
boost testing efficiency. We introduce a stepped pooled testing strategy, a
probability driven approach which significantly reduces the number of tests
required to identify infected individuals in a large population. Our
comprehensive methodology incorporates the effect of false negative and
positive rates to accurately determine not only the efficiency of pooling but
also it's accuracy. Under various plausible scenarios, we show that this
approach significantly reduces the cost of testing and also reduces the
effective false positive rate of tests when compared to a strategy of testing
every individual of a population. We also outline an optimization strategy to
obtain the pool size that maximizes the efficiency of pooling given the
diagnostic protocol parameters and local infection conditions.


  As the interactions between people increases, the impending menace of
COVID-19 outbreaks materialize, and there is an inclination to apply lockdowns.
In this context, it is essential to have easy-to-use indicators for people to
use as a reference. The basic reproduction number of confirmed positives,
$R_t$, fulfill such a role. This document proposes a data-driven approach to
nowcast $R_t$ based on previous observations' statistical behavior. As more
information arrives, the method naturally becomes more precise about the final
count of confirmed positives. Our method's strength is that it is based on the
self-reported onset of symptoms, in contrast to other methods that use the
daily report's count to infer this quantity. We show that our approach may be
the foundation for determining useful epidemy tracking indicators.


  Without contributing, defectors take more benefit from social resources than
cooperators which is the reflection of a specific character of individuals.
However, natural physical mechanisms of our society promote cooperation. Thus,
in the long run, the evolution about genetic variation is something more than
the social evolution about fitness. The loci of evolutionary paths of the
cooperation and the defection are correlated, but not a full complement of each
other. Yet, the only single specific mechanism which is operated by some rules
explains the enhancement of cooperation where the independent analysis of
defect evolutionary mechanism is ignored. Moreover, the execution of a
particular evolutionary rule through algorithm method over the long time
encounters highly sensitive influence of the model parameters. Theoretically,
biodiversity of two types relatively persists rarely. Here I describe the
evolutionary outcome in the demographic fluctuation. Using both analytical
procedure and algorithm method the article concludes that the intratype fitness
of individual species is the key factor for not only surviving, but thriving.
In consideration of the random drift, the experimental outcomes show that
dominant enhancement of cooperation over defection is qualitatively independent
of environmental scenario. Collectively, the set of the rules becomes an
evolutionary principle to cooperation enhancement.


  Populations are made up of an integer number of individuals and are subject
to stochastic birth-death processes whose rates may vary in time. Useful
quantities, like the chance of ultimate fixation, satisfy an appropriate
difference (master) equation, but closed-form solutions of these equations are
rare. Analytical insights in fields like population genetics, ecology and
evolution rely, almost exclusively, on an uncontrolled application of the
diffusion approximation (DA) which assumes the smoothness of the relevant
quantities over the set of integers. Here we combine asymptotic matching
techniques with a first-order (controlling-factor) WKB method to obtain a
theory whose range of applicability is much wider. This allows us to rederive
DA from a more general theory, to identify its limitations, and to suggest
alternative analytical solutions and scalable numerical techniques when it
fails. We carry out our analysis for the calculation of the fixation
probability in a fluctuating environment, highlighting the difference between
(on average) deleterious and beneficial mutant invasion and the intricate
distinction between weak and strong selection.


  The diffusion approximation (DA) is widely used in the analysis of stochastic
population dynamics, from population genetics to ecology and evolution. DA is
an uncontrolled approximation that assumes the smoothness of the calculated
quantity over the relevant state space and fails when this property is not
satisfied. This failure becomes severe in situations where the direction of
selection switches sign. Here we employ the WKB (large-deviations) method,
which requires only the logarithm of a given quantity to be smooth over its
state space. Combining the WKB scheme with asymptotic matching techniques, we
show how to derive the diffusion approximation in a controlled manner and how
to produce better approximations, applicable for much wider regimes of
parameters. We also introduce a scalable (independent of population size)
WKB-based numerical technique. The method is applied to a central problem in
population genetics and evolution, finding the chance of ultimate fixation in a
zero-sum, two-types competition.


  We consider random perturbations of a general eco-epidemiological model. We
prove the existence of a global random attractor, the persistence of
susceptibles preys and provide conditions for the simultaneous extinction of
infectives and predators. We also discuss the dynamics of the corresponding
random epidemiological $SI$ and predator-prey models. We obtain for this cases
a global random attractor, prove the prevalence of susceptibles/preys and
provide conditions for the extinctions of infectives/predators.


  Punctuated equilibrium is a mode of evolution in which phenetic change occurs
in rapid bursts that are separated by much longer intervals of stasis during
which mutations accumulate but no major phenotypic change occurs. Punctuated
equilibrium has been originally proposed within the framework of paleobiology,
to explain the lack of transitional forms that is typical of the fossil record.
Theoretically, punctuated equilibrium has been linked to self-organized
criticality (SOC), a model in which the size of avalanches in an evolving
system is power-law distributed, resulting in increasing rarity of major
events. We show here that, under the weak-mutation limit, a large population
would spend most of the time in stasis in the vicinity of saddle points in the
fitness landscape. The periods of stasis are punctuated by fast transitions, in
lnNe time (Ne, effective population size), when a new beneficial mutation is
fixed in the evolving population, which moves to a different saddle, or on much
rarer occasions, from a saddle to a local peak. Thus, punctuated equilibrium is
the default mode of evolution under a simple model that does not involve SOC or
other special conditions.


  Using Mickens nonstandard method, we obtain a discrete family of
nonautonomous eco-epidemiological models that include general functions
corresponding to the predation of the infected and uninfected preys. We obtain
results on the persistence and extinction of the infected preys assuming that
the bi-dimensional predator-prey subsystem that describes the dynamics in the
absence of the infection satisfies some assumptions. Some examples and
simulations are undertaken to illustrate our results.


  Knowing COVID-19 epidemiological distributions, such as the time from patient
admission to death, is directly relevant to effective primary and secondary
care planning, and moreover, the mathematical modelling of the pandemic
generally. We determine epidemiological distributions for patients hospitalised
with COVID-19 using a large dataset ($N=21{,}000-157{,}000$) from the Brazilian
Sistema de Informa\c{c}\~ao de Vigil\^ancia Epidemiol\'ogica da Gripe database.
A joint Bayesian subnational model with partial pooling is used to
simultaneously describe the 26 states and one federal district of Brazil, and
shows significant variation in the mean of the symptom-onset-to-death time,
with ranges between 11.2-17.8 days across the different states, and a mean of
15.2 days for Brazil. We find strong evidence in favour of specific probability
density function choices: for example, the gamma distribution gives the best
fit for onset-to-death and the generalised log-normal for
onset-to-hospital-admission. Our results show that epidemiological
distributions have considerable geographical variation, and provide the first
estimates of these distributions in a low and middle-income setting. At the
subnational level, variation in COVID-19 outcome timings are found to be
correlated with poverty, deprivation and segregation levels, and weaker
correlation is observed for mean age, wealth and urbanicity.


  While the SARS-CoV-2 keeps spreading world-wide, comparing its evolution
across different nations is a timely challenge of both theoretical and
practical importance. The large variety of dissimilar and country-dependent
epidemiological factors, in fact, makes extremely difficult to understand their
influence on the epidemic trends within a unique and coherent framework. We
present a geometric framework to characterize, in an integrated and
low-dimensional fashion, the epidemic plume-like trajectories traced by the
infection rate, $I$, and the fatality rate, $D$, in the $(I,D)$ plane. Our
analysis enables the definition of an epidemiometric system based on three
geometric observables rating the SARS-CoV-2 pandemic events via scales
analogous to those for the magnitude and the intensity of seismic events. Being
exquisitely geometric, our framework can be applied to classify other epidemic
data and secondary waves, raising the possibility of designing epidemic alerts
or early warning systems to enhance public and governmental responses to a
rapidly emerging outbreak.


  DNA N6-methylation (6mA) in Adenine nucleotide is a post replication
modification and is responsible for many biological functions. Experimental
methods for genome wide 6mA site detection is an expensive and manual labour
intensive process. Automated and accurate computational methods can help to
identify 6mA sites in long genomes saving significant time and money. Our study
develops a convolutional neural network based tool i6mA-CNN capable of
identifying 6mA sites in the rice genome. Our model coordinates among multiple
types of features such as PseAAC inspired customized feature vector, multiple
one hot representations and dinucleotide physicochemical properties. It
achieves area under the receiver operating characteristic curve of 0.98 with an
overall accuracy of 0.94 using 5 fold cross validation on benchmark dataset.
Finally, we evaluate our model on two other plant genome 6mA site
identification datasets besides rice. Results suggest that our proposed tool is
able to generalize its ability of 6mA site identification on plant genomes
irrespective of plant species. Web tool for this research can be found at:
https://cutt.ly/Co6KuWG. Supplementary data (benchmark dataset, independent
test dataset, comparison purpose dataset, trained model, physicochemical
property values, attention mechanism details for motif finding) are available
at https://cutt.ly/PpDdeDH.


  This article contains a series of analyses done for the SARS-CoV-2 outbreak
in Rio Grande do Sul (RS) in the south of Brazil. These analyses are focused on
the high-incidence cities such as the state capital Porto Alegre and at the
state level. We provide methodological details and estimates for the effective
reproduction number $R_t$, a joint analysis of the mobility data together with
the estimated $R_t$ as well as ICU simulations and ICU LoS (length of stay)
estimation for hospitalizations in Porto Alegre/RS.


  The emergence of self-sustaining autocatalytic networks in chemical reaction
systems has been studied as a possible mechanism for modelling how living
systems first arose. It has been known for several decades that such networks
will form within systems of polymers (under cleavage and ligation reactions)
under a simple process of random catalysis, and this process has since been
mathematically analysed. In this paper, we provide an exact expression for the
expected number of self-sustaining autocatalytic networks that will form in a
general chemical reaction system, and the expected number of these networks
that will also be uninhibited (by some molecule produced by the system). Using
these equations, we are able to describe the patterns of catalysis and
inhibition that maximise or minimise the expected number of such networks. We
apply our results to derive a general theorem concerning the trade-off between
catalysis and inhibition, and to provide some insight into the extent to which
the expected number of self-sustaining autocatalytic networks coincides with
the probability that at least one such system is present.


  Over the past two decades there has been a number of global outbreaks of
viral diseases. This has accelerated the efforts to model and forecast the
disease spreading, in order to find ways to confine the spreading regionally
and between regions. Towards this we have devised a model of geographical
spreading of viral infections due to human spatial mobility and adapted it to
the latest COVID-19 pandemic. In this the region to be modelled is overlaid
with a two-dimensional grid weighted with the population density defined cells,
in each of which a compartmental SEIRS system of delay difference equations
simulate the local dynamics (microdynamics) of the disease. The infections
between cells are stochastic and allow for the geographical spreading of the
virus over the two-dimensional space (macrodynamics). This approach allows to
separate the parameters related to the biological aspects of the disease from
the ones that represent the spatial contagious behaviour through different
kinds of mobility of people acting as virus carriers. These provide sufficient
information to trace the evolution of the pandemic in different situations. In
particular we have applied this approach to three in many ways different
countries, Mexico, Finland and Iceland and found that the model is capable of
reproducing and predicting the stochastic global path of the pandemic. This
study sheds light on how the diverse cultural and socioeconomic aspects of a
country influence the evolution of the epidemics and also the efficacy of
social distancing and other confinement measures.


  Disease-gene association through Genome-wide association study (GWAS) is an
arduous task for researchers. Investigating single nucleotide polymorphisms
(SNPs) that correlate with specific diseases needs statistical analysis of
associations. Considering the huge number of possible mutations, in addition to
its high cost, another important drawback of GWAS analysis is the large number
of false-positives. Thus, researchers search for more evidence to cross-check
their results through different sources. To provide the researchers with
alternative low-cost disease-gene association evidence, computational
approaches come into play. Since molecular networks are able to capture complex
interplay among molecules in diseases, they become one of the most extensively
used data for disease-gene association prediction. In this survey, we aim to
provide a comprehensive and an up-to-date review of network-based methods for
disease gene prediction. We also conduct an empirical analysis on 14
state-of-the-art methods. To summarize, we first elucidate the task definition
for disease gene prediction. Secondly, we categorize existing network-based
efforts into network diffusion methods, traditional machine learning methods
with handcrafted graph features and graph representation learning methods.
Thirdly, an empirical analysis is conducted to evaluate the performance of the
selected methods across seven diseases. We also provide distinguishing findings
about the discussed methods based on our empirical analysis. Finally, we
highlight potential research directions for future studies on disease gene
prediction.


  COVID-19 is a new pandemic disease that is affecting almost every country
with a negative impact on social life and economic activities. The number of
infected and deceased patients continues to increase globally. Mathematical
models can help in developing better strategies to contain a pandemic.
Considering multiple measures taken by African governments and challenging
socio-economic factors, simple models cannot fit the data. We studied the
dynamical evolution of COVID-19 in selected African countries. We derived a
time-dependent reproduction number for each country studied to offer further
insights into the spread of COVID-19 in Africa.


  The outbreaks of Coronavirus Disease 2019 (COVID-19) have impacted the world
significantly. Modeling the trend of infection and real-time forecasting of
cases can help decision making and control of the disease spread. However,
data-driven methods such as recurrent neural networks (RNN) can perform poorly
due to limited daily samples in time. In this work, we develop an integrated
spatiotemporal model based on the epidemic differential equations (SIR) and
RNN. The former after simplification and discretization is a compact model of
temporal infection trend of a region while the latter models the effect of
nearest neighboring regions. The latter captures latent spatial information.
%that is not publicly reported. We trained and tested our model on COVID-19
data in Italy, and show that it out-performs existing temporal models (fully
connected NN, SIR, ARIMA) in 1-day, 3-day, and 1-week ahead forecasting
especially in the regime of limited training data.


  The novel coronavirus disease (COVID-19) is a public health problem once
according to the World Health Organization up to June 10th, 2020, more than 7.1
million people were infected, and more than 400 thousand have died worldwide.
In the current scenario, the Brazil and the United States of America present a
high daily incidence of new cases and deaths. It is important to forecast the
number of new cases in a time window of one week, once this can help the public
health system developing strategic planning to deals with the COVID-19. In this
paper, Bayesian regression neural network, cubist regression, k-nearest
neighbors, quantile random forest, and support vector regression, are used
stand-alone, and coupled with the recent pre-processing variational mode
decomposition (VMD) employed to decompose the time series into several
intrinsic mode functions. All Artificial Intelligence techniques are evaluated
in the task of time-series forecasting with one, three, and six-days-ahead the
cumulative COVID-19 cases in five Brazilian and American states up to April
28th, 2020. Previous cumulative COVID-19 cases and exogenous variables as daily
temperature and precipitation were employed as inputs for all forecasting
models. The hybridization of VMD outperformed single forecasting models
regarding the accuracy, specifically when the horizon is six-days-ahead,
achieving better accuracy in 70% of the cases. Regarding the exogenous
variables, the importance ranking as predictor variables is past cases,
temperature, and precipitation. Due to the efficiency of evaluated models to
forecasting cumulative COVID-19 cases up to six-days-ahead, the adopted models
can be recommended as a promising models for forecasting and be used to assist
in the development of public policies to mitigate the effects of COVID-19
outbreak.


  The adaptive immune system responds to pathogens by selecting clones of cells
with specific receptors. While clonal selection in response to particular
antigens has been studied in detail, it is unknown how a lifetime of exposures
to many antigens collectively shape the immune repertoire. Here, through
mathematical modeling and statistical analyses of T cell receptor sequencing
data we demonstrate that clonal expansions during a perinatal time window leave
a long-lasting imprint on the human T cell repertoire. We demonstrate how the
empirical scaling law relating the rank of the largest clones to their size can
emerge from clonal growth during repertoire formation. We statistically
identify early founded clones and find that they are indeed highly enriched
among the largest clones. This enrichment persists even after decades of human
aging, in a way that is quantitatively predicted by a model of fluctuating
clonal selection. Our work presents a quantitative theory of human T cell
dynamics compatible with the statistical laws of repertoire organization and
provides a mechanism for how early clonal dynamics imprint the hierarchy of T
cell clone sizes with implications for pathogen defense and autoimmunity.


  Several non-pharmaceutical interventions have been proposed to control the
spread of the COVID-19 pandemic. On the large scale, these empirical solutions,
often associated with extended and complete lockdowns, attempt to minimize the
costs associated with mortality, economic losses and social factors, while
being subject to constraints such as finite hospital capacity. Here we pose the
question of how to mitigate pandemic costs subject to constraints by adopting
the language of optimal control theory. This allows us to determine top-down
policies for the nature and dynamics of social contact rates given an
age-structured model for the dynamics of the disease. Depending on the relative
weights allocated to life and socioeconomic losses, we see that the optimal
strategies range from long-term social-distancing only for the most vulnerable,
to partial lockdown to ensure not over-running hospitals, to alternating-shifts
with significant reduction in life and/or socioeconomic losses. Crucially,
commonly used strategies that involve long periods of broad lockdown are almost
never optimal, as they are highly unstable to reopening and entail high
socioeconomic costs. Using parameter estimates from data available for Germany
and the USA, we quantify these policies and use sensitivity analysis in the
relevant model parameters and initial conditions to determine the range of
robustness of our policies. Finally we also discuss how bottom-up behavioral
changes can also change the dynamics of the pandemic and show how this in
tandem with top-down control policies can mitigate pandemic costs even more
effectively.


  To reduce the biases of traditional survey-based methods, this paper proposes
an epidemic model-based approach to inference the incubation period
distribution of COVID-19 utilizing the publicly reported confirmed case number.
We construct an epidemic model, namely SEAIR, and take advantage of the dynamic
transmission process depicted by SEAIR to estimate the onset probability in
each day of exposed individuals in eight impacted countries. Based on these
estimations, the general incubation probability distribution of COVID-19 has
been revealed. The proposed method can avoid several biases of traditional
survey-based methods. However, due to the mathematical-model-based nature of
this method, the inference results are somewhat sensitive to the setting of
parameters. Therefore, this method should be practiced reasonably on the basis
of a certain understanding of the studied epidemic.


  We generalise the epidemic Renormalisation Group framework while connecting
it to a SIR model with time-dependent coefficients. We then confront the model
with COVID-19 in Denmark, Germany, Italy and France and show that the approach
works rather well in reproducing the data. We also show that a better
understanding of the time dependence of the recovery rate would require
extending the model to take into account the number of deaths whenever these
are over 15% of the total number of infected cases.


  In this study, we compare the charge transport properties of multiple (double
stranded) dsRNA sequences with corresponding dsDNA sequences. Recent studies
have presented a contradictory picture of relative charge transport
efficiencies in A-form DNA:RNA hybrids and dsDNA. Using a multiscale modelling
framework, we compute conductance of dsDNA and dsRNA using Landauer formalism
in coherent limit and Marcus-Hush theory in the incoherent limit. We find that
dsDNA conducts better than dsRNA in both the charge transport regimes. Our
analysis shows that the structural differences in the twist angle and slide of
dsDNA and dsRNA are the main reasons behind the higher conductance of dsDNA in
the incoherent hopping regime. In the coherent limit however, for the same base
pair length, the conductance of dsRNA is higher than that of dsDNA for the
morphologies where dsRNA has smaller end-to-end length relative to that of
dsDNA.


  The advent of the COVID-19 pandemic has instigated unprecedented changes in
many countries around the globe, putting a significant burden on the health
sectors, affecting the macro economic conditions, and altering social
interactions amongst the population. In response, the academic community has
produced multiple forecasting models, approaches and algorithms to best predict
the different indicators of COVID-19, such as the number of confirmed infected
cases. Yet, researchers had little to no historical information about the
pandemic at their disposal in order to inform their forecasting methods. Our
work studies the predictive performance of models at various stages of the
pandemic to better understand their fundamental uncertainty and the impact of
data availability on such forecasts. We use historical data of COVID-19
infections from 253 regions from the period of 22nd January 2020 until 22nd
June 2020 to predict, through a rolling window backtesting framework, the
cumulative number of infected cases for the next 7 and 28 days. We implement
three simple models to track the root mean squared logarithmic error in this
6-month span, a baseline model that always predicts the last known value of the
cumulative confirmed cases, a power growth model and an epidemiological model
called SEIRD. Prediction errors are substantially higher in early stages of the
pandemic, resulting from limited data. Throughout the course of the pandemic,
errors regress slowly, but steadily. The more confirmed cases a country
exhibits at any point in time, the lower the error in forecasting future
confirmed cases. We emphasize the significance of having a rigorous backtesting
framework to accurately assess the predictive power of such models at any point
in time during the outbreak which in turn can be used to assign the right level
of certainty to these forecasts and facilitate better planning.


  Coarse graining enables the investigation of molecular dynamics for larger
systems and at longer timescales than is possible at atomic resolution.
However, a coarse graining model must be formulated such that the conclusions
we draw from it are consistent with the conclusions we would draw from a model
at a finer level of detail. It has been proven that a force matching scheme
defines a thermodynamically consistent coarse-grained model for an atomistic
system in the variational limit. Wang et al. [ACS Cent. Sci. 5, 755 (2019)]
demonstrated that the existence of such a variational limit enables the use of
a supervised machine learning framework to generate a coarse-grained force
field, which can then be used for simulation in the coarse-grained space. Their
framework, however, requires the manual input of molecular features upon which
to machine learn the force field. In the present contribution, we build upon
the advance of Wang et al.and introduce a hybrid architecture for the machine
learning of coarse-grained force fields that learns their own features via a
subnetwork that leverages continuous filter convolutions on a graph neural
network architecture. We demonstrate that this framework succeeds at
reproducing the thermodynamics for small biomolecular systems. Since the
learned molecular representations are inherently transferable, the architecture
presented here sets the stage for the development of machine-learned,
coarse-grained force fields that are transferable across molecular systems.


  In this paper, we continue development of the new epidemiological model HIT,
which is suitable for analyzing and predicting the propagation of COVID-19
epidemics. This is a discrete-time model allowing a reconstruction of the
dynamics of asymptomatic virus holders using the available daily statistics on
the number of new cases. We suggest to use a new indicator, the total infection
rate, to distinguish the propagation and recession modes of the epidemic. We
check our indicator on the available data for eleven different countries and
for the whole world. Our reconstructions are very precise. In several cases, we
are able to detect the exact dates of the disastrous political decisions,
ensuring the second wave of the epidemics. It appears that for all our examples
the decisions made on the basis of the current number of new cases are wrong.
In this paper, we suggest a reasonable alternative. Our analysis shows that all
tested countries are in a dangerous zone except Sweden.


  The noise in daily infection counts of an epidemic should be super-Poissonian
due to intrinsic epidemiological and administrative clustering. Here, we use
this clustering to classify the official national SARS-CoV-2 daily infection
counts and check for infection counts that are unusually anti-clustered. We
adopt a one-parameter model of $\phi'_i$ infections per cluster, dividing any
daily count $n_i$ into $n_i/\phi'_i$ 'clusters', for 'country' $i$. We assume
that $n_i/\phi'_i$ on a given day $j$ is drawn from a Poisson distribution
whose mean is robustly estimated from the four neighbouring days, and calculate
the inferred Poisson probability $P'_{ij}$ of the observation. The $P'_{ij}$
values should be uniformly distributed. We find the value $\phi_i$ that
minimises the Kolmogorov-Smirnov distance from a uniform distribution. We
investigate the $(\phi_i, N_i)$ distribution, for total infection count $N_i$.
We find that most of the daily infection count sequences are inconsistent with
a Poissonian model. Most are found to be consistent with the $\phi_i$ model.
The 28-, 14- and 7-day least noisy sequences for several countries are best
modelled as sub-Poissonian, suggesting a distinct epidemiological family. The
28-day least noisy sequence of Algeria has a preferred model that is strongly
sub-Poissonian, with $\phi_i^{28} < 0.1$. TJ, TR, RU, BY, AL, AE, and NI have
preferred models that are also sub-Poissonian, with $\phi_i^{28} < 0.5$. A
statistically significant ($P^{\tau} < 0.05$) correlation was found between the
lack of media freedom in a country, as represented by a high Reporters sans
frontieres Press Freedom Index (PFI$^{2020}$), and the lack of statistical
noise in the country's daily counts. The $\phi_i$ model appears to be an
effective detector of suspiciously low statistical noise in the national
SARS-CoV-2 daily infection counts.


  We investigate the global dynamics of a renewal-type epidemic model with
variable susceptibility. We show that in this extended model there exists a
unique endemic equilibrium and prove that it is globally asymptotically stable
when $R_0 > 1$, i.e. when it exists. We also show that the infection-free
equilibrium, which exists always, is globally asymptotically stable for $R_0
\leq 1$.