  Neurotensin (NT) exerts naloxone-insensitive antinociceptive action through
its binding to both NTS1 and NTS2 receptors and NT analogs provide stronger
pain relief than morphine on a molecular basis. Here, we examined the
analgesic/adverse effect profile of a new NT(8-13) derivative denoted JMV2009,
in which the Pro10 residue was substituted by a silicon-containing unnatural
amino acid silaproline. We first report the synthesis and in vitro
characterization (receptor-binding affinity, functional activity and stability)
of JMV2009. We next examined its analgesic activity in a battery of acute,
tonic and chronic pain models. We finally evaluated its ability to induce
adverse effects associated with chronic opioid use, such as constipation and
analgesic tolerance or related to NTS1 activation, like hypothermia. In in
vitro assays, JMV2009 exhibited high binding affinity for both NTS1 and NTS2,
improved proteolytic resistance as well as agonistic activities similar to NT,
inducing sustained activation of p42/p44 MAPK and receptor internalization.
Intrathecal injection of JMV2009 produced dose-dependent antinociceptive
responses in the tail-flick test and almost completely abolished the
nociceptive-related behaviors induced by chemical somatic and visceral noxious
stimuli. Likewise, increasing doses of JMV2009 significantly reduced tactile
allodynia and weight bearing deficits in nerve-injured rats. Importantly,
chronic agonist treatment did not result in the development of analgesic
tolerance. Furthermore, JMV2009 did not cause constipation and was ineffective
in inducing hypothermia. These findings suggest that NT drugs can act as an
effective opioid-free medication for the management of pain or can serve as
adjuvant analgesics to reduce the opioid adverse effects.


  Some organism-specific databases about regulation in bacteria have become
larger, accelerated by high-throughput methodologies, while others are no
longer updated or accessible. Each database homogenize its datasets, giving
rise to heterogeneity across databases. Such heterogeneity mainly encompasses
different names for a gene and different network representations, generating
duplicated interactions that could bias network analyses. Abasy
(Across-bacteria systems) Atlas consolidates information from different sources
into meta-curated regulatory networks in bacteria. The high-quality networks in
Abasy Atlas enable cross-organisms analyses, such as benchmarking studies where
gold standards are required. Nevertheless, network incompleteness still casts
doubts on the conclusions of network analyses, and available sampling methods
cannot reflect the curation process. To tackle this problem, the updated
version of Abasy Atlas presented in this work provides historical snapshots of
regulatory networks. Thus, network analyses can be performed at different
completeness levels, making possible to identify potential bias and to predict
future results. We leverage the recently found constraint in the complexity of
regulatory networks to develop a novel model to quantify the total number of
regulatory interactions as a function of the genome size. This completeness
estimation is a valuable insight that may aid in the daunting task of network
curation, prediction, and validation. The new version of Abasy Atlas provides
76 networks (204,282 regulatory interactions) covering 42 bacteria (64%
Gram-positive and 36% Gram-negative) distributed in 9 species, containing 8,459
regulons and 4,335 modules.


  Due to modern transportation networks (airplanes, cruise ships, etc.) an
epidemic in a given country or city may be triggered by the arrival of external
infected agents. Posterior government quarantine policies are usually taken in
order to control the epidemic growth. We formulate a minimal epidemic evolution
model that takes into account these components. The previous and posterior
evolutions to the quarantine policy are modeled in a separate way by
considering different complexities parameters in each stage. Application of
this model to COVID-19 data in different countries is implemented. Estimations
of the infected peak time-occurrence and epidemic saturation values as well as
possible post-quarantine scenarios are analyzed over the basis of the model,
reported data, and the fraction of the population that adopts the quarantine
policy.


  Many cellular processes rely on the cell's ability to transport material to
and from the nucleus. Networks consisting of many microtubules and actin
filaments are key to this transport. Recently, the inhibition of intracellular
transport has been implicated in neurodegenerative diseases such as Alzheimer's
disease and Amyotrophic Lateral Sclerosis (ALS). Furthermore, microtubules may
contain so-called defective regions where motor protein velocity is reduced due
to accumulation of other motors and microtubule associated proteins. In this
work, we propose a new mathematical model describing the motion of motor
proteins on microtubules which incorporate a defective region. We take a
mean-field approach derived from a first principle lattice model to study motor
protein dynamics and density profiles. In particular, given a set of model
parameters we obtain a closed-form expression for the equilibrium density
profile along a given microtubule. We then verify the analytic results using
mathematical analysis on the discrete model and Monte Carlo simulations. This
work will contribute to the fundamental understanding of inhomogeneous
microtubules providing insight into microscopic interactions that may result in
the onset of neurodegenerative diseases. Our results for inhomogeneous
microtubules are consistent with prior work studying the homogeneous case.


  The number of new infections per day is a key quantity for effective epidemic
management. It can be estimated relatively directly by testing of random
population samples. Without such direct epidemiological measurement, other
approaches are required to infer whether the number of new cases is likely to
be increasing or decreasing: for example, estimating the pathogen effective
reproduction number, R, using data gathered from the clinical response to the
disease. For Covid-19 (SARS-CoV-2) such R estimation is heavily dependent on
modelling assumptions, because the available clinical case data are
opportunistic observational data subject to severe temporal confounding. Given
this difficulty it is useful to retrospectively reconstruct the time course of
infections from the least compromised available data, using minimal prior
assumptions. A Bayesian inverse problem approach applied to UK data on first
wave Covid-19 deaths and the disease duration distribution suggests that fatal
infections were in decline before full UK lockdown (24 March 2020), and that
fatal infections in Sweden started to decline only a day or two later. An
analysis of UK data using the model of Flaxman et al. (2020, Nature 584) gives
the same result under relaxation of its prior assumptions on R, suggesting an
enhanced role for non pharmaceutical interventions (NPIs) short of full lock
down in the UK context. Similar patterns appear to have occurred in the
subsequent two lockdowns. Estimates from the main UK Covid statistical
surveillance surveys, available since original publication, support these
results. Replication code for the paper is available in the supporting
information of doi/10.1111/biom.13462.


  This paper proposes a cluster-based method to analyze the evolution of
multivariate time series and applies this to the COVID-19 pandemic. On each
day, we partition countries into clusters according to both their case and
death counts. The total number of clusters and individual countries' cluster
memberships are algorithmically determined. We study the change in both
quantities over time, demonstrating a close similarity in the evolution of
cases and deaths. The changing number of clusters of the case counts precedes
that of the death counts by 32 days. On the other hand, there is an optimal
offset of 16 days with respect to the greatest consistency between cluster
groupings, determined by a new method of comparing affinity matrices. With this
offset in mind, we identify anomalous countries in the progression from
COVID-19 cases to deaths. This analysis can aid in highlighting the most and
least significant public policies in minimizing a country's COVID-19 mortality
rate.


  Alzheimer's disease (AD) is an age-specific neurodegenerative disease that
compromises cognitive functioning and impacts the quality of life of an
individual. Pathologically, AD is characterised by abnormal accumulation of
beta-amyloid (A$\beta$) and hyperphosphorylated tau protein. Despite research
advances over the last few decades, there is currently still no cure for AD.
Although, medications are available to control some behavioural symptoms and
slow the disease's progression, most prescribed medications are based on
cholinesterase inhibitors. Over the last decade, there has been increased
attention towards novel drugs, targeting alternative neurotransmitter pathways,
particularly those targeting serotonergic (5-HT) system. In this review, we
focused on 5-HT receptor (5-HTR) mediated signalling and drugs that target
these receptors. These pathways regulate key proteins and kinases such as GSK-3
that are associated with abnormal levels of A$\beta$ and tau in AD. We then
review computational studies related to 5-HT signalling pathways with the
potential for providing deeper understanding of AD pathologies. In particular,
we suggest that multiscale and multilevel modelling approaches could
potentially provide new insights into AD mechanisms, and towards discovering
novel 5-HTR based therapeutic targets.


  Effective, sensitive, and reliable diagnostic reagents are of paramount
importance for combating the ongoing coronavirus disease 2019 (COVID-19)
pandemic at a time there is no preventive vaccine nor specific drug available
for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It would be
an absolute tragedy if currently used diagnostic reagents are undermined in any
manner. Based on the genotyping of 7818 SARS-CoV-2 genome samples collected up
to May 1, 2020, we reveal that essentially all of the current COVID-19
diagnostic targets have had mutations. We further show that SARS-CoV-2 has the
most devastating mutations on the targets of various nucleocapsid (N) gene
primers and probes, which have been unfortunately used by countries around the
world to diagnose COVID-19. Our findings explain what has seriously gone wrong
with a specific diagnostic reagent made in China. To understand whether
SARS-CoV-2 genes have mutated unevenly, we have computed the mutation ratio and
mutation $h$-index of all SARS-CoV genes, indicating that the N gene is the
most non-conservative gene in the SARS-CoV-2 genome. Our findings enable
researchers to target the most conservative SARS-CoV-2 genes and proteins for
the design and development of COVID-19 diagnostic reagents, preventive
vaccines, and therapeutic medicines.


  The virus SARS-COV-2 caused disease COVID-19 has been declared a pandemic by
WHO. Currently, over 210 countries and territories have been affected. Careful,
well-designed drugs and vaccine for the total elimination of this virus seem to
be the need of the hour. In this context, the invivo mathematical modelling
studies can be extremely helpful in understanding the efficacy of the drug
interventions. These studies can also help understand the role of the crucial
inflammatory mediators and the behaviour of immune response towards this novel
coronavirus. Motivated by these facts, in this paper, we study the invivo
dynamics of Covid-19. The results obtained here are inline with some of the
clinical findings for Covid-19. This invivo modelling study involving the
crucial biomarkers of Covid-19 is the first of its kind and the results
obtained from this can be helpful to researchers, epidemiologists, clinicians
and doctors who are working in this field.


  The topology of gene expression space for a set of 12 cancer types is studied
