  The SARS-CoV-2 pandemic has caused significant mortality and morbidity
worldwide, sparing almost no community. As the disease will likely remain a
threat for years to come, an understanding of the precise influences of human
demographics and settlement, as well as the dynamic factors of climate,
susceptible depletion, and intervention, on the spread of localized epidemics
will be vital for mounting an effective response. We consider the entire set of
local epidemics in the United States; a broad selection of demographic,
population density, and climate factors; and local mobility data, tracking
social distancing interventions, to determine the key factors driving the
spread and containment of the virus. Assuming first a linear model for the rate
of exponential growth (or decay) in cases/mortality, we find that
population-weighted density, humidity, and median age dominate the dynamics of
growth and decline, once interventions are accounted for. A focus on distinct
metropolitan areas suggests that some locales benefited from the timing of a
nearly simultaneous nationwide shutdown, and/or the regional climate conditions
in mid-March; while others suffered significant outbreaks prior to
intervention. Using a first-principles model of the infection spread, we then
develop predictions for the impact of the relaxation of social distancing and
local climate conditions. A few regions, where a significant fraction of the
population was infected, show evidence that the epidemic has partially resolved
via depletion of the susceptible population (i.e., "herd immunity"), while most
regions in the United States remain overwhelmingly susceptible. These results
will be important for optimal management of intervention strategies, which can
be facilitated using our online dashboard.


  Zhang et al. (2020) used linear regression to quantify the effect of
lockdowns on the number of cases of COVID-19. We show using differential
equations from the susceptible-exposed-infected-recovered (SEIR) model and with
an example from another location not previously considered that the Zhang et
al. analysis should not be considered sound evidence that mask mandates are
sufficient to control or the primary factor controlling the spread of COVID-19.


  Raw moments are used as a way to estimate species abundance distribution. The
almost linear pattern of the log transformation of raw moments across scales
allow us to extrapolate species abundance distribution for larger areas.
However, results may produce errors. Some of these errors are due to
computational complexity, fittings of patterns, binning methods, and so on. We
provide some methods to reduce some of the errors. The main result is
introducing new techniques for evaluating a more accurate species abundance
distributions across scales through moments across scales.


  In many biological systems, chemical reactions or changes in a physical state
are assumed to occur instantaneously. For describing the dynamics of those
systems, Markov models that require exponentially distributed inter-event times
have been used widely. However, some biophysical processes such as gene
transcription and translation are known to have a significant gap between the
initiation and the completion of the processes, which renders the usual
assumption of exponential distribution untenable. In this paper, we consider
relaxing this assumption by incorporating age-dependent random time delays into
the system dynamics. We do so by constructing a measure-valued Markov process
on a more abstract state space, which allows us to keep track of the "ages" of
molecules participating in a chemical reaction.
  We study the large-volume limit of such age-structured systems. We show that,
when appropriately scaled, the stochastic system can be approximated by a
system of Partial Differential Equations (PDEs) in the large-volume limit, as
opposed to Ordinary Differential Equations (ODEs) in the classical theory. We
show how the limiting PDE system can be used for the purpose of further model
reductions and for devising efficient simulation algorithms. In order to
describe the ideas, we use a simple transcription process as a running example.
We, however, note that the methods developed in this paper apply to a wide
class of biophysical systems.


  The recent appearance of COVID-19 virus has created a global crisis due to
unavailability of any vaccine or drug that can effectively and
deterministically work against it. Naturally, different possibilities
(including herbal medicines having known therapeutic significance) have been
explored by the scientists. The systematic scientific study (beginning with in
silico study) of herbal medicines in particular and any drug in general is now
possible as the structural components (proteins) of COVID-19 are already
characterized. The main protease of COVID-19 virus is $\rm{M^{pro}}$ or
$\rm{3CL^{pro}}$ which is a key CoV enzyme and an attractive drug target as it
plays a pivotal role in mediating viral replication and transcription. In the
present study, $\rm{3CL^{pro}}$ is used to study drug:3CLpro interactions and
thus to investigate whether all or any of the main chemical constituents of
Tinospora cordifolia (e.g., berberine $\rm{(C_{20}H_{18}NO_{4})}$,
$\beta$-sitosterol $\rm{(C_{29}H_{50}O)}$, choline $\rm{(C_{5}H_{14}NO)}$,
tetrahydropalmatine $\rm{(C_{21}H_{25}NO_{4})}$ and octacosanol
$\rm{(C_{28}H_{58}O))}$ can be used as an anti-viral drug against SARS-CoV-2.
The in silico study performed using tools of network pharmacology, molecular
docking including molecular dynamics have revealed that among all considered
phytochemicals in Tinospora cordifolia, berberine can regulate $\rm{3CL^{pro}}$
protein's function due to its easy inhibition and thus can control viral
replication. The selection of Tinospora cordifolia was motivated by the fact
that the main constituents of it are known to be responsible for various
antiviral activities and the treatment of jaundice, rheumatism, diabetes, etc.


  Non-pharmaceutical interventions (NPIs) have been crucial in curbing COVID-19
in the United States (US). Consequently, relaxing NPIs through a phased
re-opening of the US amid still-high levels of COVID-19 susceptibility could
lead to new epidemic waves. This calls for a COVID-19 early warning system.
Here we evaluate multiple digital data streams as early warning indicators of
increasing or decreasing state-level US COVID-19 activity between January and
June 2020. We estimate the timing of sharp changes in each data stream using a
simple Bayesian model that calculates in near real-time the probability of
exponential growth or decay. Analysis of COVID-19-related activity on social
network microblogs, Internet searches, point-of-care medical software, and a
metapopulation mechanistic model, as well as fever anomalies captured by smart
thermometer networks, shows exponential growth roughly 2-3 weeks prior to
comparable growth in confirmed COVID-19 cases and 3-4 weeks prior to comparable
growth in COVID-19 deaths across the US over the last 6 months. We further
observe exponential decay in confirmed cases and deaths 5-6 weeks after
implementation of NPIs, as measured by anonymized and aggregated human mobility
data from mobile phones. Finally, we propose a combined indicator for
exponential growth in multiple data streams that may aid in developing an early
warning system for future COVID-19 outbreaks. These efforts represent an
initial exploratory framework, and both continued study of the predictive power
of digital indicators as well as further development of the statistical
approach are needed.


  We extend the classical SIR model of infectious disease spread to account for
time dependence in the parameters, which also include diffusivities. The
temporal dependence accounts for the changing characteristics of testing,
quarantine and treatment protocols, while diffusivity incorporates a mobile
population. This model has been applied to data on the evolution of the
COVID-19 pandemic in the US state of Michigan. For system inference, we use
recent advances; specifically our framework for Variational System
Identification (Wang et al., Comp. Meth. App. Mech. Eng., 356, 44-74, 2019;
arXiv:2001.04816 [cs.CE]) as well as Bayesian machine learning methods.


  Motivation: Uncovering the genomic causes of cancer, known as cancer driver
genes, is a fundamental task in biomedical research. Cancer driver genes drive
the development and progression of cancer, thus identifying cancer driver genes
and their regulatory mechanism is crucial to the design of cancer treatment and
intervention. Many computational methods, which take the advantages of computer
science and data science, have been developed to utilise multiple types of
genomic data to reveal cancer drivers and their regulatory mechanism behind
cancer development and progression. Due to the complexity of the mechanistic
insight of cancer genes in driving cancer and the fast development of the
field, it is necessary to have a comprehensive review about the current
computational methods for discovering different types of cancer drivers.
Results: We survey computational methods for identifying cancer drivers from
genomic data. We categorise the methods into three groups, methods for single
driver identification, methods for driver module identification, and methods
for identifying personalised cancer drivers. We also conduct a case study to
compare the performance of the current methods. We further analyse the
advantages and limitations of the current methods, and discuss the challenges
and future directions of the topic. In addition, we investigate the resources
for discovering and validating cancer drivers in order to provide a one-stop
reference of the tools to facilitate cancer driver discovery. The ultimate goal
of the paper is to help those interested in the topic to establish a solid
background to carry out further research in the field.


  Visual Molecular Dynamics (VMD) is one of the most widely used molecular
graphics software in the community of theoretical simulations. So far, however,
it still lacks a graphical user interface (GUI) for molecular manipulations
when doing some modeling tasks. For instance, translation or rotation of a
selected molecule(s) or part(s) of a molecule, which are currently only can be
achieved using tcl scripts. Here, we use tcl script develop a user-friendly GUI
for VMD, named Molcontroller, which is featured by allowing users to quickly
and conveniently perform various molecular manipulations. This GUI might be
helpful for improving the modeling efficiency of VMD users.


  Stochastic epidemic models on networks are inherently high-dimensional and
the resulting exact models are intractable numerically even for modest network
sizes. Mean-field models provide an alternative but can only capture average
quantities, thus offering little or no information about variability in the
outcome of the exact process. In this paper we conjecture and numerically prove
that it is possible to construct PDE-limits of the exact stochastic SIS
epidemics on regular and Erd\H{o}s-R\'enyi networks. To do this we first
approximate the exact stochastic process at population level by a
Birth-and-Death process (BD) (with a state space of $O(N)$ rather than
$O(2^N)$) whose coefficients are determined numerically from Gillespie
simulations of the exact epidemic on explicit networks. We numerically
demonstrate that the coefficients of the resulting BD process are
density-dependent, a crucial condition for the existence of a PDE limit.
Extensive numerical tests for Regular and Erd\H{o}s-R\'enyi networks show
excellent agreement between the outcome of simulations and the numerical
solution of the Fokker-Planck equations. Apart from a significant reduction in
dimensionality, the PDE also provides the means to derive the epidemic outbreak
threshold linking network and disease dynamics parameters, albeit in an
implicit way. Perhaps more importantly, it enables the formulation and
numerical evaluation of likelihoods for epidemic and network inference as
illustrated in a worked out example.


  Several formulations are describing the pattern of species-area relationship,
