  Between pandemics, the influenza virus exhibits periods of incremental
evolution via a process known as antigenic drift. This process gives rise to a
sequence of strains of the pathogen that are continuously replaced by newer
strains, preventing a build up of immunity in the host population. In this
paper, a parsimonious epidemic model is defined that attempts to capture the
dynamics of evolving strains within a host population. The `evolving strains'
epidemic model has many properties that lie in-between the
Susceptible-Infected-Susceptible and the Susceptible-Infected-Removed epidemic
models, due to the fact that individuals can only be infected by each strain
once, but remain susceptible to reinfection by newly emerged strains. Coupling
results are used to identify key properties, such as the time to extinction. A
range of reproduction numbers are explored to characterize the model, including
a novel quasi-stationary reproduction number that can be used to describe the
re-emergence of the pathogen into a population with `average' levels of strain
immunity, analogous to the beginning of the winter peak in influenza. Finally
the quasi-stationary distribution of the evolving strains model is explored via
simulation.


  This research was done during the DOMath program at Duke University from May
18 to July 10, 2020. At the time, Duke and other universities across the
country were wrestling with the question of how to safely welcome students back
to campus in the Fall. Because of this, our project focused on using
mathematical models to evaluate strategies to suppress the spread of the virus
on campus, specifically in dorms and in classrooms. For dorms, we show that
giving students single rooms rather than double rooms can substantially reduce
virus spread. For classrooms, we show that moving classes with size above some
cutoff online can make the basic reproduction number $R_0<1$, preventing a wide
spread epidemic. The cutoff will depend on the contagiousness of the disease in
classrooms.


  A contact-tracing strategy has been deemed necessary to contain the spread of
COVID-19 following the relaxation of lockdown measures. Using an agent-based
model, we explore one of the technology-based strategies proposed, a
contact-tracing smartphone app. The model simulates the spread of COVID-19 in a
population of agents on an urban scale. Agents are heterogeneous in their
characteristics and are linked in a multi-layered network representing the
social structure - including households, friendships, employment and schools.
We explore the interplay of various adoption rates of the contact-tracing app,
different levels of testing capacity, and behavioural factors to assess the
impact on the epidemic. Results suggest that a contact tracing app can
contribute substantially to reducing infection rates in the population when
accompanied by a sufficient testing capacity or when the testing policy
prioritises symptomatic cases. As user rate increases, prevalence of infection
decreases. With that, when symptomatic cases are not prioritised for testing, a
high rate of app users can generate an extensive increase in the demand for
testing, which, if not met with adequate supply, may render the app
counterproductive. This points to the crucial role of an efficient testing
policy and the necessity to upscale testing capacity.


  Inspired by recent work on anti-covid-19 drugs \cite{2} here we study the
Quantitative-structure property relationships(QSPR) of phytochemicals screened
against SARS-CoV-2 $3CL^{pro}$ with the help of topological indices like the
first Zagreb index $M_{1}$, second Zagreb index $M_{2}$, Randi$\acute{c}$ index
$R$, Balban index $J$ and sum-connectivity index $SCI(G)$. Our study has
raveled that the sum-connectivity index $(SCI)$ and the first Zagreb index
$(M_{1})$ are two important parameters to predict the molecular weight and the
topological polar surface area of phytochemicals respectively.


  A well-known characteristic of pandemics such as COVID-19 is the high level
of transmission heterogeneity in the infection spread: not all infected
individuals spread the disease at the same rate and some individuals
(superspreaders) are responsible for most of the infections. To quantify this
phenomenon requires the analysis of the effect of the variance and higher
moments of the infection distribution. Working in the framework of stochastic
branching processes, we derive an approximate analytical formula for the
probability of an outbreak in the high variance regime of the infection
distribution, verify it numerically and analyze its regime of validity in
various examples. We show that it is possible for an outbreak not to occur in
the high variance regime even when the basic reproduction number $R_0$ is
larger than one and discuss the implications of our results for COVID-19 and
other pandemics.


  The cell nucleus houses the chromosomes, which are linked to a soft shell of
lamin filaments. Experiments indicate that correlated chromosome dynamics and
nuclear shape fluctuations arise from motor activity. To identify the physical
mechanisms, we develop a model of an active, crosslinked Rouse chain bound to a
polymeric shell. System-sized correlated motions occur but require both motor
activity {\it and} crosslinks. Contractile motors, in particular, enhance
chromosome dynamics by driving anomalous density fluctuations. Nuclear shape
fluctuations depend on motor strength, crosslinking, and chromosome-lamina
binding. Therefore, complex chromatin dynamics and nuclear shape emerge from a
minimal, active chromosome-lamina system.


  The transmission and evolution of severe acute respiratory syndrome
coronavirus 2 (SARS-CoV-2) are of paramount importance to the controlling and
combating of coronavirus disease 2019 (COVID-19) pandemic. Currently, near
15,000 SARS-CoV-2 single mutations have been recorded, having a great
ramification to the development of diagnostics, vaccines, antibody therapies,
and drugs. However, little is known about SARS-CoV-2 evolutionary
characteristics and general trend. In this work, we present a comprehensive
genotyping analysis of existing SARS-CoV-2 mutations. We reveal that host
immune response via APOBEC and ADAR gene editing gives rise to near 65\% of
recorded mutations. Additionally, we show that children under age five and the
elderly may be at high risk from COVID-19 because of their overreacting to the
viral infection. Moreover, we uncover that populations of Oceania and Africa
react significantly more intensively to SARS-CoV-2 infection than those of
Europe and Asia, which may explain why African Americans were shown to be at
increased risk of dying from COVID-19, in addition to their high risk of
getting sick from COVID-19 caused by systemic health and social inequities.
Finally, our study indicates that for two viral genome sequences of the same
origin, their evolution order may be determined from the ratio of mutation type
C$>$T over T$>$C.


  Fenna-Matthews-Olson (FMO) bacteriochlorophylls (BChls) are molecules
responsible of the high efficiency energy transfer in the photosynthetic
process of green sulfur bacteria, controversially associated to quantum
phenomena of long lived coherence. This phenomenon is modelled using Quantum
Open Systems (QOS) without included memory effects of the surrounding
approximated as a phonon bath on thermal equilibrium. This work applies the
Hierarchical Equations of Motion method (HEOM), a non-Markovian approach, in
the modelling of the system evolution of the FMO complex to perform predictions
about the coherence time scales together with global and semi-local
entanglement during the quantum excitation.


  COVID-19 pandemic has brought the whole world to a stand-still over the last
few months. In particular the pace at which pandemic has spread has taken
everybody off-guard. The Governments across the world have responded by
imposing lock-downs, stopping/restricting travel and mandating social
distancing. On the positive side there is wide availability of information on
active cases, recoveries and deaths collected daily across regions. However,
what has been particularly challenging is to track the spread of the disease by
asymptomatic carriers termed as super-spreaders. In this paper we look at
applying Hidden Markov Model to get a better assessment of extent of spread.
The outcome of such analysis can be useful to Governments to design the
required interventions/responses in a calibrated manner. The data we have
chosen to analyze pertains to Indian scenario.


  Liquid liquid phase separation (LLPS) mediated by pi-cation bonds between
tyrosine and arginine residues are of biological importance. To understand the
interactions between proteins in the condensed phase in close analogy to
complex coacervation, we run multiple umbrella calculations between oligomers
containing tyrosine (pY) and arginine (pR). We find pR-pY complexation to be
energetically driven. Metadynamics simulations reveal that this energy of
complexation comes primarily from pi-cation bonds. On running free energy
calculation for the second binding step of complex coacervation, we find
striking similarities between this process and pi-mediated LLPS. These
calculations lead us to believe that contrary to the common notion, complex
coacervation as whole, which involves an entropic complexation followed by an
energetic aggregation is not invoked by proteins containing arginine and
tyrosine residues. Rather, the latter step in itself, in which neutral polyion
pairs aggregate together is the correct mechanism for pi-cation mediated LLPS.


  The relation between the number of passengers in the main stations and the
infection rate of COVID19 in Tokyo is empirically studied. Our analysis based
on conventional compartment model suggests: 1) Average time from the true day
of infection to the day the infections are reported is about $15$ days. 2) The
scaling relation between the density of active population and the infection
rate suggests that the increase of infection rate is linear to the active
population rather than quadratic, as that is assumed in the conventional SIR
model. 3) Notable deviations from the overall scaling relation seems to
correspond to the change of the peoples's behavior in response to the public
announcements of action regulation.


  Rooted phylogenetic networks provide a way to describe species' relationships
when evolution departs from the simple model of a tree. However, networks
inferred from genomic data can be highly tangled, making it difficult to
discern the main reticulation signals present. In this paper, we describe a
natural way to transform any rooted phylogenetic network into a simpler
canonical network, which has desirable mathematical and computational
properties, and is based only on the 'visible' nodes in the original network.
The method has been implemented and we demonstrate its application to some
examples.


  We are asking: why are the dynamics and control of Covid-19 most interesting
for mathematicians and why are mathematicians urgently needed for controlling
the pandemic?
  First we present our comments in a Bottom-up approach, i.e., following the
events from their beginning as they evolved through time. They happened
differently in different countries, and the main objective of the first part is
to compare these evolutions in a few selected countries with each other.
  The second part of the article is not "country-oriented" but
"problem-oriented". From a given problem we go Top-down to its solutions and
their applications in concrete situations. We have organized this part by the
mathematical methods that play a role in their solution. We give an overview of
the main branches of mathematics that play a role and sketch the most frequent
applications, emphasising mathematical pattern analysis in laboratory work and
statistical-mathematical models in judging the quality of tests; demographic
methods in the collection of data; different ways to model the evolution of the
pandemic mathematically; and clinical epidemiology in attempts to develop a
vaccine.


  The presence of oscillations in aggregated COVID-19 data not only raises
