  Vesicle Traffic Systems (VTSs) are the material transport mechanisms among
the compartments inside the biological cells. The compartments are viewed as
nodes that are labeled with the containing chemicals and the transport channels
are similarly viewed as labeled edges between the nodes. Understanding VTSs is
an ongoing area of research and for many cells they are partially known. For
example, there may be undiscovered edges, nodes, or their labels in a VTS of a
cell. It has been speculated that there are properties that the VTSs must
satisfy. For example, stability, i.e., every chemical that is leaving a
compartment comes back. Many synthesis questions may arise in this scenario,
where we want to complete a partially known VTS under a given property. In the
paper, we present novel encodings of the above questions into the QBF
(quantified Boolean formula) satisfiability problems. We have implemented the
encodings in a highly configurable tool and applied to a couple of
found-in-nature VTSs and several synthetic graphs. Our results demonstrate that
our method can scale up to the graphs of interest.


  Cells need to reliably sense external ligand concentrations to achieve
various biological functions such as chemotaxis or signaling. The molecular
recognition of ligands by surface receptors is degenerate in many systems
leading to crosstalk between different receptors. Crosstalk is often thought of
as a deviation from optimal specific recognition, as the binding of non-cognate
ligands can interfere with the detection of the receptor's cognate ligand,
possibly leading to a false triggering of a downstream signaling pathway. Here
we quantify the optimal precision of sensing the concentrations of multiple
ligands by a collection of promiscuous receptors. We demonstrate that crosstalk
can improve precision in concentration sensing and discrimination tasks. To
achieve superior precision, the additional information about ligand
concentrations contained in short binding events of the non-cognate ligand
should be exploited. We present a proofreading scheme to realize an approximate
estimation of multiple ligand concentrations that reaches a precision close to
the derived optimal bounds. Our results help rationalize the observed ubiquity
of receptor crosstalk in molecular sensing.


  Microbial communities routinely have several alternative stable states
observed for the same environmental parameters. Sudden and irreversible
transitions between these states make external manipulation of these systems
more complicated. To better understand the mechanisms and origins of
multistability in microbial communities, we introduce and study a model of a
microbial ecosystem colonized by multiple specialist species selected from a
fixed pool. Growth of each species can be limited by essential nutrients of two
types, e.g. carbon and nitrogen, each represented in the environment by
multiple metabolites. We demonstrate that our model has an exponentially large
number of potential stable states realized for different environmental
parameters. Using game theoretical methods adapted from the stable marriage
problem we predict all of these states based only on ranked lists of
competitive abilities of species for each of the nutrients. We show that for
every set of nutrient influxes, several mutually uninvadable stable states are
generally feasible and we distinguish them based upon their dynamic stability.
We further explore an intricate network of discontinuous transitions (regime
shifts) between these alternative states both in the course of community
assembly, or upon changes of nutrient influxes.


  Drosophila melanogaster, a small dipteran of African origin, represents one
of the best-studied model organisms. Early work in this system has uniquely
shed light on the basic principles of genetics and resulted in a versatile
collection of genetic tools that allow to uncover mechanistic links between
genotype and phenotype. Moreover, given its world-wide distribution in diverse
habitats and its moderate genome-size, Drosophila has proven very powerful for
population genetics inference and was one of the first eukaryotes whose genome
was fully sequenced. In this book-chapter, we provide a brief historical
overview of research in Drosophila and then focus on recent advances during the
genomic era. After describing different types and sources of genomic data, we
discuss mechanisms of neutral evolution including the demographic history of
Drosophila and the effects of recombination and biased gene conversion. Then,
we review recent advances in detecting genome-wide signals of selection, such
as soft and hard selective sweeps. We further provide a brief introduction to
background selection, selection of non-coding DNA and codon usage and focus on
the role of structural variants, such as transposable elements and chromosomal
inversions, during the adaptive process. Finally, we discuss how genomic data
helps to dissect neutral and adaptive evolutionary mechanisms that shape
genetic and phenotypic variation in natural populations along environmental
gradients. In summary, this book chapter serves as a starting point to
Drosophila population genomics and provides an introduction to the system and
an overview to data sources, important population genetic concepts and recent
advances in the field.


  In this work, we introduce a new type of linear classifier that is
implemented in a chemical form. We propose a novel encoding technique which
simultaneously represents multiple datasets in an array of microliter-scale
chemical mixtures. Parallel computations on these datasets are performed as
robotic liquid handling sequences, whose outputs are analyzed by
high-performance liquid chromatography. As a proof of concept, we chemically
encode several MNIST images of handwritten digits and demonstrate successful
chemical-domain classification of the digits using volumetric perceptrons. We
additionally quantify the performance of our method with a larger dataset of
binary vectors and compare the experimental measurements against predicted
results. Paired with appropriate chemical analysis tools, our approach can work
on increasingly parallel datasets. We anticipate that related approaches will
be scalable to multilayer neural networks and other more complex algorithms.
Much like recent demonstrations of archival data storage in DNA, this work
blurs the line between chemical and electrical information systems, and offers
early insight into the computational efficiency and massive parallelism which
may come with computing in chemical domains.


  Animal groups exhibit emergent properties that are a consequence of local