log-log linear, semi-log linear, among others. These patterns mainly explain
the species-area relationship for large areas, and for the small area, they
provide significant differences from real data. We consider the geometric
position of individuals of species, and base on that, we find the probability
of observing at least one individual of the species. We apply a translation of
the well-studied problem of mixed salt-water in a tank to describe the formula
of SAR. For a rectangular sample area the species-area relationship follows the
pattern, with some simplification, $S=c|A^{\beta}+a|^z$, where $S$ is the
number of species in the area of size $A$ and $a,c,z$, and $\beta$ are
constants with $z<1$ and $\beta\leq1$. We also show how the constant $z$
relates to some macroecological patterns, namely spatial aggregation,
percentage of area coverage, and the core-satellite model. We exemplify our
method using data on tropical tree species from a 50ha plot in Barro Colorado
Island (BCI), Panama, using all individuals.


  For most diseases, building large databases of labeled genetic data is an
expensive and time-demanding task. To address this, we introduce genetic
Generative Adversarial Networks (gGAN), a semi-supervised approach based on an
innovative GAN architecture to create large synthetic genetic data sets
starting with a small amount of labeled data and a large amount of unlabeled
data. Our goal is to determine the propensity of a new individual to develop
the severe form of the illness from their genetic profile alone. The proposed
model achieved satisfactory results using real genetic data from different
datasets and populations, in which the test populations may not have the same
genetic profiles. The proposed model is self-aware and capable of determining
whether a new genetic profile has enough compatibility with the data on which
the network was trained and is thus suitable for prediction. The code and
datasets used can be found at https://github.com/caio-davi/gGAN.


  Since the onset of the COVID-19 outbreak in Wuhan, China, numerous
forecasting models have been proposed to project the trajectory of coronavirus
infection cases. We propose a new discrete-time Markov chain transition matrix
model that directly incorporates stochastic behavior and for which parameter
estimation is straightforward from available data. Using such data from China's
Hubei province (for which Wuhan is the provincial capital city), the model is
shown to be flexible, robust, and accurate. As a result, it has been adopted by
the first Shanghai assistance medical team in Wuhan's Jinyintan Hospital, which
was the first designated hospital to take COVID-19 patients in the world. The
forecast has been used for preparing medical staff, intensive care unit (ICU)
beds, ventilators, and other critical care medical resources and for supporting
real-time medical management decisions. Empirical data from China's first two
months (January/February) of fighting COVID-19 was collected and used to
enhance the model by embedding NPI efficiency into the model. We applied the
model to forecast Italy, South Korea, and Iran on March 9. Later we made
forecasts for Spain, Germany, France, US on March 24. Again, the model has
performed very well, proven to be flexible, robust, and accurate for most of
these countries/regions outside China.


  Bacteria evolve resistance to antibiotics by a multitude of mechanisms. A
central, yet unsolved question is how resistance evolution affects cell growth
at different drug levels. Here we develop a fitness model that predicts growth
rates of common resistance mutants from their effects on cell metabolism. We
map metabolic effects of resistance mutations in drug-free environments and
under drug challenge; the resulting fitness trade-off defines a Pareto surface
of resistance evolution. We predict evolutionary trajectories of
dosage-dependent growth rates and resistance levels, as well as the prevalent
resistance mechanism depending on drug and nutrient levels. These predictions
are confirmed by empirical growth curves and genomic data of E. coli
populations. Our results show that resistance evolution, by coupling major
metabolic pathways, is strongly intertwined with systems biology and ecology of
microbial populations.


  The outbreak of Covid-19 led the world to an unprecedent health and
economical crisis. In an attempt to responde to this emergency researchers
worldwide are intensively studying the Covid-19 pandemic dynamics. In this
work, a SIRSi compartmental model is proposed, which is a modification of the
known classical SIR model. The proposed SIRSi model considers differences in
the immunization within a population, and the possibility of unreported or
asymptomatic cases. The model is adjusted to three major cities of S\~ao Paulo
State, in Brazil, namely, S\~ao Paulo, Santos and Campinas, providing estimates
on the duration and peaks of the outbreak.


  B cell receptor sequences diversify through mutations introduced by
purpose-built cellular machinery. A recent paper has concluded that a
"templated mutagenesis" process is a major contributor to somatic
hypermutation, and therefore immunoglobulin diversification, in mice and
humans. In this proposed process, mutations in the immunoglobulin locus are
introduced by copying short segments from other immunoglobulin genes. If true,
this would overturn decades of research on B cell diversification, and would
require a complete re-write of computational methods to analyze B cell data for
these species.
  In this paper, we re-evaluate the templated mutagenesis hypothesis. By
applying the original inferential method using potential donor templates absent
from B cell genomes, we obtain estimates of the methods's false positive rates.
We find false positive rates of templated mutagenesis in murine and human
immunoglobulin loci that are similar to or even higher than the original rate
inferences, and by considering the bases used in substitution we find evidence
that if templated mutagenesis occurs, it is at a low rate. We also show that
the statistically significant results in the original paper can easily result
from a slight misspecification of the null model.


  Coronavirus disease 2019 (COVID-19) is a continuously devastating public
health and the world economy. One of the major challenges in controlling the
COVID-19 outbreak is its asymptomatic infection and transmission, which are
elusive and defenseless in most situations. The pathogenicity and virulence of
asymptomatic COVID-19 remain mysterious. Based on the genotyping of 20656
Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome isolates,
we reveal that asymptomatic infection is linked to SARS-CoV-2 11083G>T
mutation, i.e., leucine (L) to phenylalanine (F) substitution at the residue 37
(L37F) of nonstructure protein 6 (NSP6). By analyzing the distribution of
11083G>T in various countries, we unveil that 11083G>T may correlate with the
hypotoxicity of SARS-CoV-2. Moreover, we show a global decaying tendency of the
11083G>T mutation ratio indicating that 11083G>T hinders SARS-CoV-2
transmission capacity. Sequence alignment found both NSP6 and residue 37
neighborhoods are relatively conservative over a few coronaviral species,
indicating their importance in regulating host cell autophagy to undermine
innate cellular defense against viral infection. Using machine learning and
topological data analysis, we demonstrate that mutation L37F has made NSP6
energetically less stable. The rigidity and flexibility index and several
network models suggest that mutation L37F may have compromised the NSP6
function, leading to a relatively weak SARS-CoV subtype. This assessment is a
good agreement with our genotyping of SARS-CoV-2 evolution and transmission
across various countries and regions over the past few months.


  The adaptive immune system provides a diverse set of molecules that can mount
specific responses against a multitude of pathogens. Memory is a key feature of
adaptive immunity, which allows organisms to respond more readily upon
re-infections. However, differentiation of memory cells is still one of the
least understood cell fate decisions. Here, we introduce a mathematical
framework to characterize optimal strategies to store memory to maximize the
utility of immune response over an organism's lifetime. We show that memory
production should be actively regulated to balance between affinity and
cross-reactivity of immune receptors for an effective protection against
evolving pathogens. Moreover, we predict that specificity of memory should
depend on the organism's lifespan, and shorter-lived organisms with fewer
pathogenic encounters should store more cross-reactive memory. Our framework
provides a baseline to gauge the efficacy of immune memory in light of an
organism's coevolutionary history with pathogens.


  Using the classical Susceptible-Infected-Recovered epidemiological model, an
analytical formula is derived for the number of beds occupied by Covid-19
patients. The analytical curve is fitted to data in Belgium, France, New York
City and Switzerland, with a correlation coefficient exceeding 98.8%,
suggesting that finer models are unnecessary with such macroscopic data. The
fitting is used to extract estimates of the doubling time in the ascending
phase of the epidemic, the mean recovery time and, for those who require
medical intervention, the mean hospitalization time. Large variations can be
observed among different outbreaks.


  An actively controlled Susceptible-Infected-Susceptible (actSIS) contagion
model is presented for studying epidemic dynamics with continuous-time feedback
control of infection rates. Our work is inspired by the observation that
epidemics can be controlled through decentralized disease-control strategies
such as quarantining, sheltering in place, social distancing, etc., where
individuals actively modify their contact rates with others in response to
observations of infection levels in the population. Accounting for a time lag
in observations and categorizing individuals into distinct sub-populations
based on their risk profiles, we show that the actSIS model manifests
qualitatively different features as compared with the SIS model. In a
homogeneous population of risk-averters, the endemic equilibrium is always
reduced, although the transient infection level can exhibit overshoot or
undershoot. In a homogeneous population of risk-tolerating individuals, the
system exhibits bistability, which can also lead to reduced infection. For a
heterogeneous population comprised of risk-tolerators and risk-averters, we
prove conditions on model parameters for the existence of a Hopf bifurcation
and sustained oscillations in the infected population.


  Fingerprint-based models for protein-ligand binding have demonstrated
outstanding success on benchmark datasets; however, these models may not learn
the correct binding rules. To assess this concern, we use in silico datasets
with known binding rules to develop a general framework for evaluating model
attribution. This framework identifies fragments that a model considers
necessary to achieve a particular score, sidestepping the need for a model to
be differentiable. Our results confirm that high-performing models may not
learn the correct binding rule, and suggest concrete steps that can remedy this
situation. We show that adding fragment-matched inactive molecules (decoys) to
the data reduces attribution false negatives, while attribution false positives
largely arise from the background correlation structure of molecular data.
Normalizing for these background correlations helps to reveal the true binding
logic. Our work highlights the danger of trusting attributions from
high-performing models and suggests that a closer examination of fingerprint
correlation structure and better decoy selection may help reduce
misattributions.


  The emergence of an epidemic evokes the need to monitor its spread and assess
and validate any mitigation measures enacted by governments and administrative
bodies in real time. We present here a method to observe and quantify this
spread and the response of affected populations and governing bodies and apply
it to COVID-19 as a case study. This method provides means to simultaneously
track in real time quantities such as the mortality and the recovery rates as
well as the number of new infections caused by an infected person. With
sufficient data, this method enables thorough monitoring and assessment of an
epidemic without assumptions regarding the evolution of the pandemic in the
future.


  Genome-Wide Association Studies are typically conducted using linear models