questions about the data's accuracy, it hinders understanding of the pandemic.
Spectral analysis is used to reveal additional properties of the data, and the
oscillations are replicated using sinusoidal resynthesis. The precise behavior
of the seven-day moving average is also discussed, specifically, the cause of
its jaggedness and the phase error it introduces. In comparison, other
filtering techniques and Fourier processing produce superior smoothing and have
zero phase error. Both of these are presented, and they are extended to isolate
several frequency ranges. This extracts some of the same short-term variability
that is resynthesized, and it shows that fluctuations with periods between 8
and 21 days are present in U.S. mortality data. These methods have applications
that include modeling epidemiological time-series data as well as identifying
less obvious properties of the data.


  Subjective Experience (SE) is part of the ancient mind-body problem, which
continues to be one of deepest mysteries of science. Despite major advances in
many fields, there is still no plausible causal link between SE and its
realization in the body. The core issue is the incompatibility of objective
(3rd person) public science with subjective (1st person) private experience.
Any scientific approach to SE assumes that it arose from extended evolutionary
processes and that examining evolutionary history should help us understand it.
While the core mystery remains, converging evidence from theoretical,
experimental, and computational studies yields strong constraints on SE and
some suggestions for further research. All animals confront many of the same
fitness challenges. They all need some kind of internal model to relate their
life goals and actionable sensed information to action. We understand the
evolution of the bodily aspects of human perception and emotion, but not the
SE. The first evolutionary evidence for SE appears in vertebrates and much of
its neural substrate and simulation mechanism is preserved in mammals and
humans. People exhibit the same phenomena, but there are remaining mysteries of
everyday experience that are demonstrably incompatible with current
neuroscience. In spite of this limitation, there is considerable progress on
understanding the role of SE in the success of prostheses.


  The Italian government has been one of the most responsive to COVID-19
emergency, through the adoption of quick and increasingly stringent measures to
contain the outbreak. Despite this, Italy has suffered a huge human and social
cost, especially in Lombardy. The aim of this paper is dual: i) first, to
investigate the reasons of the case fatality rate (CFR) differences across
Italian 20 regions and 107 provinces, using a multivariate OLS regression
approach; and ii) second, to build a taxonomy of provinces with similar
mortality risk of COVID-19, by using the Ward hierarchical agglomerative
clustering method. I considered health system metrics, environmental pollution,
climatic conditions, demographic variables, and three ad hoc indexes that
represent the health system saturation. The results showed that overall health
care efficiency, physician density, and average temperature helped to reduce
the CFR. By the contrary, population aged 70 and above, car and firm density,
level of air pollutants (NO2, O3, PM10, and PM2.5), relative average humidity,
COVID-19 prevalence, and all three indexes of health system saturation were
positively associated with the CFR. Population density, social vertical
integration, and altitude were not statistically significant. In particular,
the risk of dying increases with age, as 90 years old and above had a
three-fold greater risk than the 80 to 89 years old and four-fold greater risk
than 70 to 79 years old. Moreover, the cluster analysis showed that the highest
mortality risk was concentrated in the north of the country, while the lowest
risk was associated with southern provinces. Finally, since prevalence and
health system saturation indexes played the most important role in explaining
the CFR variability, a significant part of the latter may have been caused by
the massive stress of the Italian health system.


  Policy-makers require data-driven tools to assess the spread of COVID-19 and
inform the public of their risk of infection on an ongoing basis. We propose a
rigorous hybrid model-and-data-driven approach to risk scoring based on a
time-varying SIR epidemic model that ultimately yields a simplified color-coded
risk level for each community. The risk score $\Gamma_t$ that we propose is
proportional to the probability of someone currently healthy getting infected
in the next 24 hours. We show how this risk score can be estimated using
another useful metric of infection spread, $R_t$, the time-varying average
reproduction number which indicates the average number of individuals an
infected person would infect in turn. The proposed approach also allows for
quantification of uncertainty in the estimates of $R_t$ and $\Gamma_t$ in the
form of confidence intervals. Code and data from our effort have been
open-sourced and are being applied to assess and communicate the risk of
infection in the City and County of Los Angeles.


  Epidemics generally spread through a succession of waves that reflect factors
on multiple timescales. On short timescales, super-spreading events lead to
burstiness and overdispersion, while long-term persistent heterogeneity in
susceptibility is expected to lead to a reduction in the infection peak and the
herd immunity threshold (HIT). Here, we develop a general approach to encompass
both timescales, including time variations in individual social activity, and
demonstrate how to incorporate them phenomenologically into a wide class of
epidemiological models through parameterization. We derive a non-linear
dependence of the effective reproduction number Re on the susceptible
population fraction S. We show that a state of transient collective immunity
(TCI) emerges well below the HIT during early, high-paced stages of the
epidemic. However, this is a fragile state that wanes over time due to changing
levels of social activity, and so the infection peak is not an indication of
herd immunity: subsequent waves can and will emerge due to behavioral changes
in the population, driven (e.g.) by seasonal factors. Transient and long-term
levels of heterogeneity are estimated by using empirical data from the COVID-19
epidemic as well as from real-life face-to-face contact networks. These results
suggest that the hardest-hit areas, such as NYC, have achieved TCI following
the first wave of the epidemic, but likely remain below the long-term HIT.
Thus, in contrast to some previous claims, these regions can still experience
subsequent waves.


  In order to investigate the effectiveness of lockdown and social distancing
restrictions, which have been widely carried out as policy choice to curb the
ongoing COVID-19 pandemic around the world, we formulate and discuss a staged
and weighed networked system based on a classical SEAIR epidemiological model.
Five stages have been taken into consideration according to four-tier response
to Public Health Crisis, which comes from the National Contingency Plan in
China. Staggered basic reproduction number has been derived and we evaluate the
effectiveness of lockdown and social distancing policies under different
scenarios among 19 cities/regions in mainland China. Further, we estimate the
infection risk associated with the sequential release based on population
mobility between cities and the intensity of some non-pharmaceutical
interventions. Our results reveal that Level I public health emergency response
is necessary for high-risk cities, which can flatten the COVID-19 curve
effectively and quickly. Moreover, properly designed staggered-release policies
are extremely significant for the prevention and control of COVID-19,
furthermore, beneficial to economic activities and social stability and
development.


  The spread of COVID-19 during the initial phase of the first half of 2020 was
curtailed to a larger or lesser extent through measures of social distancing
imposed by most countries. In this work, we link directly, through machine
learning techniques, infection data at a country level to a single number that
signifies social distancing effectiveness. We assume that the standard SIR
model gives a reasonable description of the dynamics of spreading, and thus the
social distancing aspect can be modeled through time-dependent infection rates
that are imposed externally. We use an exponential ansatz to analyze the SIR
model, find an exact solution for the time-independent infection rate, and
derive a simple first-order differential equation for the time-dependent
infection rate as a function of the infected population. Using infected number
data from the "first wave" of the infection from eight countries, and through
physics-informed machine learning, we extract the degree of linear dependence
in social distancing that led to the specific infections. We find that in the
two extremes are Greece, with the highest decay slope on one side, and the US
on the other with a practically flat "decay". The hierarchy of slopes is
compatible with the effectiveness of the pandemic containment in each country.
Finally, we train our network with data after the end of the analyzed period,
and we make week-long predictions for the current phase of the infection that
appear to be very close to the actual infection values.


  Transmission risk of air-borne diseases in public transportation systems is a
concern. The paper proposes a modified Wells-Riley model for risk analysis in
public transportation systems to capture the passenger flow characteristics,
including spatial and temporal patterns in terms of number of boarding,
alighting passengers, and number of infectors. The model is utilized to assess
overall risk as a function of OD flows, actual operations, and factors such as
mask wearing, and ventilation. The model is integrated with a microscopic
simulation model of subway operations (SimMETRO). Using actual data from a
subway system, a case study explores the impact of different factors on
transmission risk, including mask-wearing, ventilation rates, infectiousness
levels of disease and carrier rates. In general, mask-wearing and ventilation
are effective under various demand levels, infectiousness levels, and carrier
rates. Mask-wearing is more effective in mitigating risks. Impacts from
operations and service frequency are also evaluated, emphasizing the importance
of maintaining reliable, frequent operations in lowering transmission risks.
Risk spatial patterns are also explored, highlighting locations of higher risk.


  Motivation: Protein-protein interactions (PPI) are critical to the function
of proteins in both normal and diseased cells, and many critical protein
functions are mediated by interactions.Knowledge of the nature of these
interactions is important for the construction of networks to analyse
biological data. However, only a small percentage of PPIs captured in protein
interaction databases have annotations of function available, e.g. only 4% of
PPI are functionally annotated in the IntAct database. Here, we aim to label
the function type of PPIs by extracting relationships described in PubMed
abstracts.
  Method: We create a weakly supervised dataset from the IntAct PPI database
containing interacting protein pairs with annotated function and associated
abstracts from the PubMed database. We apply a state-of-the-art deep learning
technique for biomedical natural language processing tasks, BioBERT, to build a
model - dubbed PPI-BioBERT - for identifying the function of PPIs. In order to
extract high quality PPI functions at large scale, we use an ensemble of
PPI-BioBERT models to improve uncertainty estimation and apply an interaction
type-specific threshold to counteract the effects of variations in the number
of training samples per interaction type.
  Results: We scan 18 million PubMed abstracts to automatically identify 3253
new typed PPIs, including phosphorylation and acetylation interactions, with an
overall precision of 46% (87% for acetylation) based on a human-reviewed
sample. This work demonstrates that analysis of biomedical abstracts for PPI
function extraction is a feasible approach to substantially increasing the
number of interactions annotated with function captured in online databases.


  We investigate the distribution of bubble lifetimes and bubble lengths in DNA
