  We provide a framework for experimentation at
https://mbuliga.github.io/quinegraphs/ic-vs-chem.html#icvschem with two
artificial chemistries: directed interaction combinators (dirIC, defined in
section 2) and chemlambda. We are interested if these chemistries allow for
artificial life behaviour: replication, metabolism and death.
  The main conclusion of these experiments is that graph rewrites systems which
allow conflicting rewrites are better than those which don't, as concerns their
artificial life properties. This is in contradiction with the search for good
graph rewrite systems for decentralized computing, where non-conflicting graph
rewrite systems are historically preferred.
  This continues the artificial chemistry experiments with chemlambda, lambda
calculus or interaction combinators, available from the entry page at
https://chemlambda.github.io/index.html and described in arXiv:2003.14332.


  Mathematical models of SARS-CoV-2 spread are used for guiding the design of
mitigation steps aimed at containing and decelerating the contagion, and at
identifying impending breaches of health care system surge capacity. The
challenges of having only lacunary information about daily new infections are
compounded by the geographic heterogeneity of the population. To address this
problem, we propose to account for the differences between rural and urban
settings using network-based, distributed models where the spread of the
pandemic is described in distinct local cohorts with nested SEIR models. The
setting of the model parameters takes into account the fact that SARS-CoV-2
transmission occurs mostly via human-to-human contact, and that the frequency
of contact among individuals differs between urban and rural areas, and may
change over time. Moreover, the probability that the virus spreads into an
uninfected community is associated with influx of individuals from other
communities where the infection is present. To account for these important
aspects, each node of the network is characterized by the frequency of contact
between its members and by its level of connectivity with other nodes. Census
and cell phone data can be used to set up the adjacency matrix of the network,
which can, in turn, be modified to account for different levels of mitigation
measures. In order to make the network SEIR model that we propose easy to
customize, it is formulated in terms of easily interpretable parameters that
can be estimated from available community level data. The models parameters are
estimated with Bayesian techniques using COVID-19 data for the states of Ohio
and Michigan. The network model also gives rise to a geographically distributed
computational model that explains the geographic dynamics of the contagion,
e.g., in larger cities surrounded by suburban and rural areas.


  One of the key indicators used in tracking the evolution of an infectious
disease isthe reproduction number. This quantity is usually computed using the
reportednumber of cases, but ignoring that many more individuals may be
infected (e.g.asymptomatics). We propose a statistical procedure to quantify
the impact of un-detected infectious cases on the determination of the
effective reproduction number. Our approach is stochastic, data-driven and not
relying on any compartmentalmodel. It is applied to the COVID-19 case in eight
different countries and all Italianregions, showing that the effect of
undetected cases leads to estimates of the effective reproduction numbers
larger than those obtained only with the reported cases by factors ranging from
two to ten. Our findings urge caution about deciding when and how to relax
containment measures based on the value of the reproduction number.


  The devastating impact of the COVID-19 pandemic caused by SARS coronavirus 2
(SARS CoV 2) has raised important questions about viral origin, mechanisms of
zoonotic transfer to humans, whether companion or commercial animals can act as
reservoirs for infection, and why there are large variations in SARS-CoV-2
susceptibilities across animal species. Powerful in silico modelling methods
can rapidly generate information on newly emerged pathogens to aid
countermeasure development and predict future behaviours. Here we report an in
silico structural homology modelling, protein-protein docking, and molecular
dynamics simulation study of the key infection initiating interaction between
the spike protein of SARS-Cov-2 and its target, angiotensin converting enzyme 2
(ACE2) from multiple species. Human ACE2 has the strongest binding interaction,
significantly greater than for any species proposed as source of the virus.
Binding to pangolin ACE2 was the second strongest, possibly due to the
SARS-CoV-2 spike receptor binding domain (RBD) being identical to pangolin CoV
spike RDB. Except for snake, pangolin and bat for which permissiveness has not
been tested, all those species in the upper half of the affinity range (human,
monkey, hamster, dog, ferret) have been shown to be at least moderately
permissive to SARS-CoV-2 infection, supporting a correlation between binding
affinity and permissiveness. Our data indicates that the earliest isolates of
SARS-CoV-2 were surprisingly well adapted to human ACE2, potentially explaining
its rapid transmission.


  Immunization with an inactivated virus is one of the strategies currently
being tested towards developing a SARS-CoV-2 vaccine. One of the methods used
to inactivate viruses is exposure to high doses of ionizing radiation to damage
their nucleic acids. Although gamma-rays effectively induce lesions in the RNA,
envelope proteins are also highly damaged in the process. This in turn may
alter their antigenic properties, affecting their capacity to induce an
adaptive immune response able to confer effective protection. Here, we modelled
the impact of sparsely and densely ionizing radiation on SARS-CoV-2 using the
Monte Carlo toolkit Geant4-DNA. With a realistic 3D target virus model, we
calculated the expected number of lesions in the spike and membrane proteins,
as well as in the viral RNA. We show that gamma-rays produce significant spike
protein damage, but densely ionizing charged particles induce less membrane
damage for the same level of RNA lesions, because a single ion traversal
through the nuclear envelope is sufficient to inactivate the virus. We propose
that accelerated charged particles produce inactivated viruses with little
structural damage to envelope proteins, thereby representing a new and
effective tool for developing vaccines against SARS-CoV-2 and other enveloped
viruses.


  One important decision, especially for long-lived species, is where and when
to start reproduction, i.e. where and when to recruit. For an informed decision
on recruitment, individuals may prospect potential breeding sites and gather
information on the future breeding site, hence reducing uncertainty.
Recruitment may also be delayed or advanced over a range of sexually mature
ages and this individual decision may influence survival costs of first
reproduction, senescence patterns and ultimately lifetime reproductive success.
Here we show that recruitment in the long-lived Audouin's gull is mostly
informed, meaning that mainly occurs when the individual has visited the colony
before. Individuals recruiting without gathering information were all young
first-time breeders. Interestingly, these recruits were affected by demographic
carry-over effects: first, they were much less prone to reproduce in subsequent
seasons than individuals that performed informed recruitment. More strikingly,
uninformed recruitment also affected survival, which was much lower for those
individuals breeding without gathering information on the breeding site. This
study suggests that taking stock before recruitment may be crucial for reducing
environmental uncertainty and for maximising individual fitness prospects.


  Understanding the dynamics of the COVID-19 pandemic is crucial for improved
control and social distancing strategies. To that effect, we have employed the
susceptible-exposed-infectious-recovered model, refined by contact tracing and
hospitalization data from Indian provinces Kerala, Delhi, Maharashtra, and West
Bengal, as well as from overall India. We have performed a sensitivity analysis
to identify the most crucial input parameters, and we have calibrated the model
to describe the data as best as possible. Short-term predictions reveal an
increasing and worrying trend of COVID-19 cases for all four provinces and
India as a whole, while long-term predictions also reveal the possibility of
oscillatory dynamics. Our research thus leaves the option open that COVID-19
might become a seasonal occurrence. We also simulate and discuss the impact of
media on the dynamics of the COVID-19 pandemic.


  A pandemic caused by a new coronavirus has spread worldwide, affecting
Argentina. We implement an SEIR model to analyze the disease evolution in
Buenos Aires and neighbouring cities. The model parameters are calibrated using
the number of casualties officially reported. Since infinite solutions honour
the data, we show different cases. In all of them the reproduction ratio $R_0$
decreases after early lockdown, but then raises, probably due to an increase in
contagion in highly populated slums. Therefore it is mandatory to reverse this
growing trend in $R_0$ by applying control strategies to avoid a high number of
infectious and dead individuals. The model provides an effective procedure to
estimate epidemic parameters (fatality rate, transmission probability,
infection and incubation periods) and monitor control measures during the
epidemic evolution.


  The compound 2-mercaptobenzothizaole (MBT) has been frequently detected in
wastewater and surface water and is a potential threat to both aquatic
organisms and human health (its mutagenic potential has been demonstrated).
This study investigated the degradation routes of MBT in the anode of a
microbial electrolysis cell (MEC) and the involved microbial communities. The
results indicated that graphene-modified anodes promoted the presence of more
enriched, developed, and specific communities compared to bare anodes.
Moreover, consecutive additions of the OH substituent to the benzene ring of
MBT were only detected in the reactor equipped with the graphene-treated
electrode. Both phenomena, together with the application of an external
voltage, may be related to the larger reduction of biotoxicity observed in the
MEC equipped with graphene-modified anodes (46.2 eqtox/m3 to 27.9 eqtox/m3).


  We consider here an extended SIR model, including several features of the
recent COVID-19 outbreak: in particular the infected and recovered individuals
can either be detected (+) or undetected (-) and we also integrate an intensive
care unit (ICU) capacity. Our model enables a tractable quantitative analysis
of the optimal policy for the control of the epidemic dynamics using both
lockdown and detection intervention levers. With parametric specification based
on literature on COVID-19, we investigate the sensitivities of various
quantities on the optimal strategies, taking into account the subtle trade-off
between the sanitary and the socio-economic cost of the pandemic, together with
the limited capacity level of ICU. We identify the optimal lockdown policy as
an intervention structured in 4 successive phases: First a quick and strong
lockdown intervention to stop the exponential growth of the contagion; second a
short transition phase to reduce the prevalence of the virus; third a long
period with full ICU capacity and stable virus prevalence; finally a return to
normal social interactions with disappearance of the virus. The optimal
scenario hereby avoids the second wave of infection, provided the lockdown is
released sufficiently slowly. We also provide optimal intervention measures
with increasing ICU capacity, as well as optimization over the effort on
detection of infectious and immune individuals. Whenever massive resources are
introduced to detect infected individuals, the pressure on social distancing
can be released, whereas the impact of detection of immune individuals reveals
to be more moderate.


  Beyond the genetic code, there is another layer of information encoded as