to find genetic variants associated with common diseases. In these studies,
association testing is done on a variant-by-variant basis, possibly missing out
on non-linear interaction effects between variants. Deep networks can be used
to model these interactions, but they are difficult to train and interpret on
large genetic datasets. We propose a method that uses the gradient based deep
interpretability technique named DeepLIFT to show that known diabetes genetic
risk factors can be identified using deep models along with possibly novel
associations.


  In the wake of the COVID-19 pandemic many countries implemented containment
measures to reduce disease transmission. Studies using digital data sources
show that the mobility of individuals was effectively reduced in multiple
countries. However, it remains unclear whether these reductions caused deeper
structural changes in mobility networks, and how such changes may affect
dynamic processes on the network. Here we use movement data of mobile phone
users to show that mobility in Germany has not only been reduced considerably:
Lockdown measures caused substantial and long-lasting structural changes in the
mobility network. We find that long-distance travel was reduced
disproportionately strongly. The trimming of long-range network connectivity
leads to a more local, clustered network and a moderation of the "small-world"
effect. We demonstrate that these structural changes have a considerable effect
on epidemic spreading processes by "flattening" the epidemic curve and delaying
the spread to geographically distant regions.


  Self-assembly is a ubiquitous process in synthetic and biological systems,
broadly defined as the spontaneous organization of multiple subunits (e.g.
macromolecules, particles) into ordered multi-unit structures. The vast
majority of equilibrium assembly processes give rise to two states: one
consisting of dispersed disassociated subunits, and the other, a bulk-condensed
state of unlimited size. This review focuses on the more specialized class of
{\it self-limiting assembly}, which describes equilibrium assembly processes
resulting in finite-size structures. These systems pose a generic and basic
question, how do thermodynamic processes involving non-covalent interactions
between identical subunits ``measure'' and select the size of assembled
structures? In this review, we begin with an introduction to the basic
statistical mechanical framework for assembly thermodynamics, and use this to
highlight the key physical ingredients that ensure equilibrium assembly will
terminate at finite dimensions. Then, we introduce examples of self-limiting
assembly systems, and classify them within this framework based on two broad
categories: {\it self-closing assemblies} and {\it open-boundary assemblies}.
These include well-known cases in biology and synthetic soft matter --
micellization of amphiphiles and shell/tubule formation of tapered subunits --
as well as less widely known classes of assemblies, such as short-range
attractive/long-range repulsive systems and geometrically-frustrated
assemblies. For each of these self-limiting mechanisms, we describe the
physical mechanisms that select equilibrium assembly size, as well as potential
limitations of finite-size selection. Finally, we discuss alternative
mechanisms for finite-size assemblies, and draw contrasts with the size-control
that these can achieve relative to self-limitation in equilibrium,
single-species assemblies.


  In $1993$, Holt and Lawton introduced a stochastic model of two host species
parasitized by a common parasitoid species. We introduce and analyze a
generalization of these stochastic difference equations with any number of host
species, stochastically varying parasitism rates, stochastically varying host
intrinsic fitnesses, and stochastic immigration of parasitoids. Despite the
lack of direct, host density-dependence, we show that this system is
dissipative i.e. enters a compact set in finite time for all initial
conditions. When there is a single host species, stochastic persistence and
extinction of the host is characterized using external Lyapunov exponents
corresponding to the average per-capita growth rates of the host when rare.
When a single host persists, say species $i$, a explicit expression is derived
for the average density, $P_i^*$, of the parasitoid at the stationary
distributions supporting both species. When there are multiple host species, we
prove that the host species with the largest $P_i^*$ value stochastically
persists, while the other host species are asymptotically driven to extinction.
A review of the main mathematical methods used to prove the results and future
challenges are given.


  Objective: To make an estimate of the excess deaths caused by COVID-19 in the
non-violent mortality of Peru, controlling for the effect of quarantine.
Methods: Analysis of longitudinal data from the departments of Peru using
official public information from the National Death Information System and the
Ministry of Health of Peru. The analysis is performed between January 1, 2018
and June 23, 2020 (100 days of quarantine). The daily death rate per million
inhabitants has been used. The days in which the departments were quarantined
with a limit number of accumulated cases of COVID-19 were used to estimate the
quarantine impact. Three limits were established for cases: less than 1, 10 and
100 cases. Result: In Peru, the daily death rate per million inhabitants
decreased by -1.89 (95% CI: -2.70; -1.07) on quarantine days and without
COVID-19 cases. When comparing this result with the total number of non-violent
deaths, the excess deaths during the first 100 days of quarantine is 36,230.
This estimate is 1.12 times the estimate with data from 2019 and 4.2 times the
deaths officers by COVID-19. Conclusion: Quarantine reduced nonviolent deaths;
however, they are overshadowed by the increase as a direct or indirect cause of
the pandemic. Therefore, the difference between the number of current deaths
and that of past years underestimates the real excess of deaths.


  In this paper, a susceptible-infected-removed (SIR) model has been used to
track the evolution of the spread of the COVID-19 virus in four countries of
interest. In particular, the epidemic model, that depends on some basic
characteristics, has been applied to model the time evolution of the disease in
Italy, India, South Korea and Iran. The economic, social and health
consequences of the spread of the virus have been cataclysmic. Hence, it is
essential that available mathematical models can be developed and used for the
comparison to be made between published data sets and model predictions. The
predictions estimated from the SIR model here, can be used in both the
qualitative and quantitative analysis of the spread. It gives an insight into
the spread of the virus that the published data alone cannot do by updating
them and the model on a daily basis. For example, it is possible to detect the
early onset of a spike in infections or the development of a second wave using
our modeling approach. We considered data from March to June, 2020, when
different communities are severely affected. We demonstrate predictions
depending on the model's parameters related to the spread of COVID-19 until
September 2020. By comparing the published data and model results, we conclude
that in this way, it may be possible to better reflect the success or failure
of the adequate measures implemented by governments and individuals to mitigate
and control the current pandemic.


  Analytical formulas for effective drift, diffusivity, run times, and run
lengths are derived for an intracellular transport system consisting of a cargo
attached to two cooperative but not identical molecular motors (for example,
kinesin-1 and kinesin-2) which can each attach and detach from a microtubule.
The dynamics of the motor and cargo in each phase are governed by stochastic
differential equations, and the switching rates depend on the spatial
configuration of the motor and cargo. This system is analyzed in a limit where
the detached motors have faster dynamics than the cargo, which in turn has
faster dynamics than the attached motors. The attachment and detachment rates
are also taken to be slow relative to the spatial dynamics. Through an
application of iterated stochastic averaging to this system, and the use of
renewal-reward theory to stitch together the progress within each switching
phase, we obtain explicit analytical expressions for the effective drift,
diffusivity, and processivity of the motor-cargo system. Our approach accounts
in particular for jumps in motor-cargo position that occur during attachment
and detachment events, as the cargo tracking variable makes a rapid adjustment
due to the averaged fast scales. The asymptotic formulas are in generally good
agreement with direct stochastic simulations of the detailed model based on
experimental parameters for various pairings of kinesin-1 and kinesin-2 under
assisting, hindering, or no load.


  A Profile Mixture Model is a model of protein evolution, describing sequence
data in which sites are assumed to follow many related substitution processes
on a single evolutionary tree. The processes depend in part on different amino
acid distributions, or profiles, varying over sites in aligned sequences. A
fundamental question for any stochastic model, which must be answered
positively to justify model-based inference, is whether the parameters are
identifiable from the probability distribution they determine. Here we show
that a Profile Mixture Model has identifiable parameters under circumstances in
which it is likely to be used for empirical analyses. In particular, for a tree
relating 9 or more taxa, both the tree topology and all numerical parameters
are generically identifiable when the number of profiles is less than 74.


  Transport of intracellular cargo is often mediated by teams of molecular
motors that function in a chaotic environment under varying conditions. We show
that the motors have unique steady state behavior which enables transport
modalities that are robust. Under reduced ATP concentrations, multi-motor
configurations are preferred over single motors. Higher load force drives
motors to cluster, but very high loads compel them to separate in a manner that
promotes immediate cargo movement once the load subsides. These inferences,
backed by analytical guarantees, provide unique insights into the coordination
strategies adopted by molecular motors to transport intracellular cargo.


  On the base of logic discrete equations system mathematical modeling of
COVID-19 epidemic spread was carried out in the world and in the countries with
the largest number of infected people such as the USA, Brasil, Russia and India
in the first half of 2020. It was shown that for the countries with strong
restrictive measures the spread of COVID-19 fit on a single wave with small
capacity as for a number of countries with violation of restrictive measures
the spread of the epidemic fit on a waves superposition. For countries with
large population mixing, the spread of the epidemic today also fits into a
single wave, but with a huge capacity value (for Brazil - 80 million people,
for India - 40 million people). We estimated that the epidemic spread in the
world today fits into 5 waves. The first two waves are caused by the epidemic
spread in China (the first - in Wuhan), the third - by the epidemic spread in
European countries, the fourth mainly by the epidemic spread in Russia and the
USA, the fifth wave is mainly caused by the epidemic spread in Latin America
and South Asia. It was the fifth wave that led to the spread of the coronavirus
epidemic COVID-19 entering a new phase, with an increase in the number of
infected more than 100 thousand inhabitants. For all the studied countries and
the world, for each of the superposition waves, the wave capacities and growth
indicators were calculated. The local peaks of the waves and their ending times
are determined. It was the fifth wave that led to the fact that COVID-19 spread
is entering a new phase, with the increase in the number of infected people
being more than 100 thousand inhabitants. For all the countries being examined
and for the whole world, for each of the superposition waves we calculated the
waves capacities and index of infected people growth.


  Coronavirus Disease 2019 (COVID-19) has been creating a worldwide pandemic