at physiological temperature, by performing extensive molecular dynamics
simulations with the Peyrard-Bishop-Dauxois (PBD) model, as well as an extended
version (ePBD) having a sequence-dependent stacking interaction, emphasizing
the effect of the sequences' guanine-cytosine (GC)/adenine-thymine (AT) content
on these distributions. For both models we find that base pair-dependent (GC vs
AT) thresholds for considering complementary nucleotides to be separated are
able to reproduce the observed dependence of the melting temperature on the GC
content of the DNA sequence. Using these thresholds for base pair openings, we
obtain bubble lifetime distributions for bubbles of lengths up to ten base
pairs as the GC content of the sequences is varied, which are accurately fitted
with stretched exponential functions. We find that for both models the average
bubble lifetime decreases with increasing either the bubble length or the GC
content. In addition, the obtained bubble length distributions are also fitted
by appropriate stretched exponential functions and our results show that short
bubbles have similar likelihoods for any GC content, but longer ones are
substantially more likely to occur in AT-rich sequences. We also show that the
ePBD model permits more, longer-lived, bubbles than the PBD system.


  The COVID-19 Pandemic has been described as the global challenge of our time,
an enormous human tragedy with dramatic economic impacts. This paper describes
the response and expected recovery process for Western Australia, where a rapid
and effective response was implemented. This has enabled an early transition
into an expected recovery both in health and economic terms. The positive
lessons learned from this experience are documented as they emerge in order to
support other states and nations as they address this issue globally in the
near-term and consider enduring improvements for the longer term. While the
authors have personal experience in the WA context, wider observations across
Australia and selected international benchmarks are also included. Key lessons
include the importance of good health advice in Australia's interest; timely,
synchronized and aligned action at all levels of government; a program of well
communicated, aligned health and economic measures which support all in society
allowing a very high level of appropriate community behaviour, ensuring the
health system was not overloaded; innovation in telehealth, testing, pandemic
modelling, and integrated operations which also allowed essential industries to
continue; and strong border and travel controls with highly effective isolation
preventing community spread, ultimately enabling rapid elimination of the
disease from the hospital system. In combination, these demonstrate that in the
case of Western Australia the result of first eliminating the disease from the
community, and then reopening the economy progressively at a strong pace, has
enabled a world leading outcome in both in health and economic terms. The
lessons from this experience are widely applicable, shareable both as
supporting service to other regions and through knowledge transfer.


  The de novo design of molecular structures using deep learning generative
models introduces an encouraging solution to drug discovery in the face of the
continuously increased cost of new drug development. From the generation of
original texts, images, and videos, to the scratching of novel molecular
structures, the incredible creativity of deep learning generative models
surprised us about the height machine intelligence can achieve. The purpose of
this paper is to review the latest advances in generative chemistry which
relies on generative modeling to expedite the drug discovery process. This
review starts with a brief history of artificial intelligence in drug discovery
to outline this emerging paradigm. Commonly used chemical databases, molecular
representations, and tools in cheminformatics and machine learning are covered
as the infrastructure for the generative chemistry. The detailed discussions on
utilizing cutting-edge generative architectures, including recurrent neural
network, variational autoencoder, adversarial autoencoder, and generative
adversarial network for compound generation are focused. Challenges and future
perspectives follow.


  The ongoing COVID-19 pandemic has affected most of the countries on Earth. It
has become a pandemic outbreak with more than 24 million confirmed infections
and above 840 thousand deaths worldwide. In this study, we consider a
mathematical model on COVID-19 transmission with the prosocial awareness
effect. The proposed model can have four equilibrium states based on different
parametric conditions. The local and global stability conditions for awareness
free, disease-free equilibrium is studied. Using Lyapunov function theory and
LaSalle Invariance Principle, the disease-free equilibrium is shown globally
asymptotically stable under some parametric constraints. The existence of
unique awareness free, endemic equilibrium and unique endemic equilibrium is
presented. We calibrate our proposed model parameters to fit daily cases and
deaths from Colombia. Sensitivity analysis indicates that the transmission rate
and learning factor related to awareness of susceptibles are very crucial for
reduction in disease related deaths. Finally, we assess the impact of prosocial
awareness during the outbreak and compare this strategy with popular control
measures. Results indicate that prosocial awareness has competitive potential
to flatten the curve.


  Eukaryotic transcription generally occurs in bursts of activity lasting
minutes to hours; however, state-of-the-art measurements have revealed that
many of the molecular processes that underlie bursting, such as transcription
factor binding to DNA, unfold on timescales of seconds. This temporal
disconnect lies at the heart of a broader challenge in physical biology of
predicting transcriptional outcomes and cellular decision-making from the
dynamics of underlying molecular processes. Here, we review how new dynamical
information about the processes underlying transcriptional control can be
combined with theoretical models that predict not only averaged transcriptional
dynamics, but also their variability, to formulate testable hypotheses about
the molecular mechanisms underlying transcriptional bursting and control.


  The current pandemic of SARS-CoV-2 has caused extensive damage to society.
The characterization of SARS-CoV-2 profiles has been addressed by researchers
globally with the aim of resolving this disruptive crisis. This investigation
process is indispensable for an understanding of how SARS-CoV-2 behaves in
human host cells. However, little is known about the systematic molecular
mechanisms involved in the effect of SARS-CoV-2 infection on human host cells.
Here, we have presented gene-to-gene regulatory networks in response to
SARS-CoV-2 using a Bayesian network. We examined the dynamic changes of the
SARS-CoV-2-purturbated networks established by our proposed framework for gene
network analysis, revealing that interferon signaling gradually switches to the
subsequent inflammatory-cytokine signaling cascades. Furthermore, we have
succeeded in capturing a COVID-19 patient-specific network in which
transduction of these signalings is coincidently induced. This enabled us to
explore local regulatory systems influenced by SARS-CoV-2 in host cells more
precisely at an individual level. Our panel of network analyses has provided
new insight into SARS-CoV-2 research from the perspective of cellular systems.


  Many species are unsustainable at small population densities (Allee Effect),
i.e., below a threshold named Allee threshold, the population decreases instead
of growing. In a closed local population, environmental fluctuations always
lead to extinction. Here, we show how, in spatially extended habitats,
dispersal can lead to a sustainable population in a region, provided the
amplitude of environmental fluctuations is below an extinction threshold. We
have identified two types of sustainable populations: high-density and
low-density populations (through a mean-field approximation, valid in the limit
of large dispersal length). Our results show that patches where population is
high, low or extinct, coexist when the population is close to global extinction
(even for homogeneous habitats). The extinction threshold is maximum for
characteristic dispersal distances much larger than the spatial scale of
synchrony of environmental fluctuations. The extinction threshold increases
proportionally to the square root of the dispersal rate and decreases with the
Allee threshold. The low-density population solution can allow understanding
difficulties in recovery after harvesting. This theoretical framework provides
a novel approach to address other factors, such as habitat fragmentation or
harvesting, impacting population resilience to environmental fluctuations.


  We develop an agent-based model on a network meant to capture features unique
to COVID-19 spread through a small residential college. We find that a safe
reopening requires strong policy from administrators combined with cautious
behavior from students. Strong policy includes weekly screening tests with
quick turnaround and halving the campus population. Cautious behavior from
students means wearing facemasks, socializing less, and showing up for COVID-19
testing. We also find that comprehensive testing and facemasks are the most
effective single interventions, building closures can lead to infection spikes
in other areas depending on student behavior, and faster return of test results
significantly reduces total infections.


  Controlling the spread of COVID-19 - even after a licensed vaccine is
available - requires the effective use of non-pharmaceutical interventions:
physical distancing, limits on group sizes, mask wearing, etc. To date, such
interventions have neither been uniformly nor systematically implemented in
most countries. For example, even when under strict stay-at-home orders,
numerous jurisdictions granted exceptions and/or were in close proximity to
locations with entirely different regulations in place. Here, we investigate
the impact of such geographic inconsistencies in epidemic control policies by
coupling search and mobility data to a simple mathematical model of SARS-COV2
transmission. Our results show that while stay-at-home orders decrease contacts
in most areas of the US, some specific activities and venues often see an
increase in attendance. Indeed, over the month of March 2020, between 10 and
30% of churches in the US saw increases in attendance; even as the total number
of visits to churches declined nationally. This heterogeneity, where certain
venues see substantial increases in attendance while others close, suggests
that closure can cause individuals to find an open venue, even if that requires
longer-distance travel. And, indeed, the average distance travelled to churches
in the US rose by 13% over the same period. Strikingly, our model reveals that
across a broad range of model parameters, partial measures can often be worse
than none at all where individuals not complying with policies by traveling to
neighboring areas can create epidemics when the outbreak would otherwise have
been controlled. Taken together, our data analysis and modelling results
highlight the potential unintended consequences of inconsistent epidemic
control policies and stress the importance of balancing the societal needs of a
population with the risk of an outbreak growing into a large epidemic.


  In this work, some new exact and approximate analytical solutions are
obtained for the SIR epidemic model, which is formulated in terms of
dimensionless variables and parameters. The susceptibles population (S) is in
this way explicitly related to the infectives population (I) using the Lambert
W function (both the principal and the secondary branches). A simple and
accurate relation for the fraction of the population that does not catch the
disease is also obtained. The explicit time dependences of the susceptibles,
infectives and removed populations, as well as that of the epidemic curve are
also modelled with good accuracy for any value of R0 (basic multiplication
number) using simple functions that are modified solutions of the R0 ->
infinity limiting case (logistic curve). It is also shown that for I0 << S0 the
effect of a change in the ratio I0/S0 on the population evolution curves
amounts to a time shift, their shape and relative position being unaffected.


  The formation and proliferation of protein aggregates play a central role in
