  All stem cell fate transitions, including the metabolic reprogramming of stem
cells and the somatic reprogramming of fibroblasts into pluripotent stem cells,
can be understood from a unified theoretical model of cell fates. Each cell
fate transition can be regarded as a phase transition in DNA supercoiling.
However, there has been a dearth of quantitative biophysical models to explain
and predict the behaviors of these phase transitions. The generalized Ising
model is proposed to define such phase transitions. The model predicts that,
apart from temperature-induced phase transitions, there exists DNA torsion
frequency-induced phase transitions. Major transitions in epigenetic states,
from stem cell activation to differentiation and reprogramming, can be
explained by such torsion frequency-induced phase transitions, with important
implications for regenerative medicine and medical diagnostics in the future.


  Text-based representations of chemicals and proteins can be thought of as
unstructured languages codified by humans to describe domain-specific
knowledge. Advances in natural language processing (NLP) methodologies in the
processing of spoken languages accelerated the application of NLP to elucidate
hidden knowledge in textual representations of these biochemical entities and
then use it to construct models to predict molecular properties or to design
novel molecules. This review outlines the impact made by these advances on drug
discovery and aims to further the dialogue between medicinal chemists and
computer scientists.


  The 2019 novel coronavirus (2019-nCoV) is currently causing a widespread
outbreak centered on Hubei province, China and is a major public health
concern. Taxonomically 2019-nCoV is closely related to SARS-CoV and
SARS-related bat coronaviruses, and it appears to share a common receptor with
SARS-CoV (ACE-2). Here, we perform structural modeling of the 2019-nCoV spike
glycoprotein. Our data provide support for the similar receptor utilization
between 2019-nCoV and SARS-CoV, despite a relatively low amino acid similarity
in the receptor binding module. Compared to SARS-CoV, we identify an extended
structural loop containing basic amino acids at the interface of the receptor
binding (S1) and fusion (S2) domains, which we predict to be
proteolytically-sensitive. We suggest this loop confers fusion activation and
entry properties more in line with MERS-CoV and other coronaviruses, and that
the presence of this structural loop in 2019-nCoV may affect virus stability
and transmission.


  During mammalian embryo development, reprogramming of DNA methylation plays
important roles in the erasure of parental epigenetic memory and the
establishment of na\"{i}ve pluripogent cells. Multiple enzymes that regulate
the processes of methylation and demethylation work together to shape the
pattern of genome-scale DNA methylation and guid the process of cell
differentiation. Recent availability of methylome information from single-cell
whole genome bisulfite sequencing (scBS-seq) provides an opportunity to study
DNA methylation dynamics in the whole genome in individual cells, which reveal
the heterogeneous methylation distributions of enhancers in embryo stem cells
(ESCs). In this study, we developed a computational model of enhancer
methylation inheritance to study the dynamics of genome-scale DNA methylation
reprogramming during exit from pluripotency. The model enables us to track
genome-scale DNA methylation reprogramming at single-cell level during the
embryo development process, and reproduce the DNA methylation heterogeneity
reported by scBS-seq. Model simulations show that DNA methylation heterogeneity
is an intrinsic property driven by cell division along the development process,
and the collaboration between neighboring enhancers is required for
heterogeneous methylation. Our study suggest that the mechanism of genome-scale
oscillation proposed by Rulands et al. (2018) might not necessary to the DNA
methylation during exit from pluripotency.


  The outbreak of novel coronavirus-caused pneumonia (COVID-19) in Wuhan has
attracted worldwide attention. Here, we propose a generalized SEIR model to
analyze this epidemic. Based on the public data of National Health Commission
of China from Jan. 20th to Feb. 9th, 2020, we reliably estimate key epidemic
parameters and make predictions on the inflection point and possible ending
time for 5 different regions. According to optimistic estimation, the epidemics
in Beijing and Shanghai will end soon within two weeks, while for most part of
China, including the majority of cities in Hubei province, the success of
anti-epidemic will be no later than the middle of March. The situation in Wuhan
is still very severe, at least based on public data until Feb. 15th. We expect
it will end up at the beginning of April. Moreover, by inverse inference, we
find the outbreak of COVID-19 in Mainland, Hubei province and Wuhan all can be
dated back to the end of December 2019, and the doubling time is around two
days at the early stage.


  Prediction of protein structures using computational approaches has been
explored for over two decades, paving a way for more focused research and
development of algorithms in comparative modelling, ab intio modelling and
structure refinement protocols. A tremendous success has been witnessed in
template-based modelling protocols, whereas strategies that involve
template-free modelling still lag behind, specifically for larger proteins (>
150 a.a.). Various improvements have been observed in ab initio protein
structure prediction methodologies overtime, with recent ones attributed to the
usage of deep learning approaches to construct protein backbone structure from
its amino acid sequence. This review highlights the major strategies undertaken
for template-free modelling of protein structures while discussing few tools
developed under each strategy. It will also briefly comment on the progress
observed in the field of ab initio modelling of proteins over the course of
time as seen through the evolution of CASP platform.


  In this work, we analyse the relationship between heterogeneity and
cooperation. Previous investigations suggest that this relation is nontrivial,
as some authors found that heterogeneity sustains cooperation, while others
obtained different results. Among the possible forms of heterogeneity, we focus
on the individual perception of risks and rewards related to a generic event,
that can show up in a number of social and biological systems. The modelling
approach is based on the framework of Evolutionary Game Theory. To represent
this kind of heterogeneity, we implement small and local perturbations on the
payoff matrix of simple 2-strategy games, as the Prisoner's Dilemma. So, while
usually the payoff is considered as a global and time-invariant structure, i.e.
it is the same for all individuals of a population at any time, in our model
its value is continuously affected by small variations, both in time and space
(i.e. position on a lattice). We found that such perturbations can be
beneficial or detrimental to cooperation, depending on their setting. Notably,
cooperation is strongly supported when perturbations act on the main diagonal
of the payoff matrix, whereas when they act on the off-diagonal the resulting
effect is more difficult to quantify. To conclude, the proposed model shows a
rich spectrum of possible equilibria, whose interpretation might offer insights
and enrich the description of several systems.


  Adhesive interactions between elastic structures such as graphene sheets,
carbon nanotubes, and microtubules have been shown to exhibit hysteresis due to
irrecoverable energy loss associated with bond breakage, even in static
(rate-independent) experiments. To understand this phenomenon, we start with a
minimal theory for the peeling of a thin sheet from a substrate, coupling the
local event of bond breaking to the nonlocal elastic relaxation of the sheet
and show that this can drive static adhesion hysteresis over a
bonding/debonding cycle. Using this model we quantify hysteresis in terms of
the adhesion and elasticity parameters of the system. This allows us to derive
a scaling relation that preserves hysteresis at different levels of granularity
while resolving a seeming paradox of lattice trapping in the continuum limit of
a discrete fracture process. Finally, to verify our theory, we use new
experiments to demonstrate and measure adhesion hysteresis in bundled
microtubules.


  When modelling epidemics or spread of information on online social networks,
it is crucial to include not just the density of the connections through which
infections can be transmitted, but also the variability of susceptibility.
Different people have different chance to be infected by a disease (due to age
or general health conditions), or, in case of opinions, ones are easier to be
convinced by others, or stronger at sharing their opinions. The goal of this
work is to examine the effect of multiple types of nodes on various random
graphs such as Erd\H{o}s--R\'enyi random graphs, preferential attachment random
graphs and geometric random graphs. We used two models for the dynamics: SEIR
model with vaccination and a version of voter model for exchanging opinions. In
the first case, among others, various vaccination strategies are compared to
each other, while in the second case we studied sevaral initial configurations
to find the key positions where the most effective nodes should be placed to
disseminate opinions.


  This work proposes a unified framework to leverage biological information in
network propagation-based gene prioritization algorithms. Preliminary results
on breast cancer data show significant improvements over state-of-the-art
baselines, such as the prioritization of genes that are not identified as
potential candidates by interactome-based algorithms, but that appear to be
involved in/or potentially related to breast cancer, according to a functional
analysis based on recent literature.


  The COVID-19 (formerly, 2019-nCoV) epidemic has become a global health
emergency, as such, WHO declared PHEIC. China has taken the most hit since the
outbreak of the virus, which could be dated as far back as late November by
some experts. It was not until January 23rd that the Wuhan government finally
recognized the severity of the epidemic and took a drastic measure to curtain
the virus spread by closing down all transportation connecting the outside
world. In this study, we seek to answer a few questions: How did the virus get
spread from the epicenter Wuhan city to the rest of the country? To what extent
did the measures, such as, city closure and community quarantine, help
controlling the situation? More importantly, can we forecast any significant
future development of the event had some of the conditions changed? By
collecting and visualizing publicly available data, we first show patterns and
characteristics of the epidemic development; we then employ a mathematical
model of disease transmission dynamics to evaluate the effectiveness of some
epidemic control measures, and more importantly, to offer a few tips on
preventive measures.


  Single-cell RNA sequencing (scRNA-seq) has revolutionized biological
discovery, providing an unbiased picture of cellular heterogeneity in tissues.
While scRNA-seq has been used extensively to provide insight into both healthy
systems and diseases, it has not been used for disease prediction or
diagnostics. Graph Attention Networks (GAT) have proven to be versatile for a
wide range of tasks by learning from both original features and graph
structures. Here we present a graph attention model for predicting disease
state from single-cell data on a large dataset of Multiple Sclerosis (MS)
patients. MS is a disease of the central nervous system that can be difficult
to diagnose. We train our model on single-cell data obtained from blood and
cerebrospinal fluid (CSF) for a cohort of seven MS patients and six healthy
adults (HA), resulting in 66,667 individual cells. We achieve 92 % accuracy in
predicting MS, outperforming other state-of-the-art methods such as a graph
convolutional network and a random forest classifier. Further, we use the
learned graph attention model to get insight into the features (cell types and
genes) that are important for this prediction. The graph attention model also
allow us to infer a new feature space for the cells that emphasizes the
differences between the two conditions. Finally we use the attention weights to
learn a new low-dimensional embedding that can be visualized. To the best of
our knowledge, this is the first effort to use graph attention, and deep
learning in general, to predict disease state from single-cell data. We
envision applying this method to single-cell data for other diseases.


  Population dynamics of a competitive two-species system under the influence
