  DNA bending is biologically important for genome regulation and is relevant
to a range of nanotechnological systems. Recent results suggest that sharp
bending is much easier than implied by the widely-used worm-like chain model;
many of these studies, however, remain controversial. We use a coarse-grained
model, previously fitted to DNA's basic thermodynamic and mechanical
properties, to explore strongly bent systems. We find that as the end-to-end
distance is decreased sufficiently short duplexes undergo a transition to a
state in which the bending strain is localized at a flexible kink that involves
disruption of base-pairing and stacking. This kinked state, which is not
well-described by the worm-like chain model, allows the duplex to more easily
be sharply bent. It is not completely flexible, however, due to constraints
arising from the connectivity of both DNA backbones. We also perform a detailed
comparison to recent experiments on a "molecular vice" that probes highly bent
DNA. Close agreement between simulations and experiments strengthens the
hypothesis that localised bending via kinking occurs in the molecular vice and
causes enhanced flexibility of duplex DNA. Our calculations therefore suggests
that the cost of kinking implied by this experiment is consistent with the
known thermodynamic and mechanical properties of DNA.


  DNA cyclization is a powerful technique to gain insight into the nature of
DNA bending. The worm-like chain model provides a good description of small to
moderate bending fluctuations, but some experiments on strongly-bent shorter
molecules suggest enhanced flexibility over and above that expected from the
worm-like chain. Here, we use a coarse-grained model of DNA to investigate the
thermodynamics of DNA cyclization for molecules with less than 210 base pairs.
As the molecules get shorter we find increasing deviations between our computed
equilibrium j-factor and the worm-like chain predictions of Shimada and
Yamakawa. These deviations are due to sharp kinking, first at nicks, and only
subsequently in the body of the duplex. At the shortest lengths, substantial
fraying at the ends of duplex domains is the dominant method of relaxation. We
also estimate the dynamic j-factor measured in recent FRET experiments. We find
that the dynamic j-factor is systematically larger than its equilibrium
counterpart, with the deviation larger for shorter molecules, because not all
the stress present in the fully cyclized state is present in the transition
state. These observations are important for the interpretation of recent
experiments, as only kinking within the body of the duplex is genuinely
indicative of non-worm-like chain behaviour.


  Elastic network models (ENMs) are valuable and efficient tools for
characterizing the collective internal dynamics of proteins based on the
knowledge of their native structures. The increasing evidence that the
biological functionality of RNAs is often linked to their innate internal
motions, poses the question of whether ENM approaches can be successfully
extended to this class of biomolecules. This issue is tackled here by
considering various families of elastic networks of increasing complexity
applied to a representative set of RNAs. The fluctuations predicted by the
alternative ENMs are stringently validated by comparison against extensive
molecular dynamics simulations and SHAPE experiments. We find that simulations
and experimental data are systematically best reproduced by either an all-atom
or a three-beads-per-nucleotide representation (sugar-base-phosphate), with the
latter arguably providing the best balance of accuracy and computational
complexity.


  Genotypic fitness landscapes are constructed by assessing the fitness of all
possible combinations of a given number of mutations. In the last years,
several experimental fitness landscapes have been completely resolved. As
fitness landscapes are high-dimensional, their characterization relies on
simple measures of their structure, which can be used as statistics in
empirical applications. Here we propose two new sets of measures that
explicitly capture two relevant features of fitness landscapes: epistasis and
constraints. The first set contains new measures for epistasis based on the
correlation of fitness effects of mutations. They have a natural
interpretation, capture well the interaction between mutations, can be obtained
analytically for most landscape models and can therefore be used to
discriminate between different models. The second set contains measures of
evolutionary constraints based on "chains" of forced mutations along
fitness-increasing paths. Some of these measures are non-monotonic in the
amount of epistatic interactions, but have instead a maximum for intermediate
values. We further characterize the relationships of these measures to the ones
that were previous proposed (e.g. number of peaks, roughness/slope, fraction of
non-additive components, etc). Finally, we show how these measures can help
uncovering the amount and the nature of epistatic interactions in two
experimental landscapes.


  The observed single-handedness of biological amino acids and sugars has long
been attributed to autocatalysis. However, the stability of homochiral states
in deterministic autocatalytic systems relies on cross inhibition of the two
chiral states, an unlikely scenario for early life self-replicators. Here, we
present a theory for a stochastic individual-level model of autocatalysis due
to early life self-replicators. Without chiral inhibition, the racemic state is
the global attractor of the deterministic dynamics, but intrinsic
multiplicative noise stabilizes the homochiral states, in both well-mixed and
spatially-extended systems. We conclude that autocatalysis is a viable
mechanism for homochirality, without imposing additional nonlinearities such as
chiral inhibition.


  We study the effect of intrinsic noise on the thermodynamic balance of
complex chemical networks subtending cellular metabolism and gene regulation. A
topological network property called deficiency, known to determine the
possibility of complex behavior such as multistability and oscillations, is
shown to also characterize the entropic balance. In particular, only when
deficiency is zero does the average stochastic dissipation rate equal that of
the corresponding deterministic model, where correlations are disregarded. In
fact, dissipation can be reduced by the effect of noise, as occurs in a toy
model of metabolism that we employ to illustrate our findings. This phenomenon
highlights that there is a close interplay between deficiency and the
activation of new dissipative pathways at low molecule numbers.


  Ligand diffusion through proteins is a fundamental process governing
biological signaling and enzymatic catalysis. The complex topology of protein
tunnels results in difficulties with computing ligand escape pathways by
standard molecular dynamics (MD) simulations. Here, two novel methods for
searching of ligand exit pathways and cavity exploration are proposed: memory
random acceleration MD (mRAMD), and memetic algorithms (MA). In mRAMD, finding
exit pathways is based on a non-Markovian biasing that is introduced to
optimize the unbinding force. In MA, hybrid learning protocols are exploited to
predict optimal ligand exit paths. The methods are tested on three proteins
with increasing complexity of tunnels: M2 muscarinic receptor, nitrile
hydratase, and cytochrome P450cam. In these cases, the proposed methods
outperform standard techniques that are used currently to find ligand egress
pathways. The proposed approach is general and appropriate for accelerated
transport of an object through a network of protein tunnels.


  Although the concepts of age, survival and transit time have been widely used
in many fields, including population dynamics, chemical engineering, and
hydrology, a comprehensive mathematical framework is still missing. Here we
discuss several relationships among these quantities by starting from the
evolution equation for the joint distribution of age and survival, from which
the equations for age and survival time readily follow. It also becomes
apparent how the statistical dependence between age and survival is directly
related to either the age-dependence of the loss function or the survival-time
dependence of the input function. The solution of the joint distribution
equation also allows us to obtain the relationships between the age at exit (or
death) and the survival time at input (or birth), as well as to stress the
symmetries of the various distributions under time reversal. The transit time
is then obtained as a sum of the age and survival time, and its properties are
discussed along with the general relationships between their mean values. The
special case of steady state case is analyzed in detail. Some examples,
inspired by hydrologic applications, are presented to illustrate the theory
with the specific results.


  The dynamics of a population undergoing selection is a central topic in
evolutionary biology. This question is particularly intriguing in the case
where selective forces act in opposing directions at two population scales. For
example, a fast-replicating virus strain outcompetes slower-replicating strains
at the within-host scale. However, if the fast-replicating strain causes host
morbidity and is less frequently transmitted, it can be outcompeted by
slower-replicating strains at the between-host scale. Here we consider a
stochastic ball-and-urn process which models this type of phenomenon. We prove
the weak convergence of this process under two natural scalings. The first
scaling leads to a deterministic nonlinear integro-partial differential
equation on the interval $[0,1]$ with dependence on a single parameter,
$\lambda$. We show that the fixed points of this differential equation are Beta
distributions and that their stability depends on $\lambda$ and the behavior of
the initial data around $1$. The second scaling leads to a measure-valued
Fleming-Viot process, an infinite dimensional stochastic process that is
frequently associated with a population genetics.


  Microbial colonies are experimental model systems for studying the
colonization of new territory by biological species through range expansion. We
study a generalization of the two-species Eden model, which incorporates local
frequency-dependent selection, in order to analyze how social interactions
between two species influence surface roughness of growing microbial colonies.
The model includes several classical scenarios from game theory. We then
concentrate on an expanding public goods game, where either cooperators or
defectors take over the front depending on the system parameters. We analyze in
detail the critical behavior of the nonequilibrium phase transition between
global cooperation and defection and thereby identify a new universality class
of phase transitions dealing with absorbing states. At the transition, the
number of boundaries separating sectors decays with a novel power law in time
and their superdiffusive motion crosses over from Eden scaling to a nearly
ballistic regime. In parallel, the width of the front initially obeys Eden
roughening and, at later times, passes over to selective roughening.


  Image scoring sustains cooperation in the repeated two-player prisoner's