a number of devastating neuro-degenerative diseases. Many experimental studies
indicate that the ability of existing aggregates to replicate is a key property
in generating their pathogenic effect across a range of diseases. However,
given the complexity of the process in vivo, no principled general approach
currently exists to obtain the rates of the fundamental steps that underlie
aggregate formation from measurements in living systems. In order to address
this challenge, here we present a general approach for analysing aggregation
kinetics that considers broad classes of processes that can be described by a
family of scaling solutions. Our approach is not limited only to fibrillar
aggregates, but applies to any aggregate shape. We show that the rates can
reliably be extracted by fitting of a simple logistic function, even from
experimental data in living systems, and give a very general analytical
expression that relates the scaling of the exponential rate with monomer
concentration to the microscopic details of the underlying reaction. This
approach can thus be used to infer the microscopic mechanism driving the
aggregation process from macroscopic measurements in complex systems.


  Long cell protrusions, which are effectively one-dimensional, are highly
dynamic subcellular structures. Length of many such protrusions keep
fluctuating about the mean value even in the the steady state. We develop here
a stochastic model motivated by length fluctuations of a type of appendage of
an eukaryotic cell called flagellum (also called cilium). Exploiting the
techniques developed for the calculation of level-crossing statistics of random
excursions of stochastic process, we have derived analytical expressions of
passage times for hitting various thresholds, sojourn times of random
excursions beyond the threshold and the extreme lengths attained during the
lifetime of these model flagella. We identify different parameter regimes of
this model flagellum that mimic those of the wildtype and mutants of a well
known flagellated cell. By analysing our model in these different parameter
regimes, we demonstrate how mutation can alter the level-crossing statistics
even when the steady state length remains unaffected by the same mutation.
Comparison of the theoretically predicted level crossing statistics, in
addition to mean and variance of the length, in the steady state with the
corresponding experimental data can be used in near future as stringent tests
for the validity of the models of flagellar length control. The experimental
data required for this purpose, though never reported till now, can be
collected, in principle, using a method developed very recently for flagellar
length fluctuations.


  The challenges presented by the COVID-19 epidemic have created a renewed
interest in the development of new methods to combat infectious diseases. A
prominent property of the SARS-CoV-2 transmission is the significant fraction
of asymptomatic transmission. This may influence the effectiveness of the
standard contact tracing procedure for quarantining potentially infected
individuals. However, the effects of asymptomatic transmission on the epidemic
threshold of epidemic spreading on networks are largely unknown. Here we study
the critical percolation transition in a simple epidemic network model in the
presence of a recursive contact tracing algorithm for instant quarantining. We
find that, above a certain fraction of asymptomatic transmission, standard
contact tracing loses its ability to suppress spreading below the epidemic
threshold. However, we also find that recursive contact tracing opens a
possibility to contain epidemics with a large fraction of asymptomatic or
presymptomatic transmission. In particular, we calculate the required fraction
of network nodes participating in the contact tracing for networks with
arbitrary degree distributions and for varying recursion depths and discuss the
influence of recursion depth and asymptomatic rate on the epidemic percolation
phase transition. We test and illustrate our theoretical results using
numerical simulations on infection trees and networks. We anticipate recursive
contact tracing to provide a basis for digital, app-based contact tracing tools
that extend the efficiency of contact tracing to diseases with a large fraction
of asymptomatic transmission.


  We report a confocal interferometric mid-infrared photothermal (MIP)
microscope and its application to label-free detection of biological
nanoparticles down to single virus level. We apply the interferometric
scattering principle to detect the weak photothermal effect induced by infrared
absorption in a viral particle. We validated this method with two kinds of
viruses, namely the vesicular stomatitis virus (VSV) and poxvirus. The single
virus spectra generated by MIP microscopy show high consistency in the same
group, with dominant peaks contributed by the amide I and amide II vibrations.
The ratio of these two peaks are significantly different between VSV and
poxvirus, highlighting the potential of using MIP microscopy for label-free
differentiation of viral particles. This all-optical chemical imaging method
opens a new way for spectroscopic detection of biological nanoparticles in a
label-free manner and may facilitate in predicting and controlling outbreaks of
emerging virus strains.


  The explosive growth of scientists, scientific journals, articles and
findings in recent years exponentially increases the difficulty scientists face
in navigating prior knowledge. This challenge is exacerbated by uncertainty
about the reproducibility of published findings. The availability of massive
digital archives, machine reading and extraction tools on the one hand, and
automated high-throughput experiments on the other, allow us to evaluate these
challenges at scale and identify novel opportunities for accelerating
scientific advance. Here we demonstrate a Bayesian calculus that enables the
positive prediction of robust, replicable scientific claims with findings
automatically extracted from published literature on gene interactions. We
matched these findings, filtered by science, with unfiltered gene interactions
measured by the massive LINCS L1000 high-throughput experiment to identify and
counteract sources of bias. Our calculus is built on easily extracted
publication meta-data regarding the position of a scientific claim within the
web of prior knowledge, and its breadth of support across institutions, authors
and communities, revealing that scientifically focused but socially and
institutionally independent research activity is most likely to replicate.
These findings recommend policies that go against the common practice of
channeling biomedical research funding into centralized research consortia and
institutes rather than dispersing it more broadly. Our results demonstrate that
robust scientific findings hinge upon a delicate balance of shared focus and
independence, and that this complex pattern can be computationally exploited to
decode bias and predict the replicability of published findings. These insights
provide guidance for scientists navigating the research literature and for
science funders seeking to improve it.


  Network embedding approaches are gaining momentum to analyse a large variety
of networks. Indeed, these approaches have demonstrated their efficiency for
tasks such as community detection, node classification, and link prediction.
However, very few network embedding methods have been specifically designed to
handle multiplex networks, i.e. networks composed of different layers sharing
the same set of nodes but having different types of edges. Moreover, to our
knowledge, existing approaches cannot embed multiple nodes from
multiplex-heterogeneous networks, i.e. networks composed of several layers
containing both different types of nodes and edges. In this study, we propose
MultiVERSE, an extension of the VERSE method with Random Walks with Restart on
Multiplex (RWR-M) and Multiplex-Heterogeneous (RWR-MH) networks. MultiVERSE is
a fast and scalable method to learn node embeddings from multiplex and
multiplex-heterogeneous networks. We evaluate MultiVERSE on several biological
and social networks and demonstrate its efficiency. MultiVERSE indeed
outperforms most of the other methods in the tasks of link prediction and
network reconstruction for multiplex network embedding, and is also efficient
in the task of link prediction for multiplex-heterogeneous network embedding.
Finally, we apply MultiVERSE to study rare disease-gene associations using link
prediction and clustering. MultiVERSE is freely available on github at
https://github.com/Lpiol/MultiVERSE.


  We study the post-translational escape of nascent proteins at the ribosomal
exit tunnel with the consideration of a real shape atomistic tunnel based on
the Protein Data Bank (PDB) structure of the large ribosome subunit of archeon
Haloarcula marismortui. Molecular dynamics simulations employing the Go-like
model for the proteins show that at intermediate and high temperatures,
including a presumable physiological temperature, the protein escape process at
the atomistic tunnel is quantitatively similar to that at a cylinder tunnel of
length L = 72 {\AA} and diameter d = 16 {\AA}. At low temperatures, the
atomistic tunnel, however, yields an increased probability of protein trapping
inside the tunnel while the cylinder tunnel does not cause the trapping.
All-$\beta$ proteins tend to escape faster than all-$\alpha$ proteins but this
difference is blurred on increasing the protein's chain length. A 29-residue
zinc-finger domain is shown to be severely trapped inside the tunnel. Most of
the single-domain proteins considered, however, can escape efficiently at the
physiological temperature with the escape time distribution following the
diffusion model proposed in our previous works. An extrapolation of the
simulation data to a realistic value of the friction coefficient for amino
acids indicates that the escape times of globular proteins are at the
sub-millisecond scale. It is argued that this time scale is short enough for
the smooth functioning of the ribosome by not allowing nascent proteins to jam
the ribosome tunnel.


  In this paper we explore a time-depended SEIR model, in which the dynamics of
the infection in four groups from a selected target group (population), divided
according to the infection, are modeled by a system of nonlinear ordinary
differential equations. Several basic parameters are involved in the model:
coefficients of infection rate, incubation rate, recovery rate. The
coefficients are adaptable to each specific infection, for each individual
country, and depend on the measures to limit the spread of the infection and
the effectiveness of the methods of treatment of the infected people in the
respective country. If such coefficients are known, solving the nonlinear
system is possible to be able to make some hypotheses for the development of
the epidemic. This is the reason for using Bulgarian COVID-19 data to first of
all, solve the so-called "inverse problem" and to find the parameters of the
current situation. Reverse logic is initially used to determine the parameters
of the model as a function of time, followed by computer solution of the
problem. Namely, this means predicting the future behavior of these parameters,
and finding (and as a consequence applying mass-scale measures, e.g.,
distancing, disinfection, limitation of public events), a suitable scenario for
the change in the proportion of the numbers of the four studied groups in the
future. In fact, based on these results we model the COVID-19 transmission
dynamics in Bulgaria and make a two-week forecast for the numbers of new cases
per day, active cases and recovered individuals. Such model, as we show, has
been successful for prediction analysis in the Bulgarian situation. We also
provide multiple examples of numerical experiments with visualization of the
results.


  Bacteria light harvesters are interpreted as Vicsek fractal, based upon their
morphology, even though such fractals contain only one, or perhaps two,
generations. After using fractal dimensions to describe the geometry, we
progress to use connectivity matrices to make an improved description of the
light harvesters and make connections with our previous studies.


  The detection and management of diseases become quite complicated when