of random events are analyzed and expressions for the steady-state population
mean, fluctuations, and cross-correlation of the two species are presented. It
is shown that random events cause the population mean of each specie to make
smooth transition from far above to far below of its growth rate threshold. At
the same time, the population mean of the weaker specie never reaches the
extinction point. It is also shown that, as a result of competition, the
relative population fluctuations do not die out as the growth rates of both
species are raised far above their respective thresholds. This behavior is most
remarkable at the maximum competition point where the weaker specie's
population statistics becomes completely chaotic regardless of how far its
growth rate in raised.


  Insect-borne diseases are diseases carried by insects affecting humans,
animals or plants. They have the potential to generate massive outbreaks such
as the Zika epidemic in 2015-2016 mostly distributed in the Americas, the
Pacific and Southeast Asia, and the multi-foci outbreak caused by the bacterium
{\it Xylella fastidiosa} in Europe in the 2010s. In this article, we propose
and analyze the behavior of a spatially-explicit compartmental model adapted to
pathosystems with fixed hosts and mobile vectors disseminating the disease. The
behavior of this model based on a system of partial differential equations is
complementarily characterized via a theoretical study of its equilibrium states
and a numerical study of its transitive phase using global sensitivity
analysis. The results are discussed in terms of implications concerning the
surveillance and control of the disease over a medium-to-long temporal horizon.


  The recent outbreak of COVID-19 in Mainland China is characterized by a
distinctive algebraic, sub-exponential increase of confirmed cases during the
early phase of the epidemic, contrasting an initial exponential growth expected
for an unconstrained outbreak with sufficiently large reproduction rate.
Although case counts vary significantly between affected provinces in Mainland
China, the scaling law $t^{\mu}$ is surprisingly universal, with a range of
exponents $\mu=2.1\pm0.3$. The universality of this behavior indicates that
despite social, regional, demographical, geographical, and socio-economical
heterogeneities of affected Chinese provinces, this outbreak is dominated by
fundamental mechanisms that are not captured by standard epidemiological
models. We show that the observed scaling law is a direct consequence of
containment policies that effectively deplete the susceptible population. To
this end we introduce a parsimonious model that captures both, quarantine of
symptomatic infected individuals as well as population wide isolation in
response to mitigation policies or behavioral changes. For a wide range of
parameters, the model reproduces the observed scaling law in confirmed cases
and explains the observed exponents. Quantitative fits to empirical data permit
the identification of peak times in the number of asymptomatic or
oligo-symptomatic, unidentified infected individuals, as well as estimates of
local variations in the basic reproduction number. The model implies that the
observed scaling law in confirmed cases is a direct signature of effective
contaiment strategies and/or systematic behavioral changes that affect a
substantial fraction of the susceptible population. These insights may aid the
implementation of containment strategies in potential export induced COVID-19
secondary outbreaks elsewhere or similar future outbreaks of other emergent
infectious diseases.


  Clusters appear in nature in a diversity of contexts, involving distances as
long as the cosmological ones, and down to atoms and molecules and the very
small nuclear size. They also appear in several other scenarios, in particular
in biological systems as in ants, bees, birds, fishes, gnus and rats, for
instance. Here we describe a model composed of a set of female and male
individuals that obeys simple rules that rapidly transform an uniform initial
state into a single cluster that evolves in time as a stable dynamical
structure. We show that the center of mass of the structure moves as a random
walk, and that the size of the cluster engenders a power law behavior in terms
of the number of individuals in the system. Moreover, we also examine other
possibilities, in particular the case of two distinct species that can evolve
to form one or two distinct clusters.


  Given a reaction (network) we are looking for minimal sets of species
starting from which all the species will have positive concentrations for all
positive times in the domain of existence of the solution of the induced
kinetic differential equation. We present three algorithms to solve the
problem. The first one essentially checks all the possible subsets of the sets
of species. This can obviously work for only a few dozen species because of
combinatorial explosion. The second one is based on an integer programming
reformulation of the problem. The third one walks around the state space of the
problem in a systematic way and produces all the minimal sets of the
advantageous initial species, and works also for large systems. All the
algorithms rely heavily on the concept of Volpert indices, used earlier for the
decomposition of overall reactions [Kov\'acs \textit{et al., Physical Chemistry
Chemical Physics} 2004, 6, 1236]. Relations to the permanence hypothesis,
possible economic or medical uses of the solution of the problem are analyzed,
and open problems are formulated at the end.


  The risk of dying increases exponentially with age, in humans as well as in
many other species. This increase is often attributed to the "accumulation of
damage" known to occur in many biological structures and systems. The aim of
this paper is to describe a generic model of damage accumulation and death in
which mortality increases exponentially with age. The damage-accumulation
process is modeled by a stochastic process know as a queue, and risk of dying
is a function of the accumulated damage, i.e. length of the queue. The model
has four parameters and the main characteristics of the model are: (i) damage
occurs at random times with a constant high rate; (ii) the damage is repaired
at a limited rate, and consequently damage can accumulate; (iii) the efficiency
of the repair mechanism decays linearly with age; (iv) the risk of dying is a
function of the accumulated damage. Using standard results from the
mathematical theory of queues it is shown that there is an exponential
dependence between risk of dying and age in these models, and that this
dependency holds irrespective of how the damage-accumulation process is
modeled. Furthermore, the ways in which this exponential dependence is shaped
by the model parameters are also independent of the details of the damage
accumulation process. These generic features suggest that the model could be
useful when interpreting changes in the relation between age and mortality in
real data. To examplify, historical mortality data from Sweden are interpreted
in the light of the model. The decrease in mortality seen between cohorts born
in 1905, compared to those born in 1885, can be accounted for by higher
threshold to damage. This fits well with the many advances made in public
health during the 20th century.


  Social interaction between microbes can be described at many levels of
details, ranging from the biochemistry of cell-cell interactions to the
ecological dynamics of populations. Choosing the best level to model microbial
communities without losing generality remains a challenge. Here we propose to
model cross-feeding interactions at an intermediate level between genome-scale
metabolic models of individual species and consumer-resource models of
ecosystems, which is suitable to empirical data. We applied our method to three
published examples of multi-strain Escherichia coli communities with increasing
complexity consisting of uni-, bi-, and multi-directional cross-feeding of
either substitutable metabolic byproducts or essential nutrients. The
intermediate-scale model accurately described empirical data and could quantify
exchange rates elusive by other means, such as the byproduct secretions, even
for a complex community of 14 amino acid auxotrophs. We used the three models
to study each community's limits of robustness to perturbations such as
variations in resource supply, antibiotic treatments and invasion by other
"cheaters" species. Our analysis provides a foundation to quantify
cross-feeding interactions from experimental data, and highlights the
importance of metabolic exchanges in the dynamics and stability of microbial
communities.


  How much information is transmitted when animals use emotions to communicate?
It is clear that emotions are used as communication systems in humans and other
species. The quantitative theory of emotion information presented here is based
on Shannon's mathematical theory of information in communication systems. The
theory explains myriad aspects of emotion communication and offers dozens of
new directions for research. It is superior to the "contagion" theory of
emotion spreading, which is currently dominant. One important application of
the information theory of emotion communication is that it permits the
development of emotion security systems for social networks to guard against
the widespread emotion manipulation we see online today.


  It is possible that there are post-translational circadian oscillators that
continue functioning in the absence of negative feedback transcriptional
repression in many cell types from diverse organisms. Apart from the KaiABC
system from cyanobacteria, the molecular components and interactions required
to create in-vitro ("test-tube") circadian oscillations in different cell types
are currently unknown. Inspired by the KaiABC system, I provide
"proof-of-principle" mathematical models that a protein with 2 (or more)
modification sites which selectively sequesters an effector/cofactor molecule
can function as a circadian time-keeper. The 2-site mechanism can be
implemented using two relatively simple coupled non-linear ODEs in terms of
site occupancy; the models do not require overly special fine-tuning of
parameters for generating stable limit cycle oscillations.


  Solving the recombination equation has been a long-standing challenge of
\emph{deterministic} population genetics. We review recent progress obtained by
introducing ancestral processes, as traditionally used in the context of
\emph{stochastic} models of population genetics, into the deterministic
setting. With the help of an ancestral partitioning process, which is obtained
by letting population size tend to infinity (without rescaling parameters or
time) in an ancestral recombination graph, we obtain the solution to the
recombination equation in a transparent form.


  Worldwide Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)
infection is disrupting in the economy and anxiety of people. The public
anxiety has increased the psychological burden on government and healthcare
professionals, resulting in a government worker suicide in Japan. The terrified
people are asking the government for border measures. However, are border
measures possible for this virus? By analyzing 48 almost complete virus genome
sequences, we found out that the viruses that invaded Taiwan, the United
States, and Japan were introduced independently. We identified thirteen
parsimony-informative sites and three groups (CTC, TCC, and TCT). Viruses found
outside China did not form a monophyletic clade, opposite to previous study.
These results suggest the difficulty of implementing effective border measures
against this virus.


  Due to the contribution of drug-target binding kinetics to drug efficacy,