chemical modifications on histone proteins positioned along the DNA.
Maintaining these modifications is crucial for survival and identity of cells.
How the information encoded in the histone marks gets inherited, given that
only half the parental nucleosomes are transferred to each daughter chromatin,
is a puzzle. We address this problem using ideas from Information theory and
understanding from recent biological experiments. Mapping the replication and
reconstruction of modifications to equivalent problems in communication, we ask
how well an enzyme-machinery can recover information, if they were ideal
computing machines. Studying a parameter regime where realistic enzymes can
function, our analysis predicts that, pragmatically, enzymes may implement a
threshold$-k$ filling algorithm which derives from maximum \`a posteriori
probability decoding. Simulations using our method produce modification
patterns similar to what is observed in recent experiments.


  The SIR evolutionary model predicts too sharp a decrease of the fractions of
people infected with COVID-19 in France after the start of the national
lockdown, compared to what is observed. I fit the daily hospital data: arrivals
in regular and critical care units, releases and deaths, using extended SEIR
models. These involve ratios of evolutionary timescales to branching fractions,
assumed uniform throughout a country, and the basic reproduction number, $R_0$,
before and during the national lockdown, for each region of France. The
joint-region Bayesian analysis allows precise evaluations of the time/fraction
ratios and pre-hospitalized fractions. The hospital data are well fit by the
models, except the arrivals in critical care, which decrease faster than
predicted, indicating better treatment over time. Averaged over France, the
analysis yields $R_0$= 3.4$\pm$0.1 before the lockdown and 0.65$\pm$0.04 (90%
c.l.) during the lockdown, with small regional variations. On 11 May 2020, the
Infection Fatality Rate in France was 4 $\pm$1% (90% c.l.), while the Feverish
vastly outnumber the Asymptomatic, contrary to the early phases. Without the
lockdown nor social distancing, over 2 million deaths from COVID-19 would have
occurred throughout France, while a lockdown that would have been enforced 10
days earlier would have led to less than 1000 deaths. The fraction of immunized
people reached a plateau below 1% throughout France (3% in Paris) by late April
2020 (95% c.l.), suggesting a lack of herd immunity. The widespread
availability of face masks on 11 May, when the lockdown was partially lifted,
should keep $R_0$ below unity if at least 46% of the population wear them
outside their home. Otherwise, without enhanced other social distancing, a
second wave is inevitable and cause the number of deaths to triple between
early May and October (if $R_0$=1.2) or even late June (if $R_0$=2).


  Phylogenomics is a new field which applies to tools in phylogenetics to
genome data. Due to a new technology and increasing amount of data, we face new
challenges to analyze them over a space of phylogenetic trees. Because a space
of phylogenetic trees with a fixed set of labels on leaves is not Euclidean, we
cannot simply apply tools in data science. In this paper we survey some new
developments of machine learning models using tropical geometry to analyze a
set of phylogenetic trees over a tree space.


  Endeavoring toward a transferable, predictive coarse-grained explicit-chain
model for biomolecular condensates underlain by liquid-liquid phase separation
(LLPS), we conducted multiple-chain simulations of the N-terminal intrinsically
disordered region (IDR) of DEAD-box helicase Ddx4, as a test case, to assess
the roles of electrostatic, hydrophobic, cation-$\pi$, and aromatic
interactions in amino acid sequence-dependent LLPS. We evaluated 3
residue-residue interaction schemes with a shared electrostatic potential.
Neither a common hydrophobicity scheme nor one augmented with
arginine/lysine-aromatic cation-$\pi$ interactions consistently accounted for
the experimental LLPS data on the wildtype, a charge-scrambled, an FtoA, and an
RtoK mutant of Ddx4 IDR. In contrast, interactions based on contact statistics
among folded globular protein structures reproduce the overall experimental
trend, including that the RtoK mutant has a much diminished LLPS propensity.
Consistency between simulation and LLPS experiment was also found for RtoK
mutants of P-granule protein LAF-1, underscoring that, to a degree, the
important LLPS-driving $\pi$-related interactions are embodied in classical
statistical potentials. Further elucidation will be necessary, however,
especially of phenylalanine's role in condensate assembly because experiments
on FtoA and YtoF mutants suggest that LLPS-driving phenylalanine interactions
are significantly weaker than those posited by common statistical potentials.
Protein-protein electrostatic interactions are modulated by relative
permittivity, which depends on protein concentration. Analytical theory
suggests that this dependence entails enhanced inter-protein interactions in
the condensed phase but more favorable protein-solvent interactions in the
dilute phase. The opposing trends lead to a modest overall impact on LLPS.


  Several analytical models have been used in this work to describe the
evolution of death cases arising from coronavirus (COVID-19). The Death or `D'
model is a simplified version of the SIR (susceptible-infected-recovered)
model, which assumes no recovery over time, and allows for the
transmission-dynamics equations to be solved analytically. The D-model can be
extended to describe various focuses of infection, which may account for the
original pandemic (D1), the lockdown (D2) and other effects (Dn). The evolution
of the COVID-19 pandemic in several countries (China, Spain, Italy, France, UK,
Iran, USA and Germany) shows a similar behavior in concord with the D-model
trend, characterized by a rapid increase of death cases followed by a slow
decline, which are affected by the earliness and efficiency of the lockdown
effect. These results are in agreement with more accurate calculations using
the extended SIR model with a parametrized solution and more sophisticated
Monte Carlo grid simulations, which predict similar trends and indicate a
common evolution of the pandemic with universal parameters.


  Genome organization in eukaryotes during interphase stems from the delicate
balance between non-random correlations present in the DNA polynucleotide
linear sequence and the physico/chemical reactions which shape continuously the
form and structure of DNA and chromatin inside the nucleus of the cell. It is
now clear that these mechanisms have a key role in important processes like
gene regulation, yet the detailed ways they act simultaneously and, eventually,
come to influence each other even across very different length-scales remain
largely unexplored. In this paper, we recapitulate some of the main results
concerning gene regulatory and physical mechanisms, in relation to the
information encoded in the 1D sequence and the 3D folding structure of DNA. In
particular, we stress how reciprocal crossfeeding between 1D and 3D models may
provide original insight into how these complex processes work and influence
each other.


  The current COVID-19 pandemic has proven that proper control and prevention
of infectious disease require creating and enforcing the appropriate public
policies. One critical policy imposed by the policymakers is encouraging the
population to practice social distancing (i.e. controlling the contact rate
among the population). Here we pose a mean-field game model of individuals each
choosing a dynamic strategy of making contacts, given the trade-off of gaining
utility but also risking infection from additional contacts. We compute and
compare the mean-field equilibrium (MFE) strategy, which assumes each
individual acting selfishly to maximize its own utility, to the socially
optimal strategy, which maximizes the total utility of the population. We prove
that the optimal decision of the infected is always to make more contacts than
the level at which it would be socially optimal, which reinforces the important
role of public policy to reduce contacts of the infected (e.g. quarantining,
sick paid leave). Additionally, we include cost to incentivize people to change
strategies, when computing the socially optimal strategies. We find that with
this cost, policies reducing contacts of the infected should be further
enforced after the peak of the epidemic has passed. Lastly, we compute the
price of anarchy (PoA) of this system, to understand the conditions under which
large discrepancies between the MFE and socially optimal strategies arise,
which is when intervening public policy would be most effective.


  In a pandemic like Covid-19, there are many countries of lower-earning cannot
provide a complete locked-down within the duration of the detected case. The
locked-down may result in famine throughout the region of underdeveloped
countries after the outbreak. So, a conjectural setup of an epidemic has been
studied by applying specific period of locked-down (30 days) in 5 different
scenarios. The stochastic approach to the SEIR (Susceptible, Exposed, Infected
and Recovered) model has been used to evaluate the dynamics and the effects of
locked-down. It is observed that there exist a suitable period to apply
locked-down where more susceptible escape from the infection. The effect of the
early (as soon as the infected case detected) and late (with respect to the
estimated peak of detected cases for no locked-down) implementation of the
locked-down has also been studied and found that the late implementation of
locked-down will take the least time to end the epidemic. The CFR (Case
Fatality Rate) has also been found to be varied from 7.55 to 8.02 for all the
considered scenarios.


  In this paper we conduct a simulation study of the spread of an epidemic like
COVID-19 with temporary immunity on finite spatial and non-spatial network
models. In particular, we assume that an epidemic spreads stochastically on a
scale-free network and that each infected individual in the network gains a
temporary immunity after its infectious period is over. After the temporary
immunity period is over, the individual becomes susceptible to the virus again.
When the underlying contact network is embedded in Euclidean geometry, we model
three different intervention strategies that aim to control the spread of the
epidemic: social distancing, restrictions on travel, and restrictions on
maximal number of social contacts per node. Our first finding is that on a
finite network, a long enough average immunity period leads to extinction of
the pandemic after the first peak, analogous to the concept of "herd immunity".
For each model, there is a critical average immunity length $L_c$ above which
this happens. Our second finding is that all three interventions manage to
flatten the first peak (the travel restrictions most efficiently), as well as
decrease the critical immunity length $L_c$, but elongate the epidemic.
However, when the average immunity length $L$ is shorter than $L_c$, the price
for the flattened first peak is often a high second peak: for limiting the
maximal number of contacts, the second peak can be as high as 1/3 of the first
peak, and twice as high as it would be without intervention. Thirdly,
interventions introduce oscillations into the system and the time to reach
equilibrium is, for almost all scenarios, much longer. We conclude that
network-based epidemic models can show a variety of behaviors that are not
captured by the continuous compartmental models.


  We study the reported data from the COVID-19 pandemic outbreak in January -