situation. Repurposing drugs, already shown to be free of harmful side effects,
for the treatment of COVID-19 patients is an important option in launching
novel therapeutic strategies. Therefore, reliable molecule interaction data are
a crucial basis, where drug-/protein-protein interaction networks establish
invaluable, year-long carefully curated data resources. However, these
resources have not yet been systematically exploited using high-performance
artificial intelligence approaches. Here, we combine three networks, two of
which are year-long curated, and one of which, on SARS-CoV-2-human host-virus
protein interactions, was published only most recently (30th of April 2020),
raising a novel network that puts drugs, human and virus proteins into mutual
context. We apply Variational Graph AutoEncoders (VGAEs), representing most
advanced deep learning based methodology for the analysis of data that are
subject to network constraints. Reliable simulations confirm that we operate at
utmost accuracy in terms of predicting missing links. We then predict hitherto
unknown links between drugs and human proteins against which virus proteins
preferably bind. The corresponding therapeutic agents present splendid starting
points for exploring novel host-directed therapy (HDT) options.


  COVID-19 is a devastating respiratory and inflammatory illness caused by a
new coronavirus that is rapidly spreading throughout the human population. Over
the past 6 months, severe acute respiratory syndrome coronavirus 2
(SARS-CoV-2), the virus responsible for COVID-19, has already infected over
11.6 million (25% located in United States) and killed more than 540K people
around the world. As we face one of the most challenging times in our recent
history, there is an urgent need to identify drug candidates that can attack
SARS-CoV-2 on multiple fronts. We have therefore initiated a computational
dynamics drug pipeline using molecular modeling, structure simulation, docking
and machine learning models to predict the inhibitory activity of several
million compounds against two essential SARS-CoV-2 viral proteins and their
host protein interactors; S/Ace2, Tmprss2, Cathepsins L and K, and Mpro to
prevent binding, membrane fusion and replication of the virus, respectively.
All together we generated an ensemble of structural conformations that increase
high quality docking outcomes to screen over >6 million compounds including all
FDA-approved drugs, drugs under clinical trial (>3000) and an additional >30
million selected chemotypes from fragment libraries. Our results yielded an
initial set of 350 high value compounds from both new and FDA-approved
compounds that can now be tested experimentally in appropriate biological model
systems. We anticipate that our results will initiate screening campaigns and
accelerate the discovery of COVID-19 treatments.


  The current rate of species extinction is rapidly approaching unprecedented
highs and life on Earth presently faces a sixth mass extinction event driven by
anthropogenic activity, climate change and ecological collapse. The field of
conservation genetics aims at preserving species by using their levels of
genetic diversity, usually measured as neutral genome-wide diversity, as a
barometer for evaluating population health and extinction risk. A fundamental
assumption is that higher levels of genetic diversity lead to an increase in
fitness and long-term survival of a species. Here, we argue against the
perceived importance of neutral genetic diversity for the conservation of wild
populations and species. We demonstrate that no simple general relationship
exists between neutral genetic diversity and the risk of species extinction.
Instead, a better understanding of the properties of functional genetic
diversity, demographic history, and ecological relationships, is necessary for
developing and implementing effective conservation genetic strategies.


  Predicting an accurate expected number of future COVID-19 cases is essential
to properly evaluate the effectiveness of any treatment or preventive measure.
This study aimed to identify the most appropriate mathematical model to
prospectively predict the expected number of cases without any intervention.
The total number of cases for the COVID-19 epidemic in 28 countries was
analyzed and fitted to several simple rate models including the logistic,
Gompertz, quadratic, simple square, and simple exponential growth models. The
resulting model parameters were used to extrapolate predictions for more recent
data. While the Gompertz growth models (mean R2 = 0.998) best fitted the
current data, uncertainties in the eventual case limit made future predictions
with logistic models prone to errors. Of the other models, the quadratic rate
model (mean R2 = 0.992) fitted the current data best for 25 (89 %) countries as
determined by R2 values. The simple square and quadratic models accurately
predicted the number of future total cases 37 and 36 days in advance
respectively, compared to only 15 days for the simple exponential model. The
simple exponential model significantly overpredicted the total number of future
cases while the quadratic and simple square models did not. These results
demonstrated that accurate future predictions of the case load in a given
country can be made significantly in advance without the need for complicated
models of population behavior and generate a reliable assessment of the
efficacy of current prescriptive measures against disease spread.


  The ongoing COVID-19 pandemic is challenging every part of society. From a
scientific point of view the first major task is to predict the dynamics of the
pandemic, allowing governments to allocate proper resources and measures to
fight it, as well as gauging the success of these measures by comparison with
the predictions in hindsight. The vast majority of pandemic models are based on
extensive models with large numbers of fit parameters, leading to individual
descriptions for every hot spot on the world. This makes predictions and
comparisons cumbersome, if not impossible. We here propose a different
approach, by moving away from a description over time, and instead choosing the
total number of infected people in an enclosed area as the independent
variable. Analyzing a few hot spots data, we derive an empirical formula for
the dynamics, dependent only on three variables. The final number of infections
is strictly connected to one fit parameter we call mitigation factor, which in
turn is mostly dependent only on the enclosed population. Despite its
simpleness, this description applies to every of the around 50 countries we
have analyzed, allows to separate different waves of the pandemic, provides a
figure of merit for the overall usefulness of government measures, and shows
when a pandemic is ending. Our model is robust against undetected cases, and
allows all nations, in particular those with fewer resources, to reasonably
predict the outcome of the pandemic in their country.


  This paper extends the canonical model of epidemiology, SIRD model, to allow
for time varying parameters for real-time measurement of the stance of the
COVID-19 pandemic. Time variation in model parameters is captured using the
generalized autoregressive score modelling structure designed for the typically
daily count data related to pandemic. The resulting specification permits a
flexible yet parsimonious model structure with a very low computational cost.
This is especially crucial at the onset of the pandemic when the data is scarce
and the uncertainty is abundant. Full sample results show that countries
including US, Brazil and Russia are still not able to contain the pandemic with
the US having the worst performance. Furthermore, Iran and South Korea are
likely to experience the second wave of the pandemic. A real-time exercise show
that the proposed structure delivers timely and precise information on the
current stance of the pandemic ahead of the competitors that use rolling
window. This, in turn, transforms into accurate short-term predictions of the
active cases. We further modify the model to allow for unreported cases.
Results suggest that the effects of the presence of these cases on the
estimation results diminish towards the end of sample with the increasing
number of testing.


  Diseases and other contagion phenomena in nature and society can interact
asymmetrically, such that one can benefit from the other, which in turn impairs
the first, in analogy with predator-prey systems. Here, we consider two models
for interacting disease-like dynamics with asymmetric interactions and
different associated time scales. Using rate equations for homogeneously mixed
populations, we show that the stationary prevalences and phase diagrams of each
model behave differently with respect to variations of the relative time
scales. We also characterize in detail the regime where transient oscillations
are observed, a pattern that is inherent to asymmetrical interactions but often
ignored in the literature. Our results contribute to a better understanding of
disease dynamics in particular, and interacting processes in general, and could
provide interesting insights for real-world applications, most notably, the
interplay between the dynamics of fact-checked and fake news.


  How to obtain informative representations of molecules is a crucial
prerequisite in AI-driven drug design and discovery. Recent researches abstract
molecules as graphs and employ Graph Neural Networks (GNNs) for molecular
representation learning. Nevertheless, two issues impede the usage of GNNs in
real scenarios: (1) insufficient labeled molecules for supervised training; (2)
poor generalization capability to new-synthesized molecules. To address them
both, we propose a novel framework, GROVER, which stands for Graph
Representation frOm self-superVised mEssage passing tRansformer. With carefully
designed self-supervised tasks in node-, edge- and graph-level, GROVER can
learn rich structural and semantic information of molecules from enormous
unlabelled molecular data. Rather, to encode such complex information, GROVER
integrates Message Passing Networks into the Transformer-style architecture to
deliver a class of more expressive encoders of molecules. The flexibility of
GROVER allows it to be trained efficiently on large-scale molecular dataset
without requiring any supervision, thus being immunized to the two issues
mentioned above. We pre-train GROVER with 100 million parameters on 10 million
unlabelled molecules -- the biggest GNN and the largest training dataset in
molecular representation learning. We then leverage the pre-trained GROVER for
molecular property prediction followed by task-specific fine-tuning, where we
observe a huge improvement (more than 6% on average) from current
state-of-the-art methods on 11 challenging benchmarks. The insights we gained
are that well-designed self-supervision losses and largely-expressive
pre-trained models enjoy the significant potential on performance boosting.


  The severe acute respiratory syndrome COVID-19 has been in the center of the
ongoing global health crisis in 2020. The high prevalence of mild cases
facilitates sub-notification outside hospital environments and the number of
those who are or have been infected remains largely unknown, leading to poor
estimates of the crude mortality rate of the disease. Here we use a simple
model to describe the number of accumulated deaths caused by COVID-19. The
close connection between the proposed model and an approximate solution of the
SIR model provides a system of equations whose solutions are robust estimates
of epidemiological parameters. We find that the crude mortality varies between
$10^{-4}$ and $10^{-3}$ depending on the severity of the outbreak which is
lower than previous estimates obtained from laboratory confirmed patients. We
also estimate quantities of practical interest such as the basic reproduction
number and the expected number of deaths in the asymptotic limit with and
without social distancing measures and lockdowns, which allow us to measure the
efficiency of these interventions.


  The SARS-CoV-2 pandemic has necessitated mitigation efforts around the world.