there is a high level of interest in developing methods to predict drug-target
binding kinetic parameters. During the review period, a wide range of enhanced
sampling molecular dynamics simulation-based methods has been developed for
computing drug-target binding kinetics and studying binding and unbinding
mechanisms. Here, we assess the performance of these methods considering two
benchmark systems in detail: mutant T4 lysozyme-ligand complexes and a large
set of N-HSP90-inhibitor complexes. The results indicate that some of the
simulation methods can already be usefully applied in drug discovery or lead
optimization programs but that further studies on more high-quality
experimental benchmark datasets are necessary to improve and validate
computational methods.


  Ions in channels have been imagined as hard balls in a macroscopic mechanical
model, for a very long time. Hard balls interact by collisions in such models,
randomly knocking each other on and off `binding' sites in thermal motion. But
ions have large charge, and the hard balls of classical models do not. The
electrodynamics of charge guarantee strong correlations between the movements
of ions on all time scales, even those of atomic scale thermal motion.
Correlations are present whenever Maxwell's equations apply, so they are
present in individual trajectories, not just averages. Indeed, in a series
system like an idealized narrow channel, the correlation is perfect (within the
accuracy of Maxwell's equations) because of conservation of total current (that
includes Maxwell's displacement current, the ethereal ${\varepsilon
}_0{\partial \boldsymbol{\mathrm{E}}}/{\partial t}$).The ethereal component of
current prevents spatial variation of total current in a series system. The
stochastic complexity of spatial thermal motion disappears for current }in a
series system. Total current does not depend on location in a series system
like a narrow ion channel on any time scale. Spatial variables are not needed
in the description of total current in a one dimensional channel on any time
scale, according to the Maxwell equations. Removing the spatial dependence of
total current should dramatically simplify theories and simulations, one might
imagine.


  Chemical reactions occur in energy, environmental, biological, and many other
natural systems, and the inference of the reaction networks is essential to
understand and design the chemical processes in engineering and life sciences.
Yet, revealing the reaction pathways for complex systems and processes is still
challenging due to the lack of knowledge of the involved species and reactions.
Here, we present a neural network approach that autonomously discovers reaction
pathways from the time-resolved species concentration data. The proposed
Chemical Reaction Neural Network (CRNN), by design, satisfies the fundamental
physics laws, including the Law of Mass Action and the Arrhenius Law.
Consequently, the CRNN is physically interpretable such that the reaction
pathways can be interpreted, and the kinetic parameters can be quantified
simultaneously from the weights of the neural network. The inference of the
chemical pathways is accomplished by training the CRNN with species
concentration data via stochastic gradient descent. We demonstrate the
successful implementations and the robustness of the approach in elucidating
the chemical reaction pathways of several chemical engineering and biochemical
systems. The autonomous inference by the CRNN approach precludes the need for
expert knowledge in proposing candidate networks and addresses the curse of
dimensionality in complex systems. The physical interpretability also makes the
CRNN capable of not only fitting the data for a given system but also
developing knowledge of unknown pathways that could be generalized to similar
chemical systems.


  Since the outbreak of COVID-19, many data analyses have been done. Some of
them are based on the classical epidemiological approach that assumes an
exponential growth, but a few studies report that a power-law scaling may
provide a better fit to the currently available data. Hereby, we examine the
data in China (01/20/2020--02/24/2020), and indeed find that the growth closely
follows a power-law kinetics over a significantly wide time period. The
exponents are $2.48(20)$, $2.21(6)$ and $4.26(12)$ for the number of confirmed
infections, deaths and cured cases, respectively, indicating an underlying
small-world network structure in the pandemic. While no obvious deviations from
the power-law growth can be seen yet for the number of deaths and cured cases,
negative deviations have clearly appeared in the number of infections,
particularly that for the region outside Hubei. This suggests the beginning of
the slowing-down of the virus spreading due to the huge containment effort.
Meanwhile, we find that despite the dramatic difference in magnitudes, the
growth kinetics of the infection number exhibits much similarity for Hubei
province and the region outside Hubei. On this basis, in log-log plot, we
rescale the infection number for the region outside Hubei such that it overlaps
as much as possible with the total infection number in China, from which an
approximate extrapolation yields the maximum of the pandemic around March 3,
2020, with the number of infections about $83,000$. Further, by analyzing the
kinetics of the mortality in log-log scale, we obtains a rough estimate that
near March 3, the death rate of COVID-19 would be about $4.7\%\thicksim 5.0\%$
for Hubei province and $0.7\%\thicksim1.0\%$ for the region outside Hubei. We
emphasize that our predictions may be quantitatively unreliable, since the data
analysis is purely empirical and various assumptions are used.


  Cytotoxic T lymphocytes (T) and natural killer (NK) cells are the main
cytotoxic killer cells of the human body to eliminate pathogen-infected or
tumorigenic cells (i.e. target cells). Once a NK or T cell has identified a
target cell, they form a tight contact zone, the immunological synapse (IS).
One then observes repolarization of the cell involving the rotation of the
microtubule (MT) cytoskeleton and a movement of the microtubule organizing
center (MTOC) to a position that is just underneath the plasma membrane at the
IS. Concomitantly a massive relocation of organelles attached to MTs is
observed, including the Golgi apparatus, lytic granules and mitochondria. Since
the mechanism of this relocation is still elusive we devise a theoretical model
for the molecular motor driven motion of the MT cytoskeleton confined between
membrane and nucleus. We analyze scenarios currently discussed in the
literature, the cortical sliding and the capture-shrinkage mechanisms, and
compare quantitative predictions about the spatio-temporal evolution of MTOC
position and MT cytoskeleton morphology with experiments. The model predicts
the experimentally observed biphasic nature of the process due to an interplay
between MT cytoskeleton geometry and motor forces and confirms the dominance of
the capture-shrinkage over the cortical sliding mechanism when MTOC and IS are
initially diametrically opposed. We also find that the two mechanisms act
synergistically, reducing the resources necessary for repositioning. The
localization of dyneins in the pSMAC facilitates their interaction with the
MTs. Our model also opens a way to infer details of the dynein distribution
from the experimentally observed features of the MT cytoskeleton dynamics. In a
subsequent publication, we will address the issue of general initial
configurations and situations in which the T cell established two immunological
synapses.


  Chronic Myeloid Leukaemia (CML) is a blood-derived proliferative disorder,
which is highly associated to a translocation of chromosomes 9 and 22 or the
creation of Philadelphia chromosome Ph(+) cases, inducing the synthesis of a
chimeric fusion protein, namely BCR-ABL1 (Breakpoint Cluster Region-Abelson 1
chimeric protein), which is known for driving the pathophysiology of the
disease, however variants of CML are also recognized as CML Ph(-), these
nonetheless account for a small percentage of the overall CML patients; posing
thus the question whether BCR-ABL1 fusion protein is required for the whole of
the pathophysiology of CML. Hereof, through a stochastic description, a
discrete time discrete Markov chain depicts the various protein-protein
interactions of BCR-ABL1 to better understand signalling pathways and
time-dependent evolution of these pathways, as well as to provide prospective
therapeutic protein targets to improve both the specificity of the treatment
and the life-expectancy of the patients.


  There were 27 novel coronavirus pneumonia cases found in Wuhan, China in
December 2019, named as 2019-nCoV temporarily and COVID-19 formally by WHO on
11 February, 2020. In December 2019 and January 2020, COVID-19 has spread in
large scale among the population, which brought terrible disaster to the life
and property of the Chinese people. In this paper, we will first analyze the
feature and pattern of the virus transmission, and discuss the key impact
factors and uncontrollable factors of epidemic transmission based on public
data. Then the virus transmission can be modelled and used for the inflexion
and extinction period of epidemic development so as to provide theoretical
support for the Chinese government in the decision-making of epidemic
prevention and recovery of economic production. Further, this paper
demonstrates the effectiveness of the prevention methods taken by the Chinese
government such as multi-level administrative region isolation. It is of great
importance and practical significance for the world to deal with public health
emergencies.


  We have applied a computational strategy, based on the synergy of virtual
screening, docking and molecular dynamics techniques, aimed at identifying
possible lead compounds for the non-covalent inhibition of the main protease
3CL-pro of the SARS-Cov2 Coronavirus. Based on the recently resolved 6LU7 PDB
structure, ligands were generated using a multimodal structure-based design and
then optimally docked to the 6LU7 monomer. Docking calculations show that
ligand-binding is strikingly similar in SARS-CoV and SARS-CoV2 main proteases,
irrespectively of the protonation state of the catalytic CYS-HIS dyad. The most
potent docked ligands are found to share a common binding pattern with aromatic
moieties connected by rotatable bonds in a pseudo-linear arrangement. Molecular
dynamics calculations fully confirm the stability in the 3CL-pro binding pocket
of the most potent binder identified by docking, namely a
chlorophenyl-pyridyl-carboxamide derivative.


  Highly Pathogenic Avian Influenza A H5N6 is a mutated virus of Influenza A
H5N1 and a new emerging infection that recently caused an outbreak in the
Philippines. The 2017 H5N6 outbreak resulted in a depopulation of 667,184
domestic birds. In this study, we incorporate half-saturated incidence in our
mathematical models and investigate three intervention strategies against H5N6:
isolation with treatment, vaccination and modified culling. We determine the
direction of the bifurcation when $\mathcal{R}_0 = 1$ and show that all the
models exhibit forward bifurcation. We administer optimal control and perform
numerical simulations to compare the consequences and implementation cost of
utilizing different intervention strategies in the poultry population. Despite
the challenges of applying each control strategy, we show that culling both
infected and susceptible birds is a better control strategy in prohibiting an
outbreak and avoiding further recurrence of the infection from the population
compared to confinement and vaccination.


  Environmental changes greatly influence the evolution of populations. Here,
we study the dynamics of a population of two strains, one growing slightly
faster than the other, competing for resources in a time-varying binary
environment modeled by a carrying capacity switching either randomly or
periodically between states of abundance and scarcity. The population dynamics
is characterized by demographic noise (birth and death events) coupled to a
varying environment. We elucidate the similarities and differences of the
evolution subject to a stochastically- and periodically-varying environment.
Importantly, the population size distribution is generally found to be broader
under intermediate and fast random switching than under periodic variations,
which results in markedly different asymptotic behaviors between the fixation
probability of random and periodic switching. We also determine the detailed
conditions under which the fixation probability of the slow strain is maximal.


  Chronic fetal hypoxia and infection are examples of adverse conditions during
complicated pregnancy, which impact cardiac myogenesis and increase the
lifetime risk of heart disease. However, the effects that chronic hypoxic or
inflammatory environments exert on cardiac pacemaker cells are poorly
understood. Here, we review the current evidence and novel avenues of
bench-to-bed research in this field of perinatal cardiogenesis as well as its
translational significance for early detection of future risk for
cardiovascular disease.


  Genome assembly asks to reconstruct an unknown string from many shorter
substrings of it. Even though it is one of the key problems in Bioinformatics,
it is generally lacking major theoretical advances. Its hardness stems both
from practical issues (size and errors of real data), and from the fact that
problem formulations inherently admit multiple solutions. Given these, at their
core, most state-of-the-art assemblers are based on finding non-branching paths
(unitigs) in an assembly graph. If one defines a genome assembly solution as a
closed arc-covering walk of the graph, then unitigs appear in all solutions,
being thus safe partial solutions. All all such safe walks were recently
characterized as omnitigs, leading to the first safe and complete genome
assembly algorithm. Even if omnitig finding was improved to quadratic time, it
remained open whether the crucial linear-time feature of finding unitigs can be
attained with omnitigs.
  We describe a surprising $O(m)$-time algorithm to identify all maximal
omnitigs of a graph with $n$ nodes and $m$ arcs, notwithstanding the existence
of families of graphs with $\Theta(mn)$ total maximal omnitig size. This is
based on the discovery of a family of walks (macrotigs) with the property that
all the non-trivial omnitigs are univocal extensions of subwalks of a macrotig,
with two consequences: (1) A linear-time output-sensitive algorithm enumerating
all maximal omnitigs. (2) A compact $O(m)$ representation of all maximal
omnitigs, which allows, e.g., for $O(m)$-time computation of various statistics
on them.
  Our results close a long-standing theoretical question inspired by practical
genome assemblers, originating with the use of unitigs in 1995. We envision our
results to be at the core of a reverse transfer from theory to practical and
complete genome assembly programs, as has been the case for other key
Bioinformatics problems.


  In the framework of population dynamics, the basic reproduction number R_0
is, by definition, the expected number of offspring that an individual has
during its lifetime. In constant and time periodic environments it is
calculated as the spectral radius of the so-called next-generation operator. In
continuously structured populations defined in a Banach lattice X with
concentrated states at birth one cannot define the next-generation operator in
X. In the present paper we present an approach to compute the basic
reproduction number of such models as the limit of the basic reproduction
number of a sequence of models for which R_0 can be computed as the spectral
radius of the next-generation operator. We apply these results to some
examples: the (classical) size-dependent model, a size structured cell
population model, a size structured model with diffusion in structure space
(under some particular assumptions) and a (physiological) age-structured model
with diffusion in structure space.


  We combine stochastic thermodynamics, large deviation theory, and information
theory to derive fundamental limits on the accuracy with which single cell
receptors can estimate external concentrations. As expected, if estimation is
performed by an ideal observer of the entire trajectory of receptor states,
then no energy consuming non-equilibrium receptor that can be divided into
bound and unbound states can outperform an equilibrium two-state receptor.
However, when estimation is performed by a simple observer that measures the
fraction of time the receptor is bound, we derive a fundamental limit on the
accuracy of general nonequilibrium receptors as a function of energy
consumption. We further derive and exploit explicit formulas to numerically
estimate a Pareto-optimal tradeoff between accuracy and energy. We find this
tradeoff can be achieved by nonuniform ring receptors with a number of states
that necessarily increases with energy. Our results yield a novel thermodynamic
uncertainty relation for the time a physical system spends in a pool of states,
and generalize the classic 1977 Berg-Purcell limit on cellular sensing along
multiple dimensions.


  Background: Single nucleotide variants (SNVs) are detected as different
distributions of DNA samples of distinct types of cancer patients. Even though,
it is an exacting task to select the appropriate method to identify cancer to
the greatest extent of SNVs. Results: In this paper, we proposed a biomarker
concept based on SNV patterns in different feature dimensions. Raw dataset
(2761 samples) consisting of twelve different cancers was obtained from TCGA
(The Cancer Genome Atlas). After preliminary screening of 562,321 DNA mutation
sites in the samples, the mutation sites were extracted and characterized by
cancer types in six different SNV feature dimensions. In this study, we found
that the extracted features showed similar distribution in the cluster center
of the disease type of the samples. After the initial processing of the raw
data, the sample was more focused on the subtype distribution of the cancer or
the cancer at the SNV level. We used k-nearest neighbors (KNN) to classify the
extracted features and Leave-One-Out cross verified them. The accuracy of
classifying is stable at around 97% and reached 97.43% at the highest. During
the validation phase, we found validated oncogenes in the loci of the features
with the highest importance among nine cancers. Conclusions: In summary, the
samples showed consistent patterns according to the cancer in which it belongs.
It is feasible to classify the cancer of the sample by the distribution of
different dimensions of the SNVs and has a high accuracy. And has potential
implications for the discovery of cancer-causing genes.


  We consider a system of two competing populations in two-dimensional
heterogeneous environments. The populations are assumed to move horizontally
and vertically with different probabilities, but are otherwise identical. We
regard these probabilities as dispersal strategies. We show that the
evolutionarily stable strategies are to move in one direction only. Our results
predict that it is more beneficial for the species to choose the direction with
smaller variation in the resource distribution. This finding seems to be in
agreement with the classical results of Hasting [15] and Dockery et al. [11]
for the evolution of slow dispersal, i.e. random diffusion is selected against
in spatially heterogeneous environments. These conclusions also suggest that
broader dispersal strategies should be considered regarding the movement in
heterogeneous habitats.


  Studying the conformations involved in the dimerization of cadherins is
highly relevant to understand the development of tissue and its failure, which
is associated with tumors and metastases. Experimental techniques, like X-ray
crystallography, can usually report only the most stable conformations, missing
minority states that could nonetheless be important for the recognition
mechanism. Computer simulations could be a valid complement to the experimental
approach. However, standard all-atom protein models in explicit solvent are
computationally too demanding to search thoroughly the conformational space of
multiple chains composed of several hundreds of amino acids. To reach this
goal, we resorted to a coarse-grained model in implicit solvent. The standard
problem with this kind of models is to find a realistic potential to describe
their interactions. We used coevolutionary information from cadherin
alignments, corrected by a statistical potential, to build an interaction
potential which is agnostic of the experimental conformations of the protein.
Using this model, we explored the conformational space of multi-chain systems
and validated the results comparing with experimental data. We identified
dimeric conformations that are sequence-specific and that can be useful to
rationalize the mechanism of recognition between cadherins.


  Fisher's geometric model describes biological fitness landscapes by combining
a linear map from the discrete space of genotypes to an $n$-dimensional
Euclidean phenotype space with a nonlinear, single-peaked phenotype-fitness
map. Genotypes are represented by binary sequences of length $L$, and the
phenotypic effects of mutations at different sites are represented by $L$
random vectors drawn from an isotropic Gaussian distribution. Recent work has
shown that the interplay between the genotypic and phenotypic levels gives rise
to a range of different landscape topographies that can be characterised by the
number of local fitness maxima. Extending our previous study of the mean number
of local maxima, here we focus on the distribution of the number of maxima when
the limit $L \to \infty$ is taken at finite $n$. We identify the typical scale
of the number of maxima for general $n$, and determine the full scaled
probability density and two point correlation function of maxima for the
one-dimensional case. We also elaborate on the close relation of the model to
the anti-ferromagnetic Hopfield model with $n$ random continuous pattern
vectors, and show that many of our results carry over to this setting. More
generally, we expect that our analysis can help to elucidate the fluctuation
structure of metastable states in various spin glass problems.


  Both molecular mechanical and quantum mechanical calculations play an
important role in describing the behavior and structure of molecules. In this
work, we compare for the same peptide systems the results obtained from folding
molecular dynamics simulations with previously reported results from quantum
mechanical calculations. More specifically, three molecular dynamics
simulations of 5 $\mu$s each in explicit water solvent were carried out for
three Asn-Gly-containing heptapeptides, in order to study their folding and
dynamics. Previous data, based on quantum mechanical calculations and the DFT
methods have shown that these peptides adopt $\beta$-turn structures in aqueous
solution, with type I' $\beta$-turn being the most preferred motif. The results
from our analyses indicate that for the given system the two methods diverge in
their predictions. The possibility of a force field-dependent deficiency is
examined as a possible source of the observed discrepancy.


  Lies can have a negating impact on governments, companies, and the society as
a whole. Understanding the dynamics of lying is therefore of crucial importance
across different fields of research. While lying has been studied before in
well-mixed populations, it is a fact that real interactions are rarely
well-mixed. Indeed, they are usually structured and thus best described by
networks. Here we therefore use the Monte Carlo method to study the evolution
of lying in the sender-receiver game in a one-parameter family of networks,
systematically covering complete networks, small-world networks, and
one-dimensional rings. We show that lies which benefit the sender at a cost to
the receiver, the so-called black lies, are less likely to proliferate on
networks than they do in well-mixed populations. Honesty is thus more likely to
evolve, but only when the benefit for the sender is smaller than the cost for
the receiver. Moreover, this effect is particularly strong in small-world
networks, but less so in the one-dimensional ring. For lies that favor the
receiver at a cost to the sender, the so-called altruistic white lies, we show
that honesty is also more likely to evolve than it is in well-mixed
populations. But contrary to black lies, this effect is more expressed in the
one-dimensional ring, whereas in small-world networks it is present only when
the cost to the sender is greater than the benefit for the receiver. Lastly,
for lies that benefit both the sender and the receiver, the so-called Pareto
white lies, we show that the network structure actually favors the evolution of
lying, but this only when the benefit for the sender is slightly greater than
the benefit for the receiver. In this case again the small-world topology acts
as an amplifier of the effect, while other network topologies fail to do the
same.


  We define the Sampled Moran Genealogy Process, a continuous-time Markov
process on the space of genealogies with the demography of the classical Moran
process, sampled through time. To do so, we begin by defining the Moran
Genealogy Process using a novel representation. We then extend this process to
include sampling through time. We derive exact conditional and marginal
probability distributions for the sampled process under a stationarity
assumption, and an exact expression for the likelihood of any sequence of
genealogies it generates. This leads to some interesting observations pertinent
to existing phylodynamic methods in the literature. Throughout, our proofs are
original and make use of strictly forward-in-time calculations and are exact
for all population sizes and sampling processes.


  Ross-Macdonald models are the building blocks of most vector-borne disease
models. Even for the same disease, different authors use different model
formulations, but a study of the dynamical consequences of assuming different
hypotheses is missing. In this work we present different formulations of the
basic Ross-Macdonald model together with a careful discussion of the
assumptions behind each model. The most general model presented is an agent
based model for which arbitrary distributions for latency and infectious
periods for both, host and vectors, is considered. At population level we also
developed a deterministic Volterra integral equations model for which also
arbitrary distributions in the waiting times are included. We compare the model
solutions using different distributions for the infectious and latency periods
using statistics, like the epidemic peak, or epidemic final size, to
characterize the epidemic curves. The basic reproduction number ($R_0$) for
each formulation is computed and compared with empirical estimations obtained
with the agent based models. The importance of considering realistic
distributions for the latent and infectious periods is highlighted and
discussed. We also show that seasonality is a key driver of vector-borne
disease dynamics shaping the epidemic curve and its duration.


  During cancer progression, malignant cells accumulate somatic mutations that
can lead to genetic aberrations. In particular, evolutionary events akin to
segmental duplications or deletions can alter the copy-number profile (CNP) of
a set of genes in a genome. Our aim is to compute the evolutionary distance
between two cells for which only CNPs are known. This asks for the minimum
number of segmental amplifications and deletions to turn one CNP into another.
This was recently formalized into a model where each event is assumed to alter
a copy-number by $1$ or $-1$, even though these events can affect large
portions of a chromosome. We propose a general cost framework where an event
can modify the copy-number of a gene by larger amounts. We show that any cost
scheme that allows segmental deletions of arbitrary length makes computing the
distance strongly NP-hard. We then devise a factor $2$ approximation algorithm
for the problem when copy-numbers are non-zero and provide an implementation
called \textsf{cnp2cnp}. We evaluate our approach experimentally by
reconstructing simulated cancer phylogenies from the pairwise distances
inferred by \textsf{cnp2cnp} and compare it against two other alternatives,
namely the \textsf{MEDICC} distance and the Euclidean distance. The
experimental results show that our distance yields more accurate phylogenies on
average than these alternatives if the given CNPs are error-free, but that the
\textsf{MEDICC} distance is slightly more robust against error in the data. In
all cases, our experiments show that either our approach or the \textsf{MEDICC}
approach should preferred over the Euclidean distance.


  The spatial distributions of diverse facilities are often understood in terms
of the optimization of the commute distance or the economic profit.
Incorporating more general objective functions into such optimization framework
may be useful, helping the policy decisions to meet various social and economic
demands. As an example, we consider how hospitals should be distributed to
minimize the total fatalities of tuberculosis (TB). The empirical data of Korea
shows that the fatality rate of TB in a district decreases with the areal
density of hospitals, implying their correlation and the possibility of
reducing the nationwide fatalities by adjusting the hospital distribution
across districts. Approximating the fatality rate by the probability of a
patient not to visit a hospital in her/his residential district for the
duration period of TB and evaluating the latter probability in the random-walk
framework, we obtain the fatality rate as an exponential function of the
hospital density with a characteristic constant related to each district's
effective lattice constant estimable empirically. This leads us to the optimal
hospital distribution which finds the hospital density in a district to be a
logarithmic function of the rescaled patient density. The total fatalities is
reduced by 13\% with this optimum. The current hospital density deviates from
the optimized one in different manners from district to district, which is
analyzed in the proposed model framework. The assumptions and limitations of
our study are also discussed.


  Oligonucleotide-based agents have the potential to treat or cure almost any
disease, and are one of the key therapeutic drug classes of the future.
Bioconjugated oligonucleotides, a subset of this class, are emerging from basic
research and being successfully translated to the clinic. In this review, we
first briefly describe two approaches for inhibiting specific genes using
oligonucleotides -- antisense DNA (ASO) and RNA interference (RNAi) -- followed
by a discussion on delivery to cells. We then summarize and analyze recent
developments in bioconjugated oligonucleotides including those possessing
GalNAc, cell penetrating peptides, $\alpha$-tocopherol, aptamers, antibodies,
cholesterol, squalene, fatty acids, or nucleolipids. These novel conjugates
provide a means to enhance tissue targeting, cell internalization, endosomal
escape, target binding specificity, resistance to nucleases, and more. We next
describe those bioconjugated oligonucleotides approved for patient use or in
clinical trials. Finally, we summarize the state of the field, describe current
limitations, and discuss future prospects. Biocon-jugation chemistry is at the
centerpiece of this therapeutic oligonucleotide revolution, and significant
opportunities exist for development of new modification chemistries, for
mechanistic studies at the chemical-biology interface, and for translating such
agents to the clinic.


  Recent advances in distance-based protein folding have led to a paradigm
shift in protein structure prediction. Through sufficiently precise estimation
of the inter-residue distance matrix for a protein sequence, it is now feasible
to predict the correct folds for new proteins much more accurately than ever
before. Despite the exciting progress, a dedicated visualization system that
can dynamically capture the distance-based folding process is still lacking.
Most molecular visualizers typically provide only a static view of a folded
protein conformation, but do not capture the folding process. Even among the
selected few graphical interfaces that do adopt a dynamic perspective, none of
them are distance-based. Here we present PolyFold, an interactive visual
simulator for dynamically capturing the distance-based protein folding process
through real-time rendering of a distance matrix and its compatible spatial
conformation as it folds in an intuitive and easy-to-use interface. PolyFold
integrates highly convergent stochastic optimization algorithms with on-demand
customizations and interactive manipulations to maximally satisfy the geometric
constraints imposed by a distance matrix. PolyFold is capable of simulating the
complex process of protein folding even on modest personal computers, thus
making it accessible to the general public for fostering citizen science. Open
source code of PolyFold is freely available for download at
https://github.com/Bhattacharya-Lab/PolyFold. It is implemented in
cross-platform Java and binary executables are available for macOS, Linux, and
Windows.


  Signaling pathways are responsible for the regulation of cell processes, such
as monitoring the external environment, transmitting information across
membranes, and making cell fate decisions. Given the increasing amount of
biological data available and the recent discoveries showing that many diseases
are related to the disruption of cellular signal transduction cascades, in
silico discovery of signaling pathways in cell biology has become an active
research topic in past years. However, reconstruction of signaling pathways
remains a challenge mainly because of the need for systematic approaches for
predicting causal relationships, like edge direction and activation/inhibition
among interacting proteins in the signal flow. We propose an approach for
predicting signaling pathways that integrates protein interactions, gene
expression, phenotypes, and protein complex information. Our method first finds
candidate pathways using a directed-edge-based algorithm and then defines a
graph model to include causal activation relationships among proteins, in
candidate pathways using cell cycle gene expression and phenotypes to infer
consistent pathways in yeast. Then, we incorporate protein complex coverage
information for deciding on the final predicted signaling pathways. We show
that our approach improves the predictive results of the state of the art using
different ranking metrics.


  Bimolecular binding rate constants are often used to describe the association
of large molecules, such as proteins. In this paper, we analyze a model for
such binding rates that includes the fact that pairs of molecules can bind only
in certain orientations. The model considers two spherical molecules, each with
an arbitrary number of small binding sites on their surface, and the two
molecules bind if and only if their binding sites come into contact (such
molecules are often called "patchy particles" in the biochemistry literature).
The molecules undergo translational and rotational diffusion, and the binding
sites are allowed to diffuse on their surfaces. Mathematically, the model takes
the form of a high-dimensional, anisotropic diffusion equation with mixed
boundary conditions. We apply matched asymptotic analysis to derive the
bimolecular binding rate in the limit of small, well-separated binding sites.
The resulting binding rate formula involves a factor that depends on the
electrostatic capacitance of a certain four-dimensional region embedded in five
dimensions. We compute this factor numerically by modifying a recent kinetic
Monte Carlo algorithm. We then apply a quasi chemical formalism to obtain a
simple analytical approximation for this factor and find a binding rate formula
that includes the effects of binding site competition/saturation. We verify our
results by numerical simulation.


  Within cells, vesicles and proteins are actively transported several
micrometers along the cytoskeletal filaments. The transport along microtubules
is propelled by dynein and kinesin motors, which carry the cargo in opposite
directions. Bidirectional intracellular transport is performed with great
efficiency, even under strong confinement, as for example in the axon. For this
kind of transport system, one would expect generically cluster formation. In
this work, we discuss the effect of the recently observed self-enhanced
binding-affinity along the kinesin trajectories on the MT. We introduce a
stochastic lattice-gas model, where the enhanced binding affinity is realized
via a floor-field. From Monte Carlo simulations and a mean-field analysis we
show that this mechanism can lead to self-organized symmetry-breaking and
lane-formation which indeed leads to efficient bidirectional transport in
narrow environments.


  Based on the SEIR model and the modeling of urban transportation networks, a
general-purpose simulator for the spread of epidemics in Chinese cities is
built. The Chinese public transportation system between over 340
prefectural-level cities is modeled as a multi-layer bi-partite network, with
layers representing different means of transportation (airlines, railways, sail
routes and buses), and nodes divided into two categories (central cities,
peripheral cities). At each city, an open-system SEIR model tracks the local
spread of the disease, with population in- and out-flow exchanging with the
overlying transportation network. The model accounts for (1) different
transmissivities of the epidemic on different transportation media, (2) the
transit of inbound flow at cities, (3) cross-infection on public transportation
vehicles due to path overlap, and the realistic considerations that (4) the
infected population are not entering public transportation and (5) the
recovered population are not subject to repeated infections. The model could be
used to simulate the city-level spread in China (and potentially other
countries) of an arbitrary epidemic, characterized by its basic reproduction
number, incubation period, infection period and zoonotic force, originated from
any Chinese prefectural-level city(s), during the period before effective
government interventions are implemented. Flowmaps are input into the system to
trigger inter-city dynamics, assuming different flow strength, determined from
empirical observation, within/between the bi-partite divisions of nodes. The
model is used to simulate the 2019 Coronavirus epidemic in Wuhan; it shows that
the framework is robust and reliable, and simulated results match public
city-level datasets to an extraordinary extent.


  We model the COVID-19 coronavirus epidemic in China. We use early reported
case data to predict the cumulative number of reported cases to a final size.
The key features of our model are the timing of implementation of major public
policies restricting social movement, the identification and isolation of
unreported cases, and the impact of asymptomatic infectious cases.


  Many vector-borne disease epidemic models neglect the fact that in modern
human civilization, social awareness as well as self-defence system are
overwhelming against advanced propagation of the disease. News are becoming
more effortlessly accessible through social media and mobile apps, while
apparatuses for disease prevention are inclined to be more abundant and
affordable. Here we study a simple host-vector model in which media-triggered
social awareness and seasonality in vector breeding are taken into account.
There appears a certain threshold indicating the alarming outbreak; the number
of infective human individuals above which shall actuate the self-defence
system for the susceptible subpopulation. A model where the infection rate
revolves in the likelihood of poverty, reluctancy, tiresomeness, perceiving the
disease as being easily curable, absence of medical access, and overwhelming
hungrier vectors is proposed. Further discoveries are made from undertaking
disparate time scales between human and vector population dynamics. The
resulting slow-fast system discloses notable dynamics in which solution
trajectories confine to the slow manifold and critical manifold, before finally
ending up at equilibria. How coinciding the slow manifold with the critical
manifold enhances periodic forcing is also studied. The finding on hysteresis
loops gives insights of how defining alarming outbreak critically perturbs the
basic reproductive number, which later helps keep the incidence cycle on small
magnitudes.


  Honest signals and cues have been observed as part of interspecific and
intraspecific communication among animals. Recent theories suggest that
existing signaling systems have evolved through natural selection imposed by
predators. Honest signaling in the interspecific communication can provide
insight into the evolution of anti-predation techniques. In this work, we
introduce a deterministic three-stage, two-species predator-prey model, which
modulates the impact of honest signals and cues on the interacting populations.
The model is built from a set of first principles originated from signaling and
social learning theory in which the response of predators to transmitted honest
signals or cues is determined. The predators then use the signals to decide
whether to pursue the attack or save their energy for an easier catch. Other
members from the prey population that are not familiar with signaling their
fitness observe and learn the technique. Our numerical bifurcation analysis
indicates that increasing the predator's search rate and the corresponding
assimilation efficiency gives a journey from predator-prey abundance and
scarcity, a stable transient cycle between persistence and near-extinction, a
homoclinic orbit pointing towards extinction, and ultimately, a quasi-periodic
orbit. A similar discovery is met under the increment of the prey's intrinsic
birth rate and carrying capacity. When both parameters are of sufficiently
large magnitudes, the separator between honest signal and cue takes the similar
journey from a stable equilibrium to a quasi-periodic orbit as it increases. In
the context of modeling, we conclude that under prey abundance, transmitting
error-free honest signals leads to not only a stable but also more predictable
predator-prey dynamics.


  In this paper we investigate the effect of moral suasion on ingroup
favouritism. We report a well-powered, pre-registered, two-stage 2x2
mixed-design experiment. In the first stage, groups are formed on the basis of
how participants answer to a set of questions, concerning non-morally relevant
issues in one treatment (assorting on non-moral preferences), and morally
relevant issues in another treatment (assorting on moral preferences). In the
second stage, participants choose how to split a given amount of money between
participants of their own group and participants of the other group, first in
the baseline setting and then in a setting where they are told to do what they
believe to be morally right (moral suasion). Our main results are: (i) in the
baseline, participants tend to favour their own group to a greater extent when
groups are assorted according to moral preferences, compared to when they are
assorted according to non-moral preferences; (ii) the net effect of moral
suasion is to decrease ingroup favouritism, but there is also a non-negligible
proportion of participants for whom moral suasion increases ingroup
favouritism; (iii) the effect of moral suasion is substantially stable across
group assorting and four pre-registered individual characteristics (gender,
political orientation, religiosity, pro-life vs pro-choice ethical
convictions).


  Tree shape statistics are important for investigating evolutionary mechanisms
mediating phylogenetic trees. As a step towards bridging shape statistics
between rooted and unrooted trees, we present a comparison study on two subtree
statistics known as numbers of cherries and pitchforks for the proportional to
distinguishable arrangements (PDA) and the Yule-Harding-Kingman (YHK) models.
Based on recursive formulas on the joint distribution of the number of cherries
and that of pitchforks, it is shown that cherry distributions are log-concave
for both rooted and unrooted trees under these two models. Furthermore, the
mean number of cherries and that of pitchforks for unrooted trees converge
respectively to those for rooted trees under the YHK model while there exists a
limiting gap of 1/4 for the PDA model. Finally, the total variation distances
between the cherry distributions of rooted and those of unrooted trees converge
for both models. Our results indicate that caution is required for conducting
statistical analysis for tree shapes involving both rooted and unrooted trees.


  Motivated by problems from Chemical Reaction Network Theory, we investigate
whether steady state ideals of reversible reaction networks are generated by
binomials. We take an algebraic approach considering, besides concentrations of
species, also rate constants as indeterminates. This leads us to the concept of
unconditional binomiality, meaning binomiality for all values of the rate
constants. This concept is different from conditional binomiality that applies
when rate constant values or relations among rate constants are given. We start
by representing the generators of a steady state ideal as sums of binomials,
which yields a corresponding coefficient matrix. On these grounds we propose an
efficient algorithm for detecting unconditional binomiality. That algorithm
uses exclusively elementary column and row operations on the coefficient
matrix. We prove asymptotic worst case upper bounds on the time complexity of
our algorithm. Furthermore, we experimentally compare its performance with
other existing methods.


  We introduce a game inspired by the challenges of disease management in
livestock farming and the transmission of endemic disease through a trade
network. Success in this game comes from balancing the cost of buying new stock
with the risk that it will be carrying some disease. When players follow a
simple memory-based strategy we observe a spontaneous separation into two
groups corresponding to players with relatively high, or low, levels of
infection. By modelling the dynamics of both the disease and the formation and
breaking of trade relationships, we derive the conditions for which this
separation occurs as a function of the transmission rate and the threshold
level of acceptable disease for each player. When interactions in the game are
restricted to players that neighbour each other in a small-world network,
players tend to have similar levels of infection as their neighbours. We
conclude that success in economic-epidemiological systems can originate from
misfortune and geographical circumstances as well as by innate differences in
personal attitudes towards risk.


  Although geographic isolation has been shown to play a key role in promoting
reproductive isolation, it is now believed that speciation can also happen in
sympatry and with considerable gene flow. Here we present a model of sympatric
speciation based on assortative mating that does not require a genetic
threshold for reproduction, i.e., that does not directly associate genetic
differences between individuals with reproductive incompatibilities. In the
model individuals mate with the most similar partner in their pool of potential
mates, irrespective of how dissimilar it might be. We show that assortativity
alone can lead to the formation of clusters of genetically similar individuals.
The absence of a minimal genetic similarity for mating implies the constant
generation of hybrids and brings up the old problem of species definition.
Here, we define species based on clustering of genetically similar individuals
but allowing genetic flow among different species. We show that the results
obtained with the present model are in good agreement with empirical data, in
which different species can still reproduce and generate hybrids.


  We study the effect that disturbances in the ecological landscape exert on
the spatial distribution of a population that evolves according to the nonlocal
FKPP equation. Using both numerical and analytical techniques, we explore the
three types of stationary profiles that can develop near abrupt spatial
variations in the environmental conditions vital for population growth:
sustained oscillations, decaying oscillations and exponential relaxation
towards a flat profile. We relate the wavelength and decay length of the
wrinkles to the shape of the interaction kernel. Since these spatial structures
encode details about the interaction among individuals, we discuss how
heterogeneities reveal information that would be hidden in a flat landscape.


  Massively parallel sequencing techniques have revolutionized biological and
medical sciences by providing unprecedented insight into the genomes of humans,
animals, and microbes. Modern sequencing platforms generate enormous amounts of
genomic data in the form of nucleotide sequences or reads. Aligning reads onto
reference genomes enables the identification of individual-specific genetic
variants and is an essential step of the majority of genomic analysis
pipelines. Aligned reads are essential for answering important biological
questions, such as detecting mutations driving various human diseases and
complex traits as well as identifying species present in metagenomic samples.
The read alignment problem is extremely challenging due to the large size of
analyzed datasets and numerous technological limitations of sequencing
platforms, and researchers have developed novel bioinformatics algorithms to
tackle these difficulties. Importantly, computational algorithms have evolved
and diversified in accordance with technological advances, leading to todays
diverse array of bioinformatics tools. Our review provides a survey of
algorithmic foundations and methodologies across 107 alignment methods
published between 1988 and 2020, for both short and long reads. We provide
rigorous experimental evaluation of 11 read aligners to demonstrate the effect
of these underlying algorithms on speed and efficiency of read aligners. We
separately discuss how longer read lengths produce unique advantages and
limitations to read alignment techniques. We also discuss how general alignment
algorithms have been tailored to the specific needs of various domains in
biology, including whole transcriptome, adaptive immune repertoire, and human
microbiome studies.


  In this paper, we conduct mathematical and numerical analyses to address the
following crucial questions for COVID-19: (Q1) Is it possible to contain
COVID-19? (Q2) When will be the peak and the end of the epidemic? (Q3) How do
the asymptomatic infections affect the spread of disease? (Q4) What is the
ratio of the population that needs to be infected to achieve herd immunity?
(Q5) How effective are the social distancing approaches? (Q6) What is the ratio
of the population infected in the long run? For (Q1) and (Q2), we propose a
time-dependent susceptible-infected-recovered (SIR) model that tracks 2 time
series: (i) the transmission rate at time t and (ii) the recovering rate at
time t. Such an approach is more adaptive than traditional static SIR models
and more robust than direct estimation methods. Using the data provided by
China, we show that the one-day prediction errors for the numbers of confirmed
cases are almost in 3%, and the total number of confirmed cases is precisely
predicted. Also, the turning point, defined as the day that the transmission
rate is less than the recovering rate can be accurately predicted. After that
day, the basic reproduction number $R_0$ is less than 1. For (Q3), we extend
our SIR model by considering 2 types of infected persons: detectable and
undetectable infected persons. Whether there is an outbreak in such a model is
characterized by the spectral radius of a 2 by 2 matrix that is closely related
to $R_0$. For (Q4), we show that herd immunity can be achieved after at least
1-1/$R_0$ fraction of individuals being infected. For (Q5) and (Q6), we analyze
the independent cascade (IC) model for disease propagation in a configuration
random graph. By relating the propagation probabilities in the IC model to the
transmission rates and recovering rates in the SIR model, we show 2 approaches
of social distancing that can lead to a reduction of $R_0$.


  This paper focuses on studying and understanding of stochastic dynamics in
population composition when the population is subject to rumor spreading. We
undertake the study by first developing an individual
Susceptible-Exposed-Infectious-Removed (iSEIR) model, an extension of the SEIR
model, for summarizing rumor-spreading behaviors of interacting groups in the
population. With this iSEIR model, the interacting groups may be regarded as
nodes in a multiplex network. Then various properties of the dynamic behaviors
of the interacting groups in rumor spreading can be drawn from samples of the
multiplex network. The samples are simulated based on the iSEIR model with
different settings in terms of population scale, population distribution and
transfer rate. Results from the simulation study show that effective control of
rumor spreading in the multiplex network entails an efficient management on
information flow, which may be achieved by setting appropriate immunization and
spreading thresholds in individual behavior dynamics. Under the proposed iSEIR
model we also have derived a steady-state result, named the "supersaturation
phenomenon", when the rumor spreading process becomes equilibrium, which may
help us to make the optimal or better control of information flow in the
practice.


  Motivation: The coronavirus disease 2019 (COVID-19) caused by a new type of
coronavirus has been emerging from China and led to thousands of death globally
since December 2019. Despite many groups have engaged in studying the newly
emerged virus and searching for the treatment of COVID-19, the understanding of
the COVID-19 target-ligand interactions represents a key chal-lenge. Herein, we
introduce COVID-19 Docking Server, a web server that predicts the binding modes
between COVID-19 targets and the ligands including small molecules, peptides
and anti-bodies. Results: Structures of proteins involved in the virus life
cycle were collected or constructed based on the homologs of coronavirus, and
prepared ready for docking. The meta platform provides a free and interactive
tool for the prediction of COVID-19 target-ligand interactions and following
drug discovery for COVID-19.


  We proposed a Monte-Carlo method to estimate temporal reproduction number
without complete information about symptom onsets of all cases. Province-level
analysis demonstrated the huge success of Chinese control measures on COVID-19,
that is, provinces' reproduction numbers quickly decrease to <1 by just one
week after taking actions.


  First we shortly review the different kinds of network modelling methods for
systems biology with an emphasis on the different subtypes of logical models,
which we review in more detail. Then we show the advantages of Boolean networks
models over more mechanistic modelling types like differential equation
techniques. Then follows an overlook about connections between different kinds
of models and how they can be converted to each other. We also give a short
overview about the mathematical frameworks for modelling of logical networks
and list available software packages for logical modelling. Then we give an
overview about the available standards and ontologies for storing such logical
systems biology models and their results. In the end we give a short review
about the difference between quantitative and qualitative models and describe
the mathematics that specifically deals with qualitative modelling.


  Mass spectrometry has experienced a rapid development since its first
application for protein analysis in the 1980s. While the most common use of
mass spectrometry for protein analysis is identification and quantification
workflows on peptides (digested from their parent protein), there is also a
rapidly growing use of mass spectrometry for structural proteomics. One example
is the analysis of cross-linked proteins that can give valuable structural
information, complementing the information gained by classical protein
structure determination methods, useful for integrated methods of structure
determination and modeling. For a broad and reproducible application of
cross-linking mass spectrometry a standardized representation of cross-linking
experimental results is necessary. This paper describes the developing and
release of the xlmod ontology from the HUPO-PSI. xlmod contains terms for the
description of reagents used in cross-linking experiments and of cross-linker
related chemical modifications together with their main properties relevant for
planning and performing cross-linking experiments. We also describe how xlmod
is used within the new release of HUPO-PSI-s mzIdentML data standard, for
reporting the used cross-linking reagents and results in a consistent manner.
In addition xlmod contains terms for GC-MS and LC-MS derivatization reagents
for specifying them in the upcoming mzTab-M and mzTab-L formats.


  To evaluate the effectiveness of the containment on the epidemic spreading of
the new Coronavirus disease 2019, we carry on an analysis of the time evolution
of the infection in a selected number of different Countries, by considering
well-known macroscopic growth laws, the Gompertz law, and the logistic law. We
also propose here a generalization of Gompertz law. Our data analysis permits
an evaluation of the maximum number of infected individuals. The daily data
must be compared with the obtained fits, to verify if the spreading is under
control. From our analysis it appears that the spreading reached saturation in
China, due to the strong containment policy of the national government. In
Singapore a large growth rate, recently observed, suggests the start of a new
strong spreading. For South Korea and Italy, instead, the next data on new
infections will be crucial to understand if the saturation will be reached for
lower or higher numbers of infected individuals.


  The outbreak of COVID-19 caused millions of deaths worldwide, and the number
of total infections is still rising. It is necessary to identify some
potentially effective drugs that can be used to prevent the development of
severe symptoms or even death for those infected. Fortunately, many efforts
have been made, and several effective drugs have been identified. The rapidly
increasing amount of data is of great help for training an effective and
specific deep learning model. In this study, we propose a multi-task deep
learning model for the purpose of screening commercially available and
effective inhibitors against SARS-CoV-2. First, we pretrained a model on
several heterogenous protein-ligand interaction datasets. The model achieved
competitive results on some benchmark datasets. Next, a coronavirus-specific
dataset was collected and used to fine-tune the model. Then, the fine-tuned
model was used to select commercially available drugs against SARS-CoV-2
protein targets. Overall, twenty compounds were listed as potential inhibitors.
We further explored the model interpretability and observed the predicted
important binding sites. Based on this prediction, molecular docking was also
performed to visualize the binding modes of the selected inhibitors.


  Regulating the upstream of the cytokines production could be a promising
strategy to the treatment of COVID-19. We suggest to pay more attention to the
dysregulated IFN-I production in COVID-19 and to considerate cGAS, ALK and
STING as potential therapeutic targets preventing cytokine storm. Approved
drugs like suramin and ALK inhibitors are worthy of clinical trials.


  Scientists in some fields are concerned that many, or even most, published
results are false. A high rate of false positives might arise accidentally,
from shoddy research practices. Or it might be the inevitable result of
institutional incentives that reward publication irrespective of veracity.
Recent models and discussion of scientific culture predict selection for
false-positive publications, as research labs that publish more positive
findings out-compete more diligent labs. There is widespread debate about how
scientific practices should be modified to avoid this degeneration. Some
analyses suggest that "bad science" will persist even when labs are
incentivized to undertake replication studies, and penalized for publications
that later fail to replicate. Here we develop a framework for modelling the
cultural evolution of research practices that allows labs to expend effort on
theory - enabling them, at a cost, to focus on hypotheses that are more likely
to be true on theoretical grounds. Theory restores the evolution of high effort
in laboratory practice, and it suppresses false-positive publications to a
technical minimum, even in the absence of replication. In fact, the mere
ability choose between two sets of hypotheses - one with greater chance of
being correct than the other - promotes better science than can be achieved by
having effortless access to the better set of hypotheses. Combining theory and
replication can have a synergistic effect in promoting good scientific
methodology and reducing the rate of false-positive publications. Based on our
analysis we propose four simple rules to promote good science in the face of
pressure to publish.


  Mathematical and computational modeling approaches are increasingly used as
quantitative tools in the analysis and forecasting of infectious disease
epidemics. The growing need for realism in addressing complex public health
questions is however calling for accurate models of the human contact patterns
that govern the disease transmission processes. Here we present a data-driven
approach to generate effective descriptions of population-level contact
patterns by using highly detailed macro (census) and micro (survey) data on key
socio-demographic features. We produce age-stratified contact matrices for 277
sub-national administrative regions of countries covering approximately 3.5
billion people and reflecting the high degree of cultural and societal
diversity of the focus countries. We use the derived contact matrices to model
the spread of airborne infectious diseases and show that sub-national
heterogeneities in human mixing patterns have a marked impact on epidemic
indicators such as the reproduction number and overall attack rate of epidemics
of the same etiology. The contact patterns derived here are made publicly
available as a modeling tool to study the impact of socio-economic differences
and demographic heterogeneities across populations on the epidemiology of
infectious diseases.


  We first propose a quantitative approach to detect high risk outbreaks of
independent and coinfective SIR dynamics on three empirical networks: a school,
a conference and a hospital contact network. This measurement is based on the
k-means clustering method and identifies proper samples for calculating the
mean outbreak size and the outbreak probability. Then we systematically study
the impact of different temporal correlations on high risk outbreaks over the
original and differently shuffled counterparts of each network. We observe
that, on the one hand, in the coinfection process, randomization of the
sequence of the events increases the mean outbreak size of high risk cases. On
the other hand, these correlations don't have a consistent effect on the
independent infection dynamics, and can either decrease or increase this mean.
While randomization of the daily pattern correlations has no significant effect
on the size of outbreak in either of the coinfection or independent spreading
cases. We also observer that an increase in the mean outbreak size doesn't
always coincide with an increase in the outbreak probability; therefore we
argue that merely considering the mean outbreak size of all realizations may
lead us into misestimating the outbreak risks. Our results suggest that some
sort of randomizing contacts in organization level of schools, events or
hospitals might help to suppress the spreading dynamics while the risk of an
outbreak is high.


  The condensation of biomolecules into biomolecular condensates via
liquid-liquid phase separation (LLPS) is a ubiquitous mechanism that drives
cellular organization. To enable these functions, biomolecules have evolved to
drive LLPS and facilitate partitioning into biomolecular condensates.
Determining the molecular features of proteins that encode LLPS will provide
critical insights into a plethora of biological processes. Problematically,
probing biomolecular dense phases directly is often technologically difficult
or impossible. By capitalizing on the symmetry between the conformational
behavior of biomolecules in dilute solution and dense phases, it is possible to
infer details critical to phase separation by precise measurements of the
dilute phase thus circumventing complicated characterization of dense phases.
The symmetry between dilute and dense phases is found in the size and shape of
the conformational ensemble of a biomolecule - parameters that small-angle
x-ray scattering (SAXS) is ideally suited to probe. Recent technological
advances have made it possible to accurately characterize samples of
intrinsically disordered protein regions at low enough concentration to avoid
interference from intermolecular attraction, oligomerization or aggregation,
all of which were previously roadblocks to characterizing self-assembling
proteins. Herein, we describe the pitfalls inherent to measuring such samples,
the details required for circumventing these issues and analysis methods that
place the results of SAXS measurements into the theoretical framework of LLPS.


  The phylogenetic tree of SARS-CoV-2 (nCov-19) viruses is reconstructed
according to the similarity of genome sequences. The tree topology of
Betacoronavirus is remarkably consistent with biologist's systematics. Because
the tree construction contains enough information about virus mutants, it is
suitable to study the evolutionary relationship between novel coronavirus
mutants transmitted among humans. The emergences of 14 kinds of main mutants
are studied and these strains can be classified as eight bifurcations of the
phylogenetic tree. It is found that there exist three types of virus mutations,
namely, the mutation among sub-branches of the same branch, the off-root
mutation and the root-oriented mutation between large branches of the tree.
From the point of the relation between viral mutation and host selection we
found that individuals with low immunity provide a special environment for the
positive natural selection of virus evolution. It gives a mechanism to explain
why large mutations between two distant branches generally occur in the nCov-19
phylogenetic tree. The finding is helpful to formulate strategies to control
the spread of COVID-19.


  Darwin's theory of evolution emphasized that positive selection of functional
proficiency provides the fitness that ultimately determines the structure of
life, a view that has dominated biochemical thinking of enzymes as perfectly
optimized for their specific functions. The 20th-century modern synthesis,
structural biology, and the central dogma explained the machinery of evolution,
and nearly neutral theory explained how selection competes with random fixation
dynamics that produce molecular clocks essential e.g. for dating evolutionary
histories. However, the quantitative proteomics revealed that fitness effects
not related to functional proficiency play much larger roles on long
evolutionary time scales than previously thought, with particular evidence that
some universal biophysical selection pressures act via protein expression
levels. This paper first summarizes recent progress in the 21st century towards
recovering this universal selection pressure. Then, the paper argues that
proteome cost minimization is the dominant, underlying "non-function" selection
pressure controlling most of the evolution of already functionally adapted
living systems. A theory of proteome cost minimization is described and argued
to have consequences for understanding evolutionary trade-offs, aging, cancer,
and neurodegenerative protein-misfolding diseases.


  Chlorophylls and carotenoids, key components of photosynthetic systems, are
proposed for molecular communications at the nanoscale with the mechanism of
resonance energy transfer. Both types of pigments are introduced focusing on
their exceptional properties like energy harvesting, ultra-low energy
consumption during transmissions, picoseconds delays, an ability for signal
conversions, and bio-compatibility. The theoretical considerations are
supported by two spectroscopic experiments on the photosystem II complex. The
first experiment aims at the description of the photosystem II including
calculation of the energy transfer efficiency between carotenoids and
chlorophylls. In the second one, the photochemical efficiency is estimated,
showing how effective the chlorophylls are in further energy processing. With
the experimentally determined values, communication and energetic performance
are analyzed, the probability of channel blockage is calculated and the energy
consumption per bit is estimated. Treating carotenoid molecules as
transmitters, chlorophylls as receivers, and the energy transfer between them
as a way to encode information, a throughput up to 1 Gbit/s is achievable with
a bit error rate below 10^-3, average transmission delays about 20 ps, and
energy consumption c.a. 2.0x10^-18 J/bit. These results indicate a high
potential of photosynthetic systems for nanocommunications and other related
applications, due to suitable energetic characteristics in terms of energy
harvesting abilities and low consumption for data transmissions.


  The COVID-19 epidemics, started in China in January 2020, was recognized to
have reached Italy around February 20; recent estimates show that most probably
the virus circulated in the country already in January, but was not recognized.
Data for the early dynamics of COVID-19 in Northern Italy are analyzed.


  In agricultural landscapes, the composition and spatial configuration of
cultivated and semi-natural elements strongly impact species dynamics, their
interactions and habitat connectivity. To allow for landscape structural
analysis and scenario generation, we here develop statistical tools for real
landscapes composed of geometric elements including 2D patches but also 1D
linear elements such as hedges. We design generative stochastic models that
combine a multiplex network representation and Gibbs energy terms to
characterize the distributional behavior of landscape descriptors for land-use
categories. We implement Metropolis-Hastings for this new class of models to
sample agricultural scenarios featuring parameter-controlled spatial and
temporal patterns (e.g., geometry, connectivity, crop-rotation).
Pseudolikelihood-based inference allows studying the relevance of model
components in real landscapes through statistical and functional validation,
the latter achieved by comparing commonly used landscape metrics between
observed and simulated landscapes. Models fitted to subregions of the Lower
Durance Valley (France) indicate strong deviation from random allocation, and
they realistically capture small-scale landscape patterns. In summary, our
approach of statistical modeling improves the understanding of structural and
functional aspects of agro-ecosystems, and it enables simulation-based
theoretical analysis of how landscape patterns shape biological and ecological
processes.


  Motion planning algorithms often leverage topological information about the
environment to improve planner performance. However, these methods often focus
only on the environment's connectivity while ignoring other properties such as
obstacle clearance, terrain conditions, and resource accessibility. We present
a method that augments a skeleton representing the workspace topology with such
information to guide a sampling-based motion planner to rapidly discover
regions most relevant to the problem at hand. Our approach decouples guidance
and planning, making it possible for basic planning algorithms to find desired
paths earlier in the planning process. We demonstrate the efficacy of our
approach in both robotics problems and applications in drug design. Our method
is able to produce desirable paths quickly with no change to the underlying
planner.


  We propose a time-parameterized analogy between the thermodynamic behavior of
a 3-level energy system and the progression of the HIV infection described by
the cell population evolution generated by an appropriated cellular automata
model. The development of internal energy and its fluctuations, and of the
entropy of the 3-level energy system allows the identification of an effective
temperature that uniquely characterizes the three main stages of the dynamic
process of HIV infection. Furthermore, this thermodynamical equivalence allows
obtaining the instants associated with particular time intervals of the
evolution of internal energy and entropy, which are not quantitatively
accessible by the usual dynamic models based on differential equations or
cellular automata. The maximum entropy point, which marks the threshold between
the states of positive and negative temperatures, also corresponds to the onset
of the immune system\' exhaustion and the concomitant and inexorable
progression to AIDS. Such point lies in the time interval where the energy
inversion mechanisms in the populations of non-infected and infected cells
occur. This time lag is also characterized by considerable fluctuations of the
internal energy when different samples are compared.