May 2020 in 119 countries. We observe that the time series of active cases in
individual countries (the difference of the total number of confirmed
infections and the sum of the total number of reported deaths and recovered
cases) display a strong agreement with polynomial growth and at a later
epidemic stage also with a combined polynomial growth with exponential decay.
Our results are also formulated in terms of compartment type mathematical
models of epidemics. Within these models the universal scaling characterizing
the observed regime in an advanced epidemic stage can be interpreted as an
algebraic decay of the relative reproduction number $R_0$ as $T_M/t$, where
$T_M$ is a constant and $t$ is the duration of the epidemic outbreak. We show
how our findings can be applied to improve predictions of the reported pandemic
data and estimate some epidemic parameters. Note that although the model shows
a good agreement with the reported data we do not make any claims about the
real size of the pandemics as the relation of the observed reported data to the
total number of infected in the population is still unknown.


  Feedback in sensory biochemical networks can give rise to bifurcations in
cells' behavioral response. These bifurcations share many properties with
thermodynamic critical points. Evidence suggests that biological systems may
operate near these critical points, but the functional benefit of doing so
remains poorly understood. Here we investigate a simple biochemical model with
nonlinear feedback and multicellular communication to determine if criticality
provides a functional benefit in terms of the ability to gain information about
a stochastic chemical signal. We find that when signal fluctuations are slow,
the mutual information between the signal and the intracellular readout is
maximized at criticality, because the benefit of high signal susceptibility
outweighs the detriment of high readout noise. When cells communicate,
criticality gives rise to long-range correlations in readout molecule number
among cells. Consequently, we find that communication increases the information
between a given cell's readout and the spatial average of the signal across the
population. Finally, we find that both with and without communication, the
sensory benefits of criticality compete with critical slowing down, such that
the information rate, as opposed to the information itself, is minimized at the
critical point. Our results reveal the costs and benefits of feedback-induced
criticality for multicellular sensing.


  DNA code design aims to generate a set of DNA sequences (codewords) with
minimum likelihood of undesired hybridizations among sequences and their
reverse-complement (RC) pairs (cross-hybridization). Inspired by the distinct
hybridization affinities (or stabilities) of perfect double helix constructed
by individual single-stranded DNA (ssDNA) and its RC pair, we propose a novel
similarity significance (SS) model to measure the similarity between DNA
sequences. Particularly, instead of directly measuring the similarity of two
sequences by any metric/approach, the proposed SS works in a way to evaluate
how more likely will the undesirable hybridizations occur over the desirable
hybridizations in the presence of the two measured sequences and their RC
pairs. With this SS model, we construct thermodynamically stable DNA codes
subject to several combinatorial constraints using a sorting-based algorithm.
The proposed scheme results in DNA codes with larger code sizes and wider free
energy gaps (hence better cross-hybridization performance) compared to the
existing methods.


  Background: The ongoing COVID-19 epidemic dilated rapidly throughout India.
To end the global COVID-19 pandemic major behavioral, social distancing,
contact tracing, and state interventions has been undertaken to reduce the
outbreak and avert the persistence of the coronavirus in humans in India and
worldwide. In absence of any vaccine or therapeutics, forecasting is of utmost
priority for health care planning and control the transmission of COVID-19.
  Methods: We have proposed a mathematical model that explain the transmission
dynamics of COVID-19 in India. Based on the estimated data our model predicts
the evolution of epidemics and the end of SARS-CoV-2 and aids to evaluate the
influence of different policies to control the spread of the diseases.
  Findings: With the real data for infected individuals, we find the basic
reproduction number $R_0$ for 17 states of India and overall India. A complete
figure is given to demonstrate the estimated pandemic life cycle along with the
real data or history to date. Our study reveals that the strict control
measures implemented in India substantially mitigated the disseminate of
SARS-CoV-2. Importantly, model simulations predict that 95% reduction of
outbreak on June 26, 2020 and 99% reduction of outbreak on July 26, 2020 in
India.
  Interpretation: Our model simulation demonstrates that the community-wide
elimination of SARS-CoV-2 is possible by mitigating the social distancing and
use essential precautions. Lockdown can be implemented strictly to prevent the
human-to-human transmission of COVID-19. The model-based and parameter
estimation of epidemic life cycle, and end dates, if can be done precisely, may
decrease distress and over optimism and develop the mentality for all of us for
the next stages of the outbreak evolution.


  An active loop-extrusion mechanism is regarded as the main
out--of--equilibrium mechanism responsible for the structuring of
megabase-sized domains in chromosomes. We developed a model to study the
dynamics of the chromosome fibre by solving the kinetic equations associated
with the motion of the extruder. By averaging out the position of the extruder
along the chain, we build an effective equilibrium model capable of reproducing
experimental contact maps based solely on the positions of extrusion--blocking
proteins. We assessed the quality of the effective model using numerical
simulations of chromosomal segments and comparing the results with
explicit-extruder models and experimental data.


  The novel Wuhan coronavirus known as SARS-CoV-2 has brought almost
unprecedented effects for a non-wartime setting, hitting social, economic and
health systems hard.~ Being able to bring to bear pharmaceutical interventions
to counteract its effects will represent a major turning point in the fight to
turn the tides in this ongoing battle.~ Recently, the World's most powerful
supercomputer, SUMMIT, was used to identify existing small molecule
pharmaceuticals which may have the desired activity against SARS-CoV-2 through
a high throughput virtual screening approach. In this communication, we
demonstrate how the use of Bayesian optimization can provide a valuable service
for the prioritisation of these calculations, leading to the accelerated
identification of high-performing candidates, and thus expanding the scope of
the utility of HPC systems for time critical screening


  We study model embeddability, which is a variation of the famous embedding
problem in probability theory, when apart from the requirement that the Markov
matrix is the matrix exponential of a rate matrix, we additionally ask that the
rate matrix follows the model structure. We provide a characterisation of model
embeddable Markov matrices corresponding to symmetric group-based phylogenetic
models. In particular, we provide necessary and sufficient conditions in terms
of the eigenvalues of symmetric group-based matrices. To showcase our main
result on model embeddability, we provide an application to hachimoji models,
which are eight-state models for synthetic DNA. Moreover, our main result on
model embeddability enables us to compute the volume of the set of model
embeddable Markov matrices relative to the volume of other relevant sets of
Markov matrices within the model.


  During the current Covid-19 pandemic in Italy, official data are collected