We use only reported deaths in the two weeks after the first death to determine
infection parameters, in order to make predictions of hidden variables such as
the time dependence of the number of infections. Early deaths are sporadic and
discrete so the use of network models of epidemic spread is imperative, with
the network itself a crucial random variable. Location-specific population age
distributions and population densities must be taken into account when
attempting to fit these events with parametrized models. These characteristics
render naive Bayesian model comparison impractical as the networks have to be
large enough to avoid finite-size effects. We reformulated this problem as the
statistical physics of independent location-specific `balls' attached to every
model in a six-dimensional lattice of 56448 parametrized models by elastic
springs, with model-specific `spring constants' determined by the stochasticity
of network epidemic simulations for that model. The distribution of balls then
determines all Bayes posterior expectations. Important characteristics of the
contagion are determinable: the fraction of infected patients that die
($0.017\pm 0.009$), the expected period an infected person is contagious ($22
\pm 6$ days) and the expected time between the first infection and the first
death ($25 \pm 8$ days) in the US. The rate of exponential increase in the
number of infected individuals is $0.18\pm 0.03$ per day, corresponding to 65
million infected individuals in one hundred days from a single initial
infection, which fell to 166000 with even imperfect social distancing
effectuated two weeks after the first recorded death. The fraction of compliant
socially-distancing individuals matters less than their fraction of social
contact reduction for altering the cumulative number of infections.


  Antimicrobial resistance is an emerging global health crisis that is
undermining advances in modern medicine and, if unmitigated, threatens to kill
10 million people per year worldwide by 2050. Research over the last decade has
demonstrated that the differences between genetically identical cells in the
same environment can lead to drug resistance. Fluctuations in gene expression,
modulated by gene regulatory networks, can lead to non-genetic heterogeneity
that results in the fractional killing of microbial populations causing drug
therapies to fail; this non-genetic drug resistance can enhance the probability
of acquiring genetic drug resistance mutations. Mathematical models of gene
networks can elucidate general principles underlying drug resistance, predict
the evolution of resistance, and guide drug resistance experiments in the
laboratory. Cells genetically engineered to carry synthetic gene networks
regulating drug resistance genes allow for controlled, quantitative experiments
on the role of non-genetic heterogeneity in the development of drug resistance.
In this perspective article, we emphasize the contributions that mathematical,
computational, and synthetic gene network models play in advancing our
understanding of antimicrobial resistance to discover effective therapies
against drug-resistant infections.


  Zinc, a suspected potentiator of learning and memory, is shown to affect
exocytotic release and storage in neurotransmitter-containing vesicles.
Structural and size analysis of the vesicular dense core and halo using
transmission electron microscopy was combined with single-cell amperometry to
study the vesicle size changes induced after zinc treatment and to compare
these changes to theoretical predictions based on the concept of partial
release as opposed to full quantal release. This powerful combined analytical
approach establishes the existence of an unsuspected strong link between
vesicle structure and exocytotic dynamics which can be used to explain the
mechanism of regulation of synaptic plasticity by Zn 2+ through modulation of
neurotransmitter release.


  Reactive oxygen and nitrogen species (ROS and RNS) play important roles in
various physiological processes (e.g., phagocytosis) and pathological
conditions (e.g., cancer). The primary ROS/RNS, viz., hydrogen peroxide,
peroxynitrite ion, nitric oxide, and nitrite ion, can be oxidized at different
electrode potentials and therefore detected and quantified by electroanalytical
techniques. Nanometer-sized electrochemical probes are especially suitable for
measuring ROS/RNS in single cells and cellular organelles. In this article, we
survey recent advances in localized measurements of ROS/RNS inside single cells
and discuss several methodological issues, including optimization of
nanoelectrode geometry, precise positioning of an electrochemical probe inside
a cell, and interpretation of electroanalytical data.


  The benefits, both in terms of productivity and public health, are
investigated for different levels of engagement with the test, trace and
isolate procedures in the context of a pandemic in which there is little or no
herd immunity. Simple mathematical modelling is used in the context of a
single, relatively closed workplace such as a factory or back-office where, in
normal operation, each worker has lengthy interactions with a fixed set of
colleagues. A discrete-time SEIR model on a fixed interaction graph is
simulated with parameters that are motivated by the recent COVID-19 pandemic in
the UK during a post-peak phase, including a small risk of viral infection from
outside the working environment. Two kinds of worker are assumed, transparents
who regularly test, share their results with colleagues and isolate as soon as
a contact tests positive for the disease, and opaques who do none of these.
Moreover, the simulations are constructed as a ``playable model" in which the
transparency level, disease parameters and mean interaction degree can be
varied by the user. The model is analysed in the continuum limit. All
simulations point to the double benefit of transparency in maximising
productivity and minimising overall infection rates. Based on these findings,
public policy implications are discussed on how to incentivise this mutually
beneficial behaviour in different kinds of workplace, and simple
recommendations are made.


  Networks provide a mathematically rich framework to represent social contacts
sufficient for the transmission of disease. Social networks are often highly
clustered and fail to be locally tree-like. In this paper, we study the effects
of clustering on the spread of sequential strains of a pathogen using the
generating function formulation under a complete cross-immunity coupling,
deriving conditions for the threshold of coexistence of the second strain. We
show that clustering reduces the coexistence threshold of the second strain and
its outbreak size in Poisson networks, whilst exhibiting the opposite effects
on uniform-degree models. We conclude that clustering within a population must
increase the ability of the second wave of an epidemic to spread over a
network. We apply our model to the study of multilayer clustered networks and
observe the fracturing of the residual graph at two distinct
transmissibilities.


  Over 80% of human cancers originate from the epithelium, which covers the
outer and inner surfaces of organs and blood vessels. In stratified epithelium,
the bottom layers are occupied by stem and stem-like cells that continually
divide and replenish the upper layers. In this work, we study the spread of
premalignant mutant clones and cancer initiation in stratified epithelium using
the biased voter model on stacked two-dimensional lattices. Our main result is
an estimate of the propagation speed of a premalignant mutant clone, which is
asymptotically precise in the cancer-relevant weak-selection limit. We use our
main result to study cancer initiation under a two-step mutational model of
cancer, which includes computing the distributions of the time of cancer
initiation and the size of the premalignant clone giving rise to cancer. Our
work quantifies the effect of epithelial tissue thickness on the process of
carcinogenesis, thereby contributing to an emerging understanding of the
spatial evolutionary dynamics of cancer.


  We define a new problem in comparative genomics, denoted PQ-Tree Search, that
takes as input a PQ-tree $T$ representing the known gene orders of a gene
cluster of interest, a gene-to-gene substitution scoring function $h$, integer
parameters $d_T$ and $d_S$, and a new genome $S$. The objective is to identify
in $S$ approximate new instances of the gene cluster that could vary from the
known gene orders by genome rearrangements that are constrained by $T$, by gene
substitutions that are governed by $h$, and by gene deletions and insertions
that are bounded from above by $d_T$ and $d_S$, respectively. We prove that the
PQ-Tree Search problem is NP-hard and propose a parameterized algorithm that
solves the optimization variant of PQ-Tree Search in $O^*(2^{\gamma})$ time,
where $\gamma$ is the maximum degree of a node in $T$ and $O^*$ is used to hide
factors polynomial in the input size. The algorithm is implemented as a search
tool, denoted PQFinder, and applied to search for instances of chromosomal gene
clusters in plasmids, within a dataset of 1,487 prokaryotic genomes. We report
on 29 chromosomal gene clusters that are rearranged in plasmids, where the
rearrangements are guided by the corresponding PQ-tree. One of these results,
coding for a heavy metal efflux pump, is further analysed to exemplify how
PQFinder can be harnessed to reveal interesting new structural variants of
known gene clusters. The code for the tool as well as all the data needed to
reconstruct the results are publicly available on GitHub
(github.com/GaliaZim/PQFinder).


  Protein-ligand docking is an in silico tool used to screen potential drug
compounds for their ability to bind to a given protein receptor within a
drug-discovery campaign. Experimental drug screening is expensive and time
consuming, and it is desirable to carry out large scale docking calculations in
a high-throughput manner to narrow the experimental search space. Few of the
existing computational docking tools were designed with high performance
computing in mind. Therefore, optimizations to maximize use of high-performance
computational resources available at leadership-class computing facilities
enables these facilities to be leveraged for drug discovery. Here we present
the porting, optimization, and validation of the AutoDock-GPU program for the
Summit supercomputer, and its application to initial compound screening efforts
to target proteins of the SARS-CoV-2 virus responsible for the current COVID-19
pandemic.


  The Richards family models comprise a well-known set of models with useful
parameters to describe several aspects of disease outbreaks. Some of these
models have been used to study the current Covid-19 pandemic. However, there
seems to be confusion regarding the discrimination among competing models. In
this paper a detailed application of Akaikes information approach is used to
discern among models using data from The European Union, The United States and
The United Kingdom. We argue that the epidemiological characteristics derived
from competing models should be examined to complement the selection strategy,
and the implicit properties of the models contrasted with the available data.
Detailed analytical expressions of the epidemiological characteristics implied
by the selected parametrizations are also offered.


  The well-known Linear Chain Trick (LCT) allows modelers to derive mean field
ODEs that assume gamma (Erlang) distributed passage times, by transitioning
individuals sequentially through a chain of sub-states. The time spent in these
states is the sum of $k$ exponentially distributed random variables, and is
thus gamma (Erlang) distributed. The Generalized Linear Chain Trick (GLCT)
extends this technique to the much broader phase-type family of distributions,
which includes exponential, Erlang, hypoexponential, and Coxian distributions.
Intuitively, phase-type distributions are the absorption time distributions for
continuous time Markov chains (CTMCs). Here we review CTMCs and phase-type
distributions, then illustrate how to use the GLCT to efficiently build mean
field ODE models from underlying stochastic model assumptions. We generalize
the Rosenzweig-MacArthur and SEIR models and show the benefits of using the
GLCT to compute numerical solutions. These results highlight some practical
benefits, and the intuitive nature, of using the GLCT to derive ODE models from
first principles.


  We propose a time-fractional compartmental model (SEI$_A$I$_S$HRD) comprising