interactions. Linking individual-level behaviour to coarse-grained descriptions
of animal groups has been a question of fundamental interest. Here, we present
two complementary approaches to deriving coarse-grained descriptions of
collective behaviour at so-called mesoscopic scales, which account for the
stochasticity arising from the finite sizes of animal groups. We construct
stochastic differential equations (SDEs) for a coarse-grained variable that
describes the order/consensus within a group. The first method of construction
is based on van Kampen's system-size expansion of transition rates. The second
method employs Gillespie's chemical Langevin equations. We apply these two
methods to two microscopic models from the literature, in which organisms
stochastically interact and choose between two directions/choices of foraging.
These `binary-choice' models differ only in the types of interactions between
individuals, with one assuming simple pair-wise interactions, and the other
incorporating higher-order effects. In both cases, the derived mesoscopic SDEs
have multiplicative, or state-dependent, noise. However, the different models
demonstrate the contrasting effects of noise: increasing order in the pair-wise
interaction model, whilst reducing order in the higher-order interaction model.
Although both methods yield identical SDEs for such binary-choice, or
one-dimensional, systems, the relative tractability of the chemical Langevin
approach is beneficial in generalizations to higher-dimensions. In summary,
this book chapter provides a pedagogical review of two complementary methods to
construct mesoscopic descriptions from microscopic rules and demonstrates how
resultant multiplicative noise can have counter-intuitive effects on shaping
collective behaviour.


  The Bayesian method is noted to produce spuriously high posterior
probabilities for phylogenetic trees in analysis of large datasets, but the
precise reasons for this over-confidence are unknown. In general, the
performance of Bayesian selection of misspecified models is poorly understood,
even though this is of great scientific interest since models are never true in
real data analysis. Here we characterize the asymptotic behavior of Bayesian
model selection and show that when the competing models are equally wrong,
Bayesian model selection exhibits surprising and polarized behaviors in large
datasets, supporting one model with full force while rejecting the others. If
one model is slightly less wrong than the other, the less wrong model will
eventually win when the amount of data increases, but the method may become
overconfident before it becomes reliable. We suggest that this extreme behavior
may be a major factor for the spuriously high posterior probabilities for
evolutionary trees. The philosophical implications of our results to the
application of Bayesian model selection to evaluate opposing scientific
hypotheses are yet to be explored, as are the behaviors of non-Bayesian methods
in similar situations.


  The analysis of biological networks is characterized by the definition of
precise linear constraints used to cumulatively reduce the solution space of
the computed states of a multi-omic (for instance metabolic, transcriptomic and
proteomic) model. In this paper, we attempt, for the first time, to combine
metabolic modelling and networked Cox regression, using the metabolic model of
the bacterium Helicobacter Pylori. This enables a platform both for
quantitative analysis of networked regression, but also testing the findings
from network regression (a list of significant vectors and their networked
relationships) on in vivo transcriptomic data. Data generated from the model,
using flux balance analysis to construct a Pareto front, specifically, a
trade-off of Oxygen exchange and growth rate and a trade-off of Carbon Dioxide
exchange and growth rate, is analysed and then the model is used to quantify
the success of the analysis. It was found that using the analysis,
reconstruction of the initial data was considerably more successful than a pure
noise alternative. Our methodological approach is quite general and it could be
of interest for the wider community of complex networks researchers; it is
implemented in a software tool, MoNeRe, which is freely available through the
Github platform.


  The Dirichlet Process (DP) mixture model has become a popular choice for
model-based clustering, largely because it allows the number of clusters to be
inferred. The sequential updating and greedy search (SUGS) algorithm (Wang and
Dunson, 2011) was proposed as a fast method for performing approximate Bayesian
inference in DP mixture models, by posing clustering as a Bayesian model
selection (BMS) problem and avoiding the use of computationally costly Markov
chain Monte Carlo methods. Here we consider how this approach may be extended
to permit variable selection for clustering, and also demonstrate the benefits
of Bayesian model averaging (BMA) in place of BMS. Through an array of
simulation examples and well-studied examples from cancer transcriptomics, we
show that our method performs competitively with the current state-of-the-art,
while also offering computational benefits. We apply our approach to
reverse-phase protein array (RPPA) data from The Cancer Genome Atlas (TCGA) in
order to perform a pan-cancer proteomic characterisation of 5,157 tumour
samples. We have implemented our approach, together with the original SUGS
algorithm, in an open-source R package named sugsvarsel, which accelerates
analysis by performing intensive computations in C++ and provides automated
parallel processing. The R package is freely available from:
https://github.com/ococrook/sugsvarsel


  A fundamental model of tumor growth in the presence of cytotoxic
chemotherapeutic agents is formulated. The model allows to study the role of
the Norton-Simon hypothesis in the context of dose-dense chemotherapy.
Dose-dense protocols aim at reducing the period between courses of chemotherapy
from three weeks to two weeks, in order to avoid tumor regrowth between cycles.
We address the conditions under which these protocols might be more or less
beneficial in comparison to less dense settings, depending on the sensitivity
of the tumor cells to the cytotoxic drugs. The effects of varying other
parameters of the protocol, as for example the duration of each continuous drug
infusion, are also inspected. We believe that the present model might serve as
a foundation for the development of more sophisticated models for cancer
chemotherapy.


  Methods of statistical physics have proven valuable for studying the
evolution of cooperation in social dilemma games. However, recent empirical
research shows that cooperative behavior in social dilemmas is only one kind of
a more general class of behavior, namely moral behavior, which includes
reciprocity, respecting others' property, honesty, equity, efficiency, as well
as many others. Inspired by these experimental works, we here open up the path
towards studying other forms of moral behavior with methods of statistical
physics. We argue that this is a far-reaching direction for future research
that can help us answer fundamental questions about human sociality. Why did
our societies evolve as they did? What moral principles are more likely to
emerge? What happens when different moral principles clash? Can we predict the
break out of moral conflicts in advance and contribute to their solution? These
are amongst the most important questions of our time, and methods of
statistical physics could lead to new insights and contribute towards finding
answers.


  Biorefinery sector has become a serious dispute for cleaner and sustainable
development in recent years. In the present study, pretreatment of pineapple
peel waste was carried out in high pressure reactor using various
pretreatment-enhancers. The type and concentration effect of each enhancer on
hemicellulose solubilization was systematically investigated. The binary acid
(phenol + sulfuric acid) at 180 {\deg}C was found to be superior amongst other
studied enhancers, giving 81.17% (w/v) hemicellulose solubilization in
liquid-fraction under optimized conditions. Solid residue thus obtained was
subjected to enzymatic hydrolysis that resulted into 24.50% (w/v) cellulose
breakdown. Treated solid residue was further characterized by scanning electron
microscopy and Fourier transform infrared spectroscopy to elucidate structural
changes. The pooled fractions (acid treated and enzymatically hydrolyzed) were
fermented using Clostridium acetobutylicum NRRL B 527 which resulted in butanol
production of 5.18 g/L with yield of 0.13 g butanol/g sugar consumed.
Therefore, pretreatment of pineapple peel waste evaluated in this study can be
considered as milestone in utilization of low cost feedstock, for bioenergy
production.


  Many dynamical systems arising in biology and other areas exhibit
multistationarity (two or more positive steady states with the same conserved
quantities). Although deciding multistationarity for a polynomial dynamical
system is an effective question in real algebraic geometry, it is in general
difficult to determine whether a given network can give rise to a
multistationary system, and if so, to identify witnesses to multistationarity,
that is, specific parameter values for which the system exhibits multiple
steady states. Here we investigate both problems. First, we build on work of
Conradi, Feliu, Mincheva, and Wiuf, who showed that for certain reaction
networks whose steady states admit a positive parametrization,
multistationarity is characterized by whether a certain "critical function"
changes sign. Here, we allow for more general parametrizations, which make it
much easier to determine the existence of a sign change. This is particularly
simple when the steady-state equations are linearly equivalent to binomials; we
give necessary conditions for this to happen, which hold for many networks
studied in the literature. We also give a sufficient condition for
multistationarity of networks whose steady-state equations can be replaced by
equivalent triangular-form equations. Finally, we present methods for finding
witnesses to multistationarity, which we show work well for certain structured
reaction networks, including those common to biological signaling pathways. Our
work relies on results from degree theory, on the existence of explicit
rational parametrizations of the steady states, and on the specialization of
Groebner bases.


  In systems and synthetic biology, much research has focused on the behavior
and design of single pathways, while, more recently, experimental efforts have
focused on how cross-talk (coupling two or more pathways) or inhibiting
molecular function (isolating one part of the pathway) affects systems-level
behavior. However, the theory for tackling these larger systems in general has
lagged behind. Here, we analyze how joining networks (e.g., cross-talk) or
decomposing networks (e.g., inhibition or knock-outs) affects three properties
that reaction networks may possess---identifiability (recoverability of
parameter values from data), steady-state invariants (relationships among
species concentrations at steady state, used in model selection), and
multistationarity (capacity for multiple steady states, which correspond to
multiple cell decisions). Specifically, we prove results that clarify, for a
network obtained by joining two smaller networks, how properties of the smaller
networks can be inferred from or can imply similar properties of the original
network. Our proofs use techniques from computational algebraic geometry,
including elimination theory and differential algebra.


  We analyze the optimal harvesting problem for an ecosystem of species that
experience environmental stochasticity. Our work generalizes the current
literature significantly by taking into account non-linear interactions between
species, state-dependent prices, and species injections. The key generalization
is making it possible to not only harvest, but also `seed' individuals into the
ecosystem. This is motivated by how fisheries and certain endangered species
are controlled. The harvesting problem becomes finding the optimal
harvesting-seeding strategy that maximizes the expected total income from the
harvest minus the lost income from the species injections. Our analysis shows
that new phenomena emerge due to the possibility of species injections.
  It is well-known that multidimensional harvesting problems are very hard to