by means of an entropy-like magnitude, which allows the characterization of the
regions occupied by tumor and normal samples. The comparison indicates that the
number of available states in gene expression space is much greater for tumors
than for normal tissues, suggesting the irreversibility of the progression to
the tumor phase. The entropy is nearly constant for tumors, whereas exhibits a
higher variability in normal tissues, probably due to tissue differentiation.
In addition, we show an interesting correlation between the fraction of
available states and the overlapping between the tumor and normal sample
clouds, interpreted as a way of reducing the decay rate to the tumor phase in
more ordered or structured tissues.


  SARS-CoV-2 is a severe respiratory infection that infects humans. Its
outburst entitled it as a pandemic emergence. To get a grip on this, outbreak
specific preventive and therapeutic interventions are urgently needed. It must
be said that, until now, there are no existing vaccines for coronaviruses. To
promptly and rapidly respond to pandemic events, the application of in silico
trials can be used for designing and testing medicines against SARS-CoV-2 and
speed-up the vaccine discovery pipeline, predicting any therapeutic failure and
minimizing undesired effects. Here, we present an in silico platform that
showed to be in very good agreement with the latest literature in predicting
SARS- CoV-2 dynamics and related immune system host response. Moreover, it has
been used to predict the outcome of one of the latest suggested approach to
design an effective vaccine, based on monoclonal antibody. UISS is then
potentially ready to be used as an in silico trial platform to predict the
outcome of vaccination strategy against SARS-CoV-2.


  Epidemiological models contain a set of parameters that must be adjusted
based on available observations. Once a model has been calibrated, it can be
used as a forecasting tool to make predictions and to evaluate contingency
plans. It is customary to employ only point estimators for such predictions.
However, some models may fit the same data reasonably well for a broad range of
parameter values, and this flexibility means that predictions stemming from
such models will vary widely, depending on the particular parameter values
employed within the range that give a good fit. When data are poor or
incomplete, model uncertainty widens further. A way to circumvent this problem
is to use Bayesian statistics to incorporate observations and use the full
range of parameter estimates contained in the posterior distribution to adjust
for uncertainties in model predictions. Specifically, given the epidemiological
model and a probability distribution for observations, we use the posterior
distribution of model parameters to generate all possible epidemiological
curves via the posterior predictive distribution. From the envelope of all
curves one can extract the worst-case scenario and study the impact of
implementing contingency plans according to this assessment. We apply this
approach to the potential evolution of COVID-19 in Mexico City and assess
whether contingency plans are being successful and whether the epidemiological
curve has flattened.


  The current outbreak is known as Coronavirus Disease or COVID-19 caused by
the virus SAR-COV-2 which continues to wreak havoc across the globe. The World
Health Organization (WHO) has declared the outbreak a Public Health Emergency
of International Concern. In Pakistan, the spread of the virus is on the rise
with the number of infected people and causalities rapidly increasing. In the
absence of proper vaccination and treatment, to reduce the number of infections
and casualties, the only option so far is to educate people regarding
preventive measures and to enforce countrywide lock-down. Any strategy about
the preventive measures needs to be based upon detailed analysis of the
COVID-19 outbreak and accurate scientific predictions. In this paper, we
conduct mathematical and numerical analysis to come up with reliable and
accurate predictions of the outbreak in Pakistan. The time-dependent
Susceptible-Infected-Recovered (SIR) model is used to fit the data and provide
future predictions. The turning point of the peak of the pandemic is defined as
the day when the transmission rate becomes less than the recovering rate. We
have predicted that the outbreak will reach its maximum peak occurring from
late May to 9 June with unrecovered number of Infectives in the range
20000-47000 and the cumulative number of infected cases in the range of
57500-153100. The number of Infectives will remain at the lower end in the
lock-down scenario but can rapidly double or triple if the spread of the
epidemic is not curtailed and localized. The uncertainty on single day
projection in our analysis after April 15 is found to be within 5\%.


  Discovering and isolating infected individuals is a cornerstone of epidemic
control. Because many infectious diseases spread through close contacts,
contact tracing is a key tool for case discovery and control. However, although
contact tracing has been performed widely, the mathematical understanding of
contact tracing has not been fully established and it has not been clearly
understood what determines the efficacy of contact tracing. Here, we reveal
that, compared with "forward" tracing---tracing to whom disease spreads,
"backward" tracing---tracing from whom disease spreads---is profoundly more
effective. The effectiveness of backward tracing is due to simple but
overlooked biases arising from the heterogeneity in contacts. Using simulations
on both synthetic and high-resolution empirical contact datasets, we show that
even at a small probability of detecting infected individuals, strategically
executed contact tracing can prevent a significant fraction of further
transmissions. We also show that---in terms of the number of prevented
transmissions per isolation---case isolation combined with a small amount of
contact tracing is more efficient than case isolation alone. By demonstrating
that backward contact tracing is highly effective at discovering
super-spreading events, we argue that the potential effectiveness of contact
tracing has been underestimated. Therefore, there is a critical need for
revisiting current contact tracing strategies so that they leverage all forms
of biases. Our results also have important consequences for digital contact
tracing because it will be crucial to incorporate the capability for backward
and deep tracing while adhering to the privacy-preserving requirements of these
new platforms.


  High throughput genome sequencing technologies such as RNA-Seq and Microarray
have the potential to transform clinical decision making and biomedical
research by enabling high-throughput measurements of the genome at a granular
level. However, to truly understand causes of disease and the effects of
medical interventions, this data must be integrated with phenotypic,
environmental, and behavioral data from individuals. Further, effective
knowledge discovery methods that can infer relationships between these data
types are required. In this work, we propose a pipeline for knowledge discovery
from integrated genomic and clinical data. The pipeline begins with a novel
variable selection method, and uses a probabilistic graphical model to
understand the relationships between features in the data. We demonstrate how
this pipeline can improve breast cancer outcome prediction models, and can
provide a biologically interpretable view of sequencing data.


  Computational drug design based on artificial intelligence is an emerging
research area. At the time of writing this paper, the world suffers from an
outbreak of the coronavirus SARS-CoV-2. A promising way to stop the virus
replication is via protease inhibition. We propose an evolutionary
multi-objective algorithm (EMOA) to design potential protease inhibitors for
SARS-CoV-2's main protease. Based on the SELFIES representation the EMOA
maximizes the binding of candidate ligands to the protein using the docking
tool QuickVina 2, while at the same time taking into account further objectives
like drug-likeliness or the fulfillment of filter constraints. The experimental
part analyzes the evolutionary process and discusses the inhibitor candidates.


  To model the evolution of diseases with extended latency periods and the
presence of asymptomatic patients like COVID-19, we define a simple discrete
time stochastic SIR-type epidemic model. We include both latent periods as well
as the presence of quarantine areas, to capture the evolutionary dynamics of
such diseases.


  Random effect models are popular statistical models for detecting and
correcting spurious sample correlations due to hidden confounders in
genome-wide gene expression data. In applications where some confounding
factors are known, estimating simultaneously the contribution of known and
latent variance components in random effect models is a challenge that has so
far relied on numerical gradient-based optimizers to maximize the likelihood
function. This is unsatisfactory because the resulting solution is poorly
characterized and the efficiency of the method may be suboptimal. Here we prove
analytically that maximum-likelihood latent variables can always be chosen
orthogonal to the known confounding factors, in other words, that
maximum-likelihood latent variables explain sample covariances not already
explained by known factors. Based on this result we propose a restricted
maximum-likelihood method which estimates the latent variables by maximizing
the likelihood on the restricted subspace orthogonal to the known confounding
factors, and show that this reduces to probabilistic PCA on that subspace. The
method then estimates the variance-covariance parameters by maximizing the
remaining terms in the likelihood function given the latent variables, using a
newly derived analytic solution for this problem. Compared to gradient-based
optimizers, our method attains greater or equal likelihood values, can be
computed using standard matrix operations, results in latent factors that don't
overlap with any known factors, and has a runtime reduced by several orders of
magnitude. Hence the restricted maximum-likelihood method facilitates the
application of random effect modelling strategies for learning latent variance
components to much larger gene expression datasets than possible with current
methods.


  In this work, we simulate the COVID-19 pandemic dynamics in a population
modeled as a network of groups wherein infection can propagate both via
intra-group and via inter-group interactions. Our results emphasize the
importance of diminishing the inter-group infections in the effort of
substantial flattening and delaying of the epi(demiologic) curve with
concomitant mitigation of disastrous economy and social consequences. To
exemplify with a limiting case, splitting a population into m (say, 5 or 10)
noninteracting groups while keeping intra-group interaction unchanged yields a
stretched epidemiologic curve having the maximum number of daily infections
reduced and postponed in time by the same factor $m$ (5 or 10). More generally,
our study suggests a practical approach to fight against SARS-CoV-2 virus
spread based on population splitting into groups and minimizing intermingling
between them. This strategy can be pursued by large-scale infrastructure
reorganization of activity at different levels in big logistic units (e.g.,
large productive networks, factories, enterprises, warehouses, schools,
(seasonal) harvest work). Importantly, unlike total lockdwon strategy, the
proposed approach prevents economic ruin and keeps social life at a more
bearable level than distancing everyone from anyone. Last but not least, the
declaration, for the first time in Europe, that COVID-19 epidemic ended in the
two million people Slovenia may be taken as confirming the advantage of the
strategy proposed in this paper.


  How does social distancing affect the reach of an epidemic in social
networks? We present Monte Carlo simulation results of a capacity constrained
Susceptible-Infected-Removed (SIR) model. The key modelling feature is that
individuals are limited in the number of acquaintances that they can interact
with, thereby constraining disease transmission to an infectious subnetwork of
the original social network. While increased social distancing always reduces
the spread of an infectious disease, the magnitude varies greatly depending on
the topology of the network. Our results also reveal the importance of
coordinating social distancing policies at the global level. In particular, the
public health benefits from social distancing to a group (e.g., a country) may
be completely undone if that group maintains connections with outside groups
that are not following suit.


  Most countries are suffering severely from the ongoing covid-19 pandemic