dilemma through indirect reciprocity, even though defection is the uniquely
dominant selfish behaviour in the one-shot game. Many real-world dilemma
situations, however, firstly, take place in groups and, secondly, lack the
necessary transparency to inform subjects reliably of others' individual past
actions. Instead, there is revelation of information regarding groups, which
allows for `group scoring' but not for image scoring. Here, we study how
sensitive the positive results related to image scoring are to information
based on group scoring. We combine analytic results and computer simulations to
specify the conditions for the emergence of cooperation. We show that under
pure group scoring, that is, under the complete absence of image-scoring
information, cooperation is unsustainable. Away from this extreme case,
however, the necessary degree of image scoring relative to group scoring
depends on the population size and is generally very small. We thus conclude
that the positive results based on image scoring apply to a much broader range
of informational settings that are relevant in the real world than previously
assumed.


  The adaptive immune system relies on the diversity of receptors expressed on
the surface of B and T-cells to protect the organism from a vast amount of
pathogenic threats. The proliferation and degradation dynamics of different
cell types (B cells, T cells, naive, memory) is governed by a variety of
antigenic and environmental signals, yet the observed clone sizes follow a
universal power law distribution. Guided by this reproducibility we propose
effective models of somatic evolution where cell fate depends on an effective
fitness. This fitness is determined by growth factors acting either on clones
of cells with the same receptor responding to specific antigens, or directly on
single cells with no regards for clones. We identify fluctuations in the
fitness acting specifically on clones as the essential ingredient leading to
the observed distributions. Combining our models with experiments we
characterize the scale of fluctuations in antigenic environments and we provide
tools to identify the relevant growth signals in different tissues and
organisms. Our results generalize to any evolving population in a fluctuating
environment.


  We study the dynamics of an ideal polymer chain in a crowded, viscoelastic
medium and in the presence of active forces. The motion of the centre of mass
and of individual monomers is calculated. On time scales that are comparable to
the persistence time of the active forces, monomers can move superdiffusively
while on larger time scales subdiffusive behaviour occurs. The difference
between this subdiffusion and that in absence of active forces is quantified.
We show that the polymer swells in response to active processes and determine
how this swelling depends on the viscoelastic properties of the environment.
Our results are compared to recent experiments on the motion of chromosomal
loci in bacteria.


  The molecular dynamics simulation package GROMACS runs efficiently on a wide
variety of hardware from commodity workstations to high performance computing
clusters. Hardware features are well exploited with a combination of SIMD,
multi-threading, and MPI-based SPMD/MPMD parallelism, while GPUs can be used as
accelerators to compute interactions offloaded from the CPU. Here we evaluate
which hardware produces trajectories with GROMACS 4.6 or 5.0 in the most
economical way. We have assembled and benchmarked compute nodes with various
CPU/GPU combinations to identify optimal compositions in terms of raw
trajectory production rate, performance-to-price ratio, energy efficiency, and
several other criteria. Though hardware prices are naturally subject to trends
and fluctuations, general tendencies are clearly visible. Adding any type of
GPU significantly boosts a node's simulation performance. For inexpensive
consumer-class GPUs this improvement equally reflects in the
performance-to-price ratio. Although memory issues in consumer-class GPUs could
pass unnoticed since these cards do not support ECC memory, unreliable GPUs can
be sorted out with memory checking tools. Apart from the obvious determinants
for cost-efficiency like hardware expenses and raw performance, the energy
consumption of a node is a major cost factor. Over the typical hardware
lifetime until replacement of a few years, the costs for electrical power and
cooling can become larger than the costs of the hardware itself. Taking that
into account, nodes with a well-balanced ratio of CPU and consumer-class GPU
resources produce the maximum amount of GROMACS trajectory over their lifetime.


  The goal of this paper is to prove rigorous results for the behavior of
genealogies in a one-dimensional long range biased voter model introduced by
Hallatschek and Nelson [25]. The first step, which is easily accomplished using
results of Mueller and Tribe [38], is to show that when space and time are
rescaled correctly, our biased voter model converges to a Wright-Fisher SPDE. A
simple extension of a result of Durrett and Restrepo [18] then shows that the
dual branching coalescing random walk converges to a branching Brownian motion
in which particles coalesce after an exponentially distributed amount of
intersection local time. Brunet et al. [8] have conjectured that genealogies in
models of this type are described by the Bolthausen-Sznitman coalescent, see
[39]. However, in the model we study there are no simultaneous coalescences.
Our third and most significant result concerns "tracer dynamics" in which some
of the initial particles in the biased voter model are labeled. We show that
the joint distribution of the labeled and unlabeled particles converges to the
solution of a system of stochastic partial differential equations. A new
duality equation that generalizes the one Shiga [44] developed for the
Wright-Fisher SPDE is the key to the proof of that result.


  Our main purpose is to compare classical nonself-centered, two-signal
theoretical models of the adaptive immune system with a novel, self-centered,
one-signal model developed by our research group. Our model hypothesizes that
the immune system of a fetus is capable learning the limited set of self
antigens but unable to prepare itself for the unlimited variety of nonself
antigens. We have built a computational model that simulates the development of
the adaptive immune system. For simplicity, we concentrated on humoral immunity
and its major components: T cells, B cells, antibodies, interleukins,
non-immune self cells, and foreign antigens. Our model is a microscopic one,
similar to the interacting particle models of statistical physics and
agent-based models in immunology. Furthermore, our model is stochastic: events
are considered random and modeled by a continuous time, finite state Markov
process, that is, they are controlled by finitely many independent exponential
clocks.
  We investigate under what conditions can an immune memory be created that
results in a more effective immune response to a repeated infection. The
simulations show that our self-centered model is realistic. Moreover, in case
of a primary adaptive immune reaction, it can destroy infections more
efficiently than a classical nonself-centered model.
  Predictions of our theoretical model were clinically supported by
autoimmune-related adverse events in high-dose immune checkpoint inhibitor
immunotherapy trials and also by safe and successful low-dose immune checkpoint
inhibitor combination treatment of heavily pretreated stage IV cancer patients
who had exhausted all conventional treatments.
  The MiStImm simulation tool and source codes are available at the address
https://github.com/kerepesi/MiStImm.


  Circadian clocks exhibit the robustness of period and plasticity of phase
against environmental changes such as temperature and nutrient conditions. Thus
far, however, it is unclear how both are simultaneously achieved. By
investigating distinct models of circadian clocks, we demonstrate reci- procity
between robustness and plasticity: higher robustness in the period implies
higher plasticity in the phase, where changes in period and in phase follow a
linear relationship with a negative coef- ficient. The robustness of period is
achieved by the adaptation on the limit cycle via a concentration change of a
buffer molecule, whose temporal change leads to a phase shift following a shift
of the limit-cycle orbit in phase space. Generality of reciprocity in clocks
with the adaptation mechanism is confirmed with theoretical analysis of simple
models, while biological significance is discussed.


  High-throughput experimental techniques and bioinformatics tools make it
possible to obtain reconstructions of the metabolism of microbial species.
Combined with mathematical frameworks such as flux balance analysis, which
assumes that nutrients are used so as to maximize growth, these reconstructions
enable us to predict microbial growth.
  Although such predictions are generally accurate, these approaches do not
give insights on how different nutrients are used to produce growth, and thus
are difficult to generalize to new media or to different organisms.
  Here, we propose a systems-level phenomenological model of metabolism
inspired by the virial expansion. Our model predicts biomass production given
the nutrient uptakes and a reduced set of parameters, which can be easily
determined experimentally. To validate our model, we test it against in silico
simulations and experimental measurements of growth, and find good agreement.
From a biological point of view, our model uncovers the impact that individual
nutrients and the synergistic interaction between nutrient pairs have on
growth, and suggests that we can understand the growth maximization principle
as the optimization of nutrient synergies.


  We present a single-molecule tool called the CoPro (Concentration of
Proteins) method that uses millisecond imaging with convolution analysis,
automated image segmentation and super-resolution localization microscopy to
generate robust estimates for protein concentration in different compartments
of single living cells, validated using realistic simulations of complex
multiple compartment cell types. We demonstrates its utility experimentally on
model Escherichia coli bacteria and Saccharomyces cerevisiae budding yeast
cells, and use it to address the biological question of how signals are
transduced in cells. Cells in all domains of life dynamically sense their
environment through signal transduction mechanisms, many involving gene
regulation. The glucose sensing mechanism of S. cerevisiae is a model system
for studying gene regulatory signal transduction. It uses the multi-copy
expression inhibitor of the GAL gene family, Mig1, to repress unwanted genes in
the presence of elevated extracellular glucose concentrations. We fluorescently
labelled Mig1 molecules with green fluorescent protein (GFP) via chromosomal
integration at physiological expression levels in living S. cerevisiae cells,
in addition to the RNA polymerase protein Nrd1 with the fluorescent protein
reporter mCherry. Using CoPro we make quantitative estimates of Mig1 and Nrd1
protein concentrations in the cytoplasm and nucleus compartments on a
cell-by-cell basis under physiological conditions. These estimates indicate a
4-fold shift towards higher values in concentration of diffusive Mig1 in the
nucleus if the external glucose concentration is raised, whereas equivalent
levels in the cytoplasm shift to smaller values with a relative change an order
of magnitude smaller. This compares with Nrd1 which is not involved directly in
glucose sensing, which is almost exclusively localized in the nucleus under
high and...


  In this paper we propose a wet lab algorithm for prediction of radiation fog
by DNA computing. The concept of DNA computing is essentially exploited for
generating the classifier algorithm in the wet lab. The classifier is based on
a new concept of similarity based fuzzy reasoning suitable for wet lab
implementation. This new concept of similarity based fuzzy reasoning is
different from conventional approach to fuzzy reasoning based on similarity
measure and also replaces the logical aspect of classical fuzzy reasoning by
DNA chemistry. Thus, we add a new dimension to existing forms of fuzzy
reasoning by bringing it down to nanoscale. We exploit the concept of massive
parallelism of DNA computing by designing this new classifier in the wet lab.
This newly designed classifier is very much generalized in nature and apart
from prediction of radiation fog this methodology can be applied to other types
of data also. To achieve our goal we first fuzzify the given observed
parameters in a form of synthetic DNA sequence which is called fuzzy DNA and
which handles the vague concept of human reasoning.


  MicroRNAs (miRNAs) are a class of small noncoding RNAs that can regulate many
genes by base pairing to sites in mRNAs. The functionality of miRNAs overlaps
that of short interfering RNAs (siRNAs), and many features of miRNA targeting
have been revealed experimentally by studying miRNA-mimicking siRNAs. This
review outlines the features associated with animal miRNA targeting and
describes currently available prediction tools.


  The most prominent property of life on Earth is its ability to evolve. It is
often taken for granted that self-replication--the characteristic that makes
life possible--implies evolvability, but many examples such as the lack of
evolvability in computer viruses seem to challenge this view. Is evolvability
itself a property that needs to evolve, or is it automatically present within
any chemistry that supports sequences that can evolve in principle? Here, we
study evolvability in the digital life system Avida, where self-replicating
sequences written by hand are used to seed evolutionary experiments. We use 170
self-replicators that we found in a search through 3 billion randomly generated
sequences (at three different sequence lengths) to study the evolvability of
generic rather than hand-designed self-replicators. We find that most can
evolve but some are evolutionarily sterile. From this limited data set we are
led to conclude that evolvability is a likely--but not a guaranteed-- property
of random replicators in a digital chemistry.


  Cancer cells utilize large amounts of ATP to sustain growth, relying
primarily on non-oxidative, fermentative pathways for its production. In many
types of cancers this leads, even in the presence of oxygen, to the secretion
of carbon equivalents (usually in the form of lactate) in the cell's
surroundings, a feature known as the Warburg effect. While the molecular basis
of this phenomenon are still to be elucidated, it is clear that the spilling of
energy resources contributes to creating a peculiar microenvironment for
tumors, possibly characterized by a degree of toxicity. This suggests that
mechanisms for recycling the fermentation products (e.g. a lactate shuttle) may
be active, effectively inducing a mutually beneficial metabolic coupling
between aberrant and non-aberrant cells. Here we analyze this scenario through
a large-scale in silico metabolic model of interacting human cells. By going
beyond the cell-autonomous description, we show that elementary
physico-chemical constraints indeed favor the establishment of such a coupling
under very broad conditions. The characterization we obtained by tuning the
aberrant cell's demand for ATP, amino-acids and fatty acids and/or the
imbalance in nutrient partitioning provides quantitative support to the idea
that synergistic multi-cell effects play a central role in cancer sustainment.


  Biological systems, from cells to organisms, must respond to the ever
changing environment in order to survive and function. This is not a simple
task given the often random nature of the signals they receive, as well as the
intrinsically stochastic, many body and often self-organized nature of the
processes that control their sensing and response and limited resources.
Despite a wide range of scales and functions that can be observed in the living
world, some common principles that govern the behavior of biological systems
emerge. Here I review two examples of very different biological problems:
information transmission in gene regulatory networks and diversity of adaptive
immune receptor repertoires that protect us from pathogens. I discuss the
trade-offs that physical laws impose on these systems and show that the optimal
designs of both immune repertoires and gene regulatory networks display similar
discrete tiling structures. These solutions rely on locally non-overlapping
placements of the responding elements (genes and receptors) that, overall,
cover space nearly uniformly.


  Folding of Ubiquitin (Ub) is investigated at low and neutral pH at different
temperatures using simulations of the coarse-grained Self-Organized-Polymer
model with side chains. The calculated radius of gyration, showing dramatic
variations with pH, is in excellent agreement with scattering experiments. At
$T_m$ Ub folds in a two-state manner at low and neutral pH. Clustering analysis
of the conformations sampled in equilibrium folding trajectories at $T_m$, with
multiple transitions between the folded and unfolded states, show a network of
metastable states connecting the native and unfolded states. At low and neutral
pH, Ub folds with high probability through a preferred set of conformations
resulting in a pH-dependent dominant folding pathway. Folding kinetics reveal
that Ub assembly at low pH occurs by multiple pathways involving a combination
of nucleation-collapse and diffusion collision mechanism. The mechanism by
which Ub folds is dictated by the stability of the key secondary structural
elements responsible for establishing long range contacts and collapse of Ub.
Nucleation collapse mechanism holds if the stability of these elements are
marginal, as would be the case at elevated temperatures. If the lifetimes
associated with these structured microdomains are on the order of hundreds of
$\mu sec$ then Ub folding follows the diffusion-collision mechanism with
intermediates many of which coincide with those found in equilibrium. Folding
at neutral pH is a sequential process with a populated intermediate resembling
that sampled at equilibrium. The transition state structures, obtained using a
$P_{fold}$ analysis, are homogeneous and globular with most of the secondary
and tertiary structures being native-like. Many of our findings are not only in
agreement with experiments but also provide missing details not resolvable in
standard experiments.


  In this work, we answer an open problem in the study of phylogenetic
networks. Phylogenetic trees are rooted binary trees in which all edges are
directed away from the root, whereas phylogenetic networks are rooted acyclic
digraphs. For the purpose of evolutionary model validation, biologists often
want to know whether or not a phylogenetic tree is contained in a phylogenetic
network. The tree containment problem is NP-complete even for very restricted
classes of networks such as tree-sibling phylogenetic networks. We prove that
this problem is solvable in cubic time for stable phylogenetic networks. A
linear time algorithm is also presented for the cluster containment problem.


  Stochastic reaction networks are dynamical models of biochemical reaction
systems and form a particular class of continuous-time Markov chains on
$\mathbb{N}^n$. Here we provide a fundamental characterisation that connects
structural properties of a network to its dynamical features. Specifically, we
define the notion of `stochastically complex balanced systems' in terms of the
network's stationary distribution and provide a characterisation of
stochastically complex balanced systems, parallel to that established in the
70-80ies for deterministic reaction networks. Additionally, we establish that a
network is stochastically complex balanced if and only if an associated
deterministic network is complex balanced (in the deterministic sense), thereby
proving a strong link between the theory of stochastic and deterministic
networks. Further, we prove a stochastic version of the `deficiency zero
theorem' and show that any (not only complex balanced) deficiency zero reaction
network has a product-form Poisson-like stationary distribution on all
irreducible components.
  Finally, we provide sufficient conditions for when a product-form
Poisson-like distribution on a single (or all) component(s) implies the network
is complex balanced, and explore the possibility to characterise complex
balanced systems in terms of product-form Poisson-like stationary
distributions.


  We introduce a family of models incorporating random segmental substitutions
and point mutations and demonstrate that such models reproduce algebraic length
distributions of exact matches with the slope $-4$ observed earlier in pairwise
comparisons of DNA of distantly related species. It is demonstrated that
power-law distributions of exact matches emerge when shorter sequences transfer
their DNA content to longer sequences, indicating potential mechanisms of the
increased genome complexity.


  Cellular responses often require the fast activation or repression of
specific genes, which depends on Transcription Factors (TFs) that have to
quickly find the promoters of these genes within a large genome. Transcription
Factors (TFs) search for their DNA promoter target by alternating between bulk
diffusion and sliding along the DNA, a mechanism known as facilitated
diffusion. We study a facilitated diffusion framework with switching between
three search modes: a bulk mode and two sliding modes triggered by
conformational changes between two protein conformations. In one conformation
(search mode) the TF interacts unspecifically with the DNA backbone resulting
in fast sliding. In the other conformation (recognition mode) it interacts
specifically and strongly with DNA base pairs leading to slow displacement.
From the bulk, a TF associates with the DNA at a random position that is
correlated with the previous dissociation point, which implicitly is a function
of the DNA structure. The target affinity depends on the conformation. We
derive exact expressions for the mean first passage time (MFPT) to bind to the
promoter and the conditional probability to bind before detaching when arriving
at the promoter site. We systematically explore the parameter space and compare
various search scenarios. We compare our results with experimental data for the
dimeric Lac repressor search in E.Coli bacteria. We find that a coiled DNA
conformation is absolutely necessary for a fast MFPT. With frequent spontaneous
conformational changes, a fast search time is achieved even when a TF becomes
immobilized in the recognition state due to the specific bindings. We find a
MFPT compatible with experimental data in presence of a specific TF-DNA
interaction energy that has a Gaussian distribution with a large variance.


  A crucial step in the regulation of gene expression is binding of
transcription factor (TF) proteins to regulatory sites along the DNA. But
transcription factors act at nanomolar concentrations, and noise due to random
arrival of these molecules at their binding sites can severely limit the
precision of regulation. Recent work on the optimization of information flow
through regulatory networks indicates that the lower end of the dynamic range
of concentrations is simply inaccessible, overwhelmed by the impact of this
noise. Motivated by the behavior of homeodomain proteins, such as the maternal
morphogen Bicoid in the fruit fly embryo, we suggest a scheme in which
transcription factors also act as indirect translational regulators, binding to
the mRNA of other transcription factors. Intuitively, each mRNA molecule acts
as an independent sensor of the TF concentration, and averaging over these
multiple sensors reduces the noise. We analyze information flow through this
new scheme and identify conditions under which it outperforms direct
transcriptional regulation. Our results suggest that the dual role of
homeodomain proteins is not just a historical accident, but a solution to a
crucial physics problem in the regulation of gene expression.


  Cultural transmission of reproductive success states that successful men have
more children and pass this raised fecundity to their offspring. Balaresque and
colleagues found high frequency haplotypes in a Central Asian Y chromosome
dataset, which they attribute to cultural transmission of reproductive success
by prominent historical men, including Genghis Khan. Using coalescent
simulation, we show that these high frequency haplotypes are consistent with a
neutral model, where they commonly appear simply by chance. Hence, explanations
invoking cultural transmission of reproductive success are statistically
unnecessary.


  The shapes of branching trees have been linked to disease transmission
patterns. In this paper we use the general Crump-Mode-Jagers branching process
to model an outbreak of an infectious disease under mild assumptions.
Introducing a new class of characteristic functions, we are able to derive a
formula for the limit of the frequency of the occurrences of a given shape in a
general tree. The computational challenges concerning the evaluation of this
formula are in part overcome using the Jumping Chronological Contour Process.
We apply the formula to derive the limit of the frequency of cherries,
pitchforks and double cherries in the constant rate birth-death model, and the
frequency of cherries under a non-constant death rate.


  The structure of ecological networks, in particular food webs, determines
their ability to evolve further, i.e. evolvability. The knowledge about how
food web evolvability is determined by the structures of diverse ecological
networks can guide human interventions purposefully to either promote or limit
evolvability of ecosystems. However, the focus of prior food web studies was on
stability and robustness; little is known regarding the impact of ecological
network structures on their evolvability. To correlate ecosystem structure and
evolvability, we adopt the NK model originally from evolutionary biology to
generate and assess the ruggedness of fitness landscapes of a wide spectrum of
model food webs with gradual variation in the amount of feeding loops and link
density. The variation in network structures is controlled by linkage rewiring.
Our results show that more feeding loops and lower trophic link density, i.e.
higher autonomy of species, of food webs increase the potential for the
ecosystem to generate heritable variations with improved fitness. Our findings
allow the prediction of the evolvability of actual food webs according to their
network structures, and provide guidance to enhancing or controlling the
evolvability of specific ecosystems.


  Microarray is one of the essential technologies used by the biologist to
measure genome-wide expression levels of genes in a particular organism under
some particular conditions or stimuli. As microarrays technologies have become
more prevalent, the challenges of analyzing these data for getting better
insight about biological processes have essentially increased. Due to
availability of artificial intelligence based sophisticated computational
techniques, such as artificial neural networks, fuzzy logic, genetic
algorithms, and many other nature-inspired algorithms, it is possible to
analyse microarray gene expression data in more better way. Here, we reviewed
artificial intelligence based techniques for the analysis of microarray gene
expression data. Further, challenges in the field and future work direction
have also been suggested.


  We present a comprehensive study of the interaction of human serum albumin
(HSA) with poly(acrylic acid) (PAA; number average degree of polymerization:
25) in aqueous solution. The interaction of HSA with PAA is studied in dilute
solution as the function of the concentration of added salt (20 - 100 mM) and
temperature (25 - 37$^{\circ}$C). Isothermal titration calorimetry (ITC) is
used to analyze the interaction and to determine the binding constant and
related thermodynamic data. It is found that only one PAA chain is bound per
HSA molecule. The free energy of binding $\Delta G_b$ increases with
temperature significantly. $\Delta G_b$ decreases with increasing salt
concentration and is dominated by entropic contributions due to the release of
bound counterions. Coarse-grained Langevin computer simulations treating the
counterions in an explicit manner are used study the process of binding in
detail. These simulations demonstrate that the PAA chains are bound in the
Sudlow II site of the HSA. Moreover, $\Delta G_b$ is calculated from the
simulations and found to be in very good agreement with the measured data. The
simulations demonstrate clearly that the driving force of binding is the
release of counterions in full agreement with the ITC-data.


  We consider a SIR model with temporary immunity and time dependent
transmission rate. We assume time dependent vaccination which confers the same
immunity as natural infection. We study two types of vaccination strategies: i)
optimal vaccination, in the sense that it minimizes the effort of vaccination
in the set of vaccination strategies for which, for any sufficiently small
perturbation of the disease free state, the number of infectious individuals is
monotonically decreasing; ii) Nash-equilibria strategies where all individuals
simultaneously minimize the joint risk of vaccination versus the risk of the
disease. The former case corresponds to an optimal solution for mandatory
vaccinations, while the second correspond to the equilibrium to be expected if
vaccination is fully voluntary. We are able to show the existence of both an
optimal and Nash strategies in a general setting. In general, these strategies
will not be functions but Radon measures. For specific forms of the
transmission rate, we provide explicit formulas for the optimal and the Nash
vaccination strategies.


  A key factor influencing a drug's efficacy is its residence time in the
binding pocket of the host protein. Using atomistic computer simulation to
predict this residence time and the associated dissociation process is a
desirable but extremely difficult task due to the long timescales involved.
This gets further complicated by the presence of biophysical factors such as
steric and solvation effects. In this work, we perform molecular dynamics (MD)
simulations of the unbinding of a popular prototypical hydrophobic
cavity-ligand system using a metadynamics based approach that allows direct
assessment of kinetic pathways and parameters. When constrained to move in an
axial manner, we find the unbinding time to be on the order of 4000 sec. In
accordance with previous studies, we find that the ligand must pass through a
region of sharp dewetting transition manifested by sudden and high fluctuations
in solvent density in the cavity. When we remove the steric constraints on
ligand, the unbinding happens predominantly by an alternate pathway, where the
unbinding becomes 20 times faster, and the sharp dewetting transition instead
becomes continuous. We validate the unbinding timescales from metadynamics
through a Poisson analysis, and by comparison through detailed balance to
binding timescale estimates from unbiased MD. This work demonstrates that
enhanced sampling can be used to perform explicit solvent molecular dynamics
studies at timescales previously unattainable, obtaining direct and reliable
pictures of the underlying physio-chemical factors including free energies and
rate constants.


  The search for similar genetic sequences is one of the main bioinformatics
tasks. The genetic sequences data banks are growing exponentially and the
searching techniques that use linear time are not capable to do the search in
the required time anymore. Another problem is that the clock speed of the
modern processors are not growing as it did before, instead, the processing
capacity is growing with the addiction of more processing cores and the
techniques which does not use parallel computing does not have benefits from
these extra cores. This work aims to use data indexing techniques to reduce the
searching process computation cost united with the parallelization of the
searching techniques to use the computational capacity of the multi core
processors. To verify the viability of using these two techniques
simultaneously, a software which uses parallelization techniques with inverted
indexes was developed.
  Experiments were executed to analyze the performance gain when parallelism is
utilized, the search time gain, and also the quality of the results when it
compared with others searching tools. The results of these experiments were
promising, the parallelism gain overcame the expected speedup, the searching
time was 20 times faster than the parallelized NCBI BLAST, and the searching
results showed a good quality when compared with this tool.
  The software source code is available at
https://github.com/felipealbrecht/Genoogle .


  Humans have been using symbolic representation (i.e. art) as a creative
cultural form indisputably for at least 80,000 years. A description of the
processes central to the evolution of art from sculpted earthen forms early in
human existence to paintings in museums of the modern world is absent in
scholarly literature. The present manuscript offers a memetic theory of art
evolution demonstrating typological transitions of art in the period post the
3-dimensional to 2-dimensional transition occurring in the early Aurignacian
period. The process of art evolution in the 2-dimensional form central to the
typological transitions post this period was propelled by material and/or
structure chosen to display the artistic message. Symbolic representational
tinkering whether collectively or individually creates artistic typologies of
display deemed aesthetic and non-aesthetic by the current cultural milieu.
Contemporary forms of graffiti underscore the aesthetic/nonaesthetic dichotomy.


  We study adaptation of a haploid asexual population on a fitness landscape
defined over binary genotype sequences of length $L$. We consider greedy
adaptive walks in which the population moves to the fittest among all single
mutant neighbors of the current genotype until a local fitness maximum is
reached. The landscape is of the rough mount Fuji type, which means that the
fitness value assigned to a sequence is the sum of a random and a deterministic
component. The random components are independent and identically distributed
random variables, and the deterministic component varies linearly with the
distance to a reference sequence. The deterministic fitness gradient $c$ is a
parameter that interpolates between the limits of an uncorrelated random
landscape ($c = 0$) and an effectively additive landscape ($c \to \infty$).
When the random fitness component is chosen from the Gumbel distribution,
explicit expressions for the distribution of the number of steps taken by the
greedy walk are obtained, and it is shown that the walk length varies
non-monotonically with the strength of the fitness gradient when the starting
point is sufficiently close to the reference sequence. Asymptotic results for
general distributions of the random fitness component are obtained using
extreme value theory, and it is found that the walk length attains a
non-trivial limit for $L \to \infty$, different from its values for $c=0$ and
$c = \infty$, if $c$ is scaled with $L$ in an appropriate combination.


  BioNetGen is an open-source software package for rule-based modeling of
complex biochemical systems. Version 2.2 of the software introduces numerous
new features for both model specification and simulation. Here, we report on
these additions, discussing how they facilitate the construction, simulation,
and analysis of larger and more complex models than previously possible.


  By fixing all parameters in a phylogenetic likelihood model except for one
branch length, one obtains a one-dimensional likelihood function. In this work,
we introduce a mathematical framework to characterize the shapes of such
one-dimensional phylogenetic likelihood functions. This framework is based on
analyses of algebraic structures on the space of all frequency patterns with
respect to a polynomial representation of the likelihood functions. Using this
framework, we provide conditions under which the one-dimensional phylogenetic
likelihood functions are guaranteed to have at most one stationary point, and
this point is the maximum likelihood branch length. These conditions are
satisfied by common simple models including all binary models, the Jukes-Cantor
model and the Felsenstein 1981 model.
  We then prove that for the simplest model that does not satisfy our
conditions, namely, the Kimura 2-parameter model, the one-dimensional
likelihood functions may have multiple stationary points. As a proof of
concept, we construct a non-degenerate example in which the phylogenetic
likelihood function has two local maxima and a local minimum. To construct such
examples, we derive a general method of constructing a tree and sequence data
with a specified frequency pattern at the root. We then extend the result to
prove that the space of all rescaled and translated one-dimensional
phylogenetic likelihood functions under the Kimura 2-parameter model is dense
in the space of all non-negative continuous functions on $[0, \infty)$ with
finite limits. These results indicate that one-dimensional likelihood functions
under advanced evolutionary models can be more complex than it is typically
assumed by phylogenetic inference algorithms; however, these complexities can
be effectively captured by the Kimura 2-parameter model.


  Alchemical Grid Dock (AlGDock) is open-source software designed to compute
the binding potential of mean force (BPMF) - the binding free energy between a
flexible ligand and a rigid receptor - for a small organic ligand and a
biological macromolecule. Multiple BPMFs can be used to rigorously compute
binding affinities between flexible partners. AlGDock uses replica exchange
between thermodynamic states at different temperatures and receptor-ligand
interaction strengths. Receptor-ligand interaction energies are represented by
interpolating precomputed grids. Thermodynamic states are adaptively
initialized and adjusted on-the-fly to maintain replica exchange rates. In
demonstrative calculations, when the bound ligand is treated as fully solvated,
AlGDock estimates BPMFs with a precision within 4 kT in 65% and within 8 kT for
91% of systems. It correctly identifies the native binding pose in 83% of
simulations. Performance is sometimes limited by subtle differences in the
important configuration space of sampled and targeted thermodynamic states.


  We compare the informational architecture of biological and random networks
to identify informational features that may distinguish biological networks
from random. The study presented here focuses on the Boolean network model for
regulation of the cell cycle of the fission yeast Schizosaccharomyces Pombe. We
compare calculated values of local and global information measures for the
fission yeast cell cycle to the same measures as applied to two different
classes of random networks: random and scale-free. We report patterns in local
information processing and storage that do indeed distinguish biological from
random, associated with control nodes that regulate the function of the fission
yeast cell cycle network. Conversely, we find that integrated information,
which serves as a global measure of "emergent" information processing, does not
differ from random for the case presented. We discuss implications for our
understanding of the informational architecture of the fission yeast cell cycle
network in particular, and more generally for illuminating any distinctive
physics that may be operative in life.


  An approximation to the ~4 Mbp basic genome shared by 32 strains of E. coli
representing six evolutionary groups has been derived and analyzed
computationally. A multiple-alignment of the 32 complete genome sequences was
filtered to remove mobile elements and identify the most reliable ~90% of the
aligned length of each of the resulting 496 basic-genome pairs. Patterns of
single-bp mutations (SNPs) in aligned pairs distinguish clonally inherited
regions from regions where either genome has acquired DNA fragments from
diverged genomes by homologous recombination since their last common ancestor.
Such recombinant transfer is pervasive across the basic genome, mostly between
genomes in the same evolutionary group, and generates many unique mosaic
patterns. The six least-diverged genome-pairs have one or two recombinant
transfers of length ~40-115 kbp (and few if any other transfers), each
containing one or more gene clusters known to confer strong selective advantage
in some environments. Moderately diverged genome-pairs (0.4-1% SNPs ) show
mosaic patterns of interspersed clonal and recombinant regions of varying
lengths throughout the basic genome, whereas more highly diverged pairs within
an evolutionary group or pairs between evolutionary groups having >1.3% SNPs
have few clonal matches longer than a few kbp. Many recombinant transfers
appear to incorporate fragments of the entering DNA produced by restriction
systems of the recipient cell. A simple computational model can closely fit the
data. Most recombinant transfers seem likely to be due to generalized
transduction by co-evolving populations of phages, which could efficiently
distribute variability throughout bacterial genomes.


  The first outbreak of dengue occurred in Madeira Island on 2012, featuring
one virus serotype. Aedes aegypti was the vector of the disease and it is
unlikely that it will be eliminated from the island. Therefore, a new outbreak
of dengue fever can occur and, if it happens, risk to the population increases
if two serotypes coexist. In this paper, mathematical modeling and numerical
simulations are carried out to forecast what may happen in Madeira Island in
such scenario.


  Recent work has attempted to use whole-genome sequence data from pathogens to
reconstruct the transmission trees linking infectors and infectees in
outbreaks. However, transmission trees from one outbreak do not generalize to
future outbreaks. Reconstruction of transmission trees is most useful to public
health if it leads to generalizable scientific insights about disease
transmission. In a survival analysis framework, estimation of transmission
parameters is based on sums or averages over the possible transmission trees. A
phylogeny can increase the precision of these estimates by providing partial
information about who infected whom. The leaves of the phylogeny represent
sampled pathogens, which have known hosts. The interior nodes represent common
ancestors of sampled pathogens, which have unknown hosts. Starting from
assumptions about disease biology and epidemiologic study design, we prove that
there is a one-to-one correspondence between the possible assignments of
interior node hosts and the transmission trees simultaneously consistent with
the phylogeny and the epidemiologic data on person, place, and time. We develop
algorithms to enumerate these transmission trees and show these can be used to
calculate likelihoods that incorporate both epidemiologic data and a phylogeny.
A simulation study confirms that this leads to more efficient estimates of
hazard ratios for infectiousness and baseline hazards of infectious contact,
and we use these methods to analyze data from a foot-and-mouth disease virus
outbreak in the United Kingdom in 2001. These results demonstrate the
importance of data on individuals who escape infection, which is often
overlooked. The combination of survival analysis and algorithms linking
phylogenies to transmission trees is a rigorous but flexible statistical
foundation for molecular infectious disease epidemiology.


  Biological aggregations such as insect swarms and bird flocks may arise from
a combination of social interactions and environmental cues. We focus on
nonlocal continuum equations, which are often used to model aggregations, and
yet which pose significant analytical and computational challenges. Beginning
with a particular nonlocal aggregation model [Topaz et al., Bull. Math. Bio.,
2006], we derive the minimal well-posed long-wave approximation, which is a
degenerate Cahn-Hilliard equation. Energy minimizers of this reduced, local
model retain many salient features of those of the nonlocal model, especially
for large populations and away from an aggregation's boundaries. Using the
Cahn-Hilliard model as a testbed, we investigate the degree to which an
external potential modeling food sources can be used to suppress peak
population density, which is essential for controlling locust outbreaks. A
random distribution of food sources tends to increase peak density above its
intrinsic value, while a periodic pattern of food sources can decrease it.


  Researchers working with mathematical models are often confronted by the
related problems of parameter estimation, model validation, and model
selection. These are all optimization problems, well-known to be challenging
due to non-linearity, non-convexity and multiple local optima. Furthermore, the
challenges are compounded when only partial data is available. Here, we
consider polynomial models (e.g., mass-action chemical reaction networks at
steady state) and describe a framework for their analysis based on optimization
using numerical algebraic geometry. Specifically, we use probability-one
polynomial homotopy continuation methods to compute all critical points of the
objective function, then filter to recover the global optima. Our approach
exploits the geometric structures relating models and data, and we demonstrate
its utility on examples from cell signaling, synthetic biology, and
epidemiology.


  Most theories of evolutionary diversification are based on equilibrium
assumptions: they are either based on optimality arguments involving static
fitness landscapes, or they assume that populations first evolve to an
equilibrium state before diversification occurs, as exemplified by the concept
of evolutionary branching points in adaptive dynamics theory. Recent results
indicate that adaptive dynamics may often not converge to equilibrium points
and instead generate complicated trajectories if evolution takes place in
high-dimensional phenotype spaces. Even though some analytical results on
diversification in complex phenotype spaces are available, to study this
problem in general we need to reconstruct individual-based models from the
adaptive dynamics generating the non-equilibrium dynamics. Here we first
provide a method to construct individual-based models such that they faithfully
reproduce the given adaptive dynamics attractor without diversification. We
then show that a propensity to diversify can by introduced by adding Gaussian
competition terms that generate frequency dependence while still preserving the
same adaptive dynamics. For sufficiently strong competition, the disruptive
selection generated by frequency-dependence overcomes the directional evolution
along the selection gradient and leads to diversification in phenotypic
directions that are orthogonal to the selection gradient.


  Classic experiments on the distribution of ducks around separated food
sources found consistency with the `ideal free' distribution in which the local
population is proportional to the local supply rate. Motivated by this
experiment and others, we examine the analogous problem in the microbial world:
the distribution of chemotactic bacteria around multiple nearby food sources.
In contrast to the optimization of uptake rate that may hold at the level of a
single cell in a spatially varying nutrient field, nutrient consumption by a
population of chemotactic cells will modify the nutrient field, and the uptake
rate will generally vary throughout the population. Through a simple model we
study the distribution of resource uptake in the presence of chemotaxis,
consumption, and diffusion of both bacteria and nutrients. Borrowing from the
field of theoretical economics, we explore how the Gini index can be used as a
means to quantify the inequalities of uptake. The redistributive effect of
chemotaxis can lead to a phenomenon we term `chemotactic levelling', and the
influence of these results on population fitness are briefly considered.


  Many discrete mathematics problems in phylogenetics are defined in terms of
the relative labeling of pairs of leaf-labeled trees. These relative labelings
are naturally formalized as tanglegrams, which have previously been an object
of study in coevolutionary analysis. Although there has been considerable work
on planar drawings of tanglegrams, they have not been fully explored as
combinatorial objects until recently. In this paper, we describe how many
discrete mathematical questions on trees "factor" through a problem on
tanglegrams, and how understanding that factoring can simplify analysis.
Depending on the problem, it may be useful to consider a unordered version of
tanglegrams, and/or their unrooted counterparts. For all of these definitions,
we show how the isomorphism types of tanglegrams can be understood in terms of
double cosets of the symmetric group, and we investigate their automorphisms.
Understanding tanglegrams better will isolate the distinct problems on
leaf-labeled pairs of trees and reveal natural symmetries of spaces associated
with such problems.


  Antibiotic-resistant bacterial infection becomes one of the most serious
risks to public health care today. However, discouragingly, the development of
new antibiotics has been little progressed over the last decade. There is an
urgent need of the alternative approaches to treat the antibiotic-resistant
bacteria. The novel methods, which include photothermal therapy based on gold
nano-materials and ionizing radiation such as X-rays and gamma rays, have been
reported. Studies of the effects of high-energy proton radiation on bacteria
are mainly focused on Bacillus species and its spores. The effect of proton
beams on Escherichia coli (E. coli) has been limitedly reported. The
Escherichia coli is an important biological tool to obtain the metabolic and
genetic information and also a common model microorganism for studying toxicity
and antimicrobial activity. In addition, E. coli is a common bacterium in the
intestinal tract of mammals. Herein, the morphological and physiological
changes of E. coli after proton irradiation were investigated. The diluted
solutions of the cells were used for proton beam radiation. LB agar plates were
used to count the number of colonies formed. The growing profile of the cells
was monitored by optical density at 600 nm. The morphology of the irradiated
cells was analyzed with optical microscope. Microarray analysis was performed
to examine the gene expression changes between irradiated samples and control
samples without irradiation.


  Network analysis has revealed whole-network properties hypothesized to be
general characteristics of ecosystems including pathway proliferation, and
network non-locality, homogenization, amplification, mutualism, and synergism.
Collectively these properties characterize the impact of indirect interactions
among ecosystem elements. While ecosystem networks generally trace a
thermodynamically conserved unit through the system, there appear to be several
model classes. For example, trophic (TRO) networks are built upon a food web,
usually follow energy or carbon, and are the most abundant in the literature.
Biogeochemical cycling (BGC) networks trace nutrients like nitrogen or
phosphorus and tend to have more recycling than TRO. We tested (1) the
hypothesized generality of the properties in BGC networks and (2) that they
tend to be more strongly expressed in BGC networks than in the TRO networks due
to increased recycling. We compared the properties in 22 BGC and 57 TRO
ecosystem networks from the literature using enaR. The results generally
support the hypotheses. First, five of the properties occurred in all 22 BGC
models, while network mutualism occurred in 86% of the models. Further, these
results were generally robust to a $\pm$50% uncertainty in the model
parameters. Second, the mean network statistics for the six properties were
statistically significantly greater in the BGC models than the TRO models.
These results (1) confirm the general presence of these properties in ecosystem
networks, (2) highlight the significance of model types in determining property
intensities, (3) reinforce the importance of recycling, and (4) provide a set
of indicator benchmarks for future systems comparisons. Further, this work
highlights how indirect effects distributed by network connectivity can
transform whole-ecosystem functioning, and adds to the growing domain of
network ecology.


  Phylogenetic trees are a central tool in understanding evolution. They are
typically inferred from sequence data, and capture evolutionary relationships
through time. It is essential to be able to compare trees from different data
sources (e.g. several genes from the same organisms) and different inference
methods. We propose a new metric for robust, quantitative comparison of rooted,
labeled trees. It enables clear visualizations of tree space, gives meaningful
comparisons between trees, and can detect distinct islands of tree topologies
in posterior distributions of trees. This makes it possible to select
well-supported summary trees. We demonstrate our approach on Dengue fever
phylogenies.


  The role of species interactions in controlling the interplay between the
stability of an ecosystem and its biodiversity is still not well understood.
The ability of ecological communities to recover after a small perturbation of
the species abundances (local asymptotic stability) has been well studied,
whereas the likelihood of a community to persist when the interactions are
altered (structural stability) has received much less attention. Our goal is to
understand the effects of diversity, interaction strenghts and ecological
network structure on the volume of parameter space leading to feasible
equilibria, i.e., ones in which all populations have positive abundances. We
develop a geometrical framework to study the range of conditions necessary for
feasible coexistence in both mutualistic and consumer-resource systems. Using
analytical and numerical methods, we show that feasibility is determined by
just a handful of quantities describing the interactions, yielding a nontrivial
complexity-feasibility relationship. Analyzing more than 100 empirical
networks, we show that the range of coexistence conditions in mutualistic
systems can be analytically predicted by means of a null model of random
interactions, whereas food webs are characterized by smaller coexistence
domains than those expected by chance. Finally, we characterize the geometric
shape of the feasibility domain, thereby identifying the direction of
perturbations that are more likely to cause extinctions. Interestingly, the
structure of mutualistic interactions leads to very heterogeneous responses to
perturbations, making those systems more fragile than expected by chance.


  Genome-wide proximity placement analysis of 10,598 HSGRL within the context
of the principal regulatory structures of the interphase chromatin, namely
topologically-associating domains (TADs) and specific sub-TAD structures termed
super-enhancer domains (SEDs) revealed that 0.8%-10.3% of TADs contain more
than half of HSGRL. Of the 3,127 TADs in the hESC genome, 24 (0.8%); 53 (1.7%);
259 (8.3%); and 322 (10.3%) harbor 1,110 (52.4%); 1,936 (50.9%); 1,151 (59.6%);
and 1,601 (58.3%) HSGRL sequences from four distinct families, respectively.
TADs that are enriched for HSGRL and termed rapidly-evolving in humans TADs
(revTADs) manifest distinct correlation patterns between HSGRL placements and
recombination rates. There are significant enrichment within revTAD boundaries
of hESC-enhancers, primate-specific CTCF-binding sites, human-specific
RNAPII-binding sites, hCONDELs, and H3K4me3 peaks with human-specific
enrichment at TSS in prefrontal cortex neurons (p < 0.0001 in all instances).
In hESC genome, 331 of 504 (66%) of SE-harboring TADs contain HSGRL and 68% of
SEs co-localize with HSGRL, suggesting that HSGRL rewired SE-driven GRNs within
revTADs by inserting novel and/or erasing existing regulatory sequences.
Consequently, markedly distinct features of chromatin structures evolved in
hESC compared to mouse: the SE quantity is 3-fold higher and the median SE size
is significantly larger; concomitantly, the TAD number is increased by 42%
while the median TAD size is decreased (p=9.11E-37). Present analyses revealed
a global role for HSGRL in increasing both quantity and size of SEs and
increasing the number and size reduction of TADs, which may facilitate a
convergence of TAD and SED architectures of interphase chromatin and define a
trend of increasing regulatory complexity during evolution of GRNs.


  Eukaryotes tune the transcriptional activity of their genome by altering the
nucleosome core particle through multiple chemical processes. In particular,
replacement of the canonical H2A histone with the variants macroH2A and H2A.Z
has been shown to affect DNA accessibility and nucleosome stability; however,
the processes by which this occurs remain poorly understood. Here, we elucidate
the molecular mechanisms of these variants with an extensive molecular dynamics
study of the canonical nucleosome along with three variant-containing
structures: H2A.Z, macroH2A, and an H2A mutant with macroH2A-like L1 loops.
Simulation results show that variant L1 loops play a pivotal role in
stabilizing DNA binding to the octamer through direct interactions, core
structural rearrangements, and altered allosteric networks in the nucleosome.
All variants influence dynamics; however, macroH2A-like systems have the
largest effect on energetics. In addition, we provide a comprehensive analysis
of allosteric networks in the nucleosome and demonstrate that variants take
advantage of stronger interactions between L1 loops to propagate dynamics
throughout the complex. Furthermore, we show that post-translational
modifications are enriched at key locations in these networks. Taken together,
these results provide new insights into the relationship between the structure,
dynamics, and function of the nucleosome core particle and chromatin fibers,
and how they are influenced by chromatin remodelling factors.


  We propose a dynamical model in which a network structure evolves in a
self-organized critical (SOC) manner and explain a possible origin of the
emergence of fractal and small-world networks. Our model combines a network
growth and its decay by failures of nodes. The decay mechanism reflects the
instability of large functional networks against cascading overload failures.
It is demonstrated that the dynamical system surely exhibits SOC
characteristics, such as power-law forms of the avalanche size distribution,
the cluster size distribution, and the distribution of the time interval
between intermittent avalanches. During the network evolution, fractal networks
are spontaneously generated when networks experience critical cascades of
failures that lead to a percolation transition. In contrast, networks far from
criticality have small-world structures. We also observe the crossover behavior
from fractal to small-world structure in the network evolution.


  With the wealth of high-throughput sequencing data generated by recent
large-scale consortia, predictive gene expression modelling has become an
important tool for integrative analysis of transcriptomic and epigenetic data.
However, sequencing data-sets are characteristically large, and previously
modelling frameworks are typically inefficient and unable to leverage
multi-core or distributed processing architectures. In this study, we detail an
efficient and parallelised MapReduce implementation of gene expression
modelling. We leverage the computational efficiency of this framework to
provide an integrative analysis of over fifty histone modification data-sets
across a variety of cancerous and non-cancerous cell-lines. Our results
demonstrate that the genome-wide relationships between histone modifications
and mRNA transcription are lineage, tissue and karyotype-invariant, and that
models trained on matched epigenetic/transcriptomic data from non-cancerous
cell-lines are able to predict cancerous expression with equivalent genome-wide
fidelity.


  The passage of DNA through a nanopore can be effectively decomposed into two
distinct phases, docking and actual translocation. In experiments each phase is
characterized by a distinct current signature which allows the discrimination
of the two events. However, at low voltages a clear distinction of the two
phases is lost. By using numerical simulations we clarify how the current
signature associated to the docking events depends on the applied voltage. The
simulations show that at small voltage the DNA globule enhances the pore
conductance due to an enrichment of charge carriers. At high voltage, the
globule drains substantial charge carriers from the pore region, thereby
reducing the overall conductance. The results provide a new interpretation to
the experimental data on conductance and show how docking interferes with the
translocation signal, of potential interest for sequencing applications.


  Large-scale force generation is essential for biological functions such as
cell motility, embryonic development, and muscle contraction. In these
processes, forces generated at the molecular level by motor proteins are
transmitted by disordered fiber networks, resulting in large-scale active
stresses. While these fiber networks are well characterized macroscopically,
this stress generation by microscopic active units is not well understood. Here
we theoretically study force transmission in these networks, and find that
local active forces are rectified towards isotropic contraction and strongly
amplified as fibers collectively buckle in the vicinity of the active units.
This stress amplification is reinforced by the networks' disordered nature, but
saturates for high densities of active units. Our predictions are
quantitatively consistent with experiments on reconstituted tissues and
actomyosin networks, and shed light on the role of the network microstructure
in shaping active stresses in cells and tissue.


  The process of protein search for specific binding sites on DNA is
fundamentally important since it marks the beginning of all major biological
processes. We present a theoretical investigation that probes the role of DNA
sequence symmetry, heterogeneity and chemical composition in the protein search
dynamics. Using a discrete-state stochastic approach with a first-passage
events analysis, which takes into account the most relevant physical-chemical
processes, a full analytical description of the search dynamics is obtained. It
is found that, contrary to existing views, the protein search is generally
faster on DNA with more heterogeneous sequences. In addition, the search
dynamics might be affected by the chemical composition near the target site.
The physical origins of these phenomena are discussed. Our results suggest that
biological processes might be effectively regulated by modifying chemical
composition, symmetry and heterogeneity of a genome.


  We give a review of the current state of the art of data-driven modelling of
the inositol trisphosphate receptor (IPR). After explaining that the IPR plays
a crucial role as a central regulator in calcium dynamics, several sources of
relevant experimental data are introduced. Single ion channels are best studied
by recording single-channel currents under different ligand concentrations via
the patch-clamp technique. The particular relevance of modal gating, the
spontaneous switching between different levels of channel activity that occur
even at constant ligand concentrations, is highlighted. In order to investigate
the interactions of IPRs, calcium release from small clusters of channels,
so-called calcium puffs, can be used. We then present the mathematical
framework common to all models based on single-channel data, aggregated
continuous-time Markov models, and give a short review of statistical
approaches for parameterising these models with experimental data. The process
of building a Markov model that integrates various sources of experimental data
is illustrated using two recent examples, the model by Ullah et al. and the
"Park-Drive" model by Siekmann et al., the only models that account for all
sources of data currently available. Finally, it is demonstrated that the
essential features of the Park-Drive model in different models of calcium
dynamics are preserved after reducing it to a two-state model that only
accounts for the switching between the inactive "park" and the active "drive"
mode. This highlights the fact that modal gating is the most important
mechanism of ligand regulation in the IPR. It also emphasises that data-driven
models of ion channels do not necessarily have to lead to detailed models but
can be constructed so that relevant data is selected to represent ion channels
at the appropriate level of complexity for a given application.


  It is generally recognized that a distinguishing feature of life is its
peculiar capability to avoid equilibration. The origin of this capability and
its evolution along the timeline of abiogenesis is not yet understood. We
propose to study an analog of this phenomenon that could emerge in
non-biological systems. To this end, we introduce the concept of sustainability
of transient kinetic regimes. This concept is illustrated via investigation of
cooperative effects in an extended system of compartmentalized chemical
oscillators under batch and semi-batch conditions. The computational study of a
model system shows robust enhancement of lifetimes of the decaying oscillations
which translates into the evolution of the survival function of the transient
non-equilibrium regime. This model does not rely on any form of replication.
Rather, it explores the role of a structured effective environment as a
contributor to the system-bath interactions that define non-equilibrium
regimes. We implicate the noise produced by the effective environment of a
compartmentalized oscillator as the cause of the lifetime extension.


  Cells are strongly out-of-equilibrium systems driven by continuous energy
supply. They carry out many vital functions requiring active transport of
various ingredients and organelles, some being small, others being large. The
cytoskeleton, composed of three types of filaments, determines the shape of the
cell and plays a role in cell motion. It also serves as a road network for the
so-called cytoskeletal motors. These molecules can attach to a cytoskeletal
filament, perform directed motion, possibly carrying along some cargo, and then
detach. It is a central issue to understand how intracellular transport driven
by molecular motors is regulated, in particular because its breakdown is one of
the signatures of some neuronal diseases like the Alzheimer.
  We give a survey of the current knowledge on microtubule based intracellular
transport. We first review some biological facts obtained from experiments, and
present some modeling attempts based on cellular automata. We start with
background knowledge on the original and variants of the TASEP (Totally
Asymmetric Simple Exclusion Process), before turning to more application
oriented models. After addressing microtubule based transport in general, with
a focus on in vitro experiments, and on cooperative effects in the
transportation of large cargos by multiple motors, we concentrate on axonal
transport, because of its relevance for neuronal diseases. It is a challenge to
understand how this transport is organized, given that it takes place in a
confined environment and that several types of motors moving in opposite
directions are involved. We review several features that could contribute to
the efficiency of this transport, including the role of motor-motor
interactions and of the dynamics of the underlying microtubule network.
Finally, we discuss some still open questions.


  Bond rupture mechanism enables to explain in block asymmetry of elaxation,
force-velocity curve and the path of energy dissipation in muscle


  Several recent experiments suggest that sharply bent DNA has a surprisingly
high bending flexibility, but the cause of this flexibility is poorly
understood. Although excitation of flexible defects can explain these results,
whether such excitation can occur with the level of DNA bending in these
experiments remains unclear. Intriguingly, the DNA contained preexisting nicks
in most of these experiments but whether nicks might play a role in flexibility
has never been considered in the interpretation of experimental results. Here,
using full-atom molecular dynamics simulations, we show that nicks promote DNA
basepair disruption at the nicked sites, which drastically reduces DNA bending
energy. In addition, lower temperatures suppress the nick-dependent basepair
disruption. In the absence of nicks, basepair disruption can also occur but
requires a higher level of DNA bending. Therefore, basepair disruption inside
B-form DNA can be suppressed if the DNA contains preexisting nicks. Overall,
our results suggest that the reported mechanical anomaly of sharply bent DNA is
likely dependent on preexisting nicks, therefore the intrinsic mechanisms of
sharply bent nick-free DNA remain an open question.


  It has been shown that genome spatial structures largely affect both genome
activity and DNA function. Knowing this, many researchers are currently
attempting to accurately model genome structures. Despite these increased
efforts there still exists a shortage of tools dedicated to visualizing the
genome. Creating a tool that can accurately visualize the genome can aid
researchers by highlighting structural relationships that may not be obvious
when examining the sequence information alone. Here we present a desktop
application, known as GMOL, designed to effectively visualize genome tertiary
structures at multiple scales so that researchers may better analyze their
genomic data. GMOL was developed based upon our multi-scale approach that
allows a user to zoom in and out between six separate levels within the genome.
These six scales are full genome, chromosome, loci, fiber, nucleosome, and
nucleotide. In order to store the data of the different scales, a new file
format, known as GSS, was created. With GMOL, a user can choose any unit at any
scale and scale it up or down to visualize its structure and retrieve
corresponding genome sequences from either Ensembl or a local database. Users
can also interactively manipulate and measure the whole genome structure and
extract static images and machine-readable data files in PDB format from the
multi-scale structure. By using GMOL researchers will be able to better
understand and analyze genome structure models and the impact their structural
relations have on genome activity and DNA function through GMOLs unique
features and functions, which includes the multi-scale method that can satisfy
the users requirement to not only visualize genome tertiary structure, but also
measure it.


  Detecting genomic footprints of selection is an important step in the
understanding of evolution. Accounting for linkage disequilibrium in genome
scans allows increasing the detection power, but haplotype-based methods
require individual genotypes and are not applicable on pool-sequenced samples.
We propose to take advantage of the local score approach to account for linkage
disequilibrium, accumulating (possibly small) signals from single markers over
a genomic segment, to clearly pinpoint a selection signal, avoiding windowing
methods. This method provided results similar to haplotype-based methods on two
benchmark data sets with individual genotypes. Results obtained for a divergent
selection experiment on behavior in quail, where two lines were sequenced in
pools, are precise and biologically coherent, while competing methods failed:
our approach led to the detection of signals involving genes known to act on
social responsiveness or autistic traits. This local score approach is general
and can be applied to other genome-wide analyzes such as GWAS or genome scans
for selection.


  In this paper, we investigate the effect of "heterogeneity of link weight",
heterogeneity of the frequency or amount of interactions among individuals, on
the evolution of cooperation. Based on an analysis of the evolutionary
prisoner's dilemma game on a weighted one-dimensional lattice network with
"intra-individual heterogeneity", we confirm that moderate level of link-weight
heterogeneity can facilitate cooperation. Furthermore, we identify two key
mechanisms by which link-weight heterogeneity promotes the evolution of
cooperation: mechanisms for spread and maintenance of cooperation. We also
derive the corresponding conditions under which the mechanisms can work through
evolutionary dynamics.


  Mutations can arise from the chance misincorporation of nucleotides during
DNA replication or from DNA lesions that are not repaired correctly. We
introduce a model that relates the source of mutations to their accumulation
with cell divisions, providing a framework for understanding how mutation rates
depend on sex, age and absolute time. We show that the accrual of mutations
should track cell divisions not only when mutations are replicative in origin
but also when they are non-replicative and repaired efficiently. One
implication is that the higher incidence of cancer in rapidly renewing tissues,
an observation ascribed to replication errors, could instead reflect exogenous
or endogenous mutagens. We further find that only mutations that arise from
inefficiently repaired lesions will accrue according to absolute time; thus, in
the absence of selection on mutation rates, the phylogenetic "molecular clock"
should not be expected to run steadily across species.


  Levy flights in the space of mutations model time evolution of bacterial DNA.
Parameters in the model are adjusted in order to fit observations coming from
the Long Time Evolution Experiment with E. Coli.


  In this paper, we analyze the dynamics of spreading processes taking place
over time-varying networks. A common approach to model time-varying networks is
via Markovian random graph processes. This modeling approach presents the
following limitation: Markovian random graphs can only replicate switching
patterns with exponential inter-switching times, while in real applications
these times are usually far from exponential. In this paper, we introduce a
flexible and tractable extended family of processes able to replicate, with
arbitrary accuracy, any distribution of inter-switching times. We then study
the stability of spreading processes in this extended family. We first show
that a direct analysis based on It\^o's formula provides stability conditions
in terms of the eigenvalues of a matrix whose size grows exponentially with the
number of edges. To overcome this limitation, we derive alternative stability
conditions involving the eigenvalues of a matrix whose size grows linearly with
the number of nodes. Based on our results, we also show that heuristics based
on aggregated static networks approximate the epidemic threshold more
accurately as the number of nodes grows, or the temporal volatility of the
random graph process is reduced. Finally, we illustrate our findings via
numerical simulations.


  There are several applications in computational biophysics which require the
optimization of discrete interacting states; e.g., amino acid titration states,
ligand oxidation states, or discrete rotamer angles. Such optimization can be
very time-consuming as it scales exponentially in the number of sites to be
optimized. In this paper, we describe a new polynomial-time algorithm for
optimization of discrete states in macromolecular systems. This algorithm was
adapted from image processing and uses techniques from discrete mathematics and
graph theory to restate the optimization problem in terms of "maximum
flow-minimum cut" graph analysis. The interaction energy graph, a graph in
which vertices (amino acids) and edges (interactions) are weighted with their
respective energies, is transformed into a flow network in which the value of
the minimum cut in the network equals the minimum free energy of the protein,
and the cut itself encodes the state that achieves the minimum free energy.
Because of its deterministic nature and polynomial-time performance, this
algorithm has the potential to allow for the ionization state of larger
proteins to be discovered.


  Bacteria are able to maintain a narrow distribution of cell sizes by
regulating the timing of cell divisions. In rich nutrient conditions, cells
divide much faster than their chromosomes replicate. This implies that cells
maintain multiple rounds of chromosome replication per cell division by
regulating the timing of chromosome replications. Here, we show that both cell
size and chromosome replication may be simultaneously regulated by the
long-standing initiator accumulation strategy. The strategy proposes that
initiators are produced in proportion to the volume increase and is accumulated
at each origin of replication, and chromosome replication is initiated when a
critical amount per origin has accumulated. We show that this model maps to the
incremental model of size control, which was previously shown to reproduce
experimentally observed correlations between various events in the cell cycle
and explains the exponential dependence of cell size on the growth rate of the
cell. Furthermore, we show that this model also leads to the efficient
regulation of the timing of initiation and the number of origins consistent
with existing experimental results.


  Ecological systems are dynamic and policies to manage them need to respond to
that variation. However, policy adjustments will sometimes be costly, which
means that fine-tuning a policy to track variability in the environment very
tightly will only sometimes be worthwhile. We use a classic fisheries
management question -- how to manage a stochastically varying population using
annually varying quotas in order to maximize profit -- to examine how costs of
policy adjustment change optimal management recommendations. Costs of policy
adjustment (here changes in fishing quotas through time) could take different
forms. For example, these costs may respond to the size of the change being
implemented, or there could be a fixed cost any time a quota change is made. We
show how different forms of policy costs have contrasting implications for
optimal policies. Though it is frequently assumed that costs to adjusting
policies will dampen variation in the policy, we show that certain cost
structures can actually increase variation through time. We further show that
failing to account for adjustment costs has a consistently worse economic
impact than would assuming these costs are present when they are not.


  Gene regulatory network reconstruction is a fundamental problem in
computational biology. We recently developed an algorithm, called PANDA
(Passing Attributes Between Networks for Data Assimilation), that integrates
multiple sources of 'omics data and estimates regulatory network models. This
approach was initially implemented in the C++ programming language and has
since been applied to a number of biological systems. In our current research
we are beginning to expand the algorithm to incorporate larger and most diverse
data-sets, to reconstruct networks that contain increasing numbers of elements,
and to build not only single network models, but sets of networks. In order to
accomplish these "Big Data" applications, it has become critical that we
increase the computational efficiency of the PANDA implementation. In this
paper we show how to recast PANDA's similarity equations as matrix operations.
This allows us to implement a highly readable version of the algorithm using
the MATLAB/Octave programming language. We find that the resulting M-code much
shorter (103 compared to 1128 lines) and more easily modifiable for potential
future applications. The new implementation also runs significantly faster,
with increasing efficiency as the network models increase in size. Tests
comparing the C-code and M-code versions of PANDA demonstrate that this
speed-up is on the order of 20-80 times faster for networks of similar
dimensions to those we find in current biological applications.


  In this paper I address the question - how large is a phylogenetic sample I
propose a definition of a phylogenetic effective sample size for Brownian
motion and Ornstein-Uhlenbeck processes - the regression effective sample size.
I discuss how mutual information can be used to define an effective sample size
in the non-normal process case and compare these two definitions to an already
present concept of effective sample size (the mean effective sample size).
Through a simulation study I find that the AICc is robust if one corrects for
the number of species or effective number of species. Lastly I discuss how the
concept of the phylogenetic effective sample size can be useful for
biodiversity quantification, identification of interesting clades and deciding
on the importance of phylogenetic correlations.


  One of the goals of phylogenetic research is to find the species tree
describing the evolutionary history of a set of species. But the trees derived
from geneti data with the help of tree inference methods are gene trees that
need not coincide with the species tree. This can for example happen when
so-called deep coalescence events take place. It is also known that species
trees can differ from their most likely gene trees. Therefore, as a means to
find the species tree, it has been suggested to use subtrees of the gene trees,
for example triples, and to puzzle them together in order to find the species
tree. In this paper, we will show that this approach may lead to wrong trees
regarding the minimum deep coalescence criterion (MDC). In particular, we
present an example in which the optimal MDC tree is unique, but none of its
triple subtrees fulfills the MDC criterion. In this sense, MDC is a
non-hereditary tree reconstruction method.


  We propose a stochastic individual-based model of the progression of atrial
fibrillation (AF). The model operates at patient level over a lifetime and is
based on elements of the physiology and biophysics of AF, making contact with
existing mechanistic models. The outputs of the model are times when the
patient is in normal rhythm and AF, and we carry out a population-level
analysis of the statistics of disease progression. While the model is stylised
at present and not directly predictive, future improvements are proposed to
tighten the gap between existing mechanistic models of AF, and epidemiological
data, with a view towards model-based personalised medicine.


  We introduce two new algorithms to minimise smooth difference of convex (DC)
functions that accelerate the convergence of the classical DC algorithm (DCA).
We prove that the point computed by DCA can be used to define a descent
direction for the objective function evaluated at this point. Our algorithms
are based on a combination of DCA together with a line search step that uses
this descent direction. Convergence of the algorithms is proved and the rate of
convergence is analysed under the Lojasiewicz property of the objective
function. We apply our algorithms to a class of smooth DC programs arising in
the study of biochemical reaction networks, where the objective function is
real analytic and thus satisfies the Lojasiewicz property. Numerical tests on
various biochemical models clearly show that our algorithms outperforms DCA,
being on average more than four times faster in both computational time and the
number of iterations. Numerical experiments show that the algorithms are
globally convergent to a non-equilibrium steady state of various biochemical
networks, with only chemically consistent restrictions on the network topology.


  Cells constantly need to monitor the state of the environment to detect
changes and timely respond. The detection of concentration changes of a ligand
by a set of receptors can be cast as a problem of hypothesis testing, and the
cell viewed as a Neyman-Pearson detector. Within this framework, we investigate
the role of receptor cooperativity in improving the cell's ability to detect
changes. We find that cooperativity decreases the probability of missing an
occurred change. This becomes especially beneficial when difficult detections
have to be made. Concerning the influence of cooperativity on how fast a
desired detection power is achieved, we find in general that there is an
optimal value at finite levels of cooperation, even though easy discrimination
tasks can be performed more rapidly by noncooperative receptors.


  Open-ended evolutionary dynamics remains an elusive goal for artificial
evolutionary systems. Many ideas exist in the biological literature beyond the
basic Darwinian requirements of variation, differential reproduction and
inheritance. I argue that these ideas can be seen as aspects of five
fundamental requirements for open-ended evolution: (1) robustly reproductive
individuals, (2) a medium allowing the possible existence of a practically
unlimited diversity of individuals and interactions, (3) individuals capable of
producing more complex offspring, (4) mutational pathways to other viable
individuals, and (5) drive for continued evolution. I briefly discuss
implications of this view for the design of artificial systems with greater
evolutionary potential.


  Above two dimensions, diffusion of a particle in a medium with quenched
random traps is believed to be well-described by the annealed continuous time
random walk (CTRW). We propose an approximate expression for the
first-passage-time (FPT) distribution in a given sample that enables detailed
comparison of the two problems. For a system of finite size, the number and
spatial arrangement of deep traps yield significant sample-to-sample variations
in the FPT statistics. Numerical simulations of a quenched trap model with
power-law sojourn times confirm the existence of two characteristic time scales
and a non-self-averaging FPT distribution, as predicted by theory.


  In the past years, a remarkable mapping has been found between the dynamics
of a population of M individuals undergoing random mutations and selection, and
that of a single system in contact with a thermal bath with temperature 1/M.
This correspondence holds under the somewhat restrictive condition that the
population is dominated by a single type at almost all times, punctuated by
rare successive mutations. Here we argue that such thermal dynamics will hold
more generally, specifically in systems with rugged fitness landscapes. This
includes cases with strong clonal interference, where a number of concurrent
mutants dominate the population. The problem becomes closely analogous to the
experimental situation of glasses subjected to controlled variations of
parameters such as temperature, pressure or magnetic fields. Non-trivial
suggestions from the field of glasses may be thus proposed for evolutionary
systems - including a large part of the numerical simulation procedures - that
in many cases would have been counter intuitive without this background.