pathogens contain asymptomatic phenotypes amongst their ranks, as evident
during the recent COVID-19 pandemic. Spreading of diseases has been studied
extensively under the paradigm of Susceptible - Infected - Recovered - Deceased
(SIRD) dynamics. Various game-theoretic approaches have also addressed disease
spread, many of which consider S, I, R, and D as strategies rather than as
states. Remarkably, most studies from the above approaches do not account for
the distinction between the symptomatic or asymptomatic aspect of the disease.
It is well-known that precautionary measures like washing hands, wearing masks
and social distancing significantly mitigate the spread of many contagious
diseases. Herein, we consider the adoption of such precautions as strategies
and treat S, I, R, and D as states. We also attempt to capture the differences
in epidemic spreading arising from symptomatic and asymptomatic diseases on
various network topologies. Through extensive computer simulations, we examine
that the cost of maintaining precautionary measures as well as the extent of
mass testing in a population affects the final fraction of socially responsible
individuals. We observe that the lack of mass testing could potentially lead to
a pandemic in case of asymptomatic diseases. Network topology also seems to
play an important role. We further observe that the final fraction of proactive
individuals depends on the initial fraction of both infected as well as
proactive individuals. Additionally, edge density can significantly influence
the overall outcome. Our findings are in broad agreement with the lessons
learnt from the ongoing COVID-19 pandemic.


  Patients infected with SARS-CoV-2 show a wide spectrum of clinical
manifestations ranging from mild febrile illness and cough up to acute
respiratory distress syndrome, multiple organ failure and death. Data from
patients with severe clinical manifestations compared to patients with mild
symptoms indicate that highly dysregulated exuberant inflammatory responses
correlate with severity of disease and lethality. Significantly elevated
cytokine levels, i.e. cytokine storm, seem to play a central role in severity
and lethality in COVID-19. We have previously shown that excessive cytokine
release induced by highly pathogenic avian H5N1 influenza A virus was reduced
by application of proteasome inhibitors. In the present study we present
experimental data of a central cellular pro-inflammatory signal pathways,
NF-kappaB, in the context of published clinical data from COVID-19 patients and
develop a hypothesis for a therapeutic approach aiming at the simultaneous
inhibition of whole cascades of pro-inflammatory cytokines and chemokines via
blocking the nuclear translocation of NF-kappaB by proteasome inhibitors. The
simultaneous inhibition of multiple cytokines/chemokines using clinically
approved proteasome inhibitors is expected to have a higher therapeutic
potential compared to single target approaches to prevent cascade (i.e.
triggering, synergistic, and redundant) effects of multiple induced cytokines
and may provide an additional therapeutic option to be explored for treatment
of critical stage COVID-19 patients.


  To combat the coronavirus disease 2019 (COVID-19) pandemic, the world has
vaccination, plasma therapy, herd immunity, and epidemiological interventions
as few possible options. The COVID-19 vaccine development is underway and it
may take a significant amount of time to develop the vaccine and after
development, it will take time to vaccinate the entire population, and plasma
therapy has some limitations. Herd immunity can be a plausible option to fight
COVID-19 for small countries. But for a country with huge population like
India, herd immunity is not a plausible option, because to acquire herd
immunity approximately 67% of the population has to be recovered from COVID-19
infection, which will put an extra burden on medical system of the country and
will result in a huge loss of human life. Thus epidemiological interventions
(complete lockdown, partial lockdown, quarantine, isolation, social distancing,
etc.) are some suitable strategies in India to slow down the COVID-19 spread
until the vaccine development. In this work, we have suggested the SIR model
with intervention, which incorporates the epidemiological interventions in the
classical SIR model. To model the effect of the interventions, we have
introduced \r{ho} as the intervention parameter. \r{ho} is a cumulative
quantity which covers all type of intervention. We have also discussed the
supervised machine learning approach to estimate the transmission rate
(\b{eta}) for the SIR model with intervention from the prevalence of COVID-19
data in India and some states of India. To validate our model, we present a
comparison between the actual and model-predicted number of COVID-19 cases.
Using our model, we also present predicted numbers of active and recovered
COVID-19 cases till Sept 30, 2020, for entire India and some states of India
and also estimate the 95% and 99% confidence interval for the predicted cases.


  We present a new mathematical model to explicitly capture the effects that
the three restriction measures: the lockdown date and duration, social
distancing and masks, and, schools and border closing, have in controlling the
spread of COVID-19 infections $i(r, t)$. Before restrictions were introduced,
the random spread of infections as described by the SEIR model grew
exponentially. The addition of control measures introduces a mixing of order
and disorder in the system's evolution which fall under a different
mathematical class of models that can eventually lead to critical phenomena. A
generic analytical solution is hard to obtain. We use machine learning to solve
the new equations for $i(r,t)$, the infections $i$ in any region $r$ at time
$t$ and derive predictions for the spread of infections over time as a function
of the strength of the specific measure taken and their duration. The machine
is trained in all of the COVID-19 published data for each region, county,
state, and country in the world. It utilizes optimization to learn the best-fit
values of the model's parameters from past data in each region in the world,
and it updates the predicted infections curves for any future restrictions that
may be added or relaxed anywhere. We hope this interdisciplinary effort, a new
mathematical model that predicts the impact of each measure in slowing down
infection spread combined with the solving power of machine learning, is a
useful tool in the fight against the current pandemic and potentially future
ones.


  This study aims to determine an optimal control strategy for vaccine
scheduling in COVID-19 pandemic treatment by converting widely acknowledged
infectious disease model named SEIR into an optimal control problem. The
problem is augmented by adding medication and vaccine limitations to match
real-world situations. Two version of the problem is formulated to minimize the
number of infected individuals at the same provide the optimal vaccine possible
to reduce the susceptible population to a considerably lower state. Optimal
control problems are solved using RBF-Galerkin method. These problems are
tested with a benchmarking dataset to determine required parameters. After this
step, problems are tested with recent data for New York State, USA. The results
regarding the proposed optimal control problem provides a set of evidences from
which an optimal strategy for vaccine scheduling can be chosen, when the
vaccine for COVID-19 will be available.


  Here we investigate translational regulation in bacteria by analyzing the
distribution of start codons in fully assembled genomes. We report 36 genes
(infC, rpoC, rnpA, etc.) showing a preference for non-AUG start codons in
evolutionarily diverse phyla ("non-AUG genes"). Most of the non-AUG genes are
functionally associated with translation, transcription or replication. In E.
coli, the percentage of essential genes among these 36 is significantly higher
than among all genes. Furthermore, the functional distribution of these genes
suggests that non-AUG start codons may be used to reduce gene expression during
starvation conditions, possibly through translational autoregulation or
IF3-mediated regulation.


  An accurate estimation of the Protein Space size, in light of the factors
that govern it, is a long-standing problem and of paramount importance in
evolutionary biology, since it determines the nature of protein evolvability. A
simple analysis will enable us to, firstly, reduce an unrealistic Protein Space
size of ~10^130 sequences, for a 100-residues polypeptide chain, to ~10^9
functional proteins and, secondly, estimate a robust average-mutation rate per
amino acid (x ~1.23) and infer from it, in light of the protein marginal
stability, that only a fraction of the sequence will be available at any one
time for a functional protein to evolve. Although this result does not solve
the Protein Space vastness problem, frames it in a more rational one.


  Eukarya and Bacteria are the most evolutionarily distant domains of life,
which is reflected by differences in their cellular structure and physiology.
For example, Eukarya feature membrane-bound organelles such as nuclei and
mitochondria, whereas Bacteria have none. The greater complexity of Eukarya
renders them difficult to study from both an experimental and theoretical
perspective. However, encouraged by a recent experimental result showing that
budding yeast (a unicellular eukaryote) obeys the same proportionality between
ribosomal proteome fractions and cellular growth rates as Bacteria, we derive a
set of relations describing eukaryotic growth from first principles of ribosome
biogenesis. We recover the observed ribosomal protein proportionality, and then
continue to obtain two growth-laws for the number of RNA polymerases
synthesizing ribosomal RNA per ribosome in the cell. These growth-laws, in
turn, reveal two invariants of eukaryotic growth, i.e. quantities predicted to
be conserved by Eukarya regardless of growth conditions. The invariants, which
are the first of their kind for Eukarya, clarify the coordination of
transcription and translation kinetics as required by ribosome biogenesis, and
link these kinetic parameters to cellular physiology. We demonstrate
application of the relations to the yeast S. cerevisiae and find the
predictions to be in good agreement with currently available data. We then
outline methods to quantitatively deduce several unknown kinetic and
physiological parameters. The analysis is not specific to S. cerevisiae and can
be extended to other lower (unicellular) Eukarya when data become available.
The relations may also have relevance to certain cancer cells which, like
bacteria and yeast, exhibit rapid cell proliferation and ribosome biogenesis.


  The size and quality of chemical libraries to the drug discovery pipeline are
crucial for developing new drugs or repurposing existing drugs. Existing
techniques such as combinatorial organic synthesis and High-Throughput
Screening usually make the process extraordinarily tough and complicated since
the search space of synthetically feasible drugs is exorbitantly huge. While
reinforcement learning has been mostly exploited in the literature for
generating novel compounds, the requirement of designing a reward function that
succinctly represents the learning objective could prove daunting in certain
complex domains. Generative Adversarial Network-based methods also mostly
discard the discriminator after training and could be hard to train. In this
study, we propose a framework for training a compound generator and learning a
transferable reward function based on the entropy maximization inverse
reinforcement learning paradigm. We show from our experiments that the inverse
reinforcement learning route offers a rational alternative for generating
chemical compounds in domains where reward function engineering may be less
appealing or impossible while data exhibiting the desired objective is readily
available.


  A recent experimental study found that the binding affinity between the