despite various levels of preventive measures. A common question is if and when
a country or region will reach herd immunity $h$. The classical herd immunity
level $h_C$ is defined as $h_C=1-1/R_0$, where $R_0$ is the basic reproduction
number, for covid-19 estimated to lie somewhere in the range 2.2-3.5 depending
on country and region. It is shown here that the disease-induced herd immunity
level $h_D$, after an outbreak has taken place in a country/region with a set
of preventive measures put in place, is actually substantially smaller than
$h_C$. As an illustration we show that if $R_0=2.5$ in an age-structured
community with mixing rates fitted to social activity studies, and also
categorizing individuals into three categories: low active, average active and
high active, and where preventive measures affect all mixing rates
proportionally, then the disease-induced herd immunity level is $h_D=43\%$
rather than $h_C=1-1/2.5=60\%$. Consequently, a lower fraction infected is
required for herd immunity to appear. The underlying reason is that when
immunity is induced by disease spreading, the proportion infected in groups
with high contact rates is greater than that in groups with low contact rates.
Consequently, disease-induced immunity is stronger than when immunity is
uniformly distributed in the community as in the classical herd immunity level.


  We consider a simplified model for gene regulation, where gene expression is
regulated by transcription factors (TFs), which are single proteins or protein
complexes. Proteins are in turn synthesised from expressed genes, creating a
feedback loop of regulation. This leads to a directed bipartite network in
which a link from a gene to a TF exists if the gene codes for a protein
contributing to the TF, and a link from a TF to a gene exists if the TF
regulates the expression of the gene. Both genes and TFs are modelled as binary
variables, which indicate, respectively, whether a gene is expressed or not,
and a TF is synthesised or not. We consider the scenario where for a TF to be
synthesised, all of its contributing genes must be expressed. This results in
an ``AND'' gate logic for the dynamics of TFs. By adapting percolation theory
to directed bipartite graphs, evolving according to the AND logic dynamics, we
are able to determine the necessary conditions, in the network parameter space,
under which bipartite networks can support a multiplicity of stable gene
expression patterns, under noisy conditions, as required in stable cell types.
  In particular, the analysis reveals the possibility of a bi-stability region,
where the extensive percolating cluster is or is not resilient to
perturbations. This is remarkably different from the transition observed in
standard percolation theory. Finally, we consider perturbations involving
single node removal that mimic gene knockout experiments. Results reveal the
strong dependence of the gene knockout cascade on the logic implemented in the
underlying network dynamics, highlighting in particular that avalanche sizes
cannot be easily related to gene-gene interaction networks.


  The current COVID-19 pandemic is having detrimental consequences worldwide.
The pandemic started to develop strongly by the end of January and beginning of
February 2020, first in China with subsequent rapid spread to other countries
with new epicenters of the outbreaks concentrated mainly within the 30-50
degrees North latitudinal band (e.g., South Korea, Japan, Iran, Italy, Spain).
Simultaneously, an unusual persistent anticyclonic situation prevailing at
latitudes around 40 degrees North was observed on global scale, in line with an
anomalously strong positive phase of the Arctic Oscillation. This atypical
situation could have resulted in favorable meteorological conditions for a
quicker spread of the virus over the latitude band detailed above. This
possible connection needs further attention in order to understand the
meteorological and climatological factors related to the COVID-19 outbreak, and
for anticipating the spatio-temporal distribution of possible future pandemics.


  Given the stochastic nature of gene expression, genetically identical cells
exposed to the same environmental inputs will produce different outputs. This
heterogeneity has been hypothesized to have consequences for how cells are able
to survive in changing environments. Recent work has explored the use of
information theory as a framework to understand the accuracy with which cells
can ascertain the state of their surroundings. Yet the predictive power of
these approaches is limited and has not been rigorously tested using precision
measurements. To that end, we generate a minimal model for a simple genetic
circuit in which all parameter values for the model come from independently
published data sets. We then predict the information processing capacity of the
genetic circuit for a suite of biophysical parameters such as protein copy
number and protein-DNA affinity. We compare these parameter-free predictions
with an experimental determination of protein expression distributions and the
resulting information processing capacity of E. coli cells. We find that our
minimal model captures the scaling of the cell-to-cell variability in the data
and the inferred information processing capacity of our simple genetic circuit
up to a systematic deviation.


  In the recent years, therapeutic use of antibodies has seen a huge growth,
due to their inherent proprieties and technological advances in the methods
used to study and characterize them. Effective design and engineering of
antibodies for therapeutic purposes are heavily dependent on knowledge of the
structural principles that regulate antibody-antigen interactions. Several
experimental techniques such as X-ray crystallography, cryo-electron
microscopy, NMR or mutagenesis analysis can be applied, but these are usually
expensive and time consuming. Therefore computational approaches like molecular
docking may offer a valuable alternative for the characterisation of
antibody-antigen complexes.
  Here we describe a protocol for the prediction of the 3D structure of
antibody-antigen complexes using the integrative modelling platform HADDOCK.
The protocol consists of: 1) The identification of the antibody residues
belonging to the hyper variable loops which are known to be crucial for the
binding and can be used to guide the docking; 2) The detailed steps to perform
docking with the HADDOCK 2.4 webserver following different strategies depending
on the availability of information about epitope residues.


  I employ a simple mathematical model of an epidemic process to evaluate how
four basic quantities: the reproduction number (R), the numbers of sensitive
(S) and infectious individuals(I), and total community size (N) affect
strategies to control COVID-19. Numerical simulations show that strict
suppression measures at the beginning of an epidemic can create low infectious
numbers, which thereafter can be managed by mitigation measures over longer
periods to flatten the epidemic curve. The stronger the suppression measure,
the faster it achieves the low numbers of infections that are conducive to
subsequent management. We discuss the predictions of this analysis and how it
fits into longer-term sequences of measures, including using the herd immunity
concept to leverage acquired immunity.


  Since the beginning of March 2020, the cumulative numbers of cases of
infection with the novel coronavirus SARS-CoV-2 in Germany have been reported
on a daily basis. The reports originate from national laws, according to which
positive test findings must be submitted to the Federal Health Authorities, the
Robert Koch Institute, via the local health authorities. Since an enormous
number of unreported cases can be expected, the question of how widespread the
disease has been in the population cannot be answered based on these
administrative reports. Using mathematical modeling, however, estimates can be
made. These estimates indicate that the small numbers of diagnostic tests
carried out at the beginning of the outbreak overlooked considerable parts of
the infection. In order to cover the initial phase of future waves of the
disease, wide-spread and comprehensive tests are recommended.


  Classical studies on aspiration-based dynamics suggest that a dissatisfied
individual changes strategy without taking into account the success of others.
This promotes defection spreading. The imitation-based dynamics allow
individuals to imitate successful strategies without taking into account their
own-satisfactions. In this article, we propose to study a dynamic based on
aspiration which takes into account imitation of successful strategies for
dissatisfied individuals. This helps cooperative members to resist. Individuals
compare their success to their desired satisfaction level before making a
decision to update their strategies. This mechanism helps individuals with a
minimum of self-satisfaction to maintain their strategies. If an individual is
dissatisfied, it will learn from others by choosing successful strategies. We
derive an exact expression of the fixation probability in well-mixed
populations as in structured populations in networks. As a result, we show that
selection may favor cooperation more than defection in well-mixed populations
as in populations ranged over a regular graph. We show that the best scenario
is a graph with small connectivity.


  In this paper, we use an adaptive modeling framework to model and study how
nutritional status (measured by the protein to carbohydrate ratio) may regulate
population dynamics and foraging task allocation of social insect colonies.
Mathematical analysis of our model shows that both investment to brood rearing
and brood nutrition are important for colony survival and dynamics. When
division of labor and/or nutrition are in an intermediate value range, the
model undergoes a backward bifurcation and creates multiple attractors due to
bistability. This bistability implies that there is a threshold population size
required for colony survival. When the investment in brood is large enough or
nutritional requirements are less strict the colony tends to survive, otherwise
the colony faces collapse. Our model suggests that the needs of colony survival
are shaped by the brood survival probability, which requires good nutritional
status. As a consequence, better nutritional status can lead to a better
survival rate of larvae, and thus a larger worker population.


  In a previous article [1] we have described the temporal evolution of the
Sars- Cov-2 in Italy in the time window February 24-April 1. As we can see in
[1] a generalized logistic equation captures both the peaks of the total
infected and the deaths. In this article our goal is to study the missing peak,
i.e. the currently infected one (or total currently positive). After the April
7 the large increase in the number of swabs meant that the logistical behavior
of the infected curve no longer worked. So we decided to generalize the model,
introducing new parameters. Moreover, we adopt a similar approach used in [1]
(for the estimation of deaths) in order to evaluate the recoveries. In this
way, introducing a simple conservation law, we define a model with 4
populations: total infected, currently positives, recoveries and deaths.
Therefore, we propose an alternative method to a classical SIRD model for the
evaluation of the Sars-Cov-2 epidemic. However, the method is general and thus
applicable to other diseases. Finally we study the behavior of the ratio
infected over swabs for Italy, Germany and USA, and we show as studying this
parameter we recover the generalized Logistic model used in [1] for these three
countries. We think that this trend could be useful for a future epidemic of
this coronavirus.


  Since the very first days of announcing Covid 19 as a global pandemic,
researchers around the world from different backgrounds and organizations
including national governments and the World Health Organization (WHO) have
extensively begun using Modeling and Simulation (M&S) techniques in deciding on
the optimal strategies in order to mitigate the public health and economical
effects of Covid-19. In this paper, we present our results from applying M&S on
the city of Casablanca, Morocco to comprehend (i): How it can be relatively
easy or difficult for such pandemic to penetrate to the studied region
(country, city, district, etc.). (ii): How it can spread relatively fast or
slow within that region. And finally (iii): What is the optimal strategy for
containing and mitigating the propagation of the disease while still keeping
the economy as near to normal as possible.


  A reaction-diffusion model was developed describing the spread of the