tackle. We are able to make progress, by characterizing the value function as a
viscosity solution of the associated Hamilton-Jacobi-Bellman (HJB) equations.
Moreover, we provide a verification theorem, which tells us that if a function
has certain properties, then it will be the value function. This allows us to
show heuristically, as was shown in Lungu and $\O$ksendal (Bernoulli '01), that
it is almost surely never optimal to harvest or seed from more than one
population at a time.
  We approximate the continuous-time systems by Markov chains and show that the
optimal harvesting-seeding strategies of the Markov chain approximations
converge to the correct optimal harvesting strategy. This is used to provide
numerical approximations to the optimal harvesting-seeding strategies and is a
first step towards a full understanding of the intricacies of how one should
harvest and seed interacting species. In particular, we look at three examples:
one species modeled by a Verhulst-Pearl diffusion, two competing species and a
two-species predator-prey system.


  We propose a purely probabilistic model to explain the evolution path of a
population maximum fitness. We show that after $n$ births in the population
there are about $\ln n$ upwards jumps. This is true for any mutation
probability and any fitness distribution and therefore suggests a general law
for the number of upwards jumps. Simulations of our model show that a typical
evolution path has first a steep rise followed by long plateaux. Moreover,
independent runs show parallel paths. This is consistent with what was observed
by Lenski and Travisano (1994) in their bacteria experiments.


  We propose and analyse a mathematical model for cholera considering
vaccination. We show that the model is epidemiologically and mathematically
well posed and prove the existence and uniqueness of disease-free and endemic
equilibrium points. The basic reproduction number is determined and the local
asymptotic stability of equilibria is studied. The biggest cholera outbreak of
world's history began on 27th April 2017, in Yemen. Between 27th April 2017 and
15th April 2018 there were 2275 deaths due to this epidemic. A vaccination
campaign began on 6th May 2018 and ended on 15th May 2018. We show that our
model is able to describe well this outbreak. Moreover, we prove that the
number of infected individuals would have been much lower provided the
vaccination campaign had begun earlier.


  Realistic human contact networks capable of spreading infectious disease, for
example studied in social contact surveys, exhibit both significant degree
heterogeneity and clustering, both of which greatly affect epidemic dynamics.
To understand the joint effects of these two network properties on epidemic
dynamics, the effective degree model of Lindquist et al. is reformulated with a
new moment closure to apply to highly clustered networks. A simulation study
comparing alternative ODE models and stochastic simulations is performed for
SIR (Susceptible-Infected-Removed) epidemic dynamics, including a test for the
conjectured error behaviour in Pellis et al., providing evidence that this
novel model can be a more accurate approximation to epidemic dynamics on
complex networks than existing approaches.


  In this study, we investigated the effects of specific low frequency
electromagnetic fields sequences on U937 cells, an in vitro model of human
monocyte/macrophage differentiation. U937 cells were exposed to electromagnetic
stimulation by means of the SyntheXer system using two similar sequences,
XR-BC31 and XR-BC31/F. Each sequence was a time series of twenty-nine wave
segments. Here, we report that exposure (4 days, once a day) of U937 cells to
the XR-BC31 setting, but not to the XR-BC31/F, resulted in increased expression
of the histone demethylase KDM6B along with a global reduction in histone H3
lysine 27 (H3K27) tri-methylation. Furthermore, exposure to the XR-BC31
sequence induced differentiation of U937 cells towards a macrophage-like
phenotype displaying a KDM6B dependent increase in expression and secretion of
the anti-inflammatory interleukins (ILs), IL-10 and IL-4. Importantly, all the
observed changes were highly dependent on the sequence's nature. Our results
open a new way of interpretation for the effects of low frequency
electromagnetic fields observed in vivo. Indeed, it is conceivable that a
specific low frequency electromagnetic fields treatment may cause changes in
chromatin accessibility and consequently in the expression of anti-inflammatory
mediators and in cell differentiation.


  In this article we consider a microscopic model for host-vector disease
transmission based on configuration space analysis. Using Vlasov scaling we
obtain the corresponding mesoscopic (kinetic) equations, describing the density
of susceptible and infected compartments in space. The resulting system of
equations can be seen as a generalization to a spatial SISUV model.


  Here we present Translocatome, the first dedicated database of human
translocating proteins. The core of the Translocatome database is the manually
curated data set of 213 human translocating proteins listing the source of
their experimental validation, several details of their translocation
mechanism, their local compartmentalized interactome, as well as their
involvement in signalling pathways and disease development. In addition, using
the well-established and widely used gradient boosting machine learning tool,
XGBoost, Translocatome provides translocation probability values for 13,066
human proteins identifying 1133 and 3268 high- and low-confidence translocating
proteins, respectively. The database has user-friendly search options with a
UniProt autocomplete quick search and advanced search for proteins filtered by
their localization, UniProt identifiers, translocation likelihood or data
complexity. Download options of search results, manually curated and predicted
translocating protein sets are available on its website. The update of the
database is helped by its manual curation framework and connection to the
previously published ComPPI compartmentalized protein-protein interaction
database. As shown by the application examples of merlin (NF2) and tumor
protein 63 (TP63) Translocatome allows a better comprehension of protein
translocation as a systems biology phenomenon and can be used as a
discovery-tool in the protein translocation field. The database is available
here: http://translocatome.linkgroup.hu


  A mechanism is proposed for the growth of vesicles dispersed in a liquid
solvent and a size distribution function is obtained for the vesicles, both
from the first principles calculations. This distribution function is shown to
be positively skewed and evolving in time obeying a Fokker-Planck type
equation. The critical size of the spherical vesicles is shown to grow in time
with an exponent of 0.25. A constant is suggested that characterizes how easily
a vesicle can absorb amphiphiles into its periphery.


  We present three examples of chemical reaction networks whose ordinary
differential equation scaling limit are almost identical and in all cases
stable. Nevertheless, the Markov jump processes associated to these reaction
networks display the full range of behaviors: one is stable (positive
recurrent), one is unstable (transient) and one is marginally stable (null
recurrent). We study these differences and characterize the invariant measures
by Lyapunov function techniques. In particular, we design a natural set of such
functions which scale homogeneously to infinity, taking advantage of the same
scaling behavior of the reaction rates.


  Is it possible to break the host-vector chain of transmission when there is
an influx of infectious hosts into a na\"{i}ve population and competent vector?
To address this question, a class of vector-borne disease models with an
arbitrary number of infectious stages that account for immigration of infective
individuals is formulated. The proposed model accounts for forward and backward
progression, capturing the mitigation and aggravation to and from any stages of
the infection, respectively. The model has a rich dynamic, which depends on the
patterns of infected immigrant influx into the host population and connectivity
of the transfer between infectious classes. We provide conditions under which
the answer of the initial question is positive.


  Local perturbations of individual metabolic reactions may result in different
levels of lethality, depending on their roles in metabolism and the size of
subsequent cascades induced by their failure. Moreover, essentiality of
individual metabolic reactions may show large variations within and across
species. Here we quantify their essentialities in hundreds of species by
computing the growth rate after removal of individual and pairs of reactions by
flux balance analysis. We find that about 10% of reactions are essential, i.e.,
growth stops without them, and most of the remaining reactions are redundant in
the metabolic network of each species. This large-scale and cross-species study
allows us to determine ad hoc ages of each reaction and species. We find that
when a reaction is older and contained in younger species, the reaction is more
likely to be essential. Such correlations of essentiality with the ages of
reactions and species may be attributable to the evolution of cellular
metabolism, in which alternative pathways are recruited to ensure the stability
of important reactions to various degrees across species.


  Several important biological processes are initiated by the binding of a
protein to a specific site on the DNA. The strategy adopted by a protein,
called transcription factor (TF), for searching its specific binding site on
the DNA has been investigated over several decades. In recent times the effects
obstacles, like DNA-binding proteins, on the search by TF has begun to receive
attention. RNA polymerase (RNAP) motors collectively move along a segment of
the DNA during a genomic process called transcription. This RNAP traffic is
bound to affect the diffusive scanning of the same segment of the DNA by a TF
searching for its binding site. Motivated by this phenomenon, here we develop a
kinetic model where a `particle', that represents a TF, searches for a specific
site on a one-dimensional lattice. On the same lattice another species of
particles, each representing a RNAP, hop from left to right exactly as in a
totally asymmetric simple exclusion process (TASEP) which forbids simultaneous
occupation of any site by more than one particle, irrespective of their
identities. Although the TF is allowed to attach to or detach from any lattice
site, the RNAPs can attach only to the first site at the left edge and detach
from only the last site on the right edge of the lattice. We formulate the
search as a {\it first-passage} process; the time taken to reach the target
site {\it for the first time}, starting from a well defined initial state, is
the search time. By approximate analytical calculations and Monte Carlo (MC)
computer simulations, we calculate the mean search time. We show that RNAP
traffic rectifies the diffusive motion of TF to that of a Brownian ratchet, and
the mean time of successful search can be even shorter than that required in
the absence of RNAP traffic. Moreover, we show that there is an optimal rate of
detachment that corresponds to the shortest mean search time.


  There is a long tradition of the axiomatic study of consensus methods in
phylogenetics that satisfy certain desirable properties. One
recently-introduced property is associative stability, which is desirable
because it confers a computational advantage, in that the consensus method only
needs to be computed "pairwise". In this paper, we introduce a phylogenetic
consensus method that satisfies this property, in addition to being "regular".
The method is based on the introduction of a partial order on the set of rooted
phylogenetic trees, itself based on the notion of a hierarchy-preserving map
between trees. This partial order may be of independent interest. We call the
method "lattice consensus", because it takes the unique maximal element in a
lattice of trees defined by the partial order. Aside from being associatively
stable, lattice consensus also satisfies the property of being Pareto on rooted
triples, answering in the affirmative a question of Bryant et al (2017). We
conclude the paper with an answer to another question of Bryant et al, showing
that there is no regular extension stable consensus method for binary trees.


  Recently, so-called treebased phylogenetic networks have gained considerable
interest in the literature, where a treebased network is a network that can be
constructed from a phylogenetic tree, called the base tree, by adding
additional edges. The main aim of this manuscript is to provide some sufficient
criteria for treebasedness by reducing phylogenetic networks to related graph
structures. While it is generally known that deciding whether a network is
treebased is NP-complete, one of these criteria, namely edgebasedness, can be
verified in linear time. Surprisingly, the class of edgebased networks is
closely related to a well-known family of graphs, namely the class of
generalized series parallel graphs, and we will explore this relationship in
full detail. Additionally, we introduce further classes of treebased networks
and analyze their relationships.


  Phylogenetic networks are a generalization of phylogenetic trees allowing for
the representation of non-treelike evolutionary events such as hybridization.
Typically, such networks have been analyzed based on their `level', i.e. based
on the complexity of their 2-edge-connected components. However, recently the
question of how `treelike' a phylogenetic network is has become the center of
attention in various studies. This led to the introduction of \emph{tree-based
networks}, i.e. networks that can be constructed from a phylogenetic tree,
called the \emph{base tree}, by adding additional edges. While the concept of
tree-basedness was originally introduced for rooted phylogenetic networks, it
has recently also been considered for unrooted networks. In the present study,
we compare and contrast findings obtained for unrooted \emph{binary} tree-based
networks to unrooted \emph{non-binary} networks. In particular, while it is
known that up to level 4 all unrooted binary networks are tree-based, we show
that in the case of non-binary networks, this result only holds up to level 3.


  A human respiratory syncytial virus surveillance system was implemented in
Florida in 1999, to support clinical decision-making for prophylaxis of
premature newborns. Recently, a local periodic SEIRS mathematical model was
proposed in [Stat. Optim. Inf. Comput. 6 (2018), no.1, 139--149] to describe
real data collected by Florida's system. In contrast, here we propose a
non-local fractional (non-integer) order model. A fractional optimal control
problem is then formulated and solved, having treatment as the control.
Finally, a cost-effectiveness analysis is carried out to evaluate the cost and
the effectiveness of proposed control measures during the intervention period,
showing the superiority of obtained results with respect to previous ones.


  Extreme deformations of the DNA double helix attracted a lot of attention
during the past decades. Particularly, the determination of the persistence
length of DNA with extreme local disruptions, or kinks, has become a crucial
problem in the studies of many important biological processes. In this paper we
review an approach to calculate the persistence length of the double helix by
taking into account the formation of kinks of arbitrary configuration. The
reviewed approach improves the Kratky--Porod model to determine the type and
nature of kinks that occur in the double helix, by measuring a reduction of the
persistence length of the kinkable DNA.


  Members of a social species need to make appropriate decisions about who,
how, and when to interact with others in their group. However, it has been
difficult for researchers to detect the inputs to these decisions and, in
particular, how much information individuals actually have about their social
context. We present a new method that can serve as a social assay to quantify
how patterns of aggression depend upon information about the ranks of
individuals within social dominance hierarchies. Applied to existing data on
aggression in 172 social groups across 85 species in 23 orders, it reveals
three main patterns of rank-dependent social dominance: the downward heuristic
(aggress uniformly against lower-ranked opponents), close competitors (aggress
against opponents ranked slightly below self), and bullying (aggress against
opponents ranked much lower than self). The majority of the groups (133 groups,
77%) follow a downward heuristic, but a significant minority (38 groups, 22%)
show more complex social dominance patterns (close competitors or bullying)
consistent with higher levels of social information use. These patterns are not
phylogenetically constrained and different groups within the same species can
use different patterns, suggesting that heuristics use may depend on context
and the structuring of aggression by social information should not be
considered a fixed characteristic of a species. Our approach provides new
opportunities to study the use of social information within and across species
and the evolution of social complexity and cognition.


  Stochasticity plays important roles in reaction systems. Vector fields of
probability flux and velocity characterize time-varying and steady-state
properties of these systems, including high probability paths, barriers,
checkpoints among different stable regions, as well as mechanisms of dynamic
switching among them. However, conventional fluxes on continuous space are
ill-defined and are problematic when at boundaries of the state space or when
copy numbers are small. By re-defining the derivative and divergence operators
based on the discrete nature of reactions, we introduce new formulations of
discrete fluxes. Our flux model fully accounts for the discreetness of both the
state space and the jump processes of reactions. The reactional discrete flux
satisfies the continuity equation and describes the behavior of the system
evolving along directions of reactions. The species discrete flux directly
describes the dynamic behavior in the state space of the reactants such as the
transfer of probability mass. With the relationship between these two fluxes
specified, we show how to construct time-evolving and steady-state global
flow-maps of probability flux and velocity in the directions of every species
at every microstate, and how they are related to the outflow and inflow of
probability fluxes when tracing out reaction trajectories. We also describe how
to impose proper conditions enabling exact quantification of flux and velocity
in the boundary regions, without the difficulty of enforcing artificial
reflecting conditions. We illustrate the computation of probability flux and
velocity using three model systems, namely, the birth-death process, the
bistable Schl\"ogl model, and the oscillating Schnakenberg model.


  Internet, social media, neuronal or blood vessel are organized in complex
networks. These networks are characterized by several quantities such as the
underlying graph connectivity (topology), how they grow in time, scaling laws
or by the mean time a random search can cover a network and also by shortest
path between two nodes. We present here a novel type of network property based
on a unidirectional transport mechanism occurring in the Endoplasmic Reticulum
network found in the cell cytoplasm. This mechanism is an active-waiting
transportation, where molecules have to wait a random time before being
transported from one node to the next one. We find that the consequence of this
unusual network transportation is that molecules travel together by recurrent
packets, which quite a large deviation behavior compared to classical
propagation in graphs. To conclude, this form of transportation is associated
with an efficient and robust molecular redistribution inside cells.


  The current gold standard for unambiguous identification in metabolomics
analysis is based on comparing two or more orthogonal properties from the
analysis of authentic, pure reference materials (standards) to experimental
data acquired in the same laboratory with the same analytical methods. This
represents a significant limitation for comprehensive chemical identification
of small molecules in complex samples since this process is time-consuming and
costly, and the majority of molecules are not yet represented by standards,
leading to a need for standards-free identification. To address this need, we
are advancing chemical property calculations and developing multi-attribute
scoring and matching algorithms to utilize data from multiple analytical
platforms through the utilization and creation of the in silico Chemical
Library Engine (ISiCLE) and the Multi-Attribute Matching Engine (MAME). Here,
we describe our results in a blinded analysis of synthetic chemical mixtures as
part of the U.S. Environmental Protection Agency's (EPA) Non-Targeted Analysis
Collaborative Trial (ENTACT). The blinded false negative rate (FNR), false
discovery rate (FDR), and accuracy were 57%, 77%, and 91%, respectively. For
high confidence identifications, the FDR was 35%. After unblinding of the
sample compositions, we improved our approach by optimizing the scoring
parameters used to increase confidence. The final FNR, FDR, and accuracy were
67%, 53%, and 96%, respectively. For high confidence identifications, the FDR
was 10%. This study demonstrates that standards-free small molecule
identification and multi-attribute matching methods can significantly reduce
reliance on standards.


  Evolutionary game dynamics in structured populations has been extensively
explored in past decades. However, most previous studies assume that payoffs of
individuals are fully determined by the strategic behaviors of interacting
parties and social ties between them only serve as the indicator of the
existence of interactions. This assumption neglects important information
carried by inter-personal social ties such as genetic similarity, geographic
proximity, and social closeness, which may crucially affect the outcome of
interactions. To model these situations, we present a framework of evolutionary
multiplayer games on graphs with edge diversity, where different types of edges
describe diverse social ties. Strategic behaviors together with social ties
determine the resulting payoffs of interactants. Under weak selection, we
provide a general formula to predict the success of one behavior over the
other. We apply this formula to various examples which cannot be dealt with
using previous models, including the division of labor and relationship- or
edge-dependent games. We find that labor division facilitates collective
cooperation by decomposing a many-player game into several games of smaller
sizes. The evolutionary process based on relationship-dependent games can be
approximated by interactions under a transformed and unified game. Our work
stresses the importance of social ties and provides effective methods to reduce
the calculating complexity in analyzing the evolution of realistic systems.


  Many experimental and field studies have shown that adaptation can occur very
rapidly. Two qualitatively different modes of fast adaptation have been
proposed: selective sweeps wherein large shifts in the allele frequencies occur
at a few loci and evolution via small changes in the allele frequencies at many
loci. While the first process has been thoroughly investigated within the
framework of population genetics, the latter is based on quantitative genetics
and is much less understood. Here we summarize results from our recent
theoretical studies of a quantitative genetic model of polygenic adaptation
that makes explicit reference to population genetics to bridge the gap between
the two frameworks. Our key results are that polygenic adaptation may be a
rapid process and can proceed via subtle or dramatic changes in the allele
frequency depending on the sizes of the phenotypic effects relative to a
threshold value. We also discuss how the signals of polygenic selection may be
detected in the genome. While powerful methods are available to identify
signatures of selective sweeps at loci controling quantitative traits, the
development of statistical tests for detecting small shifts of allele
frequencies at quantitative trait loci is still in its infancy.


  Linked beneficial and deleterious mutations are known to decrease the
fixation probability of a favorable mutation in large asexual populations.
While the hindering effect of strongly deleterious mutations on adaptive
evolution has been well studied, how weak deleterious mutations, either in
isolation or with superior beneficial mutations, influence the fixation of a
beneficial mutation has not been fully explored. Here, using a multitype
branching process, we obtain an accurate analytical expression for the fixation
probability when deleterious effects are weak, and exploit this result along
with the clonal interference theory to investigate the joint effect of linked
beneficial and deleterious mutations on the rate of adaptation. We find that
when the mutation rate is increased beyond the beneficial fitness effect, the
fixation probability of the beneficial mutant decreases from Haldane's
classical result towards zero. This has the consequence that above a critical
mutation rate that may depend on the population size, the adaptation rate
decreases exponentially with the mutation rate and is independent of the
population size. In addition, we find that for a range of mutation rates, both
beneficial and deleterious mutations interfere and impede the adaptation
process in large populations. We also study the evolution of mutation rates in
adapting asexual populations, and conclude that linked beneficial mutations
have a stronger influence on mutator fixation than the deleterious mutations.


  We apply tools from real algebraic geometry to the problem of
multistationarity of chemical reaction networks. A particular focus is on the
case of reaction networks whose steady states admit a monomial parametrization.
For such systems we show that in the space of total concentrations
multistationarity is scale invariant: if there is multistationarity for some
value of the total concentrations, then there is multistationarity on the
entire ray containing this value (possibly for different rate constants) -- and
vice versa. Moreover, for these networks it is possible to decide about
multistationarity independent of the rate constants by formulating
semi-algebraic conditions that involve only concentration variables. These
conditions can easily be extended to include total concentrations. Hence
quantifier elimination may give new insights into multistationarity regions in
the space of total concentrations. To demonstrate this, we show that for the
distributive phosphorylation of a protein at two binding sites
multistationarity is only possible if the total concentration of the substrate
is larger than either the total concentration of the kinase or the total
concentration of the phosphatase. This result is enabled by the chamber
decomposition of the space of total concentrations from polyhedral geometry.
Together with the corresponding sufficiency result of Bihan et al. this yields
a characterization of multistationarity up to lower dimensional regions.


  The rise and spread of antibiotic resistance causes worsening medical cost
and mortality especially for life-threatening bacteria infections, thereby
posing a major threat to global health. Prescribing behavior of physicians is
one of the important factors impacting the underlying dynamics of resistance
evolution. It remains unclear when individual prescribing decisions can lead to
the overuse of antibiotics on the population level, and whether population
optimum of antibiotic use can be reached through an adaptive social learning
process that governs the evolution of prescribing norm. Here we study a
behavior-disease interaction model, specifically incorporating a feedback loop
between prescription behavior and resistance evolution. We identify the
conditions under which antibiotic resistance can evolve as a result of the
tragedy of the commons in antibiotic overuse. Furthermore, we show that fast
social learning that adjusts prescribing behavior in prompt response to
resistance evolution can steer out cyclic oscillations of antibiotic usage
quickly towards the stable population optimum of prescribing. Our work
demonstrates that provision of prompt feedback to prescribing behavior with the
collective consequences of treatment decisions and costs that are associated
with resistance helps curb the overuse of antibiotics.


  Site-specific recombination is an enzymatic process where two sites of
precise sequence and orientation along a circle come together, are cleaved, and
the ends are recombined. Site-specific recombination on a knotted substrate
produces another knot or a two-component link depending on the relative
orientation of the sites prior to recombination. Mathematically, site-specific
recombination is modeled as coherent (knot to link) or non-coherent (knot to
knot) banding. We here survey recent developments in the study of non-coherent
bandings on knots and discuss biological implications.


  Recovery of multispecies oral biofilms is investigated following treatment by
chlorhexidine gluconate (CHX), iodine-potassium iodide (IPI) and Sodium
hypochlorite (NaOCl) both experimentally and theoretically. Experimentally,
biofilms taken from two donors were exposed to the three antibacterial
solutions (irrigants) for 10 minutes, respectively. We observe that (a) live
bacterial cell ratios decline for a week after the exposure and the trend
reverses beyond a week; after fifteen weeks, live bacterial cell ratios in
biofilms fully return to their pretreatment levels; (b) NaOCl is shown as the
strongest antibacterial agent for the oral biofilms; (c) multispecies oral
biofilms from different donors showed no difference in their susceptibility to
all the bacterial solutions. Guided by the experiment, a mathematical model for
biofilm dynamics is developed, accounting for multiple bacterial phenotypes,
quorum sensing, and growth factor proteins, to describe the nonlinear time
evolutionary behavior of the biofilms. The model captures time evolutionary
dynamics of biofilms before and after antibacterial treatment very well. It
reveals the crucial role played by quorum sensing molecules and growth factors
in biofilm recovery and verifies that the source of biofilms has a minimal to
their recovery. The model is also applied to describe the state of biofilms of
various ages treated by CHX, IPI and NaOCl, taken from different donors. Good
agreement with experimental data predicted by the model is obtained as well,
confirming its applicability to modeling biofilm dynamics in general.


  Experiments and simulations have established that dynamics in a class of
living and abiotic systems that are far from equilibrium exhibit super
diffusive behavior at long times, which in some cases (for example evolving
tumor) is preceded by slow glass-like dynamics. By using the evolution of a
collection of tumor cells, driven by mechanical forces and subject to cell
birth and apoptosis, as a case study we show theoretically that on short time
scales the mean square displacement is sub-diffusive due to jamming, whereas at
long times it is super diffusive. The results obtained using stochastic
quantization method, which is needed because of the absence of
fluctuation-dissipation theorem (FDT), show that the super-diffusive behavior
is universal and impervious to the nature of cell-cell interactions.
Surprisingly, the theory also quantitatively accounts for the non-trivial
dynamics observed in simulations of a model soap foam characterized by creation
and destruction of spherical bubbles, which suggests that the two
non-equilibrium systems belong to the same universality class. The theoretical
prediction for the super diffusion exponent is in excellent agreement with
simulations for collective motion of tumor cells and dynamics associated with
soap bubbles.


  The fundamental equations of various disciplines often seem to share the same
basic structure. Natural selection increases information in the same way that
Bayesian updating increases information. Thermodynamics and the forms of common
probability distributions express maximum increase in entropy, which appears
mathematically as loss of information. Physical mechanics follows paths of
change that maximize Fisher information. The information expressions typically
have analogous interpretations as the Newtonian balance between force and
acceleration, representing a partition between direct causes of change and
opposing changes in the frame of reference. This web of vague analogies hints
at a deeper common mathematical structure. I suggest that the Price equation
expresses that underlying universal structure. The abstract Price equation
describes dynamics as the change between two sets. One component of dynamics
expresses the change in the frequency of things, holding constant the values
associated with things. The other component of dynamics expresses the change in
the values of things, holding constant the frequency of things. The separation
of frequency from value generalizes Shannon's separation of the frequency of
symbols from the meaning of symbols in information theory. The Price equation's
generalized separation of frequency and value reveals a few simple invariances
that define universal geometric aspects of change. For example, the
conservation of total frequency, although a trivial invariance by itself,
creates a powerful constraint on the geometry of change. That constraint plus a
few others seem to explain the common structural forms of the equations in
different disciplines. From that abstract perspective, interpretations such as
selection, information, entropy, force, acceleration, and physical work arise
from the same underlying geometry expressed by the Price equation.


  The Internet of Bio-Nano-Things is a new paradigm that can bring novel
remotely controlled actuation and sensing techniques inside the human body.
Towards precise bionano sensing techniques in the brain, we investigate the
challenges of modelling spatial distribution of astrocyte networks in
developing a mathematical framework that lay the groundwork for future
early-detection techniques of neurodegenerative disease. In this paper, we
investigate the effect of the $\beta$-amyloid plaques in astrocytes with the
Alzheimer's disease. We developed a computation model of healthy and
Alzheimer's diseases astrocytes networks from the state of the art models and
results that account for the intracellular pathways, IP$_3$ dynamics, gap
junctions, voltage-gated calcium channels and astrocytes volumes. We also
implemented different types of astrocytes network topologies including shortcut
networks, regular degree networks, Erd\"os R\'enyi networks and link radius
networks. A proposed multi-scale stochastic computational model captures the
relationship between the intracellular and intercellular scales. Lastly, we
designed and evaluated a single-hop communication system with frequency
modulation using metrics such as propagation extend, molecular delay and
channel gain. The results show that the more unstable but at the same time
lower level oscillations of Alzheimer's astrocyte networks can create a
multi-scale effect on communication between astrocytes with increased molecular
delay and lower channel gain compared to healthy astrocytes, with an elevated
impact on Erd\"os R\'enyi networks and link radius networks topologies.


  Mechanism of vibronic amplification of transport of excitons was discussed in
relation to quantum photosynthesis. Vibrones (some modes of vibrations of
molecules) are observed experimentally in photosynthetic systems. In the
present paper we discuss a model of vibronic amplification of quantum transfer
where generation of vibrones as a coherent vibrational mode is described by an
analog of semiclassical theory of laser. We consider two models --- a model of
nonequilibrium three level system with vibronic mode, and some variant of a
model of laser without inversion. We conjecture that dark states discussed in
relation to quantum photosynthesis might be related to mechanism of vibronic
"laser" without inversion which amplifies the transfer of excitons. We prove
that in presence of vibronic mode transfer rate of excitons increases and
compute dependence of the transfer rate on parameters of the model.


  The optimal design of DNA origami systems that assemble rapidly and robustly
is hampered by the lack of a model for self-assembly that is sufficiently
detailed yet computationally tractable. Here, we propose a model for DNA
origami that strikes a balance between these two criteria by representing these
systems on a lattice at the level of binding domains. The free energy of
hybridization between individual binding domains is estimated with a
nearest-neighbour model. Double helical segments are treated as rigid rods, but
we allow flexibility at points where the backbone of one of the strands is
interrupted, which provides a reasonably realistic representation of partially
and fully assembled states. Particular attention is paid to the constraints
imposed by the double helical twist, as they determine where strand crossovers
between adjacent helices can occur. To improve the efficiency of sampling
configuration space, we develop Monte Carlo methods for sampling scaffold
conformations in near-assembled states, and we carry out simulations in the
grand canonical ensemble, enabling us to avoid considering states with unbound
staples. We demonstrate that our model can quickly sample assembled
configurations of a small origami design previously studied with the oxDNA
model, as well as a design with staples that span longer segments of the
scaffold. The sampling ability of our method should allow for good statistics
to be obtained when studying the assembly pathways, and is suited to
investigating in particular the effects of design and assembly conditions on
these pathways and their resulting final assembled structures.


  Precision medicine aims for personalized prognosis and therapeutics by
utilizing recent genome-scale high-throughput profiling techniques, including
next-generation sequencing (NGS). However, translating NGS data faces several
challenges. First, NGS count data are often overdispersed, requiring
appropriate modeling. Second, compared to the number of involved molecules and
system complexity, the number of available samples for studying complex
disease, such as cancer, is often limited, especially considering disease
heterogeneity. The key question is whether we may integrate available data from
all different sources or domains to achieve reproducible disease prognosis
based on NGS count data. In this paper, we develop a Bayesian Multi-Domain
Learning (BMDL) model that derives domain-dependent latent representations of
overdispersed count data based on hierarchical negative binomial factorization
for accurate cancer subtyping even if the number of samples for a specific
cancer type is small. Experimental results from both our simulated and NGS
datasets from The Cancer Genome Atlas (TCGA) demonstrate the promising
potential of BMDL for effective multi-domain learning without "negative
transfer" effects often seen in existing multi-task learning and transfer
learning methods.


  The impact of information dissemination on epidemic control is essentially
subject to individual behaviors. Unlike information-driven behaviors,
vaccination is determined by many cost-related factors, whose correlation with
the information dissemination should be better understood. To this end, we
propose an evolutionary vaccination game model in multiplex networks by
integrating an information-epidemic spreading process into the vaccination
dynamics, and explore how information dissemination influences vaccination. The
spreading process is described by a two-layer coupled
susceptible-alert-infected-susceptible (SAIS) model, where the strength
coefficient between two layers is defined to characterize the tendency and
intensity of information dissemination. We find that information dissemination
can increase the epidemic threshold, however, more information transmission
cannot promote vaccination. Specifically, increasing information dissemination
even leads to a decline of the vaccination equilibrium and raises the final
infection density. Moreover, we study the impact of strength coefficient and
individual sensitivity on social cost, and unveil the role of information
dissemination in controlling the epidemic with numerical simulations.


  Classical molecular dynamics (MD) simulations will be able to reach sampling
in the second timescale within five years, producing petabytes of simulation
data at current force field accuracy. Notwithstanding this, MD will still be in
the regime of low-throughput, high-latency predictions with average accuracy.
We envisage that machine learning (ML) will be able to solve both the accuracy
and time-to-prediction problem by learning predictive models using expensive
simulation data. The synergies between classical, quantum simulations and ML
methods, such as artificial neural networks, have the potential to drastically
reshape the way we make predictions in computational structural biology and
drug discovery.


  Head moves are a type of rearrangement moves for phylogenetic networks. They
have mostly been studied as part of more encompassing types of moves, such as
rSPR moves. Here, we study head moves as a type of moves on themselves. We show
that the tiers ($k>0$) of phylogenetic network space are connected by local
head moves. Then we show tail moves and head moves are closely related:
sequences of tail moves can be converted to sequences of head moves and vice
versa, changing the length by at most a constant factor. Because the tiers of
network space are connected by rSPR moves, this gives a second proof of the
connectivity of these tiers. Furthermore, we show that these tiers have small
diameter by reproving the connectivity a third time. As the head move
neighbourhood is in general small, this makes head moves a good candidate for
local search heuristics. Finally we prove that finding the shortest sequence of
head moves between two networks is NP-hard.


  The idea of incompetence as a learning or adaptation function was introduced
in the context of evolutionary games as a fixed parameter. However, live
organisms usually perform different nonlinear adaptation functions such as a
power law or exponential fitness growth. Here, we examine how the functional
form of the learning process may affect the social competition between
different behavioral types. Further, we extend our results for the evolutionary
games where fluctuations in the environment affect the behavioral adaptation of
competing species and demonstrate importance of the starting level of
incompetence for survival. Hence, we define a new concept of learning
advantages that becomes crucial when environments are constantly changing and
requiring rapid adaptation from species. This may lead to the evolutionarily
weak phase when even evolutionary stable populations become vulnerable to
invasions.


  The mobility of membrane proteins controls many biological functions. The
application of the model of Saffman and Delbr\"uck to the diffusion of membrane
proteins does not account for all the experimental measurements. These
discrepancies have triggered a lot of studies on the role of the mechanical
factors in the mobility. After a short review of the Saffman and Delbr\"uck
model and of some key experiments, we explore the various ways to incorporate
the effects of the different mechanical factors. Our approach focuses on the
coupling of the protein to the membrane, which is the central element in the
modelling. We present a general, polaron-like model, its recent application to
the mobility of a curvature sensitive protein, and its various extensions to
other couplings that may be relevant in future experiments.


  We study the most probable trajectories of the concentration evolution for
the transcription factor activator in a genetic regulation system, with
non-Gaussian stable L\'evy noise in the synthesis reaction rate taking into
account. We calculate the most probable trajectory by spatially maximizing the
probability density of the system path, i.e., the solution of the associated
nonlocal Fokker-Planck equation. We especially examine those most probable
trajectories from low concentration state to high concentration state (i.e.,
the likely transcription regime) for certain parameters, in order to gain
insights into the transcription processes and the tipping time for the
transcription likely to occur. This enables us: (i) to visualize the progress
of concentration evolution (i.e., observe whether the system enters the
transcription regime within a given time period); (ii) to predict or avoid
certain transcriptions via selecting specific noise parameters in particular
regions in the parameter space. Moreover, we have found some peculiar or
counter-intuitive phenomena in this gene model system, including (a) a smaller
noise intensity may trigger the transcription process, while a larger noise
intensity can not, under the same asymmetric L\'evy noise. This phenomenon does
not occur in the case of symmetric L\'evy noise; (b) the symmetric L\'evy
motion always induces transition to high concentration, but certain asymmetric
L\'evy motions do not trigger the switch to transcription. These findings
provide insights for further experimental research, in order to achieve or to
avoid specific gene transcriptions, with possible relevance for medical
advances.


  Cell fate determination, the process through which cells commit to
differentiated states is commonly mediated by gene regulatory motifs with
mutually exclusive expression states. The classical deterministic picture for
cell fate determination includes bistability and hysteresis, which enables the
persistence of the acquired cellular state after withdrawal of the stimulus,
ensuring a robust cellular response. However, the stochasticity inherent to
gene expression dynamics is not compatible with hysteresis, since the
stationary solution of the governing Chemical Master Equation does not depend
on the initial conditions. In this work, we provide a quantitative description
of a transient hysteresis phenomenon that reconciles experimental evidence of
hysteretic behaviour in gene regulatory networks with their inherent
stochasticity. Under sufficiently slow dynamics, the dependency of the
non-stationary solutions on the initial state of the cells can lead to what we
denote here as transient hysteresis. To quantify this phenomenon, we provide an
estimate of the convergence rate to the equilibrium. We also introduce the
equation of a natural landscape capturing the evolution of the system that,
unlike traditional cell fate potential landscapes, is compatible with the
notion of coexistence at the microscopic level.


  We studied acetylhistidine (AcH), bare or microsolvated with a zinc cation by
simulations in isolation. First, a global search for minima of the potential
energy surface combining both, empirical and first-principles methods, is
performed individually for either one of five possible protonation states.
Comparing the most stable structures between tautomeric forms of negatively
charged AcH shows a clear preference for conformers with the neutral imidazole
ring protonated at the N-epsilon-2 atom. When adding a zinc cation to the
system, the situation is reversed and N-delta-1-protonated structures are
energetically more favorable. Obtained minima structures then served as basis
for a benchmark study to examine the goodness of commonly applied levels of
theory, i.e. force fields, semi-empirical methods, density-functional
approximations (DFA), and wavefunction-based methods with respect to high-level
coupled-cluster calculations, i.e. the DLPNO-CCSD(T) method. All tested force
fields and semi-empirical methods show a poor performance in reproducing the
energy hierarchies of conformers, in particular of systems involving the zinc
cation. Meta-GGA, hybrid, double hybrid DFAs, and the MP2 method are able to
describe the energetics of the reference method within chemical accuracy, i.e.
with a mean absolute error of less than 1kcal/mol. Best performance is found
for the double hybrid DFA B3LYP+XYG3 with a mean absolute error of 0.7 kcal/mol
and a maximum error of 1.8 kcal/mol. While MP2 performs similarly as
B3LYP+XYG3, computational costs, i.e. timings, are increased by a factor of 4
in comparison due to the large basis sets required for accurate results.


  We study simple stochastic scenarios, based on birth-and-death Markovian
processes, that describe populations with Allee effect, to account for the role
of demographic stochasticity. In the mean-field deterministic limit we recover
well-known deterministic evolution equations widely employed in population
ecology. The mean-time to extinction is in general obtained by the
Wentzel-Kramers-Brillouin (WKB) approximation for populations with strong and
weak Allee effects. An exact solution for the mean time to extinction can be
found via a recursive equation for special cases of the stochastic dynamics. We
study the conditions for the validity of the WKB solution and analyze the
boundary between the weak and strong Allee effect by comparing exact solutions
with numerical simulations.


  Activation of most primary sensory neurons results in transduction currents
that are carried by cations. One notable exception is the vertebrate olfactory
receptor neuron (ORN), where the transduction current is carried largely by the
anion $Cl^-$. However, it remains unclear why ORNs use an anionic current for
signal amplification. We have sought to provide clarification on this topic by
studying the so far neglected dynamics of $Na^+$, $Ca^{2+}$, $K^+$ and $Cl^-$
in the small space of olfactory cilia during an odorant response. Using
computational modeling and simulations we compared the outcomes of signal
amplification based on either $Cl^-$ or $Na^+$ currents. We found that
amplification produced by $Na^+$ influx instead of a $Cl^-$ efflux is
problematic due to several reasons: First, the $Na^+$ current amplitude varies
greatly depending on mucosal ion concentration changes. Second, a $Na^+$
current leads to a large increase in the ciliary $Na^+$ concentration during an
odorant response. This increase inhibits and even reverses $Ca^{2+}$ clearance
by $Na^+/Ca^{2+}/K^+$ exchange, which is essential for response termination.
Finally, a $Na^+$ current increases the ciliary osmotic pressure, which could
cause swelling to damage the cilia. By contrast, a transduction pathway based
on $Cl^-$ efflux circumvents these problems and renders the odorant response
robust and reliable.


  Understanding wildlife-livestock interactions is crucial for the design and
management of protected areas that aim to conserve large mammal communities
undergoing conflicts with humans worldwide. An example of the need to quantify
the strength and direction of species interactions is the conservation of big
cats in newly established protected areas in China. Currently, free-ranging
livestock degrade the food and habitat of the endangered Amur tiger and Amur
leopard in the forest landscapes of Northeast China, but quantitative
assessments of how livestock affect the use of habitat by the major ungulate
prey of these predators are very limited. Here, we examined livestock-ungulate
interactions using large-scale camera-trap data in the newly established Tiger
and Leopard National Park in Northeast China, which borders Russia. We used
N-mixture models, two-species occupancy models and activity pattern overlap to
understand the effects of cattle grazing on three ungulate species (wild boar,
roe deer and sika deer) at a fine spatiotemporal scale. Our results showed that
incorporating the biotic interactions with cattle had significant negative
effects on encounters with three ungulates; sika deer were particularly
displaced as more cattle encroached on forest habitat, as they exhibited low
levels of co-occurrence with cattle in terms of habitat use. These results,
combined with spatiotemporal overlap, suggested fine-scale avoidance
behaviours, and they can help to refine strategies for the conservation of
tigers, leopards and their prey in human-dominated transboundary landscapes.
Progressively controlling cattle and the impact of cattle on biodiversity while
simultaneously addressing the economic needs of local communities should be key
priority actions for the Chinese government.


  The endangered Amur tiger is confined primarily to a narrow area along the
border with Russia in Northeast China. Little is known about the foraging
strategies of this small subpopulation in Hunchun Nature Reserve on the Chinese
side of the border; at this location, the prey base and land use patterns are
distinctly different from those in the larger population of the Sikhote-Alin
Mountains of Russia. Using dietary analysis of scats and camera-trapping data
from Hunchun Nature Reserve, we assessed spatiotemporal overlap of tigers and
their prey and identified prey selection patterns to enhance understanding of
the ecological requirements of tigers in Northeast China. Results indicated
that wild prey constituted 94.9% of the total biomass consumed by tigers;
domestic livestock represented 5.1% of the diet. Two species, wild boar and
sika deer , collectively represented 83% of the biomass consumed by tigers.
Despite lower spatial overlap of tigers and wild boar compared to tigers and
sika deer, tigers preferentially preyed on boar, likely facilitated by high
temporal overlap in activity patterns. Tigers exhibit significant spatial
overlap with sika deer, likely favoring a high level of tiger predation on this
large-sized ungulate. However, tigers did not preferred roe deer (Capreolus
pygargus) and showed a low spatial overlap with roe deer. Overall, our results
suggest that tiger prey selection is determined by prey body size and also
overlap in tiger and prey use of time or space. Also, we suggest that
strategies designed to minimize livestock forays into forested lands may be
important for decreasing the livestock depredation by tigers. This study offers
a framework to simultaneously integrate food habit analysis with the
distribution of predators and prey through time and space to provide a
comprehensive understanding of foraging strategies of large carnivores.


  Achieving a comprehensive understanding of the behaviour of proteins is
greatly facilitated by the knowledge of their structures, thermodynamics and
dynamics. All this information can be provided in an effective manner in terms
of structural ensembles. A structural ensemble can be obtained by determining
the structures, populations and interconversion rates for all the main states
that a protein can occupy. To achieve this goal, integrative methods that
combine experimental and computational approaches provide powerful tools. Here
we focus on cryo-electron microscopy, which has become over recent years an
invaluable resource to bridge the gap from order to disorder in structural
biology. In this review, we provide a perspective of the current challenges and
opportunities in determining protein structural ensembles using integrative
approaches that can combine cryo-electron microscopy data with other available
sources of information, along with an overview of the tools available to the
community.


  We develope a theoretical framework, based on exclusion process, that is
motivated by a biological phenomenon called transcript slippage (TS). In this
model a discrete lattice represents a DNA strand while each of the particles
that hop on it unidirectionally, from site to site, represents a RNA polymerase
(RNAP). While walking like a molecular motor along a DNA track in a
step-by-step manner, a RNAP simultaneously synthesizes a RNA chain; in each
forward step it elongates the nascent RNA molecule by one unit, using the DNA
track also as the template. At some special "slippery" position on the DNA,
which we represent as a defect on the lattice, a RNAP can lose its grip on the
nascent RNA and the latter's consequent slippage results in a final product
that is either longer or shorter than the corresponding DNA template. We
develope an exclusion model for RNAP traffic where the kinetics of the system
at the defect site captures key features of TS events. We demonstrate the
interplay of the crowding of RNAPs and TS. A RNAP has to wait at the defect
site for longer period in a more congested RNAP traffic, thereby increasing the
likelihood of its suffering a larger number of TS events. The qualitative
trends of some of our results for a simple special case of our model are
consistent with experimental observations. The general theoretical framework
presented here will be useful for guiding future experimental queries and for
analysis of the experimental data with more detailed versions of the same
model.


  Rapidly evolving pathogens like influenza viruses can persist by accumulating
antigenic novelty fast enough to evade the adaptive immunity of the host
population, yet without continuous accumulation of genetic diversity. This
dynamical state is often compared to the Red Queen evolving as fast as it can
just to maintain its foothold in the host population: Accumulation of antigenic
novelty is balanced by the build-up of host immunity. Such Red Queen States
(RQS) of continuous adaptation in large rapidly mutating populations are well
understood in terms of Traveling Wave (TW) theories of population genetics.
Here we shall make explicit the mapping of the established Multi-strain
Susceptible-Infected-Recovered (SIR) model onto the TW theory and demonstrate
that a pathogen can persist in RQS if cross-immunity is long-ranged and its
population size is large populations allowing for rapid adaptation. We then
investigate the stability of this state focusing on the rate of extinction and
the rate of "speciation" defined as antigenic divergence of viral strains
beyond the range of cross-inhibition. RQS states are transient, but in a
certain range of evolutionary parameters can exist for the time long compared
to the typical time to the most recent common ancestor ($T_{MRCA}$). In this
range the steady TW is unstable and the antigenic advance of the lead strains
relative to the typical co-circulating viruses tends to oscillate. This results
in large fluctuations in prevalence that facilitate extinction. We shall
demonstrate that the rate of TW fission into antigenically uncoupled viral
populations is related to fluctuations of $T_{MRCA}$ and construct a "phase
diagram" identifying different regimes of viral phylodynamics as a function of
evolutionary parameters.


  Nanoparticles coated with hydrophilic polymers often show a reduction in
unspecific interactions with the biological environment, which improves their
biocompatibility. The molecular determinants of this reduction are not very
well understood yet, and their knowledge may help improving nanoparticle
design. Here we address, using molecular dynamics simulations, the interactions
of human serum albumin, the most abundant serum protein, with two promising
hydrophilic polymers used for the coating of therapeutic nanoparticles,
poly(ethylene-glycol) and poly-sarcosine. By simulating the protein immersed in
a polymer-water mixture, we show that the two polymers have a very similar
affinity for the protein surface, both in terms of the amount of polymer
adsorbed and also in terms of the type of amino acids mainly involved in the
interactions. We further analyze the kinetics of adsorption and how it affects
the polymer conformations. Minor differences between the polymers are observed
in the thickness of the adsorption layer, that are related to the different
degree of flexibility of the two molecules. In comparison poly-alanine, an
isomer of poly-sarcosine known to self-aggregate and induce protein
aggregation, shows a significantly larger affinity for the protein surface than
PEG and PSar, which we show to be related not to a different patterns of
interactions with the protein surface, but to the different way the polymer
interacts with water.


  Spatial heterogeneity can have dramatic effects on the biochemical networks
that drive cell regulation and decision-making. For this reason, a number of
methods have been developed to model spatial heterogeneity and incorporated
into widely used modeling platforms. Unfortunately, the standard approaches for
specifying and simulating chemical reaction networks become untenable when
dealing with multi-state, multi-component systems that are characterized by
combinatorial complexity. To address this issue, we developed MCell-R, a
framework that extends the particle-based spatial Monte Carlo simulator, MCell,
with the rule-based model specification and simulation capabilities provided by
BioNetGen and NFsim. The BioNetGen syntax enables the specification of
biomolecules as structured objects whose components can have different internal
states that represent such features as covalent modification and conformation
and which can bind components of other molecules to form molecular complexes.
The network-free simulation algorithm used by NFsim enables efficient
simulation of rule-based models even when the size of the network implied by
the biochemical rules is too large to enumerate explicitly, which frequently
occurs in detailed models of biochemical signaling. The result is a framework
that can efficiently simulate systems characterized by combinatorial complexity
at the level of spatially-resolved individual molecules over biologically
relevant time and length scales.


  The problem of controlling the mean and the variance of a species of interest
in a simple gene expression is addressed. It is shown that the protein mean
level can be globally and robustly tracked to any desired value using a simple
PI controller that satisfies certain sufficient conditions. Controlling both
the mean and variance however requires an additional control input, e.g. the
mRNA degradation rate, and local robust tracking of mean and variance is proved
to be achievable using multivariable PI control, provided that the reference
point satisfies necessary conditions imposed by the system. Even more
importantly, it is shown that there exist PI controllers that locally, robustly
and simultaneously stabilize all the equilibrium points inside the admissible
region. The results are then extended to the mean control of a gene expression
with protein dimerization. It is shown that the moment closure problem can be
circumvented without invoking any moment closure technique. Local stabilization
and convergence of the average dimer population to any desired reference value
is ensured using a pure integral control law. Explicit bounds on the controller
gain are provided and shown to be valid for any reference value. As a
byproduct, an explicit upper-bound of the variance of the monomer species,
acting on the system as unknown input due to the moment openness, is obtained.
The results are illustrated by simulation.


  There is a general consensus among archaeologists that replacement of
Neanderthals by anatomically modern humans in Europe occurred around 40K to 35K
YBP. However, the causal mechanism for this replacement continues to be
debated. Searching for specific fitness advantages in the archaeological record
has proven difficult, as these may be obscured, absent, or subject to
interpretation. Proposed models have therefore featured either fitness
advantages in favor of anatomically modern humans, or invoked neutral drift
under various preconditions. To bridge this gap, we rigorously compare the
system-level properties of fitness- and drift-based explanations of Neanderthal
replacement. Our stochastic simulations and analytical predictions show that,
although both fitness and drift can produce fixation, they present important
differences in 1) required initial conditions, 2) reliability, 3) time to
replacement, and 4) path to replacement (population histories). These results
present useful opportunities for comparison with archaeological and genetic
data. We find far greater agreement between the available empirical evidence
and the system-level properties of replacement by differential fitness, rather
than by neutral drift.


  Data from seven data storage tags recovered from Atlantic salmon marked as
smolts were analyzed for depth movements and patterns of deep diving during the
marine migration. The salmon mostly stayed at the surface and showed diurnal
activity especially from autumn until spring. During the first months at sea
the salmon stayed at shallower depths (< 100 m). The salmon took short deep
dives (> 100 m), that were rare or absent during the first summer at sea but
increased in frequency and duration especially in late winter. The maximum
depth of the dives varied from 419 to 1187 m. Most of dives were short, (< 5
hours) but could last up to 33 hours. The duration of dives increased in late
winter until spring and the overall depth and maximum depth per dive increased
exponentially over time. The initiation of the dives was more common in
evenings and at night, suggesting nocturnal diving. We hypothesized that deep
diving is related to feeding of salmon as mesopelagic fish can be important
food for salmon during winter.


  In this paper we present a framework for investigating coloured noise in
reaction-diffusion systems. We start by considering a deterministic
reaction-diffusion equation and show how external forcing can cause temporally
correlated or coloured noise. Here, the main source of external noise is
considered to be fluctuations in the parameter values representing the inflow
of particles to the system. First, we determine which reaction systems, driven
by extrinsic noise, can admit only one steady state, so that effects, such as
stochastic switching, are precluded from our analysis. To analyse the steady
state behaviour of reaction systems, even if the parameter values are changing,
necessitates a parameter-free approach, which has been central to algebraic
analysis in chemical reaction network theory. To identify suitable models we
use tools from real algebraic geometry that link the network structure to its
dynamical properties. We then make a connection to internal noise models and
show how power spectral methods can be used to predict stochastically driven
patterns in systems with coloured noise. In simple cases we show that the power
spectrum of the coloured noise process and the power spectrum of the
reaction-diffusion system modelled with white noise multiply to give the power
spectrum of the coloured noise reaction-diffusion system.


  We consider a stochastic model for evolution of group-structured populations
in which interactions between group members correspond to the Prisoner's
Dilemma or the Hawk-Dove game. Selection operates at two organization levels:
individuals compete with peer group members based on individual payoff, while
groups also compete with other groups based on average payoff of group members.
In the Prisoner's Dilemma, this creates a tension between the two levels of
selection, as defectors are favored at the individual level, whereas groups
with at least some cooperators outperform groups of defectors at the
between-group level. In the limit of infinite group size and infinite number of
groups, we derive a non-local PDE that describes the probability distribution
of group compositions in the population. For special families of payoff
matrices, we characterize the long-time behavior of solutions of our equation,
finding a threshold level of between-group selection required to sustain
density steady states and the survival of cooperation. When all-cooperator
groups are most fit, the average and most abundant group composition at steady
state range from featuring all-defector groups when individual-level selection
dominates to featuring all-cooperator groups when group-level selection
dominates. When the most fit groups have a mix of cooperators and defectors,
then the average and most abundant group compositions always feature a smaller
fraction of cooperators than required for the optimal mix, even in the limit
where group-level selection is infinitely stronger than individual-level
selection. In such cases, the conflict between the two levels of selection
cannot be decoupled, and cooperation cannot be sustained at all in the case
when between-group competition favors an even mix of cooperators and defectors.


  Changes in morphological characteristics of cells from two cultured cancer
cell lines, HeLa and A549, induced by photodynamic treatment with Radachlorin
photosensitizer have been monitored using digital holographic microscopy. The
observed dose-dependent post-treatment dynamics of phase shift variations
demonstrated several scenarios of cell death. In particular the phase shift
increase at low doses can be associated with apoptosis while its decrease at
high doses can be associated with necrosis. Two cell types were shown to be
differently responsive to treatment at the same doses. Although the sequence of
death scenarios with increasing irradiation dose was demonstrated to be the
same, each specific scenario was realized at substantially different doses.
Results obtained by holographic microscopy were confirmed by confocal
fluorescence microscopy with the commonly used test assay.


  Stochastic fluctuations in the copy number of gene products have perceivable
effects on the functioning of gene regulatory networks (GRN). The Master
equation (ME) provides a theoretical basis for studying such effects. However,
solving the ME can be a task that ranges from simple to difficult to impossible
using conventional methods. Therefore, discovering new techniques for solving
the ME is an important part of research on stochastic GRN. In this paper, we
present a novel approach to obtaining the generating function (GF), which
contains the same information as the ME, for a one gene system that includes
multi-step post-transcription and post-translation processes. The novelty of
the approach lies in the separation of the mRNAs from proteins. The GF for the
mRNAs is obtained using a formalism involving operators and vector states.
Using the same formalism, the GF for the proteins is solved for a particular
path taken by all mRNAs in the time-copy number plane; then, the GF is summed
over all possible paths. We prove a theorem that shows the summation of all
paths to be equivalent to an equation similar to the ME for the mRNAs. On a
system with six gene products in total and randomly selected initial
conditions, we confirm the validity of our results by comparing them with
Gillespie simulations.


  Population structure can be modelled by evolutionary graphs, which can have a
substantial, but very subtle influence on the fate of the arising mutants.
Individuals are located on the nodes of these graphs, competing with each other
to eventually take over the graph via the links. Many applications for this
framework can be envisioned, from the ecology of river systems and cancer
initiation in colonic crypts to biotechnological search for optimal mutations.
In all these applications, it is not only important where and when novel
variants arise and how likely it is that they ultimately take over, but also
how long this process takes. More concretely, how is the probability to take
over the population related to the associated time? We study this problem for
all possible undirected and unweighted graphs up to a certain size. To move
beyond the graph size where an exhaustive search is possible, we devise a
genetic algorithm to find graphs with either high or low fixation probability
and either short or long fixation time and study their structure in detail
searching for common themes. Our work unravels structural properties that
maximize or minimize fixation probability and time, which allows us to
contribute to a first map of the universe of evolutionary graphs.


  Microbial genome web portals have a broad range of capabilities that address
a number of information-finding and analysis needs for scientists. This article
compares the capabilities of the major microbial genome web portals to aid
researchers in determining which portal(s) are best suited to solving their
information-finding and analytical needs. We assessed both the bioinformatics
tools and the data content of BioCyc, KEGG, Ensembl Bacteria, KBase, IMG, and
PATRIC. For each portal, our assessment compared and tallied the available
capabilities. The strengths of BioCyc include its genomic and metabolic tools,
multi-search capabilities, table-based analysis tools, regulatory network tools
and data, omics data analysis tools, breadth of data content, and large amount
of curated data. The strengths of KEGG include its genomic and metabolic tools.
The strengths of Ensembl Bacteria include its genomic tools and large number of
genomes. The strengths of KBase include its genomic tools and metabolic models.
The strengths of IMG include its genomic tools, multi-search capabilities,
large number of genomes, table-based analysis tools, and breadth of data
content. The strengths of PATRIC include its large number of genomes,
table-based analysis tools, metabolic models, and breadth of data content.


  Next-generation sequencing technology enables routine detection of bacterial
pathogens for clinical diagnostics and genetic research. Whole genome
sequencing has been of importance in the epidemiologic analysis of bacterial
pathogens. However, few whole genome sequencing-based genotyping pipelines are
available for practical applications. Here, we present the whole genome
sequencing-based single nucleotide polymorphism (SNP) genotyping method and
apply to the evolutionary analysis of methicillin-resistant Staphylococcus
aureus. The SNP genotyping method calls genome variants using next-generation
sequencing reads of whole genomes and calculates the pair-wise Jaccard
distances of the genome variants. The method may reveal the high-resolution
whole genome SNP profiles and the structural variants of different isolates of
methicillin-resistant S. aureus (MRSA) and methicillin-susceptible S. aureus
(MSSA) strains. The phylogenetic analysis of whole genomes and particular
regions may monitor and track the evolution and the transmission dynamic of
bacterial pathogens. The computer programs of the whole genome sequencing-based
SNP genotyping method are available to the public at
https://github.com/cyinbox/NGS.


  We consider an asexually reproducing population on a finite type space whose
evolution is driven by exponential birth, death and competition rates, as well
as the possibility of mutation at a birth event. On the individual-based level
this population can be modelled as a measure-valued Markov process. Multiple
variations of this system have been studied in the simultaneous limit of large
populations and rare mutations, where the regime is chosen such that mutations
are separated. We consider the deterministic system, resulting from the large
population limit, and then let the mutation probability tend to zero. This
corresponds to a much higher frequency of mutations, where multiple microscopic
types are present at the same time. The limiting process resembles an adaptive
walk or flight and jumps between different equilibria of coexisting types. The
graph structure on the type space, determined by the possibilities to mutate,
plays an important role in defining this jump process. In a variation of the
above model, where the radius in which mutants can be spread is limited, we
study the possibility of crossing valleys in the fitness landscape and derive
different kinds of limiting walks.


  Dense packing of hydrophobic residues in the cores of globular proteins
determines their stability. Recently, we have shown that protein cores possess
packing fraction $\phi \approx 0.56$, which is the same as dense, random
packing of amino acid-shaped particles. In this article, we compare the
structural properties of protein cores and jammed packings of amino acid-shaped
particles in much greater depth by measuring their local and connected void
regions. We find that the distributions of surface Voronoi cell volumes and
local porosities obey similar statistics in both systems. We also measure the
probability that accessible, connected void regions percolate as a function of
the size of a spherical probe particle and show that both systems possess the
same critical probe size. By measuring the critical exponent $\tau$ that
characterizes the size distribution of connected void clusters at the onset of
percolation, we show that void percolation in packings of amino acid-shaped
particles and protein cores belong to the same universality class, which is
different from that for void percolation in jammed sphere packings. We propose
that the connected void regions of proteins are a defining feature of proteins
and can be used to differentiate experimentally observed proteins from decoy
structures that are generated using computational protein design software. This
work emphasizes that jammed packings of amino acid-shaped particles can serve
as structural and mechanical analogs of protein cores, and could therefore be
useful in modeling the response of protein cores to cavity-expanding and
-reducing mutations.


  The activity of biological cells is primarily based on chemical reactions and
typically modeled as a reaction-diffusion system. Cells are, however, highly
crowded with macromolecules, including a variety of molecular machines such as
enzymes. The working cycles of these machines are often coupled with their
internal motion (conformational changes). In the crowded environment of a cell,
motion interference between neighboring molecules is not negligible, and this
interference can affect the reaction dynamics through machine operation. To
simulate such a situation, we propose a reaction-diffusion model consisting of
particles whose shape depends on an internal state variable, for crowds of
nano- to micro-machines. The interference between nearby particles is naturally
introduced through excluded volume repulsion. In the simulations, we observed
segregation and flow-like patterns enhanced by crowding out of relevant
molecules, as well as molecular synchronization waves and phase transitions.
The presented model is simple and extensible for diverse molecular machinery,
and may serve as a framework to study the interplay between the mechanical
stress/strain network and the chemical reaction network in the cell.
Applications to more macroscopic systems, e.g., crowds of cells, are also
discussed.


  TF-MoDISco (Transcription Factor Motif Discovery from Importance Scores) is
an algorithm for identifying motifs from basepair-level importance scores
computed on genomic sequence data. This technical note focuses on version
v0.5.6.5. The implementation is available at
https://github.com/kundajelab/tfmodisco/tree/v0.5.6.5