cellular receptor human angiotensin converting enzyme 2 (ACE2) and
receptor-binding domain (RBD) in spike (S) protein of novel severe acute
respiratory syndrome coronavirus 2 (SARS-CoV-2) is more than 10-fold higher
than that of the original severe acute respiratory syndrome coronavirus
(SARS-CoV). However, main-chain structures of the SARS-CoV-2 RBD are almost the
same with that of the SARS-CoV RBD. Understanding physical mechanism
responsible for the outstanding affinity between the SARS-CoV-2 S and ACE2 is
the "urgent challenge" for developing blockers, vaccines and therapeutic
antibodies against the coronavirus disease 2019 (COVID-19) pandemic.
Considering the mechanisms of hydrophobic interaction, hydration shell, surface
tension, and the shielding effect of water molecules, this study reveals a
hydrophobic-interaction-based mechanism by means of which SARS-CoV-2 S and ACE2
bind together in an aqueous environment. The hydrophobic interaction between
the SARS-CoV-2 S and ACE2 protein is found to be significantly greater than
that between SARS-CoV S and ACE2. At the docking site, the hydrophobic portions
of the hydrophilic side chains of SARS-CoV-2 S are found to be involved in the
hydrophobic interaction between SARS-CoV-2 S and ACE2. We propose a method to
design live attenuated viruses by mutating several key amino acid residues of
the spike protein to decrease the hydrophobic surface areas at the docking
site. Mutation of a small amount of residues can greatly reduce the hydrophobic
binding of the coronavirus to the receptor, which may be significant reduce
infectivity and transmissibility of the virus.


  In the cellular phenomena of cytoplasmic streaming, molecular motors carrying
cargo along a network of microtubules entrain the surrounding fluid. The
piconewton forces produced by individual motors are sufficient to deform long
microtubules, as are the collective fluid flows generated by many moving
motors. Studies of streaming during oocyte development in the fruit fly
$D.~melanogaster$ have shown a transition from a spatially-disordered
cytoskeleton, supporting flows with only short-ranged correlations, to an
ordered state with a cell-spanning vortical flow. To test the hypothesis that
this transition is driven by fluid-structure interactions we study a
discrete-filament model and a coarse-grained continuum theory for motors moving
on a deformable cytoskeleton, both of which are shown to exhibit a
$swirling~instability$ to spontaneous large-scale rotational motion, as
observed.


  Prediction of protein-ligand complexes for flexible proteins remains still a
challenging problem in computational structural biology and drug design. Here
we present two novel deep neural network approaches with significant
improvement in efficiency and accuracy of binding mode prediction on a large
and diverse set of protein systems compared to standard docking. Whereas the
first graph convolutional network is used for re-ranking poses the second
approach aims to generate and rank poses independent of standard docking
approaches. This novel approach relies on the prediction of distance matrices
between ligand atoms and protein C_alpha atoms thus incorporating side-chain
flexibility implicitly.


  Improved monitoring and associated inferential tools to efficiently identify
declining bird populations, particularly of rare or sparsely distributed
species, is key to informed conservation and management across large
spatio-temporal regions. We assess abundance trends for 106 bird species in a
network of eight national park forests located within the northeast USA from
2006-2019 using a novel hierarchical model. We develop a multi-species,
multi-region removal sampling model that shares information across species and
parks to enable inference on rare species and sparsely sampled parks and to
evaluate the effects of local forest structure. Trends in bird abundance over
time varied widely across parks, but species showed similar trends within
parks. Three parks (Acadia, Marsh-Billings-Rockefeller, and Morristown)
decreased in bird abundance across all species, while three parks (Saratoga,
Roosevelt-Vanderbilt, and Weir-Farm) increased in abundance. Bird abundance
peaked at medium levels of basal area and high levels of percent forest and
forest regeneration, with percent forest having the largest effect. Variation
in these effects across parks could be a result of differences in forest
structural stage and diversity. Our novel hierarchical model enables estimates
of abundance at the network, park, guild, and species levels. We found large
variation in abundance trends across parks but not across bird guilds,
suggesting that local forest condition may have a broad and consistent effect
on the entire bird community within a given park. Management should target the
three parks with overall decreasing trends in bird abundance to further
identify what specific factors are driving observed declines across the bird
community. Understanding how bird communities respond to local forest structure
and other stressors is crucial for informed and lasting management.


  Predicting the properties of a molecule from its structure is a challenging
task. Recently, deep learning methods have improved the state of the art for
this task because of their ability to learn useful features from the given
data. By treating molecule structure as graphs, where atoms and bonds are
modeled as nodes and edges, graph neural networks (GNNs) have been widely used
to predict molecular properties. However, the design and development of GNNs
for a given data set rely on labor-intensive design and tuning of the network
architectures. Neural architecture search (NAS) is a promising approach to
discover high-performing neural network architectures automatically. To that
end, we develop an NAS approach to automate the design and development of GNNs
for molecular property prediction. Specifically, we focus on automated
development of message-passing neural networks (MPNNs) to predict the molecular
properties of small molecules in quantum mechanics and physical chemistry data
sets from the MoleculeNet benchmark. We demonstrate the superiority of the
automatically discovered MPNNs by comparing them with manually designed GNNs
from the MoleculeNet benchmark. We study the relative importance of the choices
in the MPNN search space, demonstrating that customizing the architecture is
critical to enhancing performance in molecular property prediction and that the
proposed approach can perform customization automatically with minimal manual
effort.


  Recently, deep neural network (DNN)-based drug-target interaction (DTI)
models were highlighted for their high accuracy with affordable computational
costs. Yet, the models' insufficient generalization remains a challenging
problem in the practice of in-silico drug discovery. We propose two key
strategies to enhance generalization in the DTI model. The first is to predict
the atom-atom pairwise interactions via physics-informed equations
parameterized with neural networks and provides the total binding affinity of a
protein-ligand complex as their sum. We further improved the model
generalization by augmenting a broader range of binding poses and ligands to
training data. We validated our model, PIGNet, in the comparative assessment of
scoring functions (CASF) 2016, demonstrating the outperforming docking and
screening powers than previous methods. Our physics-informing strategy also
enables the interpretation of predicted affinities by visualizing the
contribution of ligand substructures, providing insights for further ligand
optimization.


  Accurate modeling of lockdown effects on SARS-CoV-2 epidemic evolution is a
key issue in order e.g. to inform health-care decisions on emergency
management. The compartmental and spatial models so far proposed use parametric
descriptions of the contact rate, often assuming a time-invariant effect of the
lockdown. In this paper we show that these assumptions may lead to erroneous
evaluations on the ongoing pandemic. Thus, we develop a new class of
nonparametric compartmental models able to describe how the impact of the
lockdown varies in time. Exploiting regularization theory, hospitalized data
are mapped into an infinite-dimensional space, hence obtaining a function which
takes into account also how social distancing measures and people's growing
awareness of infection's risk evolves as time progresses. This permits to
reconstruct a continuous-time profile of SARS-CoV-2 reproduction number with a
resolution never reached before. When applied to data collected in Lombardy,
the most affected Italian region, our model illustrates how people behaviour
changed during the restrictions and its importance to contain the epidemic.
Results also indicate that, at the end of the lockdown, around 12% of people in
Lombardy and 5% in Italy was affected by SARS-CoV-2. Then, we discuss how the
situation evolved after the end of the lockdown showing that the reproduction
number is dangerously increasing in the last weeks due to holiday relax
especially in the younger population and increased migrants arrival, reaching
values larger than one on August 1, 2020. Since several countries still observe
a growing epidemic, including Italy, and all could be subject to a second wave
after the summer, the proposed reproduction number tracking methodology can be
of great help to health care authorities to prevent another SARS-CoV-2
diffusion or to assess the impact of lockdown restrictions to contain the
spread.


  Genomics, especially multi-omics, has made precision medicine feasible. The
completion and publicly accessible multi-omics resource with clinical outcome,
such as The Cancer Genome Atlas (TCGA) is a great test bed for developing
computational methods that integrate multi-omics data to predict patient cancer
phenotypes. We have been utilizing TCGA multi-omics data to predict cancer
patient survival, using a variety of approaches, including prior-biological
knowledge (such as pathways), and more recently, deep-learning methods. Over
time, we have developed methods such as Cox-nnet, DeepProg, and two-stage
Cox-nnet, to address the challenges due to multi-omics and multi-modality.
Despite the limited sample size (hundreds to thousands) in the training
datasets as well as the heterogeneity nature of human populations, these
methods have shown significance and robustness at predicting patient survival
in independent population cohorts. In the following, we would describe in
detail these methodologies, the modeling results, and important biological
insights revealed by these methods.


  Modeling the effects of mutations on the binding affinity plays a crucial
role in protein engineering and drug design. In this study, we develop a novel
deep learning based framework, named GraphPPI, to predict the binding affinity
changes upon mutations based on the features provided by a graph neural network
(GNN). In particular, GraphPPI first employs a well-designed pre-training
scheme to enforce the GNN to capture the features that are predictive of the
effects of mutations on binding affinity in an unsupervised manner and then
integrates these graphical features with gradient-boosting trees to perform the
prediction. Experiments showed that, without any annotated signals, GraphPPI
can capture meaningful patterns of the protein structures. Also, GraphPPI
achieved new state-of-the-art performance in predicting the binding affinity
changes upon both single- and multi-point mutations on five benchmark datasets.
In-depth analyses also showed GraphPPI can accurately estimate the effects of
mutations on the binding affinity between SARS-CoV-2 and its neutralizing
antibodies. These results have established GraphPPI as a powerful and useful
computational tool in the studies of protein design.


  COVID-19 outbreaks have proven to be very difficult to isolate and extinguish
before they spread out. An important reason behind this might be that
epidemiological barriers consisting in stopping symptomatic people are likely
to fail because of the contagion time before onset, mild cases and/or
asymptomatics carriers. Motivated by these special COVID-19 features, we study
a scheme for containing an outbreak in a city that consists in adding an extra
firewall block between the outbreak and the rest of the city. We implement a
coupled compartment model with stochastic noise to simulate a localized
outbreak that is partially isolated and analyze its evolution with and without
firewall for different plausible model parameters. We explore how further
improvements could be achieved if the epidemic evolution would trigger policy
changes for the flux and/or lock-down in the different blocks. Our results show
that a substantial improvement is obtained by merely adding an extra block
between the outbreak and the bulk of the city.


  Record Dynamics (RD) deals with complex systems evolving through a sequence