COVID-19 virus considering the mean daily movement of susceptible, exposed and
asymptomatic individuals. The model was calibrated using data on the confirmed
infection and death from France as well as their initial spatial distribution.
First, the system of partial differential equations is studied, then the basic
reproduction number, R0 is derived. Second, numerical simulations, based on a
combination of level-set and finite differences, shown the spatial spread of
COVID-19 from March 16 to June 16. Finally, scenarios of unlockdown are
compared according to variation of distancing, or partially spatial lockdown.


  We here propose to model active and cumulative cases data from COVID-19 by a
continuous effective model based on a modified diffusion equation under
Lifshitz scaling with a dynamic diffusion coefficient. The proposed model is
rich enough to capture different aspects of a complex virus diffusion as
humanity has been recently facing. The model being continuous it is bound to be
solved analytically and/or numerically. So, we investigate two possible models
where the diffusion coefficient associated with possible types of contamination
are captured by some specific profiles. The active cases curves here derived
were able to successfully describe the pandemic behavior of Germany and Spain.
Moreover, we also predict some scenarios for the evolution of COVID-19 in
Brazil. Furthermore, we depicted the cumulative cases curves of COVID-19,
reproducing the spreading of the pandemic between the cities of S\~ao Paulo and
S\~ao Jos\'e dos Campos, Brazil. The scenarios also unveil how the lockdown
measures can flatten the contamination curves. We can find the best profile of
the diffusion coefficient that better fit the real data of pandemic.


  We investigate adaptive strategies to robustly and optimally control the
COVID-19 pandemic via social distancing measures based on the example of
Germany. Our goal is to minimize the number of fatalities over the course of
two years without inducing excessive social costs. We consider a tailored model
of the German COVID-19 outbreak with different parameter sets to design and
validate our approach. Our analysis reveals that an open-loop optimal control
policy can significantly decrease the number of fatalities when compared to
simpler policies under the assumption of exact model knowledge. In a more
realistic scenario with uncertain data and model mismatch, a feedback strategy
that updates the policy weekly using model predictive control (MPC) leads to a
reliable performance, even when applied to a validation model with deviant
parameters. On top of that, we propose a robust MPC-based feedback policy using
interval arithmetic that adapts the social distancing measures cautiously and
safely, thus leading to a minimum number of fatalities even if measurements are
inaccurate and the infection rates cannot be precisely specified by social
distancing. Our theoretical findings support various recent studies by showing
that 1) adaptive feedback strategies are required to reliably contain the
COVID-19 outbreak, 2) well-designed policies can significantly reduce the
number of fatalities compared to simpler ones while keeping the amount of
social distancing measures on the same level, and 3) imposing stronger social
distancing measures early on is more effective and cheaper in the long run than
opening up too soon and restoring stricter measures at a later time.


  We analyze risk factors correlated with the initial transmission growth rate
of the recent COVID-19 pandemic in different countries. The number of cases
follows in its early stages an almost exponential expansion; we chose as a
starting point in each country the first day $d_i$ with 30 cases and we fitted
for 12 days, capturing thus the early exponential growth. We looked then for
linear correlations of the exponents $\alpha$ with other variables, for a
sample of 126 countries. We find a positive correlation, {\it i.e. faster
spread of COVID-19}, with high confidence level with the following variables,
with respective $p$-value: low Temperature ($4\cdot10^{-7}$), high ratio of old
vs.~working-age people ($3\cdot10^{-6}$), life expectancy ($8\cdot10^{-6}$),
number of international tourists ($1\cdot10^{-5}$), earlier epidemic starting
date $d_i$ ($2\cdot10^{-5}$), high level of physical contact in greeting habits
($6 \cdot 10^{-5}$), lung cancer prevalence ($6 \cdot 10^{-5}$), obesity in
males ($1 \cdot 10^{-4}$), share of population in urban areas
($2\cdot10^{-4}$), cancer prevalence ($3 \cdot 10^{-4}$), alcohol consumption
($0.0019$), daily smoking prevalence ($0.0036$), UV index ($0.004$, 73
countries). We also find a correlation with low Vitamin D levels
($0.002-0.006$, smaller sample, $\sim 50$ countries, to be confirmed on a
larger sample). There is highly significant correlation also with blood types:
positive correlation with types RH- ($3\cdot10^{-5}$) and A+ ($3\cdot10^{-3}$),
negative correlation with B+ ($2\cdot10^{-4}$). Several of the above variables
are intercorrelated and likely to have common interpretations. We performed a
Principal Component Analysis, in order to find their significant independent
linear combinations. We also analyzed a possible bias: countries with low
GDP-per capita might have less testing and we discuss correlation with the
above variables.


  This study focuses on the task of supervised prediction of aging-related
genes from -omics data. Unlike gene expression methods for this task that
capture aging-specific information but ignore interactions between genes (i.e.,
their protein products), or protein-protein interaction (PPI) network methods
for this task that account for PPIs but the PPIs are context-unspecific, we
recently integrated the two data types into an aging-specific PPI subnetwork,
which yielded more accurate aging-related gene predictions. However, a dynamic
aging-specific subnetwork did not improve prediction performance compared to a
static aging-specific subnetwork, despite the aging process being dynamic. This
could be because the dynamic subnetwork was inferred using a naive Induced
subgraph approach. Instead, we recently inferred a dynamic aging-specific
subnetwork using a methodologically more advanced notion of network propagation
(NP), which improved upon Induced dynamic aging-specific subnetwork in a
different task, that of unsupervised analyses of the aging process. Here, we
evaluate whether our existing NP-based dynamic subnetwork will improve upon the
dynamic as well as static subnetwork constructed by the Induced approach in the
considered task of supervised prediction of aging-related genes. The existing
NP-based subnetwork is unweighted, i.e., it gives equal importance to each of
the aging-specific PPIs. Because accounting for aging-specific edge weights
might be important, we additionally propose a weighted NP-based dynamic
aging-specific subnetwork. We demonstrate that a predictive machine learning
model trained and tested on the weighted subnetwork yields higher accuracy when
predicting aging-related genes than predictive models run on the existing
unweighted dynamic or static subnetworks, regardless of whether the existing
subnetworks were inferred using NP or the Induced approach.


  Multiple species in the ecosystem are believed to compete cyclically for