of the susceptible, exposed, infected (asymptomatic and symptomatic),
hospitalized, recovered and dead population for the Covid-19 pandemic. We study
the properties and dynamics of the proposed model. The conditions under which
the disease-free and endemic equilibrium points are asymptotically stable are
discussed. Furthermore, we study the sensitivity of the parameters and use the
data from Tennessee state (as a case study) to discuss identifiability of the
parameters of the model. The non-negative parameters in the model are obtained
by solving inverse problems with empirical data from California, Florida,
Georgia, Maryland, Tennessee, Texas, Washington and Wisconsin. The basic
reproduction number is seen to be slightly above the critical value of one
suggesting that stricter measures such as the use of face-masks, social
distancing, contact tracing, and even longer stay-at-home orders need to be
enforced in order to mitigate the spread of the virus. As stay-at-home orders
are rescinded in some of these states, we see that the number of cases began to
increase almost immediately and may continue to rise until the end of the year
2020 unless stricter measures are taken.


  We consider the development of hyperbolic transport models for the
propagation in space of an epidemic phenomenon described by a classical
compartmental dynamics. The model is based on a kinetic description at discrete
velocities of the spatial movement and interactions of a population of
susceptible, infected and recovered individuals. Thanks to this, the unphysical
feature of instantaneous diffusive effects, which is typical of parabolic
models, is removed. In particular, we formally show how such reaction-diffusion
models are recovered in an appropriate diffusive limit. The kinetic transport
model is therefore considered within a spatial network, characterizing
different places such as villages, cities, countries, etc. The transmission
conditions in the nodes are analyzed and defined. Finally, the model is solved
numerically on the network through a finite-volume IMEX method able to maintain
the consistency with the diffusive limit without restrictions due to the
scaling parameters. Several numerical tests for simple epidemic network
structures are reported and confirm the ability of the model to correctly
describe the spread of an epidemic.


  Socio-economic disparities quite often have a central role in the unfolding
of large-scale catastrophic events. One of the most concerning aspects of the
ongoing COVID-19 pandemics is that it disproportionately affects people from
Black and African American backgrounds creating an unexpected infection gap.
Interestingly, the abnormal impact on these ethnic groups seem to be almost
uncorrelated with other risk factors, including co-morbidity, poverty, level of
education, access to healthcare, residential segregation, and response to
cures. A proposed explanation for the observed incidence gap is that people
from African American backgrounds are more often employed in low-income service
jobs, and are thus more exposed to infection through face-to-face contacts, but
the lack of direct data has not allowed to draw strong conclusions in this
sense so far. Here we introduce the concept of dynamic segregation, that is the
extent to which a given group of people is internally clustered or exposed to
other groups, as a result of mobility and commuting habits. By analysing census
and mobility data on more than 120 major US cities, we found that the dynamic
segregation of African American communities is significantly associated with
the weekly excess COVID-19 incidence and mortality in those communities. The
results confirm that knowing where people commute to, rather than where they
live, is much more relevant for disease modelling.


  With the unfolding of the COVID-19 pandemic, mathematical modeling of
epidemics has been perceived and used as a central element in understanding,
predicting, and governing the pandemic event. However, soon it became clear
that long term predictions were extremely challenging to address. Moreover, it
is still unclear which metric shall be used for a global description of the
evolution of the outbreaks. Yet a robust modeling of pandemic dynamics and a
consistent choice of the transmission metric is crucial for an in-depth
understanding of the macroscopic phenomenology and better-informed mitigation
strategies. In this study, we propose a Markovian stochastic framework designed
to describe the evolution of entropy during the COVID-19 pandemic and the
instantaneous reproductive ratio. We then introduce and use entropy-based
metrics of global transmission to measure the impact and temporal evolution of
a pandemic event. In the formulation of the model, the temporal evolution of
the outbreak is modeled by the master equation of a nonlinear Markov process
for a statistically averaged individual, leading to a clear physical
interpretation. We also provide a full Bayesian inversion scheme for
calibration. The time evolution of the entropy rate, the absolute change in the
system entropy, and the instantaneous reproductive ratio are natural and
transparent outputs of this framework. The framework has the appealing property
of being applicable to any compartmental epidemic model. As an illustration, we
apply the proposed approach to a simple modification of the
Susceptible-Exposed-Infected-Removed (SEIR) model. Applying the model to the
Hubei region, South Korean, Italian, Spanish, German, and French COVID-19
data-sets, we discover a significant difference in the absolute change of
entropy but highly regular trends for both the entropy evolution and the
instantaneous reproductive ratio.


  We study a multi-type SIR epidemic process among a heterogeneous population
that interacts through a network. When we base social contact on a random graph
with given vertex degrees, we give limit theorems on the fraction of infected
individuals. For a given social distancing individual strategies, we establish
the epidemic reproduction number $R_0$ which can be used to identify network
vulnerability and inform vaccination policies. In the second part of the paper
we study the equilibrium of the social distancing game, in which individuals
choose their social distancing level according to an anticipated global
infection rate, which then must equal the actual infection rate following their
choices. We give conditions for the existence and uniqueness of equilibrium.
For the case of random regular graphs, we show that voluntary social distancing
will always be socially sub-optimal.


  Within the context of SEIR models, we consider a lockdown that is both
imposed and lifted at an early stage of an epidemic. We show that, in these
models, although such a lockdown may delay deaths, it eventually does not avert
a significant number of fatalities. Therefore, in these models, the efficacy of
a lockdown cannot be gauged by simply comparing figures for the deaths at the
end of the lockdown with the projected figure for deaths by the same date
without the lockdown. We provide a simple but robust heuristic argument to
explain why this conclusion should generalize to more elaborate compartmental
models. We qualitatively discuss some important effects of a lockdown, which go
beyond the scope of simple models, but could cause it to increase or decrease
an epidemic's final toll. Given the significance of these effects in India, and
the limitations of currently available data, we conclude that simple
epidemiological models cannot be used to reliably quantify the impact of the
Indian lockdown on fatalities caused by the COVID-19 pandemic.


  Visual analytics techniques are useful tools to support decision-making and
cope with increasing data, which is particularly important when monitoring
natural or artificial phenomena. When monitoring disease progression, visual
analytics approaches help decision-makers choose to understand or even prevent
dissemination paths. In this paper, we propose a new visual analytics tool for
monitoring COVID-19 dissemination. We use k-nearest neighbors of cities to
mimic neighboring cities and analyze COVID-19 dissemination based on the
comparison of a city under consideration and its neighborhood. Moreover, such
analysis is performed based on periods, which facilitates the assessment of
isolation policies. We validate our tool by analyzing the progression of
COVID-19 in neighboring cities of S\~ao Paulo state, Brazil.


  We analyzed 5,484 close contacts of COVID-19 cases from Italy, all of them
tested for SARS-CoV-2 infection. We found an infection fatality ratio of 2.2%
(95%CI 1.69-2.81%) and identified male sex, age >70 years, cardiovascular
comorbidities, and infection early in the epidemics as risk factors for death.


  Redlining is the discriminatory practice whereby institutions avoided
investment in certain neighborhoods due to their demographics. Here we explore
the lasting impacts of redlining on the spread of COVID-19 in New York City
(NYC). Using data available through the Home Mortgage Disclosure Act, we
construct a redlining index for each NYC census tract via a multi-level
logistical model. We compare this redlining index with the COVID-19 statistics
for each NYC Zip Code Tabulation Area. Accurate mappings of the pandemic would
aid the identification of the most vulnerable areas and permit the most
effective allocation of medical resources, while reducing ethnic health
disparities.


  On June 24th, Governor Cooper announced that North Carolina will not be
moving into Phase 3 of its reopening process at least until July 17th. Given
the recent increases in daily positive cases and hospitalizations, this
decision was not surprising. However, given the political and economic
pressures which are forcing the state to reopen, it is not clear what actions
will help North Carolina to avoid the worst. We use a compartmentalized model
to study the effects of social distancing measures and testing capacity
combined with contact tracing on the evolution of the pandemic in North
Carolina until the end of the year. We find that going back to restrictions
that were in place during Phase 1 will slow down the spread but if the state
wants to continue to reopen or at least remain in Phase 2 or Phase 3 it needs
to significantly expand its testing and contact tracing capacity. Even under
our best-case scenario of high contact tracing effectiveness, the number of
contact tracers the state currently employs is inadequate.


  Graph Neural Networks (GNN) have been extensively used to extract meaningful
representations from graph structured data and to perform predictive tasks such
as node classification and link prediction. In recent years, there has been a
lot of work incorporating edge features along with node features for prediction
tasks. One of the main difficulties in using edge features is that they are
often handcrafted, hard to get, specific to a particular domain, and may
contain redundant information. In this work, we present a framework for
creating new edge features, applicable to any domain, via a combination of
self-supervised and unsupervised learning. In addition to this, we use
Forman-Ricci curvature as an additional edge feature to encapsulate the local
geometry of the graph. We then encode our edge features via a Set Transformer
and combine them with node features extracted from popular GNN architectures
for node classification in an end-to-end training scheme. We validate our work
on three biological datasets comprising of single-cell RNA sequencing data of
neurological disease, \textit{in vitro} SARS-CoV-2 infection, and human
COVID-19 patients. We demonstrate that our method achieves better performance
on node classification tasks over baseline Graph Attention Network (GAT) and
Graph Convolutional Network (GCN) models. Furthermore, given the attention
mechanism on edge and node features, we are able to interpret the cell types
and genes that determine the course and severity of COVID-19, contributing to a
growing list of potential disease biomarkers and therapeutic targets.


  The COVID-19 pandemic poses challenges for continuing economic activity while
reducing health risks. While these challenges can be mitigated through testing,
testing budget is often limited. Here we study how institutions, such as
nursing homes, should utilize a fixed test budget for early detection of an
outbreak. Using an extended network-SEIR model, we show that given a certain
budget of tests, it is generally better to test smaller subgroups of the
population frequently than to test larger groups but less frequently. The
numerical results are consistent with an analytical expression we derive for
the size of the outbreak at detection in an exponential spread model. Our work
provides a simple guideline for institutions: distribute your total tests over
several batches instead of using them all at once. We expect that in the
appropriate scenarios, this easy-to-implement policy recommendation will lead
to earlier detection and better mitigation of local COVID-19 outbreaks.


  Control theory has seen recently impactful applications in network science,