with medical swabs following a pure convenience criterion which, at least in an
early phase, has privileged the exam of patients showing evident symptoms.
However, there are evidences of a very high proportion of asymptomatic patients
(e. g. Aguilar et al., 2020; Chugthai et al, 2020; Li, et al., 2020; Mizumoto
et al., 2020a, 2020b and Yelin et al., 2020). In this situation, in order to
estimate the real number of infected (and to estimate the lethality rate), it
should be necessary to run a properly designed sample survey through which it
would be possible to calculate the probability of inclusion and hence draw
sound probabilistic inference. Some researchers proposed estimates of the total
prevalence based on various approaches, including epidemiologic models, time
series and the analysis of data collected in countries that faced the epidemic
in earlier time (Brogi et al., 2020). In this paper, we propose to estimate the
prevalence of Covid-19 in Italy by reweighting the available official data
published by the Istituto Superiore di Sanit\`a so as to obtain a more
representative sample of the Italian population. Reweighting is a procedure
commonly used to artificially modify the sample composition so as to obtain a
distribution which is more similar to the population (Valliant et al., 2018).
In this paper, we will use post-stratification of the official data, in order
to derive the weights necessary for reweighting them using age and gender as
post-stratification variables thus obtaining more reliable estimation of
prevalence and lethality.


  As the COVID-19 outbreak is developing the two most frequently reported
statistics seem to be the raw confirmed case and case fatalities counts.
Focusing on Italy, one of the hardest hit countries, we look at how these two
values could be put in perspective to reflect the dynamics of the virus spread.
In particular, we find that merely considering the confirmed case counts would
be very misleading. The number of daily tests grows, while the daily fraction
of confirmed cases to total tests has a change point. It (depending on region)
generally increases with strong fluctuations till (around, depending on region)
15th-22nd March and then decreases linearly after. Combined with the increasing
trend of daily performed tests, the raw confirmed case counts are not
representative of the situation and are confounded with the sampling effort.
This we observe when regressing on time the logged fraction of positive tests
and for comparison the logged raw confirmed count. Hence, calibrating model
parameters for this virus's dynamics should not be done based only on confirmed
case counts (without rescaling by the number of tests), but take also
fatalities and hospitalization count under consideration as variables not prone
to be distorted by testing efforts. Furthermore, reporting statistics on the
national level does not say much about the dynamics of the disease, which are
taking place at the regional level. These findings are based on the official
data of total death counts up to 15th April 2020 released by ISTAT and up to
10th May 2020 for the number of cases. In this work we do not fit models but we
rather investigate whether this task is possible at all. This work also informs
about a new tool to collect and harmonize official statistics coming from
different sources in the form of a package for the R statistical environment
and presents the COVID-19 Data Hub.


  We give a counterexample to the conjecture of Martin and Thatte that two
balanced rooted binary leaf-labelled trees on $n$ leaves have a maximum
agreement subtree (MAST) of size at least $n^{\frac{1}{2}}$. In particular, we
show that for any $c>0$, there exist two balanced rooted binary leaf-labelled
trees on $n$ leaves such that any MAST for these two trees has size less than
$c n^{\frac{1}{2}}$. We also improve the lower bound of the size of such a MAST
to $n^{\frac{1}{6}}$.


  We analyse the temporal and regional structure in mortality rates related to
COVID-19 infections. We relate the fatality date of each deceased patient to
the corresponding day of registration of the infection, leading to a nowcasting
model which allows us to estimate the number of present-day infections that
will, at a later date, prove to be fatal. The numbers are broken down to the
district level in Germany. Given that death counts generally provide more
reliable information on the spread of the disease compared to infection counts,
which inevitably depend on testing strategy and capacity, the proposed model
and the presented results allow to obtain reliable insight into the current
state of the pandemic in Germany.


  Proteins play a key role in facilitating the infectiousness of the 2019 novel
coronavirus. A specific spike protein enables this virus to bind to human
cells, and a thorough understanding of its 3-dimensional structure is therefore
critical for developing effective therapeutic interventions. However, its
structure may continue to evolve over time as a result of mutations. In this
paper, we use a data science perspective to study the potential structural
impacts due to ongoing mutations in its amino acid sequence. To do so, we
identify a key segment of the protein and apply a sequential Monte Carlo
sampling method to detect possible changes to the space of low-energy
conformations for different amino acid sequences. Such computational approaches
can further our understanding of this protein structure and complement
laboratory efforts.


  Fire has been an integral part of the Earth for millennia. Several recent
wildfires have exhibited an unprecedented spatial and temporal extent and their
control is beyond national firefighting capabilities. Prescribed or controlled
burning treatments are debated as a potential measure for ameliorating the
spread and intensity of wildfires. Machine learning analysis using random
forests was performed in a spatio-temporal data set comprising a large number
of savanna fires across 22 years. Results indicate that fire return interval
was not an important predictor of fire spread rate or fire intensity, having a
feature importance of 3.5%, among eight other predictor variables. Manipulating
burn seasonality showed a feature importance of 6% or less regarding fire
spread rate or fire intensity. While manipulated fire return interval and
seasonality moderated both fire spread rate and intensity, their overall
effects were low in comparison with meteorological (hydrological and climatic)
variables. The variables with the highest feature importance regarding fire
spread rate resulted in fuel moisture with 21%, relative humidity with 15%,
wind speed with 14%, and last years rainfall with 14%. The variables with the
highest feature importance regarding fire intensity included fuel load with
21.5%, fuel moisture with 16.5%, relative humidity with 12.5%, air temperature
with 12.5%, and rainfall with 12.5%. Predicting fire spread rate and intensity
has been a poor endeavour thus far and we show that more data of the variables
already monitored would not result in higher predictive accuracy.


  The COVID-19 epidemic that emerged in Wuhan China at the end of 2019 hit
Italy particularly hard, yielding the implementation of strict national
lockdown rules (Phase 1). There is now a hot ongoing debate in Italy and abroad
on what the best strategy is to restart a country to exit a national lockdown
(Phase 2). Previous studies have focused on modelling possible restarting
scenarios at the national level, overlooking the fact that Italy, as other
nations around the world, is divided in administrative regions who can
independently oversee their own share of the Italian National Health Service.
In this study, we show that regionalism, and heterogeneity between regions, is
essential to understand the spread of the epidemic and, more importantly, to
design effective post Lock-Down strategies to control the disease. To achieve
this, we model Italy as a network of regions and parameterize the model of each
region on real data spanning almost two months from the initial outbreak. Using
the model, we confirm the effectiveness at the regional level of the national
lockdown strategy implemented so far by the Italian government to mitigate the
spread of the disease and show its efficacy at the regional level. We also
propose that differentiated, albeit coordinated, regional interventions can be
effective in Phase 2 to restart the country and avoid future recurrence of the
epidemic, while avoiding saturation of the regional health systems and
mitigating impact on costs. Our study and methodology can be easily extended to
other levels of granularity (provinces or counties in the same region or states
in other federal countries, etc.) to support policy- and decision-makers.


  Network-based intervention strategies can be effective and cost-efficient
approaches to curtailing harmful contagions in myriad settings. As studied,
these strategies are often impractical to implement, as they typically assume
complete knowledge of the network structure, which is unusual in practice. In
this paper, we investigate how different immunization strategies perform under
realistic conditions where the strategies are informed by partially-observed
network data. Our results suggest that global immunization strategies, like
degree immunization, are optimal in most cases; the exception is at very high
levels of missing data, where stochastic strategies, like acquaintance
immunization, begin to outstrip them in minimizing outbreaks. Stochastic
strategies are more robust in some cases due to the different ways in which
they can be affected by missing data. In fact, one of our proposed variants of
acquaintance immunization leverages a logistically-realistic ongoing
survey-intervention process as a form of targeted data-recovery to improve with
increasing levels of missing data. These results support the effectiveness of
targeted immunization as a general practice. They also highlight the risks of
considering networks as idealized mathematical objects: overestimating the
accuracy of network data and foregoing the rewards of additional inquiry.


  In this article we propose a compartmental model for the dynamics of
Coronavirus Disease 2019 (COVID-19). We take into account the presence of
asymptomatic infections and the main policies that have been adopted so far to
contain the epidemic: isolation (or social distancing) of a portion of the
population, quarantine for confirmed cases and testing. We model isolation by
separating the population in two groups: one composed by key-workers that keep
working during the pandemic and have a usual contact rate, and a second group
consisting of people that are enforced/recommended to stay at home. We refer to
quarantine as strict isolation, and it is applied to confirmed infected cases.
In the proposed model, the proportion of people in isolation, the level of
contact reduction and the testing rate are control parameters that can vary in
time, representing policies that evolve in different stages. We obtain an
explicit expression for the basic reproduction number $\mathcal{R}_0$ in terms
of the parameters of the disease and of the control policies. In this way we
can quantify the effect that isolation and testing have in the evolution of the
epidemic. We present a series of simulations to illustrate different realistic
scenarios. From the expression of $\mathcal{R}_0$ and the simulations we
conclude that isolation (social distancing) and testing among asymptomatic
cases are fundamental actions to control the epidemic, {and the stricter these
measures are and the sooner they are implemented,} the more lives can be saved.
Additionally, we show that people that remain in isolation significantly reduce
their probability of contagion, so risk groups should be recommended to
maintain a low contact rate during the course of the epidemic.


  Most epidemic models are spatially aggregate and the index which is most used
for planning and policy numbers, the r number, typically refers to a single
system of interest. Even if r numbers are calculated for each of adjacent
areas, regions or countries for example, there is no interaction between them.
Here we aim to offer a fine-grained geography: models of epidemics in spatially
disaggregated systems with interactions. This offers the possibility of new
insights into the dynamics of epidemics and of policies aimed at mitigation and
control.


  Predicting accurate protein-ligand binding affinity is important in drug
discovery but remains a challenge even with computationally expensive
biophysics-based energy scoring methods and state-of-the-art deep learning
approaches. Despite the recent advances in the deep convolutional and graph
neural network based approaches, the model performance depends on the input
data representation and suffers from distinct limitations. It is natural to
combine complementary features and their inference from the individual models
for better predictions. We present fusion models to benefit from different
feature representations of two neural network models to improve the binding
affinity prediction. We demonstrate effectiveness of the proposed approach by
performing experiments with the PDBBind 2016 dataset and its docking pose
complexes. The results show that the proposed approach improves the overall
prediction compared to the individual neural network models with greater
computational efficiency than related biophysics based energy scoring
functions. We also discuss the benefit of the proposed fusion inference with
several example complexes. The software is made available as open source at
https://github.com/llnl/fast.


  The origins of the large differences observed to the rates with which the
diverse particles are conveyed along axonal microtubules are still a matter of
debate in the literature. There is evidence that certain neurodegenerative
diseases may be triggered by disturbances to the related transport processes.
Motivated by this, we employ a model to investigate the mobility properties of
certain cargoes which dynamics are coupled with that of molecular motors on
crowded microtubules. For certain initial and boundary conditions, we use the
method of characteristics to resolve perturbatively the pair of equations of
the Burgers type resulting from a mean-field approach to the original
microscopic stochastic model. Extensions to the non-perturbative limits are
explored numerically. In this context, we were able to figure out conditions
under which cargos average velocities may differ up to orders of magnitude just
by changing the number of motors on the considered track. We then discuss
possibilities to connect these theoretical predictions with available
experimental data about axon transport.


  Susceptible-Invective-Recovered (SIR) mathematical models are in high demand
due to the COVID-19 pandemic. They are used in their standard formulation, or
through the many variants, trying to fit and hopefully predict the number of
new cases for the next days or weeks, in any place, city, or country. Such is
key knowledge for the authorities to prepare for the health systems demand or
to apply restrictions to slow down the infectives curve. Even when the model
can be easily solved ---by the use of specialized software or by programming
the numerical solution of the differential equations that represent the
model---, the prediction is a non-easy task, because the behavioral change of
people is reflected in a continuous change of the parameters. A relevant
question is what we can use of one city to another; if what happened in Madrid
could have been applied to New York and then, if what we have learned from this
city would be of use for S\~ao Paulo. With this idea in mind, we present an
analysis of a spreading-rate related measure of COVID-19 as a function of
population density and population size for all US counties, as long as for
Brazilian cities and German cities. Contrary to what is the common hypothesis
in epidemics modeling, we observe a higher {\em per-capita} contact rate for
higher city's population density and population size. Also, we find that the
population size has a more explanatory effect than the population density. A
contact rate scaling theory is proposed to explain the results.


  This work introduces a geometric framework and a novel network architecture
for creating correspondences between samples of different conditions. Under
this formalism, the latent space is a fiber bundle stratified into a base space
encoding conditions, and a fiber space encoding the variations within
conditions. Furthermore, this latent space is endowed with a natural pull-back
metric. The correspondences between conditions are obtained by minimizing an
energy functional, resulting in diffeomorphism flows between fibers.
  We illustrate this approach using MNIST and Olivetti and benchmark its
performances on the task of batch correction, which is the problem of
integrating multiple biological datasets together.


  Coronavirus disease 2019 (COVID-19) has impacted almost every part of human
life worldwide, posing a massive threat to human health. There is no specific
drug for COVID-19, highlighting the urgent need for the development of
effective therapeutics. To identify potentially repurposable drugs, we employed
a systematic approach to mine candidates from U.S. FDA-approved drugs and
preclinical small-molecule compounds by integrating the gene expression
perturbation data for chemicals from the Library of Integrated Network-Based
Cellular Signatures project with a publicly available single-cell RNA
sequencing dataset from mild and severe COVID-19 patients. We identified 281
FDA-approved drugs that have the potential to be effective against SARS-CoV-2
infection, 16 of which are currently undergoing clinical trials to evaluate
their efficacy against COVID-19. We experimentally tested the inhibitory
effects of tyrphostin-AG-1478 and brefeldin-a on the replication of the
single-stranded ribonucleic acid (ssRNA) virus influenza A virus. In
conclusion, we have identified a list of repurposable anti-SARS-CoV-2 drugs
using a systems biology approach.


  Virtual screening of phytochemicals was performed through molecular docking,
simulation, in silico ADMET and drug-likeness prediction to identify the
potential hits that can inhibit the effects of SARS-CoV-2. Considering the
published literature on medicinal importance, total 154 phytochemicals with
analogous structure from limonoids and triterpenoids were selected to search
potential inhibitors for the five therapeutic protein targets of SARS-CoV-2,
i.e., 3CLpro (main protease), PLpro (papain-like protease), SGp-RBD (spike
glycoprotein-receptor binding domain), RdRp (RNA dependent RNA polymerase) and
ACE2 (angiotensin-converting enzyme 2). The in silico computational results
revealed that the phytochemicals such as glycyrrhizic acid, limonin,
7-deacetyl-7-benzoylgedunin, maslinic acid, corosolic acid, obacunone and
ursolic acid were found to be effective against the target proteins of
SARS-CoV-2. The protein-ligand interaction study revealed that these
phytochemicals bind with the amino acid residues at the active site of the
target proteins. Therefore, the core structure of these potential hits can be
used for further lead optimization to design drugs for SARS-CoV-2. Also, the
medicinal plants containing these phytochemicals like licorice, neem, tulsi,
citrus and olives can be used to formulate suitable therapeutic approaches in
traditional medicines.


  By equipping a previously reported dynamic causal model of COVID-19 with an
isolation state, we modelled the effects of self-isolation consequent on
tracking and tracing. Specifically, we included a quarantine or isolation state
occupied by people who believe they might be infected but are asymptomatic, and
only leave if they test negative. We recovered maximum posteriori estimates of
the model parameters using time series of new cases, daily deaths, and tests
for the UK. These parameters were used to simulate the trajectory of the
outbreak in the UK over an 18-month period. Several clear-cut conclusions
emerged from these simulations. For example, under plausible (graded)
relaxations of social distancing, a rebound of infections within weeks is
unlikely. The emergence of a later second wave depends almost exclusively on
the rate at which we lose immunity, inherited from the first wave. There exists
no testing strategy that can attenuate mortality rates, other than by deferring
or delaying a second wave. A sufficiently powerful tracking and tracing
policy--implemented at the time of writing (10th May 2020)--will defer any
second wave beyond a time horizon of 18 months. Crucially, this deferment is
within current testing capabilities (requiring an efficacy of tracing and
tracking of about 20% of asymptomatic infected cases, with less than 50,000
tests per day). These conclusions are based upon a dynamic causal model for
which we provide some construct and face validation, using a comparative
analysis of the United Kingdom and Germany, supplemented with recent
serological studies.


  An epidemiological model is developed for the spread of COVID-19 in South
Africa. A variant of the classical compartmental SEIR model, called the SEIQRDP
model, is used. As South Africa is still in the early phases of the global
COVID-19 pandemic with the confirmed infectious cases not having peaked, the
SEIQRDP model is first parameterized on data for Germany, Italy, and South
Korea - countries for which the number of infectious cases are well past their
peaks. Good fits are achieved with reasonable predictions of where the number
of COVID-19 confirmed cases, deaths, and recovered cases will end up and by
when. South African data for the period from 23 March to 8 May 2020 is then
used to obtain SEIQRDP model parameters. It is found that the model fits the
initial disease progression well, but that the long-term predictive capability
of the model is rather poor. The South African SEIQRDP model is subsequently
recalculated with the basic reproduction number constrained to reported values.
The resulting model fits the data well, and long-term predictions appear to be
reasonable. The South African SEIQRDP model predicts that the peak in the
number of confirmed infectious individuals will occur at the end of October
2020, and that the total number of deaths will range from about 10,000 to
90,000, with a nominal value of about 22,000. All of these predictions are
heavily dependent on the disease control measures in place, and the adherence
to these measures. These predictions are further shown to be particularly
sensitive to parameters used to determine the basic reproduction number. The
future aim is to use a feedback control approach together with the South
African SEIQRDP model to determine the epidemiological impact of varying
lockdown levels proposed by the South African Government.


  The World Health Organization (WHO) made the assessment that COVID-19 can be
characterized as a pandemic on March 11, 2020. To the COVID-19 outbreak, there
is no vaccination and no treatment. The only way to control the COVID-19
outbreak is sustained physical distancing. In this work, a simple control law
was proposed to keep the infected individuals during the COVID-19 outbreak
below the desired number. The proposed control law keeps the value of infected
individuals controlled only adjusting the social distancing level. The
stability analysis of the proposed control law is done and the uncertainties in
the parameters were considered. A version of the proposed controller to daily
update was developed. This is a very simple approach to the developed control
law and can be calculated in a spreadsheet. In the end, numerical simulations
were done to show the behavior of the number of infected individuals during an
epidemic disease when the proposed control law is used to adjust the social
distancing level.


  The human gastrointestinal tract is an environment that hosts an ecosystem of
microorganisms essential to human health. Vital biological processes emerge
from fundamental inter- and intra-species molecular interactions that influence
the assembly and composition of the gut microbiota ecology. Here we quantify
the complexity of the ecological relationships within the human infant gut
microbiota ecosystem as a function of the information contained in the
nonlinear associations of a sequence of increasingly-specified maximum entropy
representations of the system. Our paradigm frames the ecological state, in
terms of the presence or absence of individual microbial ecological units that
are identified by amplicon sequence variants (ASV) in the gut microenvironment,
as a function of both the ecological states of its neighboring units and, in a
departure from standard graphical model representations, the associations among
the units within its neighborhood. We characterize the order of the system
based on the relative quantity of statistical information encoded by high-order
statistical associations of the infant gut microbiota.


  We are currently living in a state of uncertainty due to the pandemic caused
by the Sars-CoV-2 virus. There are several factors involved in the epidemic
spreading such as the individual characteristics of each city/country. The true
shape of the epidemic dynamics is a large, complex system such as most of the
social systems. In this context, Complex networks are a great candidate to
analyze these systems due to their ability to tackle structural and dynamical
properties. Therefore this study presents a new approach to model the COVID-19
epidemic using a multi-layer complex network, where nodes represent people,
edges are social contacts, and layers represent different social activities.
The model improves the traditional SIR and it is applied to study the Brazilian
epidemic by analyzing possible future actions and their consequences. The
network is characterized using statistics of infection, death, and
hospitalization time. To simulate isolation, social distancing, or
precautionary measures we remove layers and/or reduce the intensity of social
contacts. Results show that even taking various optimistic assumptions, the
current isolation levels in Brazil still may lead to a critical scenario for
the healthcare system and a considerable death toll (average of 149,000). If
all activities return to normal, the epidemic growth may suffer a steep
increase, and the demand for ICU beds may surpass 3 times the country's
capacity. This would surely lead to a catastrophic scenario, as our estimation
reaches an average of 212,000 deaths even considering that all cases are
effectively treated. The increase of isolation (up to a lockdown) shows to be
the best option to keep the situation under the healthcare system capacity,
aside from ensuring a faster decrease of new case occurrences (months of
difference), and a significantly smaller death toll (average of 87,000).


  Incidence vs Cumulative Cases (ICC) curves are introduced and shown to
provide a simple framework for parameter identification in the case of the most
elementary epidemiological model, consisting of susceptible, infected, and
removed compartments. This novel methodology is used to estimate the basic
reproduction ratio of recent outbreaks, including the ongoing COVID-19
epidemic.


  Ca$^{2+}$ plays an important role in cell signal transduction. Its
intracellular propagation is the most basic process of Ca$^{2+}$ signaling,
such as calcium wave and double messenger system. In this work, with both
numerical simulation and mean field ansatz, the 3-dimensional probability
distribution of Ca$^{2+}$, which is read out by phosphorylation, is studied in
two scenarios with boundaries. The coverage of distribution of Ca$^{2+}$ is
found at an order of magnitude of $\mu$m, which is consistent with experimental
observed calcium spike and wave. Our results suggest that the double messenger
system may occur in the ER-PM junction to acquire great efficiency. The buffer
effect of kinase is also discussed by calculating the average position of
phosphorylations and free Ca$^{2+}$. The results are helpful to understand the
mechanism of Ca$^{2+}$ signaling.


  We consider the pandemic spreading of COVID-19 in India after the outbreak of
the coronavirus in Wuhan city, China. We estimate the transmission rate of the
initial infecting individuals of COVID-19 in India by using the officially
reported data at the early stage of the epidemic with the help of Susceptible
(S), Exposed (E), Infected (I), and Removed (R) population model, the so-called
SEIR dynamical model. Numerical analysis and model verification are performed
to calibrate the system parameters with official public information about the
number of people infected, and then to evaluate several COVID -19 scenarios
potentially applicable to India. Our findings provide an estimation of disease
occurrence in the near future and also demonstrate the importance of
governmental and individual efforts to control the effects and time of the
pandemic-related critical situations. We also give special emphasis to
individual reactions in the containment process.


  COVID-19 is a global pandemic caused by the SARS-CoV-2 coronavirus. T cells
play a key role in the adaptive antiviral immune response by killing infected
cells and facilitating the selection of virus-specific antibodies. However
neither the dynamics and cross-reactivity of the SARS-CoV-2-specific T cell
response nor the diversity of resulting immune memory are well understood. In
this study we use longitudinal high-throughput T cell receptor (TCR) sequencing
to track changes in the T cell repertoire following two mild cases of COVID-19.
In both donors we identified CD4+ and CD8+ T cell clones with transient clonal
expansion after infection. The antigen specificity of CD8+ TCR sequences to
SARS-CoV-2 epitopes was confirmed by both MHC tetramer binding and presence in
large database of SARS-CoV-2 epitope-specific TCRs. We describe characteristic
motifs in TCR sequences of COVID-19-reactive clones and show preferential
occurence of these motifs in publicly available large dataset of repertoires
from COVID-19 patients. We show that in both donors the majority of
infection-reactive clonotypes acquire memory phenotypes. Certain T cell clones
were detected in the memory fraction at the pre-infection timepoint, suggesting
participation of pre-existing cross-reactive memory T cells in the immune
response to SARS-CoV-2.


  Stochastic discrete-time SIS and SIR models of endemic diseases are
introduced and analyzed. For the deterministic, mean-field model, the basic
reproductive number $R_0$ determines their global dynamics. If $R_0\le 1$, then
the frequency of infected individuals asymptotically converges to zero. If
$R_0>1$, then the infectious class uniformly persists for all time; conditions
for a globally stable, endemic equilibrium are given. In contrast, the
infection goes extinct in finite time with probability one in the stochastic
models for all $R_0$ values. To understand the length of the transient prior to
extinction as well as the behavior of the transients, the quasi-stationary
distributions and the associated mean time to extinction are analyzed using
large deviation methods. When $R_0>1$, these mean times to extinction are shown
to increase exponentially with the population size $N$. Moreover, as $N$
approaches $\infty$, the quasi-stationary distributions are supported by a
compact set bounded away from extinction; sufficient conditions for convergence
to a Dirac measure at the endemic equilibrium of the deterministic model are
also given. In contrast, when $R_0<1$, the mean times to extinction are bounded
above $1/(1-\alpha)$ where $\alpha<1$ is the geometric rate of decrease of the
infection when rare; as $N$ approaches $\infty$, the quasi-stationary
distributions converge to a Dirac measure at the disease-free equilibrium for
the deterministic model. For several special cases, explicit formulas for
approximating the quasi-stationary distribution and the associated mean
extinction are given. These formulas illustrate how for arbitrarily small $R_0$
values, the mean time to extinction can be arbitrarily large, and how for
arbitrarily large $R_0$ values, the mean time to extinction can be arbitrarily
large.


  To better predict the dynamics of spread of COVID-19 epidemics, it is
important not only to investigate the network of local and long-range
contagious contacts, but also to understand the temporal dynamics of
infectiousness and detectable symptoms. Here we present a model of infection
spread in a well-mixed group of individuals, which usually corresponds to a
node in large-scale epidemiological networks. The model uses delay equations
that take into account the duration of infection and is based on
experimentally-derived time courses of viral load, virus shedding, severity and
detectability of symptoms. We show that because of an early onset of
infectiousness, which is reported to be synchronous or even precede the onset
of detectable symptoms, the tracing and immediate testing of everyone who came
in contact with the detected infected individual reduces the spread of
epidemics, hospital load, and fatality rate. We hope that this more precise
node dynamics could be incorporated into complex large-scale epidemiological
models to improve the accuracy and credibility of predictions.


  In this paper, we deal with the study of the impact of nationwide measures
COVID-19 anti-pandemic. We drive two processes to analyze COVID-19 data
considering measures. We associate level of nationwide measure with value of
parameters related to the contact rate of the model. Then a parametric solve,
with respect to those parameters of measures, shows different possibilities of
the evolution of the pandemic. Two machine learning tools are used to forecast
the evolution of the pandemic. Finally, we show comparison between
deterministic and two machine learning tools.


  Recent computational advances in the accurate prediction of protein
three-dimensional (3D) structures from amino acid sequences now present a
unique opportunity to decipher the interrelationships between proteins. This
task entails--but is not equivalent to--a problem of 3D structure comparison
and classification. Historically, protein domain classification has been a
largely manual and subjective activity, relying upon various heuristics.
Databases such as CATH represent significant steps towards a more systematic
(and automatable) approach, yet there still remains much room for the
development of more scalable and quantitative classification methods, grounded
in machine learning. We suspect that re-examining these relationships via a
Deep Learning (DL) approach may entail a large-scale restructuring of
classification schemes, improved with respect to the interpretability of
distant relationships between proteins. Here, we describe our training of DL
models on protein domain structures (and their associated physicochemical
properties) in order to evaluate classification properties at CATH's
"homologous superfamily" (SF) level. To achieve this, we have devised and
applied an extension of image-classification methods and image segmentation
techniques, utilizing a convolutional autoencoder model architecture. Our DL
architecture allows models to learn structural features that, in a sense,
'define' different homologous SFs. We evaluate and quantify pairwise
'distances' between SFs by building one model per SF and comparing the loss
functions of the models. Hierarchical clustering on these distance matrices
provides a new view of protein interrelationships--a view that extends beyond
simple structural/geometric similarity, and towards the realm of
structure/function properties.


  Bayesian analysis of publicly available time series of cases and fatalities
in different geographical regions of India during April 2020 is reported. It is
found that the initial apparent rapid growthin infections could be partly due
to confounding factors such as initial rapid ramp-up of disease surveillance. A
brief discussion is given of the fallacies which arise if this possibility is
neglected. The growth after April 10 is consistent with a time independent but
region dependent exponential. From this, R0 is extracted using both known cases
and fatalities. The two estimates are seen to agree in many cases; for these
CFR is reported. It is seen that CFR and R0 increase together. Some public
health implications of this observation are discussed, including a target
doubling interval if medical facilities are to remain adequate.


  Sequences of nucleotides (for DNA and RNA) or amino acids (for proteins) are
central objects in biology. Among the most important computational problems is
that of sequence alignment, i.e. arranging sequences from different organisms
in such a way to identify similar regions, to detect evolutionary relationships
between sequences, and to predict biomolecular structure and function. This is
typically addressed through profile models, which capture
position-specificities like conservation in sequences, but assume an
independent evolution of different positions. Over the last years, it has been
well established that coevolution of different amino-acid positions is
essential for maintaining three-dimensional structure and function. Modeling
approaches based on inverse statistical physics can catch the coevolution
signal in sequence ensembles; and they are now widely used in predicting
protein structure, protein-protein interactions, and mutational landscapes.
Here, we present DCAlign, an efficient alignment algorithm based on an
approximate message-passing strategy, which is able to overcome the limitations
of profile models, to include coevolution among positions in a general way, and
to be therefore universally applicable to protein- and RNA-sequence alignment
without the need of using complementary structural information. The potential
of DCAlign is carefully explored using well-controlled simulated data, as well
as real protein and RNA sequences.


  We introduce a version of the Minority Game where the total number of
available choices is $D>2$, but the agents only have two available choices to
switch. For all agents at an instant in any given choice, therefore, the other
choice is distributed between the remaining $D-1$ options. This brings in the
added complexity in reaching a state with the maximum resource utilization, in
the sense that the game is essentially a set of MG that are coupled and played
in parallel. We show that a stochastic strategy, used in the MG, works well
here too. We discuss the limits in which the model reduces to other known
models. Finally, we study an application of the model in the context of
population movement between various states within a country during an ongoing
epidemic. We show that the total infected population in the country could be as
low as that achieved with a complete stoppage of inter-region movements for a
prolonged period, provided that the agents instead follow the above mentioned
stochastic strategy for their movement decisions between their two choices. The
objective for an agent is to stay in the lower infected state between their two
choices. We further show that it is the agents moving once between any two
states, following the stochastic strategy, who are less likely to be infected
than those not having (or not opting for) such a movement choice, when the risk
of getting infected during the travel is not considered. This shows the
incentive for the moving agents to follow the stochastic strategy.


  We present a phenomenological procedure of dealing with the COVID--19 data
provided by government health agencies of eleven different countries. Instead
of using the (exact or approximate) solutions to the SIR (or other) model(s) to
fit the data by adjusting the time--independent parameters included in those
models, we introduce dynamical parameters whose time--dependence may be
phenomenologically obtained by adequately extrapolating a chosen subset of the
daily provided data. This phenomenological approach works extremely well to
properly adjust the number of infected (and removed) individuals in time, for
the countries we consider. Besides, it can handle the sub--epidemic events that
some countries may experience. In this way, we obtain the evolution of the
pandemic without using any a priori model based on differential equations.


  In this paper, we present a model to predict the spread of the Covid-19
epidemic and apply it to the specific case of Italy. We started from a simple
Susceptible, Infected, Recovered (SIR) model and we added the condition that,
after a certain time, the basic reproduction number $R_0$ exponentially decays
in time, as empirically suggested by world data. Using this model, we were able
to reproduce the real behavior of the epidemic with an average error of 5\%.
Moreover, we illustrate possible future scenarios, associated to different
intervals of $R_0$. This model has been used since the beginning of March 2020,
predicting the Italian peak of the epidemic in April 2020 with about 100.000
detected active cases. The real peak of the epidemic happened on the 20th of
April 2020, with 108.000 active cases. This result shows that the model had
predictive power for the italian case.


  Many proteins that undergo sequential enzymatic modification in the Golgi
cisternae are displayed at the plasma membrane as cell identity markers. The
modified proteins, called glycans, represent a molecular code. The fidelity of
this glycan code is measured by how accurately the glycan synthesis machinery
realises the desired target glycan distribution for a particular cell type and
niche. In this paper, we quantitatively analyse the tradeoffs between the
number of cisternae and the number and specificity of enzymes, in order to
synthesize a prescribed target glycan distribution of a certain complexity. We
find that to synthesize complex distributions, such as those observed in real
cells, one needs to have multiple cisternae and precise enzyme partitioning in
the Golgi. Additionally, for fixed number of enzymes and cisternae, there is an
optimal level of specificity of enzymes that achieves the target distribution
with high fidelity. Our results show how the complexity of the target glycan
distribution places functional constraints on the Golgi cisternal number and
enzyme specificity.


  An ultrametric model of epidemic spread of infections based on the classical
SIR model is proposed. Ultrametrics on a set of individuals based on theire
hierarchical clustering relativly to the average time of infectious contact is
introduced. The general equations of the ultrametric SIR model are written down
and their particular implementation using $p$-adic parameterization is
presented. A numerical analysis of the $p$-adic SIR model and a comparison of
its behavior with the classical SIR model are performed. The concept of
hierarchical isolation and the scenario of its management in order to reduce
the level of epidemic spread is considered.


  The ongoing COVID-19 pandemic has overwhelmingly demonstrated the need to
accurately evaluate the effects of implementing new or altering existing
nonpharmaceutical interventions. Since these interventions applied at the
societal level cannot be evaluated through traditional experimental means,
public health officials and other decision makers must rely on statistical and
mathematical epidemiological models. Nonpharmaceutical interventions are
typically focused on contacts between members of a population, and yet most
epidemiological models rely on homogeneous mixing which has repeatedly been
shown to be an unrealistic representation of contact patterns. An alternative
approach is individual based models (IBMs), but these are often time intensive
and computationally expensive to implement, requiring a high degree of
expertise and computational resources. More often, decision makers need to know
the effects of potential public policy decisions in a very short time window
using limited resources. This paper presents an estimation algorithm for an IBM
designed to evaluate nonpharmaceutical interventions. By utilizing recursive
relationships, our method can quickly compute the expected epidemiological
outcomes even for large populations based on any arbitrary contact network. We
utilize our methods to evaluate the effects of relaxing current social
distancing measures in Iowa, USA, at various times and to various degrees.
\verb!R! code for our method is provided in the supplementary material, thereby
allowing others to utilize our approach for other regions.


  We revisit the modeling of the diauxic growth of a pure microorganism on two
distinct sugars which was first described by Monod. Most available models are
deterministic and make the assumption that all cells of the microbial ecosystem
behave homogeneously with respect to both sugars, all consuming the first one
and then switching to the second when the first is exhausted. We propose here a
stochastic model which describes what is called "metabolic heterogeneity". It
allows to consider small populations as in microfluidics as well as large
populations where billions of individuals coexist in the medium in a batch or
chemostat. We highlight the link between the stochastic model and the
deterministic behavior in real large cultures using a large population
approximation. Then the influence of model parameter values on model dynamics
is studied, notably with respect to the lag-phase observed in real systems
depending on the sugars on which the microorganism grows. It is shown that both
metabolic parameters as well as initial conditions play a crucial role on
system dynamics.


  OBJECTIVES: to describe the first wave of the COVID-19 pandemic with a focus
on undetected cases and to evaluate different post-lockdown scenarios. DESIGN:
the study introduces a SEIR compartmental model, taking into account the
region-specific fraction of undetected cases, the effects of mobility
restrictions, and the personal protective measures adopted, such as wearing a
mask and washing hands frequently. SETTING AND PARTICIPANTS: the model is
experimentally validated with data of all the Italian regions, some European
countries, and the US. MAIN OUTCOME MEASURES: the accuracy of the model results
is measured through the mean absolute percentage error (MAPE) and Lewis
criteria; fitting parameters are in good agreement with previous literature.
RESULTS: the epidemic curves for different countries and the amount of
undetected and asymptomatic cases are estimated, which are likely to represent
the main source of infections in the near future. The model is applied to the
Hubei case study, which is the first place to relax mobility restrictions.
Results show different possible scenarios. Mobility and the adoption of
personal protective measures greatly influence the dynamics of the infection,
determining either a huge and rapid secondary epidemic peak or a more delayed
and manageable one. CONCLUSIONS: mathematical models can provide useful
insights for healthcare decision makers to determine the best strategy in case
of future outbreaks.


  The optimal design of compounds through manipulating properties at the
molecular level is often the key to considerable scientific advances and
improved process systems performance. This paper highlights key trends,
challenges, and opportunities underpinning the Computer-Aided Molecular Design
(CAMD) problems. A brief review of knowledge-driven property estimation methods
and solution techniques, as well as corresponding CAMD tools and applications,
are first presented. In view of the computational challenges plaguing
knowledge-based methods and techniques, we survey the current state-of-the-art
applications of deep learning to molecular design as a fertile approach towards
overcoming computational limitations and navigating uncharted territories of
the chemical space. The main focus of the survey is given to deep generative
modeling of molecules under various deep learning architectures and different
molecular representations. Further, the importance of benchmarking and
empirical rigor in building deep learning models is spotlighted. The review
article also presents a detailed discussion of the current perspectives and
challenges of knowledge-based and data-driven CAMD and identifies key areas for
future research directions. Special emphasis is on the fertile avenue of hybrid
modeling paradigm, in which deep learning approaches are exploited while
leveraging the accumulated wealth of knowledge-driven CAMD methods and tools.


  We study the potential scenarios from a
Susceptible-Infected-Recovered-Asymptomatic-Symptomatic-Dead (SIRASD) model. As
a novelty, we consider populations that differ in their degree of compliance
with social distancing policies following socioeconomic attributes that are
observed in emerging and developing countries. Considering epidemiological
parameters estimated from data of the propagation of SARS-CoV-2 in Brazil --
where there is a significant stake of the population making their living in the
informal economy and thus prone to not follow self-isolation -- we assert that
if the confinement measures are lifted too soon, namely as much as one week of
consecutive declining numbers of new cases, it is very likely the appearance of
a second peak. Our approach should be valid for any country where the number of
people involved in the informal economy is a large proportion of the total
labor force. In summary, our results point out the crucial relevance of target
policies for supporting people in the informal economy to properly comply with
preventive measures during the pandemic.


  In Diffusive Molecular Communication (DMC), information is transmitted by
diffusing molecules. Synaptic signaling, as a natural implementation of this
paradigm, encompasses functional components that, once understood, can
facilitate the development of synthetic DMC systems. To unleash this potential,
however, a thorough understanding of the synaptic communication channel based
on biophysical principles is needed. Since synaptic transmission critically
depends also on non-neural cells, such understanding requires the consideration
of the so-called tripartite synapse. In this paper, we develop a comprehensive
channel model of the tripartite synapse encompassing a three-dimensional,
finite-size spatial model of the synaptic cleft, molecule uptake at the
presynaptic neuron and at glial cells, reversible binding to individual
receptors at the postsynaptic neuron, and spillover to the extrasynaptic space.
Based on this model, we derive analytical time domain expressions for the
channel impulse response (CIR) of the synaptic DMC system and for the number of
molecules taken up at the presynaptic neuron and at glial cells, respectively.
These expressions provide insight into the impact of macroscopic physical
channel parameters on the decay rate of the CIR and the reuptake rate, and
reveal fundamental limits for synaptic signal transmission induced by chemical
reaction kinetics and the channel geometry. Adapted to realistic parameters,
our model produces plausible results when compared to experimental and
simulation studies and we provide results from particle-based computer
simulations to further validate the analytical model. The proposed
comprehensive channel model admits a wide range of synaptic configurations
making it suitable for the investigation of many practically relevant
questions, such as the impact of glial cell uptake and spillover on signal
transmission in the tripartite synapse.


  The control of biological systems presents interesting applications such as
cell reprogramming or drug target identification. A common type of control
strategy consists in a set of interventions that, by fixing the values of some
variables, force the system to evolve to a desired state. This work presents a
new approach for finding control strategies in biological systems modeled by
Boolean networks. In this context, we explore the properties of trap spaces,
subspaces of the state space which the dynamics cannot leave. Trap spaces for
biological networks can often be efficiently computed, and provide useful
approximations of attraction basins. Our approach provides control strategies
for a target phenotype that are based on interventions that allow the control
to be eventually released. Moreover, our method can incorporate information
about the attractors to find new control strategies that would escape usual
percolation-based methods. We show the applicability of our approach to two
cell fate decision models.


  An epidemiological compartmental model was used to simulate social distancing
strategies to contain the COVID-19 pandemic and prevent a second wave in Sao
Paulo, Brazil. Optimization using genetic algorithm was used to determine the
optimal solutions. Our results suggest the best-case strategy for Sao Paulo is
to maintain or increase the current magnitude of social distancing for at least
60 more days and increase the current levels of personal protection behaviors
by a minimum of 10% (e.g., wearing facemasks, proper hand hygiene and avoid
agglomeration). Followed by a long-term oscillatory level of social distancing
with a stepping-down approach every 80 days over a period of two years with
continued protective behavior.


  This paper has a twofold contribution. The first is a data driven approach
for predicting the Covid 19 pandemic dynamics, based on data from more advanced
countries. The second is to report and discuss the results obtained with this
approach for Brazilian states, as of May 4th, 2020. We start by presenting
preliminary results obtained by training an LSTM SAE network, which are
somewhat disappointing. Then, our main approach consists in an initial
clustering of the world regions for which data is available and where the
pandemic is at an advanced stage, based on a set of manually engineered
features representing a country response to the early spread of the pandemic. A
Modified Auto-Encoder network is then trained from these clusters and learns to
predict future data for Brazilian states. These predictions are used to
estimate important statistics about the disease, such as peaks. Finally, curve
fitting is carried out on the predictions in order to find the distribution
that best fits the outputs of the MAE, and to refine the estimates of the peaks
of the pandemic. Results indicate that the pandemic is still growing in Brazil,
with most states peaks of infection estimated between the 25th of April and the
19th of May 2020. Predicted numbers reach a total of 240 thousand infected
Brazilians, distributed among the different states, with S\~ao Paulo leading
with almost 65 thousands estimated, confirmed cases. The estimated end of the
pandemics (with 97 percent of cases reaching an outcome) starts as of May 28th
for some states and rests through August 14th, 2020.


  In this thesis we develop minimal models of the relationship between
motility, growth, and evolution of cancer cells. We utilise simple simulations
of a population of individual cells in space to examine how changes in
mechanical properties of invasive cells and their surroundings can affect the
speed of cell migration. We also find that the growth rate of large lesions
depends weakly on the migration speed of escaping cells, and has stronger and
more complex dependencies on the rates of other stochastic processes in the
model, namely the rate at which cells transition to being motile and the
reverse rate at which cells cease to be motile.
  To examine how the rates of growth and evolution of an ensemble of cancerous
lesions depends on their geometry and underlying fitness landscape, we develop
an analytical framework in which the spatial structure is coarse grained and
the cancer treated as a continuously growing system with stochastic migration
events. Both approaches conclude that the whole ensemble can undergo
migration-driven exponential growth regardless of the dependence of size on
time of individual lesions, and that the relationship between growth rate and
rate of migration is determined by the geometrical constraints of individual
lesions. We also find that linear fitness landscapes result in
faster-than-exponential growth of the ensemble, and we can determine the
expected number of driver mutations present in several important cases of the
model.
  Finally, we study data from a clinical study of the effectiveness of a new
low-dose combined chemotherapy. This enables us to test some important
hypotheses about the growth rate of pancreatic cancers and the speed with which
evolution occurs in reality. Despite this, we find that the frequency of
resistant mutants is far too high to be explained without resorting to novel
mechanisms of cross-resistance to multiple drugs.


  Many animals live in societies where individuals frequently interact socially
with each other. The social structures of these systems can be studied in depth
by means of network analysis. A large number of studies on animal social
networks in many species have in recent years been carried out in the
biological research field of animal behaviour and have provided new insights
into behaviour, ecology, and social evolution. This line of research is
currently not so well connected to the field of complex systems as could be
expected. The purpose of this paper is to provide an introduction to animal
social networks for complex systems scientists and highlight areas of synergy.
We believe that an increased integration of animal social networks with the
interdisciplinary field of complex systems and networks would be beneficial for
various reasons. Increased collaboration between researchers in this field and
biologists studying animal social systems could be valuable in solving
challenges that are of importance to animal social network research.
Furthermore, animal social networks provide the opportunity to investigate
hypotheses about complex systems across a range of natural real-world social
systems. In this paper, we describe what animal social networks are and main
research themes where they are studied; we give an overview of the methods
commonly used to study animal social networks; we highlight challenges in the
study of animal social networks where complex systems expertise may be
particularly valuable; and we consider aspects of animal social networks that
may be of particular interest to complex systems researchers. We hope that this
will help to facilitate further interdisciplinary collaborations involving
animal social networks, and further integration of these networks into the
field of complex systems.


  PyRoss is an open-source Python library that offers an integrated platform
for inference, prediction and optimisation of NPIs in age- and
contact-structured epidemiological compartment models. This report outlines the
rationale and functionality of the PyRoss library, with various illustrations
and examples focusing on well-mixed, age-structured populations. The PyRoss
library supports arbitrary structured models formulated stochastically (as
master equations) or deterministically (as ODEs) and allows mid-run
transitioning from one to the other. By supporting additional compartmental
subdivision ad libitum, PyRoss can emulate time-since-infection models and
allows medical stages such as hospitalization or quarantine to be modelled and
forecast. The PyRoss library enables fitting to epidemiological data, as
available, using Bayesian parameter inference, so that competing models can be
weighed by their evidence. PyRoss allows fully Bayesian forecasts of the impact
of idealized NPIs by convolving uncertainties arising from epidemiological
data, model choice, parameters, and intrinsic stochasticity. Algorithms to
optimize time-dependent NPI scenarios against user-defined cost functions are
included. PyRoss's current age-structured compartment framework for well-mixed
populations will in future reports be extended to include compartments
structured by location, occupation, use of travel networks and other attributes
relevant to assessing disease spread and the impact of NPIs. We argue that such
compartment models, by allowing social data of arbitrary granularity to be
combined with Bayesian parameter estimation for poorly-known disease variables,
could enable more powerful and robust prediction than other approaches to
detailed epidemic modelling. We invite others to use the PyRoss library for
research to address today's COVID-19 crisis, and to plan for future pandemics.


  A number of models in mathematical epidemiology have been developed to
account for control measures such as vaccination or quarantine. However,
COVID-19 has brought unprecedented social distancing measures, with a challenge
on how to include these in a manner that can explain the data but avoid
overfitting in parameter inference. We here develop a simple time-dependent
model, where social distancing effects are introduced analogous to
coarse-grained models of gene expression control in systems biology. We apply
our approach to understand drastic differences in COVID-19 infection and
fatality counts, observed between Hubei (Wuhan) and other Mainland China
provinces. We find that these unintuitive data may be explained through an
interplay of differences in transmissibility, effective protection, and
detection efficiencies between Hubei and other provinces. More generally, our
results demonstrate that regional differences may drastically shape infection
outbursts. The obtained results demonstrate the applicability of our developed
method to extract key infection parameters directly from publically available
data so that it can be globally applied to outbreaks of COVID-19 in a number of
countries. Overall, we show that applications of uncommon strategies, such as
methods and approaches from molecular systems biology research to mathematical
epidemiology, may significantly advance our understanding of COVID-19 and other
infectious diseases.


  The COVID-19 pandemic has plagued the world for months. The U.S. has taken
measures to counter it. On a daily basis, newly confirmed cases have been
reported. In the early days, these numbers showed an increasing trend.
Recently, the numbers have been generally flattened out. This report tries to
estimate the hidden number of currently alive infections in the population by
using the confirmed cases. A major result indicates an existing infections
estimate at about 10-50 times the daily confirmed new cases, with the stringent
social distancing policy tipping to the upper end of this range. It clarifies
the relationship between the infection rate and the test rate to put the
epidemic under control, which says that the test rate shall keep up at the same
pace as infection rate to prevent an outbreak. This relationship is meaningful
in the wake of business re-opening in the U.S. and the world. The report also
reveals the connections of all the measures taken to the epidemic spread. A
stratified sampling method is proposed to add to the current tool kits of
epidemic control. Again, this report is a summary of some straight observations
and thoughts, not through a thorough study backed with field data. The results
appear obvious and suitable for general education to interested policymakers
and the public.


  The goal of this work is to study the optimal controls for the COVID-19
epidemic in Brazil. We consider an age-structured SEIRQ model with quarantine
compartment, where the controls are the quarantine entrance parameters. We then
compare the optimal controls for different quarantine lengths and distribution
of the total control cost by assessing their respective reductions in deaths in
comparison to the same period without quarantine. The best strategy provides a
calendar of when to relax the isolation measures for each age group. Finally,
we analyse how a delay in the beginning of the quarantine affects this calendar
by changing the initial conditions.


  The Coronavirus Disease (COVID-19) pandemic caused by the SARS-coronavirus 2
(SARS-CoV-2) urgently calls for the design of drugs directed against this new
virus. Given its essential role in proteolytic processing, the main protease
Mpro has been identified as an attractive candidate for drugs against
SARS-CoV-2 and similar coronaviruses. Recent high-throughput screening studies
have identified a set of existing, small-molecule drugs as potent Mpro
inhibitors. Amongst these, Ebselen (2-Phenyl-1,2-benzoselenazol-3-one), a
glutathione peroxidase mimetic seleno-organic compound, is particularly
attractive. Recent experiments suggest that its effectiveness is higher than
that of other molecules that also act at the enzyme catalytic site. By relying
on extensive simulations with all-atom models, in this study we examine at a
molecular level the potential of Ebselen to decrease Mpro catalytic activity.
Our results indicate that Ebselen exhibits a distinct affinity for the
catalytic site cavity of Mpro. In addition, our molecular models reveal a
second, previously unkown binding site for Ebselen in the dimerization region
localized between the II and III domains of the protein. A detailed analysis of
the free energy of binding indicates that the affinity of Ebselen to this
second binding site is in fact significantly larger than that to the catalytic
site. A strain analysis indicates that Ebselen bound between the II-III domains
exerts a pronounced allosteric effect that regulates catalytic site access
through surface loop interactions, and induces a displacement and
reconfiguration of water hotspots, including the catalytic water, which could
interfere with normal enzymatic function. Taken together, these findings
provide a framework for the future design of more potent and specific Mpro
inhibitors, based on the Ebselen scaffold, that could lead to new therapeutic
strategies for COVID-19.