of metastable stages. These are macroscopically distinguishable and appear
stationary, except for the sudden and rapid changes, called quakes, which
induce the transitions from one stage to the next. This phenomenology is well
known in physics as "physical aging", but from the vantage point of RD the
evolution of a class of systems of physical, biological and cultural origin is
rooted in a hierarchically structured configuration space and can therefore be
analyzed by similar statistical tools. This colloquium paper strives to present
in a coherent fashion methods and ideas that have gradually evolved over time.
To this end, it first describes the differences and similarities between RD and
two widespread paradigms of complex dynamics, Self Organized Criticality and
Continuous Time Random Walks. It then outlines the Poissonian nature of records
events in white noise time series, and connects it to the statistics of quakes
in metastable hierarchical systems, arguing that the relaxation effects of
quakes can generally be described by power laws unrelated to criticality.
Several different applications of RD have been developed over the years. Some
of these are described, showinghe basic RD hypothesis, the log time homogeneity
of quake dynamics, can be empirically verified in a given context. The
discussion summarizes the paper and briefly mentions applications not discussed
in detail. Finally, the outlook points to possible improvements and to new
areas of research where RG could be of use.


  Motivated by analogies between the spreading of human-to-human infections and
of chemical processes, we develop a comprehensive model that accounts both for
infection and for transport. In this analogy, the three different populations
of infection models correspond to three chemical species. Areal densities
emerge as the key variables, thus capturing the effect of spatial density. We
derive expressions for the kinetics of the infection rates and for the
important parameter R0, that include areal density and its spatial
distribution. Coupled with mobility the model allows the study of various
effects. We first present results for a batch reactor, the chemical process
equivalent of the SIR model. Because density makes R0 a decreasing function of
the process extent, the infection curves are different and smaller than for the
standard SIR model. We show that the effect of the initial conditions is
limited to the onset of the epidemic. We derive effective infection curves for
a number of cases, including a back-and-forth commute between regions of low
and high R0 environments. We then consider spatially distributed systems. We
show that diffusion leads to traveling waves, which in 1-D geometries propagate
at a constant speed and with a constant shape, both of which are sole functions
of R0. The infection curves are slightly different than for the batch problem,
as diffusion mitigates the infection intensity, thus leading to an effective
lower R0. The dimensional wave speed is found to be proportional to the product
of the square root of the diffusivity and of an increasing function of R0,
confirming the importance of restricting mobility in arresting the propagation
of infection. We examine the interaction of infection waves under various
conditions and scenarios, and extend the wave propagation analysis to 2-D
heterogeneous systems.


  Spike (S) glycoproteins mediate the coronavirus entry into the host cell. The
S1 subunit of S-proteins contains the receptor-binding domain (RBD) that is
able to recognize different host receptors, highlighting its remarkable
capacity to adapt to their hosts along the viral evolution. While RBD in spike
proteins is determinant for the virus-receptor interaction, the active residues
lie at the receptor-binding motif (RBM), a region located in RBD that plays a
fundamental role binding the outer surface of their receptors. Here, we address
the hypothesis that SARS-CoV and SARS-CoV-2 strains able to use
angiotensin-converting enzyme 2 (ACE2) proteins have adapted their RBM along
the viral evolution to explore specific conformational topology driven by the
residues YGF to infect host cells. We also speculate that this YGF-based
mechanism can act as a protein signature located at the RBM to distinguish
coronaviruses able to use ACE2 as a cell entry receptor.


  Phylogenetic trees are important tools in the study of evolutionary
relationships between species. Measures such as the index of Sackin, Colless,
and Total Cophenetic have been extensively used to quantify tree balance, one
key property of phylogenies. Recently a new proposal has been introduced, based
on the spectrum of the Laplacian matrix associated with the tree. In this work,
we calculate the Laplacian matrix analytically for two extreme cases,
corresponding to fully balanced and fully unbalanced trees. For maximally
balanced trees no closed form for the Laplacian matrix was derived, but we
present an algorithm to construct it. We show that Laplacian matrices of fully
balanced trees display self-similar patterns that result in highly degenerated
eigenvalues. Degeneracy is the main signature of this topology, since it is
totally absent in fully unbalanced trees. We also establish some analytical and
numerical results about the largest eigenvalue of Laplacian matrices for these
topologies.


  The study of species organization and their clustering by genetic or
phenotypic similarity is carried out with the tools of phylogenetic trees. An
important structural property of phylogenetic trees is the balance, which
measures how taxa are distributed among clades. Tree balance can be measured
using indices such as the Sackin ($S$) and the Total Cophenetic ($\Phi$), which
are based on the distance between nodes of the tree and its root. Here, we
propose a new metric for tree balance, $\bar{d}$, the Area per Pair (APP) of
the tree, which is a re-scaled version of the so called tree area. We compute
$\bar{d}$ for the rooted caterpillar and maximally balanced trees and we also
obtain exact formulas for its expected value and variance under the Yule model.
The variance of APP for Yule trees has the remarkable property of converging to
an asymptotic constant value for large trees. We compare the Sackin, Total
Cophenetic and APP indices for hundreds of empirical phylogenies and show that
APP represents the observed distribution of tree balances better than the two
other metrics.


  Trust and trustworthiness form the basis for continued social and economic
interactions, and they are also fundamental for cooperation, fairness, honesty,
and indeed for many other forms of prosocial and moral behavior. However, trust
entails risks, and building a trustworthy reputation requires effort. So how
did trust and trustworthiness evolve, and under which conditions do they
thrive? To find answers, we operationalize trust and trustworthiness using the
trust game with the trustor's investment and the trustee's return of the
investment as the two key parameters. We study this game on different networks,
including the complete network, random and scale-free networks, and in the
well-mixed limit. We show that in all but one case the network structure has
little effect on the evolution of trust and trustworthiness. Specifically, for
well-mixed populations, lattices, random and scale-free networks, we find that
trust never evolves, while trustworthiness evolves with some probability
depending on the game parameters and the updating dynamics. Only for the
scale-free network with degree non-normalized dynamics, we find parameter
values for which trust evolves but trustworthiness does not, as well as values
for which both trust and trustworthiness evolve. We conclude with a discussion
about mechanisms that could lead to the evolution of trust and outline
directions for future work.


  Coronavirus outbreak is one of the most challenging pandemics for the entire
human population of the planet Earth. Techniques such as the isolation of
infected persons and maintaining social distancing are the only preventive
measures against the epidemic COVID-19. The actual estimation of the number of
infected persons with limited data is an indeterminate problem faced by data
scientists. There are a large number of techniques in the existing literature,
including reproduction number, the case fatality rate, etc., for predicting the
duration of an epidemic and infectious population. This paper presents a case
study of different techniques for analysing, modeling, and representation of
data associated with an epidemic such as COVID-19. We further propose an
algorithm for estimating infection transmission states in a particular area.
This work also presents an algorithm for estimating end-time of an epidemic
from Susceptible Infectious and Recovered model. Finally, this paper presents
empirical and data analysis to study the impact of transmission probability,
rate of contact, infectious, and susceptible on the epidemic spread.


  Disgust is a basic emotion that serves to avoid contaminants, and is central
to the behavioural immune system. While disgust-motivated avoidance occurs in
bonobos and chimpanzees, humans show uniquely high levels of contamination
sensitivity. Current theory postulates that human disgust is primarily a
genetic adaptation, side-lining social transmission and denying parental
modelling. Here, we test whether this strong view is warranted by simulating
100000 years of cultural and biological evolution in various scenarios. We
modelled disgust as a trait that governed the extent to which individuals
forewent potentially contaminated nutrition. This lowered their risk of
gastrointestinal illness, at the cost of increasing starvation risk and
birth-interval. As environmental levels and costs of contamination rose, so did
relative fitness of high-disgust individuals. Evolutionary shifts in disgust as
polygenic trait occurred as a consequence of germ-cell mutations, but were most
prominent in populations with high initial genetic variance. Crucially,
cultural transmission between generations operated at a higher rate, even if
parental modelling was eliminated. This study serves not only as evidence of
cultural evolution shaping the behavioural immune system, but also as an
illustration of emerging theories that paint affective and cognitive mechanisms
as socially transmitted rather than biologically determined functions.


  COVID-19 presents a great threat to public health worldwide and the
infectious agent SARS-CoV-2 is currently the target of much research aiming at
inhibition. The virus' main protease is a dimeric enzyme that has only recently
begun to be thoroughly described, opening the door for virtual screening more
broadly. Here, a PAIN-filtered flavonoid database was screened against four
sites of the protease: a free (normal) conformation of the Substrate Binding
Site (NSBS), an induced-fit state of the SBS (ISBS), a Dimerization Site (DS)
and a Cryptic Site (CS). The mean binding energies of the top five ligands from
each site were -9.52, -11.512, -7.042 and -10.348 kcal/mol for the NSBS, the
ISBS, the DS and the CS, respectively. For the DS and CS, these top five
compounds were selected as candidates to bind their respective site. In the
case of SBS, the top 30 ligands with the lowest binding energies from NSBS and
ISBS were contrasted and the ones present in both lists were selected as the
final candidates. The final list was: Dorsilurin E (FL3FQUNP0001), Euchrenone
a11 (FL2FALNP0014), Kurziflavolactone C (FL2FA9NC0016), Licorice glycoside E
(FL2F1AGSN001) and Taxifolin 3'- (6"-phenyl- acetylglucoside) (FL4DACGS0020)
for the SBS; Sanggenol O (FL2FALNP0020), CHEMBL2171573, Kanzonol E
(FL3F1ANP0001), CHEMBL2171584 and Abyssynoflavanone VI (FL2FACNP0014) for DS
and CHEMBL2171598, CHEMBL2171577, Denticulaflavanol (FL5FAANR0001),
Kurzichalcolactone (FL1CA9NC0001) and CHEMBL2171578 for CS. Virtual screening
integrated several confirmation methods, including cross-docking assays and
positive and negative controls. All 15 compounds are currently subjected to
molecular dynamics so as to theoretically validate their binding to the
protease.


  The SARS-CoV-2 pandemic peaked very early in comparison to the thresholds