especially in connections with applications in network medicine. A key topic of
research is that of finding minimal external interventions that offer control
over the dynamics of a given network, a problem known as network
controllability. We propose in this article a new solution for this problem
based on genetic algorithms. We tailor our solution for applications in
computational drug repurposing, seeking to maximise its use of FDA-approved
drug targets in a given disease-specific protein-protein interaction network.
We show how our algorithm identifies a number of potentially efficient drugs
for breast, ovarian, and pancreatic cancer. We demonstrate our algorithm on
several benchmark networks from cancer medicine, social networks, electronic
circuits, and several random networks with their edges distributed according to
the Erd\H{o}s-R\'{e}nyi, the small-world, and the scale-free properties.
Overall, we show that our new algorithm is more efficient in identifying
relevant drug targets in a disease network, advancing the computational
solutions needed for new therapeutic and drug repurposing approaches.


  We examine the question of whether it is possible for a diagonalizing system,
to learn to use environmental reward and punishment as an information, in order
to appropriately adapt. More specifically, we study the possiblity of such a
system, to learn to use diagonalization on the basis of a rewarding function.
Relevant phenomena regarding memory are also investigated.


  Across the world, scholars are racing to predict the spread of the novel
coronavirus, COVID-19. Such predictions are often pursued by numerically
simulating epidemics with a large number of plausible combinations of relevant
parameters. It is essential that any forecast of the epidemic trajectory
derived from the resulting ensemble of simulated curves is presented with
confidence intervals that communicate the uncertainty associated with the
forecast. Here we argue that the state-of-the-art approach for summarizing
ensemble statistics does not capture crucial epidemiological information. In
particular, the current approach systematically suppresses information about
the projected trajectory peaks. The fundamental problem is that each time step
is treated separately in the statistical analysis. We suggest using curve-based
descriptive statistics to summarize trajectory ensembles. The results presented
allow researchers to report more representative confidence intervals, resulting
in more realistic projections of epidemic trajectories and -- in turn -- enable
better decision making in the face of the current and future pandemics.


  By coupling peptides with DNA tags (i.e., 'barcodes'), it is now possible to
harness high-throughput sequencing (HTS) technologies to enable highly
multiplexed peptide-based assays, which have a variety of potential
applications including broad characterization of the epitopes recognized by
antibodies. While the processing of HTS data, in general, is already well
supported, there are very few software tools that have been developed for
working with data generated in these highly-multiplexed peptide assays. In
order to fill this gap, we present PepSIRF (Peptide-based Serological Immune
Response Framework), which is a flexible and comprehensive software package
designed specifically for the analysis of HTS data from highly-multiplexed
peptide-based assays.


  The advent of DNA and RNA sequencing has revolutionized the study of genomics
and molecular biology. Next generation sequencing (NGS) technologies like
Illumina, Ion Torrent, SOLiD sequencing etc. have brought about a quick and
cheap way to sequence genomes. Recently, third generation sequencing (TGS)
technologies like PacBio and Oxford Nanopore Technology (ONT) have also been
developed. Different technologies use different underlying methods for
sequencing and are prone to different error rates. Though many tools exist for
error correction of sequencing data from NGS and TGS methods, no standard
method is available yet to evaluate the accuracy and effectiveness of these
error-correction tools. In this study, we present a Software Package for Error
Correction Tool Assessment on nuCLEic acid sequences (SPECTACLE) providing
comprehensive algorithms to evaluate error-correction methods for DNA and RNA
sequencing, for NGS and TGS platforms. We also present a compilation of
sequencing datasets for Illumina, PacBio and ONT platforms that present
challenging scenarios for error-correction tools. Using these datasets and
SPECTACLE, we evaluate the performance of 23 different error-correction tools
and present unique and helpful insights into their strengths and weaknesses. We
hope that our methodology will standardize the evaluation of DNA and RNA
error-correction tools in the future.


  Forecasting the effect of COVID-19 is essential to design policies that may
prepare us to handle the pandemic. Many methods have already been proposed,
particularly, to forecast reported cases and deaths at country-level and
state-level. Many of these methods are based on traditional epidemiological
model which rely on simulations or Bayesian inference to simultaneously learn
many parameters at a time. This makes them prone to over-fitting and slow
execution. We propose an extension to our model SIkJ$\alpha$ to forecast deaths
and show that it can consider the effect of many complexities of the epidemic
process and yet be simplified to a few parameters that are learned using fast
linear regressions. We also present an evaluation of our method against seven
approaches currently being used by the CDC, based on their two weeks forecast
at various times during the pandemic. We demonstrate that our method achieves
better root mean squared error compared to these seven approaches during
majority of the evaluation period. Further, on a 2 core desktop machine, our
approach takes only 3.18s to tune hyper-parameters, learn parameters and
generate 100 days of forecasts of reported cases and deaths for all the states
in the US. The total execution time for 184 countries is 11.83s and for all the
US counties ($>$ 3000) is 101.03s.


  Computational prediction of RNA structures is an important problem in
computational structural biology. Studies of RNA structure formation often
assume that the process starts from a fully synthesized sequence. Experimental
evidence, however, has shown that RNA folds concurrently with its elongation.
We investigate RNA secondary structure formation, including pseudoknots, that
takes into account the cotranscriptional effects. We propose a
single-nucleotide resolution kinetic model of the folding process of RNA
molecules, where the polymerase-driven elongation of an RNA strand by a new
nucleotide is included as a primitive operation, together with a stochastic
simulation method that implements this folding concurrently with the
transcriptional synthesis. Numerical case studies show that our
cotranscriptional RNA folding model can predict the formation of conformations
that are favored in actual biological systems. Our new computational tool can
thus provide quantitative predictions and offer useful insights into the
kinetics of RNA folding.


  Homeostasis refers to a phenomenon whereby the output $x_o$ of a system is
approximately constant on variation of an input $\mathcal{I}$. Homeostasis
occurs frequently in biochemical networks and in other networks of interacting
elements where mathematical models are based on differential equations
associated to the network. These networks can be abstracted as digraphs
$\mathcal{G}$ with a distinguished input node $\iota$, a different
distinguished output node $o$, and a number of regulatory nodes
$\rho_1,\ldots,\rho_n$. In these models the input-output map $x_o(\mathcal{I})$
is defined by a stable equilibrium $X_0$ at $\mathcal{I}_0$. Stability implies
that there is a stable equilibrium $X(\mathcal{I})$ for each $\mathcal{I}$ near
$\mathcal{I}_0$ and infinitesimal homeostasis occurs at $\mathcal{I}_0$ when
$(dx_o/d\mathcal{I})(\mathcal{I}_0) = 0$. We show that there is an
$(n+1)\times(n+1)$ homeostasis matrix $H(\mathcal{I})$ for which
$dx_o/d\mathcal{I} = 0$ if and only if $\det(H) = 0$. We note that the entries
in $H$ are linearized couplings and $\det(H)$ is a homogeneous polynomial of
degree $n+1$ in these entries. We use combinatorial matrix theory to factor the
polynomial $\det(H)$ and thereby determine a menu of different types of
possible homeostasis associated with each digraph $\mathcal{G}$. Specifically,
we prove that each factor corresponds to a subnetwork of $\mathcal{G}$. The
factors divide into two combinatorially defined classes: structural and
appendage. Structural factors correspond to feedforward motifs and appendage
factors correspond to feedback motifs. Finally, we discover an algorithm for
determining the homeostasis subnetwork motif corresponding to each factor of
$\det(H)$ without performing numerical simulations on model equations. The
algorithm allows us to classify low degree factors of $\det(H)$.


  Despite its historical and biological stability, the sex ratio at birth (SRB)
has risen in parts of the world in the last several decades. The resultant
demographic consequences, mostly on sex imbalance, are well documented,
typically including "missing girls/women" and "marriage squeeze." However, the
SRB-induced impact on demographic dynamics, particularly its underlying
mechanism, has not been explored in depth. We aim to investigate the impact of
the SRB rise on the size, structure, and growth of a population, particularly
emphasizing on population aging. We provide a simple framework, derived from
classical stable population models, to analyze how the SRB rise can reduce the
population size and make the population old. We demonstrate that the cohorts
born with a higher SRB are smaller in size than those with a lower SRB. As the
affected cohorts are born into the population, their smaller size will reduce
the total population size, thereby lifting the fraction of old people that were
born with the original SRB and have the same size as before. The resultant
population aging speed increases as the cohorts with the new SRB take an
increasing share of the population. This study adds that, in addition to
fertility and mortality, the SRB can be a driving factor of population
dynamics, especially when it moves far above normal biological levels.


  Networks are abundant in the life sciences. Outstanding challenges include
how to characterize similarities between networks, and in extension how to
integrate information across networks. Yet, network alignment remains a core
algorithmic problem. Here, we present a novel learning algorithm called
evolutionary heat diffusion-based network alignment (EDNA) to address this
challenge. EDNA uses the diffusion signal as a proxy for computing node
similarities between networks. Comparing EDNA with state-of-the-art algorithms
on a popular protein-protein interaction network dataset, using four different
evaluation metrics, we achieve (i) the most accurate alignments, (ii) increased
robustness against noise, and (iii) superior scaling capacity. The EDNA
algorithm is versatile in that other available network alignments/embeddings
can be used as an initial baseline alignment, and then EDNA works as a wrapper
around them by running the evolutionary diffusion on top of them. In
conclusion, EDNA outperforms state-of-the-art methods for network alignment,
thus setting the stage for large-scale comparison and integration of networks.


  Genome-Wide Association Studies (GWAS) explain only a small fraction of