survival and thus maintain balance in nature. Stochasticity has also an
inevitable role in this dynamics. Considering these attributes of nature, the
stochastic dynamics of the rock-paper-scissor model based on the idea of cyclic
dominance becomes an effective tool to capture different aspects of ecosystem.
The evolutionary dynamics of this model crucially depends on different
interactions representing different natural habits. In this framework we
explore the role of mortality of individual organism in the collective survival
of a species. For this purpose a new parameter called `natural death' is
introduced. It is meant for bringing about the decease of an individual
irrespective of any intra- and interspecific interaction. We perform Monte
Carlo simulation followed by the stability analysis of different fixed points
of defined rate equations and observe that the natural death rate is
surprisingly one of the most significant factors in deciding whether an
ecosystem would come up with a coexistence or a single species survival.


  We introduce a new algorithm called {\sc Rec-Gen} for reconstructing the
genealogy or \textit{pedigree} of an extant population purely from its genetic
data. We justify our approach by giving a mathematical proof of the
effectiveness of {\sc Rec-Gen} when applied to pedigrees from an idealized
generative model that replicates some of the features of real-world pedigrees.
Our algorithm is iterative and provides an accurate reconstruction of a large
fraction of the pedigree while having relatively low \emph{sample complexity},
measured in terms of the length of the genetic sequences of the population. We
propose our approach as a prototype for further investigation of the pedigree
reconstruction problem toward the goal of applications to real-world examples.
As such, our results have some conceptual bearing on the increasingly important
issue of genomic privacy.


  We address a novel approach for stochastic individual-based modelling of a
single species population. Individuals are distinguished by their remaining
lifetimes, which are regulated by the interplay between the inexorable running
of time and the individual's nourishment history. A food-limited environment
induces intraspecific competition and henceforth the carrying capacity of the
medium may be finite, often emulating the qualitative features of logistic
growth. Inherently non-logistic behavior is also obtained by suitable change of
the few parameters involved, composing a wide variety of dynamical features.
Some analytical results are obtained. Beyond the rich phenomenology observed,
we expect that possible modifications of our model may account for an even
broader scope of collective population growth phenomena.


  Cell type identification from single-cell transcriptomic data is a common
goal of single-cell RNA sequencing (scRNAseq) data analysis. Neural networks
have been employed to identify cell types from scRNAseq data with high
performance. However, it requires a large mount of individual cells with
accurate and unbiased annotated types to build the identification models.
Unfortunately, labeling the scRNAseq data is cumbersome and time-consuming as
it involves manual inspection of marker genes. To overcome this challenge, we
propose a semi-supervised learning model to use unlabeled scRNAseq cells and
limited amount of labeled scRNAseq cells to implement cell identification.
Firstly, we transform the scRNAseq cells to "gene sentences", which is inspired
by similarities between natural language system and gene system. Then genes in
these sentences are represented as gene embeddings to reduce data sparsity.
With these embeddings, we implement a semi-supervised learning model based on
recurrent convolutional neural networks (RCNN), which includes a shared
network, a supervised network and an unsupervised network. The proposed model
is evaluated on macosko2015, a large scale single-cell transcriptomic dataset
with ground truth of individual cell types. It is observed that the proposed
model is able to achieve encouraging performance by learning on very limited
amount of labeled scRNAseq cells together with a large number of unlabeled
scRNAseq cells.


  The demographic factors have a substantial impact on the overall casualties
caused by the COVID-19. In this study, the spatial association between the key
demographic variables and COVID-19 cases and deaths were analyzed using the
spatial regression models. Total 13 (for COVID-19 case factor) and 8 (for
COVID-19 death factor) key variables were considered for the modelling. Total
five spatial regression models such as Geographically weighted regression
(GWR), Spatial Error Model (SEM), Spatial Lag Model (SLM), Spatial Error_Lag
model (SEM_SLM), and Ordinary Least Square (OLS) were performed for the spatial
modelling and mapping of model estimates. The local R2 values, which suggesting
the influences of the selected demographic variables on overall casualties
caused by COVID-19, was found highest in Italy and the UK. The moderate local
R2 was observed for France, Belgium, Netherlands, Ireland, Denmark, Norway,
Sweden, Poland, Slovakia, and Romania. The lowest local R2 value for COVID-19
cases was accounted for Latvia and Lithuania. Among the 13 variables, the
highest local R2 was calculated for total population (R2 = 0.92), followed by
death crude death rate (R2 = 0.9), long time illness (R2 = 0.84), population
with age >80 (R2 = 0.59), employment (R2 = 0.46), life expectancy at 65 (R2 =
0.34), crude birth rate (R2 = 0.31), life expectancy (R2 = 0.31), Population
with age 65-80 (R2 = 0.29), Population with age 15-24 (R2 = 0.27), Population
with age 25-49 (R2 = 0.27), Population with age 0-14 (R2 = 0.23), and
Population with age 50-65 (R2 = 0.23), respectively.


  Covid-19 has caused hundred of thousands of deaths and an economic damage
amounting to trillions of dollars, creating a desire for the rapid development
of vaccine. Once available, vaccine is gradually produced, evoking the question
on how to distribute it best. While official vaccination guidelines largely
focus on the question to whom vaccines should be provided first (e.g. to risk
groups), here we propose a strategy for their distribution in time and space,
which sequentially prioritizes regions with a high local infection growth rate.
To demonstrate this strategy, we develop a simple statistical model describing
the time-evolution of infection patterns and their response to vaccination, for
infectious diseases like Covid-19. For inhomogeneous infection patterns,
locally well-mixed populations and basic reproduction numbers $R_0\sim 1.5-4$
the proposed strategy at least halves the number of deaths in our simulations
compared to the standard practice of distributing vaccines proportionally to
the population density. For $R_0\sim 1$ we still find a significant increase of
the survival rate. The proposed vaccine distribution strategy can be further
tested in detailed modelling works and could excite discussions on the
importance of the spatiotemporal distribution of vaccines for official
guidelines.


  Hepatocellular carcinoma (HCC) is the most common type of primary liver
cancer and the fourth most common cause of cancer-related death worldwide.
Understanding the underlying gene mutations in HCC provides great prognostic
value for treatment planning and targeted therapy. Radiogenomics has revealed
an association between non-invasive imaging features and molecular genomics.
However, imaging feature identification is laborious and error-prone. In this
paper, we propose an end-to-end deep learning framework for mutation prediction
in APOB, COL11A1 and ATRX genes using multiphasic CT scans. Considering
intra-tumour heterogeneity (ITH) in HCC, multi-region sampling technology is
implemented to generate the dataset for experiments. Experimental results
demonstrate the effectiveness of the proposed model.


  Gene mutation prediction in hepatocellular carcinoma (HCC) is of great
diagnostic and prognostic value for personalized treatments and precision
medicine. In this paper, we tackle this problem with multi-instance multi-label
learning to address the difficulties on label correlations, label
representations, etc. Furthermore, an effective oversampling strategy is
applied for data imbalance. Experimental results have shown the superiority of
the proposed approach.


  Antimicrobial peptides (AMPs) are anti-infectives that have potential as a
novel and untapped class of biotherapeutics. Modes of action of antimicrobial
peptides imply interaction with cell envelope. Comprehensive understanding of
peculiarities of interactions of antimicrobial peptides with cell envelope is
necessary to perform the task-oriented design of new biotherapeutics, against
which for microbes it is hard to work out resistance. In order to enable a de
novo design with low costs and in high throughput, in silico predictive models
have to be required. To develop the performant predictive model, comprehensive
knowledge on mechanisms of action of AMPs has to be possessed. The last
knowledge will allow us to encode amino acid sequences expressively and to get
success to the choosing of the accurate classifier of AMPs. A shared protective
layer of microbial cells is inner, plasmatic membrane. The interaction of AMP
with a biological membrane (native and/or artificial) is the most
comprehensively studied. We provide a review of mechanisms and results of
interaction of AMP with the cell membrane, relying on the survey of
physicochemical, aggregative and structural features of AMPs. Potency and
mechanism of action of AMP have presented in the terms of amino acid
compositions and distributions of the polar and apolar residues along the
chain, that is in such physicochemical features of peptides as the
hydrophobicity, hydrophilicity, and amphiphilicity. Many different approaches
were used to classify AMPs. The survey of the knowledge on sequences,
structures, and modes of actions of AMP, allows concluding that, only the
physicochemical features of AMPs give the capability to perform the unambiguous
classification. Comprehensive knowledge of physicochemical features of AMP is
necessary to develop task-oriented methods of design of peptide-based
antibiotics de novo.


  COVID-19 has created a global pandemic with high morbidity and mortality in
2020. Novel coronavirus (nCoV), also known as Severe Acute Respiratory Syndrome
Coronavirus 2 (SARS-CoV2), is responsible for this deadly disease.
International Committee on Taxonomy of Viruses (ICTV) has declared that nCoV is
highly genetically similar to SARS-CoV epidemic in 2003 (89% similarity).
Limited number of clinically validated Human-nCoV protein interaction data is
available in the literature. With this hypothesis, the present work focuses on
developing a computational model for nCoV-Human protein interaction network,
using the experimentally validated SARS-CoV-Human protein interactions.
Initially, level-1 and level-2 human spreader proteins are identified in
SARS-CoV-Human interaction network, using Susceptible-Infected-Susceptible
(SIS) model. These proteins are considered as potential human targets for nCoV
bait proteins. A gene-ontology based fuzzy affinity function has been used to
construct the nCoV-Human protein interaction network at 99.98% specificity
threshold. This also identifies the level-1 human spreaders for COVID-19 in
human protein-interaction network. Level-2 human spreaders are subsequently
identified using the SIS model. The derived host-pathogen interaction network
is finally validated using 7 potential FDA listed drugs for COVID-19 with
significant overlap between the known drug target proteins and the identified
spreader proteins.


  We report the performance of our newly introduced Ensemble Docking with
Enhanced sampling of pocket Shape (EDES) protocol coupled to a template-based
algorithm to generate near-native ligand conformations in the 2019 iteration of
the Grand Challenge organized by the D3R consortium. Using either AutoDock4.2
or HADDOCK2.2 docking programs (each software in two variants of the protocol)
our method generated native-like poses among the top 5 submitted for evaluation
for most of the 20 targets with similar performances. The protein selected for
GC4 was the human beta-site amyloid precursor protein cleaving enzyme 1
(BACE-1), a transmembrane aspartic-acid protease. We identified at least one
pose whose heavy-atoms RMSD was less than 2.5 {\AA} from the native
conformation for 16 (80%) and 17 (85%) of the twenty targets using AutoDock and
HADDOCK, respectively. Dissecting the possible sources of errors revealed that:
i) our EDES protocol (with minor modifications) was able to sample
sub-{\aa}ngstrom conformations for all 20 protein targets, reproducing the
correct conformation of the binding site within ~1 {\AA} RMSD; ii) as already
shown by some of us in GC3, even in the presence of near-native protein
structures, a proper selection of ligand conformers is crucial for the success
of ensemble-docking calculations. Importantly, our approach performed best
among the protocols exploiting only structural information of the apo protein
to generate conformations of the receptor for ensemble-docking calculations.


  The sudden outbreak of the Coronavirus disease (COVID-19) swept across the
world in early 2020, triggering the lockdowns of several billion people across
many countries, including China, Spain, India, the U.K., Italy, France,
Germany, and most states of the U.S. The transmission of the virus accelerated
rapidly with the most confirmed cases in the U.S., and New York City became an
epicenter of the pandemic by the end of March. In response to this national and
global emergency, the NSF Spatiotemporal Innovation Center brought together a
taskforce of international researchers and assembled implemented strategies to
rapidly respond to this crisis, for supporting research, saving lives, and
protecting the health of global citizens. This perspective paper presents our
collective view on the global health emergency and our effort in collecting,
analyzing, and sharing relevant data on global policy and government responses,
geospatial indicators of the outbreak and evolving forecasts; in developing
research capabilities and mitigation measures with global scientists, promoting
collaborative research on outbreak dynamics, and reflecting on the dynamic
responses from human societies.


  In this article, we consider a dynamic epidemiology model for the spread of
the COVID-19 infection. Starting from the classical SEIR model, the model is
modified so as to better describe characteristic features of the underlying
pathogen and its infectious modes. In line with the large number of secondary
infections not related to contact with documented infectious individuals, the
model includes a cohort of asymptomatic or oligosymptomatic infectious
individuals, not accounted for in the data of new daily counts of infections. A
Bayesian particle filtering algorithm is used to update dynamically the
relevant cohort and simultaneously estimate the transmission rate as the new
data on the number of new infections and disease related death become
available. The underlying assumption of the model is that the infectivity rate
is dynamically changing during the epidemics, either because of a mutation of
the pathogen or in response to mitigation and containment measures. The
sequential Bayesian framework naturally provides a quantification of the
uncertainty in the estimate of the model parameters, including the reproduction
number, and of the size of the different cohorts. Moreover, we introduce a
dimensionless quantity, which is the equilibrium ratio between asymptomatic and
symptomatic cohort sizes, and propose a simple formula to estimate the
quantity. This ratio leads naturally to another dimensionless quantity that
plays the role of the basic reproduction number $R_0$ of the model. When we
apply the model and particle filter algorithm to COVID-19 infection data from
several counties in Northeastern Ohio and Southeastern Michigan we found the
proposed reproduction number $R_0$ to have a consistent dynamic behavior within
both states, thus proving to be a reliable summary of the success of the
mitigation measures.


  A three dimensional mesoscopic model is applied to study the properties of
short DNA chains in a confining environment. The cylindrical channel is
represented by a hard-wall repulsive potential incorporated in the system
Hamiltonian. The macroscopic helical parameters are computed performing
statistical averages over the ensemble of microscopic base pair fluctuations.
The average molecule elongation, measured by the end-to-end distance, is
derived as a function of the channel potential parameters both for a
homogeneous and a heterogeneous chain. The overall results suggest that the
mesoscopic model, with the channel potential term, yields consistent
quantitative estimates for the stretching and twisting of short chains.


  Numerous living systems are hierarchically organised, whereby replicating
components are grouped into reproducing collectives -- e.g., organelles are
grouped into cells, and cells are grouped into multicellular organisms. In such
systems, evolution can operate at two levels: evolution among collectives,
which tends to promote selfless cooperation among components within collectives
(called altruism), and evolution within collectives, which tends to promote
cheating among components within collectives. The balance between within- and
among-collective evolution thus exerts profound impacts on the fitness of these
systems. Here, we investigate how this balance depends on the size of a
collective (denoted by $N$) and the mutation rate of components ($m$) through
mathematical analyses and computer simulations of multiple population genetics
models. We first confirm a previous result that increasing $N$ or $m$
accelerates within-collective evolution relative to among-collective evolution,
thus promoting the evolution of cheating. Moreover, we show that when within-
and among-collective evolution exactly balance each other out, the following
scaling relation generally holds: $Nm^{\alpha}$ is a constant, where scaling
exponent $\alpha$ depends on multiple parameters, such as the strength of
selection and whether altruism is a binary or quantitative trait. This relation
indicates that although $N$ and $m$ have quantitatively distinct impacts on the
balance between within- and among-collective evolution, their impacts become
identical if $m$ is scaled with a proper exponent. Our results thus provide a
novel insight into conditions under which cheating or altruism evolves in
hierarchically-organised replicating systems.


  Nationwide lockdown for COVID-19 created an urgent demand for public
transportation among migrant workers stranded at different parts of India to
return to their native places. Arranging transportation could spike the number
of COVID-19 infected cases. Hence, this paper investigates the potential surge
in confirmed and active cases of COVID-19 infection and assesses the train and
bus fleet size required for the repatriating migrant workers. The expected to
repatriate migrant worker population was obtained by forecasting the 2011
census data and comparing it with the information reported in the news media. A
modified susceptible-exposed-infected-removed (SEIR) model was proposed to
estimate the surge in confirmed and active cases of COVID-19 patients in
India's selected states with high outflux of migrants. The developed model
considered combinations of different levels of the daily arrival rate of
migrant workers, total migrant workers in need of transportation, and the
origin of the trip dependent symptomatic cases on arrival. Reducing the daily
arrival rate of migrant workers for states with very high outflux of migrants
(i.e., Uttar Pradesh and Bihar) can help to lower the surge in confirmed and
active cases. Nevertheless, it could create a disparity in the number of days
needed to transport all repatriating migrant workers to the home states. Hence,
travel arrangements for about 100,000 migrant workers per day to Uttar Pradesh
and Bihar, about 50,000 per day to Rajasthan and Madhya Pradesh, 20,000 per day
to Maharashtra and less than 20,000 per day to other states of India was
recommended.


  A major difficulty in the analysis of propagation of the coronavirus is that
many infected individuals show no symptoms of Covid-19. This implies a lack of
information on the total counts of infected individuals and of recovered and
immunized individuals. In this paper, we consider parametric time varying
Markov processes of Coronavirus propagation and show how to estimate the model
parameters and approximate the unobserved counts from daily numbers of infected
and detected individuals and total daily death counts. This model-based
approach is illustrated in an application to French data.


  We examine the spatial modeling of the outbreak of COVID-19 in two regions:
the autonomous community of Andalusia in Spain and the mainland of Greece. We
start with a 0D compartmental epidemiological model consisting of Susceptible,
Exposed, Asymptomatic, (symptomatically) Infected, Hospitalized, Recovered, and
deceased populations. We emphasize the importance of the viral latent period
and the key role of an asymptomatic population. We optimize model parameters
for both regions by comparing predictions to the cumulative number of infected
and total number of deaths via minimizing the $\ell^2$ norm of the difference
between predictions and observed data. We consider the sensitivity of model
predictions on reasonable variations of model parameters and initial
conditions, addressing issues of parameter identifiability. We model both
pre-quarantine and post-quarantine evolution of the epidemic by a
time-dependent change of the viral transmission rates that arises in response
to containment measures. Subsequently, a spatially distributed version of the
0D model in the form of reaction-diffusion equations is developed. We consider
that, after an initial localized seeding of the infection, its spread is
governed by the diffusion (and 0D model "reactions") of the asymptomatic and
symptomatically infected populations, which decrease with the imposed
restrictive measures. We inserted the maps of the two regions, and we imported
population-density data into COMSOL, which was subsequently used to solve
numerically the model PDEs. Upon discussing how to adapt the 0D model to this
spatial setting, we show that these models bear significant potential towards
capturing both the well-mixed, 0D description and the spatial expansion of the
pandemic in the two regions. Veins of potential refinement of the model
assumptions towards future work are also explored.


  In this study, we develop the mathematical model to understand the coupling
between the spreading dynamics of infectious diseases and the mobility dynamics
through urban transportation systems. We first describe the mobility dynamics
of the urban population as the process of leaving from home, traveling to and
from the activity locations, and engaging in activities. We then embed the
susceptible-exposed-infectious-recovered (SEIR) process over the mobility
dynamics and develops the spatial SEIR model with travel contagion
(Trans-SEIR), which explicitly accounts for contagions both during travel and
during daily activities. We investigate the theoretical properties of the
proposed model and show how activity contagion and travel contagion contribute
to the average number of secondary infections. In the numerical experiments, we
explore how the urban transportation system may alter the fundamental dynamics
of the infectious disease, change the number of secondary infections, promote
the synchronization of the disease across the city, and affect the peak of the
disease outbreaks. The Trans-SEIR model is further applied to the understand
the disease dynamics during the COVID-19 outbreak in New York City, where we
show how the activity and travel contagion may be distributed and how effective
travel control can be implemented with only limited resources. The Trans-SEIR
model along with the findings in our study may have significant contributions
to improving our understanding of the coupling between urban transportation and
disease dynamics, the development of quarantine and control measures of disease
system, and promoting the idea of disease-resilient urban transportation
networks.


  Generation of mechanical oscillation is ubiquitous to wide variety of
intracellular processes. We show that catchbonding behaviour of motor proteins
provides a generic mechanism of generating spontaneous oscillations in
motor-cytoskeletal filament complexes. We obtain the phase diagram to
characterize how this novel catch bond mediated mechanism can give rise to
bistability and sustained limit cycle oscillations and results in very
distinctive stability behaviour, including bistable and non-linearly stabilised
in motor-microtubule complexes in biologically relevant regimes. Hitherto, it
was thought that the primary functional role of the biological catchbond was to
improve surface adhesion of bacteria and cell when subjected to external forces
or flow field. Instead our theoretical study shows that the imprint of this
catch bond mediated physical mechanism would have ramifications for whole gamut
of intracellular processes ranging from oscillations in mitotic spindle
oscillations to activity in muscle fibres.


  The fundamental models of epidemiology describe the progression of an
infectious disease through a population using compartmentalized differential
equations, but do not incorporate population-level heterogeneity in infection
susceptibility. We show that variation strongly influences the rate of
infection, while the infection process simultaneously sculpts the
susceptibility distribution. These joint dynamics influence the force of
infection and are, in turn, influenced by the shape of the initial variability.
Intriguingly, we find that certain susceptibility distributions (the
exponential and the gamma) are unchanged through the course of the outbreak,
and lead naturally to power-law behavior in the force of infection; other
distributions often tend towards these "eigen-distributions" through the
process of contagion. The power-law behavior fundamentally alters predictions
of the long-term infection rate, and suggests that first-order epidemic models
that are parameterized in the exponential-like phase may systematically and
significantly over-estimate the final severity of the outbreak.


  We study a discrete Susceptible-Infected-Recovered (SIR) model for the spread
of infectious disease on a homogeneous tree and the limit behavior of the model
in the case when the tree vertex degree tends to infinity. We obtain the
distribution of the time it takes for a susceptible vertex to get infected in
terms of a solution of a non-linear integral equation under broad assumptions
on the model parameters. Namely, infection rates are assumed to be
time-dependent, and recovery times are given by random variables with a fairly
arbitrary distribution. We then study the behavior of the model in the limit
when the tree vertex degree tends to infinity, and infection rates are
appropriately scaled. We show that in this limit the integral equation of the
discrete model implies an equation for the susceptible population compartment.
This is a master equation in the sense that both the infectious and the
recovered compartments can be explicitly expressed in terms of its solution. In
other words, the master equation implies a continuous SIR model for the joint
time evolution of all three population compartments.


  Mathematical approaches from dynamical systems theory are used in a range of
fields. This includes biology where they are used to describe processes such as
protein-protein interaction and gene regulatory networks. As such networks
increase in size and complexity, detailed dynamical models become cumbersome,
making them difficult to explore and decipher. This necessitates the
application of simplifying and coarse graining techniques in order to derive
explanatory insight. Here we demonstrate that Zwanzig-Mori projection methods
can be used to arbitrarily reduce the dimensionality of dynamical networks
while retaining their dynamical properties. We show that a systematic expansion
around the quasi-steady state approximation allows an explicit solution for
memory functions without prior knowledge of the dynamics. The approach not only
preserves the same steady states but also replicates the transients of the
original system. The method also correctly predicts the dynamics of multistable
systems as well as networks producing sustained and damped oscillations.
Applying the approach to a gene regulatory network from the vertebrate neural
tube, a well characterised developmental transcriptional network, identifies
features of the regulatory network responsible dfor its characteristic
transient behaviour. Taken together, our analysis shows that this method is
broadly applicable to multistable dynamical systems and offers a powerful and
efficient approach for understanding their behaviour.


  Many microbes live in habitats below their optimum temperature. Retention of
metabolic heat by aggregation or insulation would boost growth. Generation of
excess metabolic heat may also provide benefit. A cell that makes excess
metabolic heat pays the cost of production, whereas the benefit may be shared
by neighbors within a zone of local heat capture. Metabolic heat as a shareable
public good raises interesting questions about conflict and cooperation of heat
production and capture. Metabolic heat may also be deployed as a weapon.
Species with greater thermotolerance gain by raising local temperature to
outcompete less thermotolerant taxa. Metabolic heat may provide defense against
bacteriophage attack, by analogy with fever in vertebrates. This article
outlines the theory of metabolic heat in microbial conflict and cooperation,
presenting several predictions for future study.


  Cells are known to sense and respond to their mechanical microenvironment in
profound ways. Various evidence has implicated the adhesome, a body of adhesion
associated proteins, in several cell fate decisions, including differentiation
and proliferation. Furthermore, recent work has demonstrated that substrate
topography can modulate the epigenetic state of fibroblasts, facilitate cell
reprogramming, and temper inflammation. However, the effectors responsible for
such phenomena remain poorly understood. Here we explore the effect of
biomaterial design, adhesion, and intracellular forces of the epigenetic and
transcriptional regulation of cell state during macrophage activation and
somatic cell reprogramming.


  During an infectious disease outbreak, biases in the data and complexities of
the underlying dynamics pose significant challenges in mathematically modelling
the outbreak and designing policy. Motivated by the ongoing response to
COVID-19, we provide a toolkit of statistical and mathematical models beyond
the simple SIR-type differential equation models for analysing the early stages
of an outbreak and assessing interventions. In particular, we focus on
parameter estimation in the presence of known biases in the data, and the
effect of non-pharmaceutical interventions in enclosed subpopulations, such as
households and care homes. We illustrate these methods by applying them to the
COVID-19 pandemic.


  We demonstrate that the epidemic renormalisation group approach to pandemics
provides an effective and simple way to investigate the dynamics of disease
transmission and spreading across different regions of the world. The framework
also allows for reliable projections on the impact of travel limitations and
social distancing measures on global epidemic spread. We test and calibrate it
on reported cases while unveiling the mechanism that governs the delay in the
relative peaks of newly infected cases among different regions of the globe. We
discover that social distancing measures are more effective than travel
limitations across borders in delaying the epidemic peak. We further provide
the link to compartmental models such as the simplistic and time-honoured
SIR-like models. We also show how to generalise the framework to account for
the interactions across several regions of the world, replacing or
complementing large scale simulations.


  Geometry and topology are the main factors that determine the functional
properties of proteins. In this work, we show how to use the Gauss linking
integral (GLN) in the form of a matrix diagram - for a pair of a loop and a
tail - to study both the geometry and topology of proteins with closed loops
e.g. lassos. We show that the GLN method is a significantly faster technique to
detect entanglement in lasso proteins in comparison with other methods. Based
on the GLN technique, we conduct comprehensive analysis of all proteins
deposited in the PDB and compare it to the statistical properties of the
polymers. We found that there are significantly more lassos with negative
crossings than those with positive ones in proteins, the average value of
maxGLN (maximal GLN between loop and pieces of tail) depends logarithmically on
the length of a tail similarly as in the polymers. Next, we show the how high
and low GLN values correlate with the internal exibility of proteins, and how
the GLN in the form of a matrix diagram can be used to study folding and
unfolding routes. Finally, we discuss how the GLN method can be applied to
study entanglement between two structures none of which are closed loops. Since
this approach is much faster than other linking invariants, the next step will
be evaluation of lassos in much longer molecules such as RNA or loops in a
single chromosome.


  Sedentism was a decisive moment in the history of humankind. In a review
article Kay and Kaplan quantified land use for early human settlements and
found that sedentism and the emergence of farming go hand in hand. For these
settlements two primary land use categories, farming and living, can be
identified, whereas for hunter gatherer societies no distinct differences can
be made. It is natural to search for this in the behavior of two different
groups, settlers and farmers. The development of two distinct zones and the two
groups lead us to the hypothesis that the emergence of settlements is the
result of diffusion-driven Turing instability. In this short communication we
further specify this and show that this results in a regular settlement
arrangement as can still be seen today in agricultural regions.


  In December 2019, the first patients in Wuhan, China were diagnosed with a
primary atypical pneumonia, which showed to be unknown and contagious. Since
then, known as COVID-19 disease, the responsible viral pathogen, SARS-CoV-2,
has spread around the world in a pandemic. Decisions on how to deal with the
crisis are often based on simulations of the pandemic spread of the virus. The
results of some of these, as well as their methodology and possibilities for
improvement, will be described in more detail in this paper in order to inform
beyond the current public health dogma called "flatten-the-curve". There are
several ways to model an epidemic in order to simulate the spread of diseases.
Depending on the timeliness, scope and quality of the associated real data,
these multivariable models differ in the value of used parameters, but also in
the selection of considered influencing factors. It was exemplarily shown that
epidemics in their course are simulated more realistically by models that
assume subexponential growth. Furthermore, various simulations of the COVID-19
pandemic were presented in an European perspective, compared against each other
and discussed in more detail. It is difficult to estimate how credible the
simulations of the pandemic models currently are, so it remains to be seen
whether the spread of the pandemic can be effectively reduced by the measures
taken. Whether a model works well in reality is largely determined by the
quality and scope of its underlying data. Past studies have shown that
countermeasures are able to reduce reproduction numbers or transmission rates
in epidemics. In addition to that, the presented modelling study provides a
good framework for the creation of subexponential-growth-models for assessing
the spread of COVID-19.


  Comprehensive simulations of the entire biochemistry of cells have great
potential to help physicians treat disease and help engineers design biological
machines. But such simulations must model networks of millions of molecular
species and reactions.
  The Stochastic Simulation Algorithm (SSA) is widely used for simulating
biochemistry, especially systems with species populations small enough that
discreteness and stochasticity play important roles. However, existing serial
SSA methods are prohibitively slow for comprehensive networks, and existing
parallel SSA methods, which use periodic synchronization, sacrifice accuracy.
  To enable fast, accurate, and scalable simulations of biochemistry, we
present an exact parallel algorithm for SSA that partitions a biochemical
network into many SSA processes that simulate in parallel. Our parallel SSA
algorithm exactly coordinates the interactions among these SSA processes and
the species state they share by structuring the algorithm as a parallel
discrete event simulation (DES) application and using an optimistic parallel
DES simulator to synchronize the interactions. We anticipate that our method
will enable unprecedented biochemical simulations.


  We present an early version of a
Susceptible-Exposed-Infected-Recovered-Deceased (SEIRD) mathematical model
based on partial differential equations coupled with a heterogeneous diffusion
model. The model describes the spatio-temporal spread of the COVID-19 pandemic,
and aims to capture dynamics also based on human habits and geographical
features. To test the model, we compare the outputs generated by a
finite-element solver with measured data over the Italian region of Lombardy,
which has been heavily impacted by this crisis between February and April 2020.
Our results show a strong qualitative agreement between the simulated forecast
of the spatio-temporal COVID-19 spread in Lombardy and epidemiological data
collected at the municipality level. Additional simulations exploring
alternative scenarios for the relaxation of lockdown restrictions suggest that
reopening strategies should account for local population densities and the
specific dynamics of the contagion. Thus, we argue that data-driven simulations
of our model could ultimately inform health authorities to design effective
pandemic-arresting measures and anticipate the geographical allocation of
crucial medical resources.


  Now that various countries are or will soon be moving towards relaxing
shelter-in-place rules, it is important that people use a face covering, to
avoid an exponential resurgence of the spreading of the coronavirus disease
(COVID-19). Here we report a pre-registered online experiment (N=2,459) using a
heterogenous, although not representative, sample of people living in the USA,
where we test the relative effect of messages highlighting that the coronavirus
is a threat to "you" vs "your family" vs "your community" vs "your country" on
self-reported intentions to wear a face covering. Results show that focusing on
"your community" promotes intentions to wear a face covering relative to the
baseline; the trend is the same when comparing "your community" to the other
conditions, but not significant. We also conducted pre-registered analyses of
gender differences on intentions to wear a face covering. We find that men less
than women intend to wear a face covering, but this difference almost
disappears in counties where wearing a face covering is mandatory. We also find
that men less than women believe that they will be seriously affected by the
coronavirus, and this partly mediates gender differences in intentions to wear
a face covering (this is particularly ironic because official statistics
actually show that men are affected by the COVID-19 more seriously than women).
Finally, we also find gender differences in self-reported negative emotions
felt when wearing a face covering. Men more than women agree that wearing a
face covering is shameful, not cool, a sign of weakness, and a stigma; and
these gender differences also mediate gender differences in intentions to wear
a face covering.


  Accurately estimating the effectiveness of stay-at-home orders (SHOs) on
reducing social contact and disease spread is crucial for mitigating pandemics.
Leveraging individual-level location data for 10 million smartphones, we
observe that by April 30th---when nine in ten Americans were under a
SHO---daily movement had fallen 70% from pre-COVID levels. One-quarter of this
decline is causally attributable to SHOs, with wide demographic differences in
compliance, most notably by political affiliation. Likely Trump voters reduce
movement by 9% following a local SHO, compared to a 21% reduction among their
Clinton-voting neighbors, who face similar exposure risks and identical
government orders. Linking social distancing behavior with an epidemic model,
we estimate that reductions in movement have causally reduced SARS-CoV-2
transmission rates by 49%.


  Accurate modeling of opinion dynamics has the potential to help us understand
polarization and what makes effective political discourse possible or
impossible. Here, we use physics-based methods to model the evolution of
political opinions within a continuously distributed population. We utilize a
network-free system of determining political influence and a local-attraction,
distal-repulsion dynamic for reaction to perceived content. Our approach allows
for the incorporation of intergroup bias such that messages from trusted
in-group sources enjoy greater leeway than out-group ones. We are able to
extrapolate these nonlinear microscopic dynamics to macroscopic population
distributions by using probabilistic functions representing biased
environments. The framework we put forward can reproduce real-world political
distributions and experimentally observed dynamics, and is amenable to further
refinement as more data becomes available.


  Strong social distancing restrictions have been crucial to controlling the
COVID-19 outbreak thus far, and the next question is when and how to relax
these restrictions. A sequential timing of relaxing restrictions across groups
is explored in order to identify policies that simultaneously reduce health
risks and economic stagnation relative to current policies. The goal will be to
mitigate health risks, particularly among the most fragile sub-populations,
while also managing the deleterious effect of restrictions on economic
activity. The results of this paper show that a properly constructed sequential
release of age-defined subgroups from strict social distancing protocols can
lead to lower overall fatality rates than the simultaneous release of all
individuals after a lockdown. The optimal release policy, in terms of
minimizing overall death rate, must be sequential in nature, and it is
important to properly time each step of the staggered release. This model
allows for testing of various timing choices for staggered release policies,
which can provide insights that may be helpful in the design, testing, and
planning of disease management policies for the ongoing COVID-19 pandemic and
future outbreaks.


  We consider a minimalist model for the Sterile Insect Technique (SIT),
assuming that residual fertility can occur in the sterile male
population.Taking into account that we are able to get regular measurements
from the biological system along the control duration, such as the size of the
wild insect population, we study different control strategies that involve
either continuous or periodic impulsive releases. We show that a combination of
open-loop control with constant large releases and closed-loop nonlinear
control, i.e. when releases are adjusted according to the wild population size
estimates, leads to the best strategy in terms both of number of releases and
total quantity of sterile males to be released.Last but not least, we show that
SIT can be successful only if the residual fertility is less than a threshold
value that depends on the wild population biological parameters. However, even
for small values, the residual fertility induces the use of such large releases
that SIT alone is not always reasonable from a practical point of view and thus
requires to be combined with other control tools. We provide applications
against a mosquito species, \textit{Aedes albopictus}, and a fruit fly,
\textit{Bactrocera dorsalis}, and discuss the possibility of using SIT when
residual fertility, among the sterile males, can occur.


  SkeleDock is a scaffold docking algorithm which uses the structure of a
protein-ligand complex as a template to model the binding mode of a chemically
similar system. This algorithm was evaluated in the D3R Grand Challenge 4 pose
prediction challenge, where it achieved competitive performance. Furthermore,
we show that, if crystallized fragments of the target ligand are available,
SkeleDock can outperform rDock docking software at predicting the binding mode.
This article also addresses the capacity of this algorithm to model macrocycles
and deal with scaffold hopping. SkeleDock can be accessed at
https://playmolecule.org/SkeleDock/.


  Smartphone-based electronic contact tracing is currently considered an
essential tool towards easing lockdowns, curfews, and shelter-in-place orders
issued by most governments around the world in response to the 2020 novel
coronavirus (SARS-CoV-2) crisis. While the focus on developing smartphone-based
contact tracing applications or apps has been on privacy concerns stemming from
the use of such apps, an important question that has not received sufficient
attention is: How reliable will such smartphone-based electronic contact
tracing be?
  This is a technical question related to how two smartphones reliably register
their mutual proximity. Here, we examine in detail the technical prerequisites
required for effective smartphone-based contact tracing. The underlying
mechanism that any contact tracing app relies on is called Neighbor Discovery
(ND), which involves smartphones transmitting and scanning for Bluetooth
signals to record their mutual presence whenever they are in close proximity.
The hardware support and the software protocols used for ND in smartphones,
however, were not designed for reliable contact tracing. In this paper, we
quantitatively evaluate how reliably can smartphones do contact tracing. Our
results point towards the design of a wearable solution for contact tracing
that can overcome the shortcomings of a smartphone-based solution to provide
more reliable and accurate contact tracing. To the best of our knowledge, this
is the first study that quantifies, both, the suitability and also the
drawbacks of smartphone-based contact tracing. Further, our results can be used
to parameterize a ND protocol to maximize the reliability of any contact
tracing app that uses it.


  We present and analyze a model aiming at recovering as dynamical outcomes of
tree-grass interactions the wide range of vegetation physiognomies observable
in the savanna biome along rainfall gradients at regional/continental scales.
The model is based on two ordinary differential equations (ODE), for woody and
grass biomass. It is parameterized from literature and retains mathematical
tractability, since we restricted it to the main processes, notably tree-grass
asym-metric interactions (either facilitative or competitive) and the
grass-fire feedback. We used a fully qualitative analysis to derive all
possible long term dynamics and express them in a bifur-cation diagram in
relation to mean annual rainfall and fire frequency. We delineated domains of
monostability (forest, grassland, savanna), of bistability (e.g.
forest-grassland or forest-savanna) and even tristability. Notably, we
highlighted regions in which two savanna equilibria may be jointly stable
(possibly in addition to forest or grassland). We verified that common
knowledge about decreasing woody biomass with increasing fire frequency is
recovered for all levels of rainfall, contrary to previous attempts using
analogous ODE frameworks. Thus, this framework appears able to render more
realistic and diversified outcomes than often thought of. Our model can help
figure out the ongoing dynamics of savanna vegetation in large territories for
which local data are sparse or absent. To explore the bifurcation diagram with
different combinations of the model parameters, we have developed a
user-friendly R-Shiny application freely available at :
https://gitlab.com/cirad-apps/tree-grass.


  Our understanding of risk preferences can be sharpened by considering their
evolutionary basis. The existing literature has focused on two sources of risk:
idiosyncratic risk and aggregate risk. We introduce a new source of risk,
heritable risk, in which there is a positive correlation between the fitness of
a newborn agent and the fitness of her parent. Heritable risk was plausibly
common in our evolutionary past and it leads to a strictly higher growth rate
than the other sources of risk. We show that the presence of heritable risk in
the evolutionary past may explain the tendency of people to exhibit skewness
loving today.


  The premise of this paper is the following value proposition: Models are good
when they describe a system of phenomena, and they are better when they can
predict the effect of interventions upon the system.
  We introduce a formalism by which dynamical system models can be obtained by
specifying basic constituent processes in the form of petri nets, mediated by
the Law of Mass Action. We prove the universality of this procedure with
respect to ordinary first order differential equations in rational functions.
In addition, we introduce a simple graphical procedure for calculating
dynamical systems from petri nets in our formalism, which we use to recover
several well-known dynamical systems. For proof-of-concept, we use our method
to obtain a differential equation model for the HES1 gene transcription network
that predicts post-intervention behaviours of the network.


  We present results of different approaches to model the evolution of the
COVID-19 epidemic in Argentina, with a special focus on the megacity conformed
by the city of Buenos Aires and its metropolitan area, including a total of 41
districts with over 13 million inhabitants. We first highlight the relevance of
interpreting the early stage of the epidemic in light of incoming infectious
travelers from abroad. Next, we critically evaluate certain proposed solutions
to contain the epidemic based on instantaneous modifications of the
reproductive number. Finally, we build increasingly complex and realistic
models, ranging from simple homogeneous models used to estimate local
reproduction numbers, to fully coupled inhomogeneous (deterministic or
stochastic) models incorporating mobility estimates from cell phone location
data. The models are capable of producing forecasts highly consistent with the
official number of cases with minimal parameter fitting and fine-tuning. We
discuss the strengths and limitations of the proposed models, focusing on the
validity of different necessary first approximations, and caution future
modeling efforts to exercise great care in the interpretation of long-term
forecasts, and in the adoption of non-pharmaceutical interventions backed by
numerical simulations.


  Notch signaling is an evolutionary conserved cell-cell communication pathway.
Besides regulating cell-fate decisions at an individual cell level, Notch
signaling coordinates the emergent spatiotemporal patterning in a tissue
through ligand-receptor interactions among transmembrane molecules of
neighboring cells, as seen in embryonic development, angiogenesis, or wound
healing. Due to its ubiquitous nature, Notch signaling is also implicated in
several aspects of cancer progression, including tumor angiogenesis, stemness
of cancer cells and cellular invasion. Here, we review experimental and
computational models to help understand the operating principles of cell
patterning driven by Notch signaling. First, we discuss the basic mechanisms of
spatial patterning via canonical lateral inhibition and lateral induction
mechanisms, including examples from angiogenesis, inner ear development and
cancer metastasis. Next, we analyze additional layers of complexity in the
Notch pathway such as the effects of varying cell sizes and shapes,
ligand-receptor binding within the same cell, variable binding affinity of
different ligand/receptor subtypes, and filopodia. Finally, we discuss some
recent evidence of mechanosensitivity in the Notch pathway in driving
collective epithelial cell migration and cardiovascular morphogenesis.


  Countries are recording health information on the global spread of COVID-19
using different methods, sometimes changing the rules after a few days. They
are all publishing the number of new individuals infected, cured and dead,
along with some supplementary data. These figures are often recorded in a
non-uniform manner and do not match the standard definitions of these
variables. However, in this paper we show that the Kaplan-Meier curves
calculated with them could provide useful information about the dynamics of the
disease in different countries. Our aim is to present a robust and simple model
to show certain characteristics of the evolution of the dynamic process,
showing that the differences of evolution among the countries is reflected in
the corresponding Kaplan-Meier-type curves. We compare the curves obtained for
the most affected countries so far, proposing possible interpretations of the
properties that distinguish them.