predicted by an analysis of prior lockdown regimes. The most convenient
explanation is that some, external factor changed the value of the basic
reproduction number, $r_{\rm 0}$; and there certainly are arguments for this.
Other factors could, nonetheless, have played a role. This research attempts to
reconcile the observed peak with the thresholds predicted by lockdown regimes
similar to the one in force at the time. It contemplates the effect of two,
different, hypothetical errors in the data: The first is that the true level of
infection has been underestimated by a multiplicative factor, while the second
is that of an imperceptible, pre-existing, immune fraction of the population.
While it is shown that it certainly is possible to manufacture the perception
of an early peak as extreme as the one observed, solely by way of these two
phenomena, the values need to be fairly high. The phenomena would not, by any
measure, be insignificant. It also remains an inescapable fact that the early
peak in infections coincided with a fairly profound change in $r_{\rm 0}$; in
all the contemplated scenarios of data-deficiency.


  Imported COVID-19 cases, if unchecked, can jeopardize the effort of domestic
containment. We aim to find out what sustainable border control options for
different entities (e.g., countries, states) exist during the reopening phases,
given their own choice of domestic control measures and new technologies such
as contact tracing. We propose a SUIHR model, which represents an extension to
the discrete time SIR models. The model focuses on studying the spreading of
virus predominantly by asymptomatic and pre-symptomatic patients. Imported risk
and (1-tier) contact tracing are both built into the model. Under plausible
parameter assumptions, we seek sustainable border control policies, in
combination with sufficient internal measures, which allow entities to confine
the virus without the need to revert back to more restrictive life styles or to
rely on herd immunity. When the base reproduction number of COVID-19 exceeds
2.5, even 100% effective contact tracing alone is not enough to contain the
spreading. For an entity that has completely eliminated the virus domestically,
and resumes "normal", very strict pre-departure screening and test and
isolation upon arrival combined with effective contact tracing can only delay
another outbreak by 6 months. However, if the total net imported cases are
non-increasing, and the entity employs a confining domestic control policy,
then the total new cases can be contained even without border control.


  Inferring the processes underlying the emergence of observed patterns is a
key challenge in theoretical ecology. Much effort has been made in the past
decades to collect extensive and detailed information about the spatial
distribution of tropical rainforests, as demonstrated, e.g., in the 50 ha
tropical forest plot on Barro Colorado Island, Panama. These kind of plots have
been crucial to shed light on diverse qualitative features, emerging both at
the single-species or the community level, like the spatial aggregation or
clustering at short scales. Here, we build on the progress made in the study of
the density correlation functions applied to biological systems, focusing on
the importance of accurately defining the borders of the set of trees, and
removing the induced biases. We also pinpoint the importance of combining the
study of correlations with the scale dependence of fluctuations in density,
which are linked to the well known empirical Taylor's power law. Density
correlations and fluctuations, in conjunction, provide an unique opportunity to
interpret the behaviors and possibly to allow comparisons between data and
models. We also study such quantities in models of spatial patterns and, in
particular, we find that a spatially explicit neutral model generates patterns
with many qualitative features in common with the empirical ones.


  This paper studies the mathematical properties of collectively canalizing
Boolean functions, a class of functions that has arisen from applications in
systems biology. Boolean networks are an increasingly popular modeling
framework for regulatory networks, and the class of functions studied here
captures a key feature of biological network dynamics, namely that a subset of
one or more variables, under certain conditions, can dominate the value of a
Boolean function, to the exclusion of all others. These functions have rich
mathematical properties to be explored. The paper shows how the number and type
of such sets influence a function's behavior and define a new measure for the
canalizing strength of any Boolean function. We further connect the concept of
collective canalization with the well-studied concept of the average
sensitivity of a Boolean function. The relationship between Boolean functions
and the dynamics of the networks they form is important in a wide range of
applications beyond biology, such as computer science, and has been studied
with statistical and simulation-based methods. But the rich relationship
between structure and dynamics remains largely unexplored, and this paper is
intended as a contribution to its mathematical foundation.


  The ability of organisms to accurately sense their environment and respond
accordingly is critical for evolutionary success. However, exactly how the
sensory ability influences fitness is a topic of active research, while the
necessity of a time delay between when unreliable environmental cues are sensed
and when organisms can mount a response has yet to be explored at any length.
Accounting for this delay in phenotype response in models of population growth,
we find that a critical error probability can exist under certain environmental
conditions: an organism with a sensory system with any error probability less
than the critical value can achieve the same long-term growth rate as an
organism with a perfect sensing system. We also observe a trade off between the
evolutionary value of sensory information and robustness to error, mediated by
the rate at which the phenotype distribution relaxes to steady-state. The
existence of the critical error probability could have several important
evolutionary consequences, primarily that sensory systems operating at the
non-zero critical error probability may be evolutionarily optimal.


  Cooperation in evolutionary biology means paying a cost, c, to enjoy
benefits, b. A defector is one who does not pay any cost but enjoys the
benefits of cooperators. Human societies, especially, have evolved a strategy
to discourage defection, punishment. Costly punishment is a type of punishment
where an agent in a biological network or some cooperative scheme pays a cost
to ensure another agent incurs some cost. A previous study shows how parameters
like diversity in neighbors across games, density of connectivity, and costly
punishment influence the evolution of cooperation in non-regular networks. In
this study, evolution in regular networks due to the influence of costly
punishment is also considered, specifically spatial lattices. This study
compares observations between non-regular and regular networks as these
parameters change and brings a clearer understanding of interactions that occur
in these networks. The models, results and analysis bring a new understanding
to game theory and punishment. The results show that costly punishment never
arises as a sole evolutionary strategy. However, in evolutionary games where
costly punishers could evolve more favorable strategies, the initial presence
of costly punishers would bring about high average payoffs in all types of
regular networks and heterogeneous networks. In regular networks, every node
has the same degree, k. Although punishment is conventionally thought to be
anti-progressive, in the presence of diversity in neighbors, it magnifies the
payoff of a group for all heterogeneous networks. In regular networks however,
diversity in neighbors is not required for costly punishment to boost average
payoff. This suggests an answer to the question on why costly punishment has
been favored by natural selection, which is particularly obvious in human
populations.


  We consider age-structured models with an imposed refractory period between
births. These models can be used to formulate alternative population control
strategies to China's one-child policy. By allowing any number of births, but
with an imposed delay between births, we show how the total population can be
decreased and how a relatively older age distribution can be generated. This
delay represents a more "continuous" form of population management for which
the strict one-child policy is a limiting case. Such a policy approach could be
more easily accepted by society. Our analyses provide an initial framework for
studying demographics and how social constraints influence population
structure.


  This is work in progress. We make it accessible hoping that people might find
the idea useful. We propose a discrete, recursive 5-compartment model for the
spread of epidemics, which we call {\em SEPIR-model}. Under mild assumptions
which typically are fulfilled for the Covid-19 pandemic it can be used to
reproduce the development of an epidemic from a small number of parameters
closely related to the data. We demonstrate this at the development in Germany
and Switzerland. It also allows model predictions assuming nearly constant
reproduction numbers. Thus it might be a useful tool for shedding light on
which interventions might be most effective in the future. In future work we
will discuss other aspects of the model and more countries.


  We explore the use of Array-RQMC, a randomized quasi-Monte Carlo method
designed for the simulation of Markov chains, to reduce the variance when
simulating stochastic biological or chemical reaction networks with
$\tau$-leaping. The task is to estimate the expectation of a function of
molecule copy numbers at a given future time $T$ by the sample average over $n$
sample paths, and the goal is to reduce the variance of this sample-average
estimator. We find that when the method is properly applied, variance
reductions by factors in the thousands can be obtained. These factors are much
larger than those observed previously by other authors who tried RQMC methods
for the same examples. Array-RQMC simulates an array of realizations of the
Markov chain and requires a sorting function to reorder these chains according
to their states, after each step. The choice of sorting function is a key
ingredient for the efficiency of the method, although in our experiments,
Array-RQMC was never worse than ordinary Monte Carlo, regardless of the sorting
method. The expected number of reactions of each type per step also has an
impact on the efficiency gain.


  Boolean network (BN) is a simple model widely used to study complex dynamic
behaviour of biological systems. Nonetheless, it might be difficult to gather
enough data to precisely capture the behavior of a biological system into a set
of Boolean functions. These issues can be dealt with to some extent using
parametrised Boolean networks (ParBNs), as it allows to leave some update
functions unspecified. In this paper, we attack the control problem for ParBNs
with asynchronous semantics. While there is an extensive work on controlling
BNs without parameters, the problem of control for ParBNs has not been in fact
addressed yet. The goal of control is to ensure the stabilisation of a system
in a given state using as few interventions as possible. There are many ways to
control BN dynamics. Here, we consider the one-step approach in which the
system is instantaneously perturbed out of its actual state. A naive approach
to handle control of ParBNs is using parameter scan and solve the control
problem for each parameter valuation separately using known techniques for
non-parametrised BNs. This approach is however highly inefficient as the
parameter space of ParBNs grows doubly-exponentially in the worst case. In this
paper, we propose a novel semi-symbolic algorithm for the one-step control
problem of ParBNs, that builds on a symbolic data structures to avoid scanning
individual parameters. We evaluate the performance of our approach on real
biological models.