heritability for most complex human phenotypes. Genomic heritability estimates
the variance explained by the SNPs on the whole genome using mixed models and
accounts for the many small contributions of SNPs in the explanation of a
phenotype.
  This paper approaches heritability from a machine learning perspective, and
examines the close link between mixed models and ridge regression. Our
contribution is twofold. First, we propose estimating genomic heritability
using a predictive approach via ridge regression and Generalized Cross
Validation (GCV). We show that this is consistent with classical mixed model
based estimation. Second, we derive simple formulae that express prediction
accuracy as a function of the ratio n/p, where n is the population size and p
the total number of SNPs. These formulae clearly show that a high heritability
does not imply an accurate prediction when p>n.
  Both the estimation of heritability via GCV and the prediction accuracy
formulae are validated using simulated data and real data from UK Biobank.


  A configuration model for epidemic spread, based on a scale-free network, is
introduced. We obtain the analytical solutions describing both unstable and
stable dynamics of the epidemic spreading, and demonstrate how these regimes
can interchange during the epidemic. We apply the model to the COVID-19 case
and demonstrate the predictive features of our model.


  The rheology of biological cells is not only determined by their cytoskeletal
networks, but also by the molecular motors that crosslink and contract them.
Recently it has been found that the assemblies of myosin II molecular motors in
non-muscle cells are mixtures of fast and slow motor variants. Using computer
simulations and an analytical mean field theory of a crossbridge model for
myosin II motors, here we show that such motor ensembles effectively behave as
active Maxwell elements. We calculate storage and loss moduli as a function of
the model parameters and show that the rheological properties cross over from
viscous to elastic as one increases the ratio of slow to fast motors. This
suggests that cells tune their mechanical properties by regulating the
composition of their myosin assemblies.


  Multiple studies have suggested the spread of COVID-19 is affected by factors
such as climate, BCG vaccinations, pollution and blood type. We perform a joint
study of these factors using the death growth rates of 40 regions worldwide
with both machine learning and Bayesian methods. We find weak, non-significant
(< 3$\sigma$) evidence for temperature and relative humidity as factors in the
spread of COVID-19 but little or no evidence for BCG vaccination prevalence or
$\text{PM}_{2.5}$ pollution. The only variable detected at a statistically
significant level (>3$\sigma$) is the rate of positive COVID-19 tests, with
higher positive rates correlating with higher daily growth of deaths.


  The high fidelity of DNA polymerase is critical for the faithful replication
of genomic DNA. Several approaches were proposed to quantify the fidelity of
DNA polymerase. Direct measurements of the error frequency of the replication
products definitely give the true fidelity but turn out very hard to implement.
Two biochemical kinetic approaches, the steady-state assay and the
transient-state assay, were then suggested and widely adopted. In these assays,
the error frequency is indirectly estimated by using the steady-state or the
transient-state kinetic theory combined with the measured kinetic rates.
However, whether these indirectly estimated fidelities are equivalent to the
true fidelity has never been clarified theoretically, and in particular there
are different strategies to quantify the proofreading efficiency of DNAP but
often lead to inconsistent results. The reason for all these confusions is that
it's mathematically challenging to formulate a rigorous and general theory of
the true fidelity. Recently we have succeeded to establish such a theoretical
framework. In this paper, we develop this theory to make a comprehensive
examination on the theoretical foundation of the kinetic assays and the
relation between fidelities obtained by different methods. We conclude that
while the steady-state assay and the transient-state assay can always measure
the true fidelity of exonuclease-deficient DNA polymerases, they only do so for
exonuclease-efficient DNA polymerases conditionally (the proper way to use
these assays to quantify the proofreading efficiency is also suggested). We
thus propose a new kinetic approach, the single-molecule assay, which
indirectly but precisely characterizes the true fidelity of either
exonuclease-deficient or exonuclease-efficient DNA polymerases.


  The SARS-CoV-2 pandemic has transformed the way that the world functions.
Health issues and population safety have driven countries' economies to a
critical state; therefore, sustaining economic activity while keeping workers
safe has become a worldwide goal. In this paper, we present a novel safety
protocol based on rapid antibody testing (ABT). Using discrete event
simulation, we evaluated its performance on the cumulative number of infected
workers, effective reproductive number ($R_e$), and active work force within a
company. Using a synthetic experiment, we showed that ABT twice a week (ABT 3)
performed the best, detecting 5.7% of infected workers, compared to 16.9% when
no ABT was applied. $R_e$ was reduced from 1.75 to 0.84, with a slight decrease
in the active workers within the firm. A sensitivity analysis on the duration
of the shedding period and sensitivity of ABT was performed and led to the same
qualitative conclusions. We applied this protocol in a Chilean winery: the
estimation of the initial $R_e$ of 1.3 was reduced to 0.7 when the ABT 3
protocol was implemented, with a 27% decrease in the number of infected
workers. Although ABT is not approved for COVID-19 diagnosis, our study shows
that upgraded safety standards can already be implemented in workspaces.


  regulation largely unexplored, in part due to methodological limitations.
Indeed, we review evidence demonstrating that commonly used methods, such as
transcriptomics, are inadequate because the variability in mRNAs coding for
ribosomal proteins (RP) does not necessarily correspond to RP variability. Thus
protein remodeling of ribosomes should be investigated by methods that allow
direct quantification of RPs, ideally of isolated ribosomes. We review such
methods, focusing on mass spectrometry and emphasizing method-specific biases
and approaches to control these biases. We argue that using multiple
complementary methods can help reduce the danger of interpreting reproducible
systematic biases as evidence for ribosome remodeling.


  Background: The COVID-19 pandemic clock is ticking and the survival of many
of mankind's modern institutions and or survival of many individuals is at
stake. There is a need for treatments to significantly reduce the morbidity and
mortality of COVID-19. Hence, we delved deep into the SARS-CoV-2 genome, which
is the virus that has caused COVID-19. SARS-CoV-2 is from the same family as
SARS-CoV in which three small viral RNAs (svRNA) were recently identified;
those svRNAs play a significant role in the virus pathogenesis in mice.
  Contribution: In this paper, we report potential orthologs of those three
svRNAs in the SARS-CoV-2 genome. Instead of off-the-shelf search and alignment
algorithms, which failed to discover the orthologs, we used a special alignment
scoring that does not penalize C/T and A/G mismatches as much as the other
mutations. RNA bases C and U both can bind to G; similarly, A and G both can
bind to U, hence, our scoring. We also validate this hypothesis using a novel,
independent computational experiment. To validate our results, we confirmed the
discovered orthologs are fully conserved in all the tested publicly available
genomes of various strains of SARS-CoV-2; the loci at which the SARS-CoV-2
orthologs occur are close to the loci at which SARS-CoV svRNAs occur. We also
report potential targets for these svRNAs. We hypothesize that the discovered
orthologs play a role in pathogenesis of SARS-CoV-2, and therefore,
antagomir-mediated inhibition of these SARS-CoV-2 svRNAs inhibits COVID-19.


  The novel Coronavirus (COVID-19) incidence in India is currently experiencing
exponential rise but with apparent spatial variation in growth rate and
doubling time rate. We classify the states into five clusters with low to the
high-risk category and study how the different states moved from one cluster to
the other since the onset of the first case on $30^{th}$ January 2020 till the
end of unlock 1 that is $30^{th}$ June 2020. We have implemented a new
clustering technique called the incrementalKMN (Prasad, R. K., Sarmah, R.,
Chakraborty, S.(2019))


  This work is an attempt to transfer to biology the methods developed in
physics for formulating and solving the kinetic equations in which the kernel
of the integral operator in spatial coordinates is slowly decreasing with
increasing distance and belongs to the class of Levy distributions. An
algorithm is proposed for the reconstruction of the step-length probability
density function (PDF) on a moderate number of trajectories of biological
objects (migrants) and for the derivation of the Green's function of the
corresponding integro-differential kinetic equation for the density of migrants
in the entire space-time range, including the construction of an approximate
self-similar solution. A wide class of time-dependent superdiffusion processes
with a model power-law step-length PDF is considered, which corresponds to
"Levy walks with rests" for given values of the migrant's constant velocity and
the average time T of the migrant's stay between runs. The algorithm is tested
within the framework of a synthetic diagnostics, consisting in the generation
of artificial experimental data for trajectories of migrants and the subsequent
reconstruction of the parameters of the step-length PDF and T. For different
volumes of synthetic data, to obtain a general idea of the distributions under
study (non-parametric case) and to evaluate the accuracy of recovering the
parameters of the PDF (in the case of a parametric representation), the method
of balanced identification is used. The approximate self-similar solution for
the parameters of step-length PDF and T is shown to provide reasonable accuracy
of the space-time evolution of migrant's density.


  Molecular data from tumor profiles is high dimensional. Tumor profiles can be
characterized by tens of thousands of gene expression features. Due to the size
of the gene expression feature set machine learning methods are exposed to
noisy variables and complexity. Tumor types present heterogeneity and can be
subdivided in tumor subtypes. In many cases tumor data does not include tumor
subtype labeling thus unsupervised learning methods are necessary for tumor
subtype discovery. This work aims to learn meaningful and low dimensional
representations of tumor samples and find tumor subtype clusters while keeping
biological signatures without using tumor labels. The proposed method named
Latent Kernel Feature Selection (LKFS) is an unsupervised approach for gene
selection in tumor gene expression profiles. By using Autoencoders a low
dimensional and denoised latent space is learned as a target representation to
guide a Multiple Kernel Learning model that selects a subset of genes. By using
the selected genes a clustering method is used to group samples. In order to
evaluate the performance of the proposed unsupervised feature selection method
the obtained features and clusters are analyzed by clinical significance. The
proposed method has been applied on three tumor datasets which are Brain, Renal
and Lung, each one composed by two tumor subtypes. When compared with benchmark
unsupervised feature selection methods the results obtained by the proposed
method reveal lower redundancy in the selected features and a better clustering
performance.