  Long-distance regulatory interactions between enhancers and their target
genes are commonplace in higher eukaryotes. Interposed boundaries or insulators
are able to block these long distance regulatory interactions. The mechanistic
basis for insulator activity and how it relates to enhancer
action-at-a-distance remains unclear. Here we explore the idea that topological
loops could simultaneously account for regulatory interactions of distal
enhancers and the insulating activity of boundary elements. We show that while
loop formation is not in itself sufficient to explain action at a distance,
incorporating transient non-specific and moderate attractive interactions
between the chromatin fibers strongly enhances long-distance regulatory
interactions and is sufficient to generate a euchromatin-like state. Under
these same conditions, the subdivision of the loop into two topologically
independent loops by insulators inhibits inter-domain interactions. The
underlying cause of this effect is a suppression of crossings in the contact
map at intermediate distances. Thus our model simultaneously accounts for
regulatory interactions at a distance and the insulator activity of boundary
elements. This unified model of the regulatory roles of chromatin loops makes
several testable predictions that could be confronted with \emph{in vitro}
experiments, as well as genomic chromatin conformation capture and fluorescent
microscopic approaches.


  An algorithm is given in this paper for the computation of dynamically
equivalent weakly reversible realizations with the maximal number of reactions,
for chemical reaction networks (CRNs) with mass action kinetics.


  No existing algorithm can start with arbitrary RNA sequences and return the
precise, three-dimensional structures that ensures their biological function.
This chapter outlines current algorithms for automated RNA structure prediction
(including our own FARNA-FARFAR), highlights their successes, and dissects
their limitations, using a tetraloop and the sarcin/ricin motif as examples.
The barriers to future advances are considered in light of three particular
challenges: improving computational sampling, reducing reliance on
experimentally solved structures, and avoiding coarse-grained representations
of atomic-level interactions. To help meet these challenges and better
understand the current state of the field, we propose an ongoing community-wide
CASP-style experiment for evaluating the performance of current structure
prediction algorithms.


  The Weibull function is widely used to describe skew distributions observed
in nature. However, the origin of this ubiquity is not always obvious to
explain. In the present paper, we consider the well-known Galton-Watson
branching process describing simple replicative systems. The shape of the
resulting distribution, about which little has been known, is found essentially
indistinguishable from the Weibull form in a wide range of the branching
parameter; this can be seen from the exact series expansion for the cumulative
distribution, which takes a universal form. We also find that the branching
process can be mapped into a process of aggregation of clusters. In the
branching and aggregation process, the number of events considered for
branching and aggregation grows cumulatively in time, whereas, for the binomial
distribution, an independent event occurs at each time with a given success
probability.


  A complete macromolecule modeling package must be able to solve the simplest
structure prediction problems. Despite recent successes in high resolution
structure modeling and design, the Rosetta software suite fares poorly on
deceptively small protein and RNA puzzles, some as small as four residues. To
illustrate these problems, this manuscript presents extensive Rosetta results
for four well-defined test cases: the 20-residue mini-protein Trp cage, an even
smaller disulfide-stabilized conotoxin, the reactive loop of a serine protease
inhibitor, and a UUCG RNA tetraloop. In contrast to previous Rosetta studies,
several lines of evidence indicate that conformational sampling is not the
major bottleneck in modeling these small systems. Instead, approximations and
omissions in the Rosetta all-atom energy function currently preclude
discriminating experimentally observed conformations from de novo models at
atomic resolution. These molecular "puzzles" should serve as useful model
systems for developers wishing to make foundational improvements to this
powerful modeling suite.


  A simple stochastic model which describes microtubule dynamics and explicitly
takes into account the relevant biochemical processes is presented. The model
incorporates binding and unbinding of monomers and random phosphate release
inside the polymer. It is shown that this theoretical approach provides a
microscopic picture of the dynamic instability phenomena of microtubules. The
cap size, the concentration dependence of the catastrophe times and the delay
before observing catastrophes following a dilution can be quantitatively
predicted by this approach in a direct and simple way. Furthermore, the model
can be solved analytically to a large extend, thus offering a valuable starting
point for more refined studies of microtubules dynamics.


  We present a general framework to study the thermodynamic denaturation of
double-stranded DNA under superhelical stress. We report calculations of
position- and size-dependent opening probabilities for bubbles along the
sequence. Our results are obtained from transfer-matrix solutions of the
Zimm-Bragg model for unconstrained DNA and of a self-consistent linearization
of the Benham model for superhelical DNA. The numerical efficiency of our
method allows for the analysis of entire genomes and of random sequences of
corresponding length ($10^6-10^9$ base pairs). We show that, at physiological
conditions, opening in superhelical DNA is strongly cooperative with average
bubble sizes of $10^2-10^3$ base pairs (bp), and orders of magnitude higher
than in unconstrained DNA. In heterogeneous sequences, the average degree of
base-pair opening is self-averaging, while bubble localization and statistics
are dominated by sequence disorder. Compared to random sequences with identical
GC-content, genomic DNA has a significantly increased probability to open large
bubbles under superhelical stress. These bubbles are frequently located
directly upstream of transcription start sites.


  Volcano plot displays unstandardized signal (e.g. log-fold-change) against
noise-adjusted/standardized signal (e.g. t-statistic or -log10(p-value) from
the t test). We review the basic and an interactive use of the volcano plot,
and its crucial role in understanding the regularized t-statistic. The joint
filtering gene selection criterion based on regularized statistics has a curved
discriminant line in the volcano plot, as compared to the two perpendicular
lines for the "double filtering" criterion. This review attempts to provide an
unifying framework for discussions on alternative measures of differential
expression, improved methods for estimating variance, and visual display of a
microarray analysis result. We also discuss the possibility to apply volcano
plots to other fields beyond microarray.


  Structure formation and constant reorganization of the actin cytoskeleton are
key requirements for the function of living cells. Here we show that a minimal
reconstituted system consisting of actin filaments, crosslinking molecules and
molecular-motor filaments exhibits a generic mechanism of structure formation,
characterized by a broad distribution of cluster sizes. We demonstrate that the
growth of the structures depends on the intricate balance between
crosslinker-induced stabilization and simultaneous destabilization by molecular
motors, a mechanism analogous to nucleation and growth in passive systems. We
also show that the intricate interplay between force generation, coarsening and
connectivity is responsible for the highly dynamic process of structure
formation in this heterogeneous active gel, and that these competing mechanisms
result in anomalous transport, reminiscent of intracellular dynamics.


  Time-dependent density functional theory (TD-DFT) and correlated ab initio
methods have been applied to the electronically excited states of vitamin B12
(cyanocobalamin or CNCbl). Different experimental techniques have been used to
probe the excited states of CNCbl, revealing many issues that remain poorly
understood from an electronic structure point of view. Due to its efficient
scaling with size, TD-DFT emerges as one of the most practical tools that can
be used to predict the electronic properties of these fairly complex molecules.
However, the description of excited states is strongly dependent on the type of
functional used in the calculations. In the present contribution, the choice of
a proper functional for vitamin B12 was evaluated in terms of its agreement
with both experimental results and correlated ab initio calculations. Three
different functionals, i.e. B3LYP, BP86, and LC-BLYP, were tested. In addition,
the effect of relative contributions of DFT and HF to the exchange-correlation
functional was investigated as a function of the range-separation parameter,
{\mu}. The issues related to the underestimation of charge transfer (CT)
excitation energies by TD-DFT was validated by Tozer's L diagnostic, which
measures the spatial overlap between occupied and virtual orbitals involved in
the particular excitation. The nature of low-lying excited states was also
analyzed based on a comparison of TD-DFT and ab initio results. Based on an
extensive comparision against experimental results and ab initio benchmark
calculations, the BP86 functional was found to be the most appropriate in
describing the electronic properties of CNCbl. Finally, an analysis of
electronic transitions and a new re-assignment of some excitations are
discussed.


  Lasso peptides constitute a class of bioactive peptides sharing a knotted
structure where the C-terminal tail of the peptide is threaded through and
trapped within an N-terminalmacrolactamring. The structural characterization of
lasso structures and differentiation from their unthreaded topoisomers is not
trivial and generally requires the use of complementary biochemical and
spectroscopic methods. Here we investigated two antimicrobial peptides
belonging to the class II lasso peptide family and their corresponding
unthreaded topoisomers: microcin J25 (MccJ25), which is known to yield
two-peptide product ions specific of the lasso structure under collisioninduced
dissociation (CID), and capistruin, for which CID does not permit to
unambiguously assign the lasso structure. The two pairs of topoisomers were
analyzed by electrospray ionization Fourier transform ion cyclotron resonance
mass spectrometry (ESI-FTICR MS) upon CID, infrared multiple photon
dissociation (IRMPD), and electron capture dissociation (ECD). CID and
ECDspectra clearly permitted to differentiate MccJ25 from its non-lasso
topoisomer MccJ25-Icm, while for capistruin, only ECD was informative and
showed different extent of hydrogen migration (formation of c\bullet/z from
c/z\bullet) for the threaded and unthreaded topoisomers. The ECD spectra of the
triply-charged MccJ25 and MccJ25-lcm showed a series of radical b-type product
ions {\eth}b0In{\TH}. We proposed that these ions are specific of
cyclic-branched peptides and result from a dual c/z\bullet and y/b
dissociation, in the ring and in the tail, respectively. This work shows the
potentiality of ECD for structural characterization of peptide topoisomers, as
well as the effect of conformation on hydrogen migration subsequent to electron
capture.


  Nanospray and collisionally-induced dissociation are used to evaluate the
presence and absence of interstrand co-chelation of zinc ions in dimers of
metallothionein. As was reported in a previous publication from this
laboratory, co-chelation stabilizes the dimer to collisional activation, and
facilitates asymmetrical zinc ion transfers during fragmentation. In the case
of metallothionein, dimers of the holoprotein are found to share zinc ions,
while dimers of metallothionein, in which one domain has been denatured, do
not. Zinc ions are silent to most physicochemical probes, e.g., NMR and
Mossbauer spectroscopies, and the capability of mass spectrometry to provide
information on zinc complexes has widespread potential application in
biochemistry.


  Epidemics seldom occur as isolated phenomena. Typically, two or more viral
agents spread within the same host population and may interact dynamically with
each other. We present a general model where two viral agents interact via an
immunity mechanism as they propagate simultaneously on two networks connecting
the same set of nodes. Exploiting a correspondence between the propagation
dynamics and a dynamical process performing progressive network generation, we
develop an analytic approach that accurately captures the dynamical interaction
between epidemics on overlay networks. The formalism allows for overlay
networks with arbitrary joint degree distribution and overlap. To illustrate
the versatility of our approach, we consider a hypothetical delayed
intervention scenario in which an immunizing agent is disseminated in a host
population to hinder the propagation of an undesirable agent (e.g. the spread
of preventive information in the context of an emerging infectious disease).


  The dynamics of classical hard particles in a quasi one-dimensional channel
were studied since the 1960s and used for explaining processes in chemistry,
physics and biology and in applications. Here we show that in a previously
un-described file made of anomalous, independent, particles (with jumping times
taken from, {\psi}_{\alpha} (t) t^(-1-{\alpha}), 0<{\alpha}<1), particles form
clusters. At steady state, the percentage of particles in clusters is about,
\surd(1-{\alpha}^3), only for anomalous {\alpha}, characterizing the
criticality of a phase transition. The asymptotic mean square displacement per
particle in the file is about, log^2(t). We show numerically that this exciting
phenomenon of a phase transition is very stable, and relate it with the
mysterious phenomenon of rafts in biological membranes, and with regulation of
biological channels.


  Using analytical techniques and Langevin dynamics simulations, we investigate
the dynamics of polymer translocation into a narrow channel of width $R$
embedded in two dimensions, driven by a force proportional to the number of
monomers in the channel. Such a setup mimics typical experimental situations in
nano/micro-fluidics. During the the translocation process if the monomers in
the channel can sufficiently quickly assume steady state motion, we observe the
scaling $\tau\sim N/F$ of the translocation time $\tau$ with the driving force
$F$ per bead and the number $N$ of monomers per chain. With smaller channel
width $R$, steady state motion cannot be achieved, effecting a non-universal
dependence of $\tau$ on $N$ and $F$. From the simulations we also deduce the
waiting time distributions under various conditions for the single segment
passage through the channel entrance. For different chain lengths but the same
driving force, the curves of the waiting time as a function of the
translocation coordinate $s$ feature a maximum located at identical
$s_{\mathrm{max}}$, while with increasing the driving force or the channel
width the value of $s_{\mathrm{max}}$ decreases.


  We present a novel optimal allocation model for perennial plants, in which
assimilates are not allocated directly to vegetative or reproductive parts but
instead go first to a storage compartment from where they are then optimally
redistributed. We do not restrict considerations purely to periods favourable
for photosynthesis, as it was done in published models of perennial species,
but analyse the whole life period of a perennial plant. As a result, we obtain
the general scheme of perennial plant development, for which annual and
monocarpic strategies are special cases.
  We not only re-derive predictions from several previous optimal allocation
models, but also obtain more information about plants' strategies during
transitions between favourable and unfavourable seasons. One of the model's
predictions is that a plant can begin to re-establish vegetative tissues from
storage, some time before the beginning of favourable conditions, which in turn
allows for better production potential when conditions become better. By means
of numerical examples we show that annual plants with single or multiple
reproduction periods, monocarps, evergreen perennials and polycarpic perennials
can be studied successfully with the help of our unified model.
  Finally, we build a bridge between optimal allocation models and models
describing trade-offs between size and the number of seeds: a modelled plant
can control the distribution of not only allocated carbohydrates but also seed
size. We provide sufficient conditions for the optimality of producing the
smallest and largest seeds possible.


  We present an extension of Felsenstein's algorithm to indel models defined on
entire sequences, without the need to condition on one multiple alignment. The
algorithm makes use of a generalization from probabilistic substitution
matrices to weighted finite-state transducers. Our approach may equivalently be
viewed as a probabilistic formulation of progressive multiple sequence
alignment, using partial-order graphs to represent ensemble profiles of
ancestral sequences. We present a hierarchical stochastic approximation
technique which makes this algorithm tractable for alignment analyses of
reasonable size.


  Spatial evolutionary games are studied with myopic players whose payoff
interest, as a personal character, is tuned from selfishness to other-regarding
preference via fraternity. The players are located on a square lattice and
collect income from symmetric two-person two-strategy (called cooperation and
defection) games with their nearest neighbors. During the elementary steps of
evolution a randomly chosen player modifies her strategy in order to maximize
stochastically her utility function composed from her own and the co-players'
income with weight factors $1-Q$ and Q. These models are studied within a wide
range of payoff parameters using Monte Carlo simulations for noisy strategy
updates and by spatial stability analysis in the low noise limit. For fraternal
players ($Q=1/2$) the system evolves into ordered arrangements of strategies in
the low noise limit in a way providing optimum payoff for the whole society.
Dominance of defectors, representing the "tragedy of the commons", is found
within the regions of prisoner's dilemma and stag hunt game for selfish players
(Q=0). Due to the symmetry in the effective utility function the system
exhibits similar behavior even for Q=1 that can be interpreted as the "lovers'
dilemma".


  Regulatory dynamics in biology is often described by continuous rate
equations for continuously varying chemical concentrations. Binary
discretization of state space and time leads to Boolean dynamics. In the
latter, the dynamics has been called unstable if flip perturbations lead to
damage spreading. Here we find that this stability classification strongly
differs from the stability properties of the original continuous dynamics under
small perturbations of the state vector. In particular, random networks of
nodes with large sensitivity yield stable dynamics under small perturbations.


  Considering the recent experimental discovery of Green et al that present day
non-Africans have 1 to 4% of their nuclear DNA of Neanderthal origin, we
propose here a model which is able to quantify the interbreeding events between
Africans and Neanderthals at the time they coexisted in the Middle East. The
model consists of a solvable system of deterministic ordinary differential
equations containing as a stochastic ingredient a realization of the neutral
Wright-Fisher drift process. By simulating the stochastic part of the model we
are able to apply it to the interbreeding of African and Neanderthal
subpopulations and estimate the only parameter of the model, which is the
number of individuals per generation exchanged between subpopulations. Our
results indicate that the amount of Neanderthal DNA in non-Africans can be
explained with maximum probability by the exchange of a single pair of
individuals between the subpopulations at each 77 generations, but larger
exchange frequencies are also allowed with sizeable probability. The results
are compatible with a total interbreeding population of order 10,000
individuals and with all living humans being descendents of Africans both for
mitochondrial DNA and Y chromosome.


  Understanding why we age is a long-lived open problem in evolutionary
biology. Aging is prejudicial to the individual and evolutionary forces should
prevent it, but many species show signs of senescence as individuals age. Here,
I will propose a model for aging based on assumptions that are compatible with
evolutionary theory: i) competition is between individuals; ii) there is some
degree of locality, so quite often competition will between parents and their
progeny; iii) optimal conditions are not stationary, mutation helps each
species to keep competitive. When conditions change, a senescent species can
drive immortal competitors to extinction. This counter-intuitive result arises
from the pruning caused by the death of elder individuals. When there is change
and mutation, each generation is slightly better adapted to the new conditions,
but some older individuals survive by random chance. Senescence can eliminate
those from the genetic pool. Even though individual selection forces always win
over group selection ones, it is not exactly the individual that is selected,
but its lineage. While senescence damages the individuals and has an
evolutionary cost, it has a benefit of its own. It allows each lineage to adapt
faster to changing conditions. We age because the world changes.


  An algorithm is given in this paper for the computation of dynamically
equivalent weakly reversible realizations with the maximal number of reactions,
for chemical reaction networks (CRNs) with mass action kinetics. The original
problem statement can be traced back at least 30 years ago. The algorithm uses
standard linear and mixed integer linear programming, and it is based on
elementary graph theory and important former results on the dense realizations
of CRNs. The proposed method is also capable of determining if no dynamically
equivalent weakly reversible structure exists for a given reaction network with
a previously fixed complex set.


  Human mobility and activity patterns mediate contagion on many levels,
including the spatial spread of infectious diseases, diffusion of rumors, and
emergence of consensus. These patterns however are often dominated by specific
locations and recurrent flows and poorly modeled by the random diffusive
dynamics generally used to study them. Here we develop a theoretical framework
to analyze contagion within a network of locations where individuals recall
their geographic origins. We find a phase transition between a regime in which
the contagion affects a large fraction of the system and one in which only a
small fraction is affected. This transition cannot be uncovered by continuous
deterministic models due to the stochastic features of the contagion process
and defines an invasion threshold that depends on mobility parameters,
providing guidance for controlling contagion spread by constraining mobility
processes. We recover the threshold behavior by analyzing diffusion processes
mediated by real human commuting data.


  A general theory of top-down cascades in complex networks is described which
explains two similar types of perturbation amplifications in the complex
networks of business supply chains (the `bullwhip effect') and ecological food
webs (trophic cascades). The dependence of the strength of the effects on the
interaction strength and covariance in the dynamics as well as the graph
structure allows both explanation and prediction of widely recognized effects
in each type of system.


  Stability of evolutionary dynamics of non-repeated Prisoner's Dilemma game
with non-uniform interaction rates [1], via benefit and cost dilemma is studied
. Moreover, the stability condition (b+c/b-c)2 < r1r3 is derived in case of
coexistence between cooperators and defectors .If r1,r3 -> infinity cooperation
is the dominant strategy and defectors can no longer exploit cooperators.


  Gene-knockout experiments on single-cell organisms have established that
expression of a substantial fraction of genes is not needed for optimal growth.
This problem acquired a new dimension with the recent discovery that
environmental and genetic perturbations of the bacterium Escherichia coli are
followed by the temporary activation of a large number of latent metabolic
pathways, which suggests the hypothesis that temporarily activated reactions
impact growth and hence facilitate adaptation in the presence of perturbations.
Here we test this hypothesis computationally and find, surprisingly, that the
availability of latent pathways consistently offers no growth advantage, and
tends in fact to inhibit growth after genetic perturbations. This is shown to
be true even for latent pathways with a known function in alternate conditions,
thus extending the significance of this adverse effect beyond apparently
nonessential genes. These findings raise the possibility that latent pathway
activation is in fact derivative of another, potentially suboptimal, adaptive
response.


  It is well known that non-enzymatic template directed molecular replicators X
+ nO ---> 2X exhibit parabolic growth d[X]/dt = k [X]^{1/2}. Here, we analyze
the dependence of the effective replication rate constant k on hybridization
energies, temperature, strand length, and sequence composition. First we derive
analytical criteria for the replication rate k based on simple thermodynamic
arguments. Second we present a Brownian dynamics model for oligonucleotides
that allows us to simulate their diffusion and hybridization behavior. The
simulation is used to generate and analyze the effect of strand length,
temperature, and to some extent sequence composition, on the hybridization
rates and the resulting optimal overall rate constant k. Combining the two
approaches allows us to semi-analytically depict a fitness landscape for
template directed replicators. The results indicate a clear replication
advantage for longer strands at low temperatures.


  Single-nucleotide-resolution chemical mapping for structured RNA is being
rapidly advanced by new chemistries, faster readouts, and coupling to
computational algorithms. Recent tests have shown that selective 2'-hydroxyl
acylation by primer extension (SHAPE) can give near-zero error rates (0-2%) in
modeling the helices of RNA secondary structure. Here, we benchmark the method
using six molecules for which crystallographic data are available: tRNA(phe)
and 5S rRNA from Escherichia coli, the P4-P6 domain of the Tetrahymena group I
ribozyme, and ligand-bound domains from riboswitches for adenine, cyclic
di-GMP, and glycine. SHAPE-directed modeling of these highly structured RNAs
gave an overall false negative rate (FNR) of 17% and a false discovery rate
(FDR) of 21%, with at least one helix prediction error in five of the six
cases. Extensive variations of data processing, normalization, and modeling
parameters did not significantly mitigate modeling errors. Only one varation,
filtering out data collected with deoxyinosine triphosphate during primer
extension, gave a modest improvement (FNR = 12%, and FDR = 14%). The residual
structure modeling errors are explained by the insufficient information content
of these RNAs' SHAPE data, as evaluated by a nonparametric bootstrapping
analysis. Beyond these benchmark cases, bootstrapping suggests a low level of
confidence (<50%) in the majority of helices in a previously proposed
SHAPE-directed model for the HIV-1 RNA genome. Thus, SHAPE-directed RNA
modeling is not always unambiguous, and helix-by-helix confidence estimates, as
described herein, may be critical for interpreting results from this powerful
methodology.


  An integro-differential equation on a tree graph is used to model the
evolution and spatial distribution of a population of organisms in a river
network. Individual organisms become mobile at a constant rate, and disperse
according to an advection-diffusion process with coefficients that are constant
on the edges of the graph. Appropriate boundary conditions are imposed at the
outlet and upstream nodes of the river network. The local rates of population
growth/decay and that by which the organisms become mobile, are assumed
constant in time and space. Imminent extinction of the population is understood
as the situation whereby the zero solution to the integro-differential equation
is stable. Lower and upper bounds for the eigenvalues of the dispersion
operator, and related Sturm-Liouville problems are found, and therefore
sufficient conditions for imminent extinction are given in terms of the
physical variables of the problem.


  The widespread relevance of complex networks is a valuable tool in the
analysis of a broad range of systems. There is a demand for tools which enable
the extraction of meaningful information and allow the comparison between
different systems. We present a novel measure of similarity between nodes in
different networks as a generalization of the concept of self-similarity. A
similarity matrix is assembled as the distance between feature vectors that
contain the in and out paths of all lengths for each node. Hence, nodes
operating in a similar flow environment are considered similar regardless of
network membership. We demonstrate that this method has the potential to be
influential in tasks such as assigning identity or function to uncharacterized
nodes. In addition an innovative application of graph partitioning to the raw
results extends the concept to the comparison of networks in terms of their
underlying role-structure.


  A non-linear Penner type interaction is introduced and studied in the random
matrix model of homo-RNA. The asymptotics in length of the partition function
is discussed for small and large $N$ (size of matrix). The interaction doubles
the coupling ($v$) between the bases and the dependence of the combinatoric
factor on ($v,N$) is found. For small $N$, the effect of interaction changes
the power law exponents for the secondary and tertiary structures. The specific
heat shows different analytical behavior in the two regions of $N$, with a
peculiar double peak in its second derivative for N=1 at low temperature.
Tapping the model indicates the presence of multiple solutions.


  The key findings of classical population genetics are derived using a
framework based on information theory using the entropies of the allele
frequency distribution as a basis. The common results for drift, mutation,
selection, and gene flow will be rewritten both in terms of information
theoretic measurements and used to draw the classic conclusions for balance
conditions and common features of one locus dynamics. Linkage disequilibrium
will also be discussed including the relationship between mutual information
and r^2 and a simple model of hitchhiking.


  When beneficial mutations are relatively common, competition between multiple
unfixed mutations can reduce the rate of fixation in well-mixed asexual
populations. We introduce a one dimensional model with a steady accumulation of
beneficial mutations. We find a transition between periodic selection and
multiple-mutation regimes. In the multiple-mutation regime, the increase of
fitness along the lattice bears a striking similarity to surface growth
phenomena, with power law growth and saturation of the interface width. We also
find significant differences compared to the well-mixed model. In our lattice
model, the transition between regimes happens at a much lower mutation rate due
to slower fixation times in one dimension. Also the rate of fixation is reduced
with increasing mutation rate due to the more intense competition, and it
saturates with large population size.


  The use of fossil evidence to calibrate divergence time estimation has a long
history. More recently Bayesian MCMC has become the dominant method of
divergence time estimation and fossil evidence has been re-interpreted as the
specification of prior distributions on the divergence times of calibration
nodes. These so-called "soft calibrations" have become widely used but the
statistical properties of calibrated tree priors in a Bayesian setting has not
been carefully investigated. Here we clarify that calibration densities, such
as those defined in BEAST 1.5, do not represent the marginal prior distribution
of the calibration node. We illustrate this with a number of analytical results
on small trees. We also describe an alternative construction for a calibrated
Yule prior on trees that allows direct specification of the marginal prior
distribution of the calibrated divergence time, with or without the restriction
of monophyly. This method requires the computation of the Yule prior
conditional on the height of the divergence being calibrated. Unfortunately, a
practical solution for multiple calibrations remains elusive. Our results
suggest that direct estimation of the prior induced by specifying multiple
calibration densities should be a prerequisite of any divergence time dating
analysis.


  Unlike their model membrane counterparts, biological membranes are richly
decorated with a heterogeneous assembly of membrane proteins. These proteins
are so tightly packed that their excluded area interactions can alter the free
energy landscape controlling the conformational transitions suffered by such
proteins. For membrane channels, this effect can alter the critical membrane
tension at which they undergo a transition from a closed to an open state, and
therefore influence protein function \emph{in vivo}. Despite their obvious
importance, crowding phenomena in membranes are much less well studied than in
the cytoplasm.
  Using statistical mechanics results for hard disk liquids, we show that
crowding induces an entropic tension in the membrane, which influences
transitions that alter the projected area and circumference of a membrane
protein. As a specific case study in this effect, we consider the impact of
crowding on the gating properties of bacterial mechanosensitive membrane
channels, which are thought to confer osmoprotection when these cells are
subjected to osmotic shock. We find that crowding can alter the gating energies
by more than $2\;k_BT$ in physiological conditions, a substantial fraction of
the total gating energies in some cases.
  Given the ubiquity of membrane crowding, the nonspecific nature of excluded
volume interactions, and the fact that the function of many membrane proteins
involve significant conformational changes, this specific case study highlights
a general aspect in the function of membrane proteins.


  Theories to explain the prevalence of sex and recombination have long been a
central theme of evolutionary biology. Yet despite decades of attention
dedicated to the evolution of sex and recombination, the widespread pattern of
sex-differences in the recombination rate is not well understood and has
received relatively little theoretical attention. Here, we argue that female
meiotic drivers - alleles that increase in frequency by exploiting the
asymmetric cell division of oogenesis - present a potent selective pressure
favoring the modification of the female recombination rate. Because
recombination plays a central role in shaping patterns of variation within and
among dyads, modifiers of the female recombination rate can function as potent
suppressors or enhancers of female meiotic drive. We show that when female
recombination modifiers are unlinked to female drivers, recombination modifiers
that suppress harmful female drive can spread. By contrast, a recombination
modifier tightly linked to a driver can increase in frequency by enhancing
female drive. Our results predict that rapidly evolving female recombination
rates, particularly around centromeres, should be a common outcome of meiotic
drive. We discuss how selection to modify the efficacy of meiotic drive may
contribute to commonly observed patterns of sex-differences in recombination.


  Compartmentalization into biochemically distinct organelles constantly
exchanging material is one of the hallmarks of eukaryotic cells. In the most
naive picture of inter-organelle transport driven by concentration gradients,
concentration differences between organelles should relax. We determine the
conditions under which cooperative transport, i.e. based on molecular
recognition, allows for the existence and maintenance of distinct organelle
identities. Cooperative transport is also shown to control the flux of material
transiting through a compartmentalized system, dramatically increasing the
transit time under high incoming flux. By including chemical processing of the
transported species, we show that this property provides a strong functional
advantage to a system responsible for protein maturation and sorting.


  We investigate a model for spatial epidemics explicitly taking into account
bi-directional movements between base and destination locations on individual
mobility networks. We provide a systematic analysis of generic dynamical
features of the model on regular and complex metapopulation network topologies
and show that significant dynamical differences exist to ordinary
reaction-diffusion and effective force of infection models. On a lattice we
calculate an expression for the velocity of the propagating epidemic front and
find that in contrast to the diffusive systems, our model predicts a saturation
of the velocity with increasing traveling rate. Furthermore, we show that a
fully stochastic system exhibits a novel threshold for attack ratio of an
outbreak absent in diffusion and force of infection models. These insights not
only capture natural features of human mobility relevant for the geographical
epidemic spread, they may serve as a starting point for modeling important
dynamical processes in human and animal epidemiology, population ecology,
biology and evolution.


  A new understanding of the consequences of how ecosystem elements are
interconnected is emerging from the development and application of Ecological
Network Analysis. The relative importance of indirect effects is central to
this understanding, and the ratio of indirect flow to direct flow (I/D) is one
indicator of their importance. Two methods have been proposed for calculating
this indicator. The unit approach shows what would happen if each system member
had a unit input or output, while the realized technique determines the ratio
using the observed system inputs or outputs. When using the unit method, the
input oriented and output oriented ratios can be different, potentially leading
to conflicting results. However, we show that the input and output oriented I/D
ratios are identical using the realized method when the system is at steady
state. This work is a step in the maturation of Ecological Network Analysis
that will let it be more readily testable empirically and ultimately more
useful for environmental assessment and management.


  Connectivity patterns of ecological elements are often the core concern of
ecologists working at multiple levels of organization (e.g., populations,
ecosystems, and landscapes) because these patterns often reflect the forces
shaping the system's development as well as constraining their operation. One
reason these patterns of direct connections are critical is that they establish
the pathways through which elements influence each other indirectly. Here, we
tested a hypothesized consequence of connectivity in ecosystems: the
homogenization of resource distributions in flow networks. Specifically, we
tested the generality of the systems ecology hypothesis of resource
homogenization in 50 empirically derived trophic ecosystem models representing
35 distinct ecosystems. We applied Ecological Network Analysis to calculate
resource homogenization for these models. We further evaluated the robustness
of our results in two ways. First, we verified the close correspondence between
the input- and output-oriented homogenization values to ensure that our results
were not biased by our decision to focus on the output orientation. Second, we
conducted a Monte Carlo based uncertainty analysis to determine the robustness
of our results to +/-5% error introduced into the original flow matrices for
each model. Our results show that resource homogenization occurs universally in
the 50 ecosystem models tested. We confirm that our results are not biased by
using the output-oriented homogenization values because there is a significant
linear regression between the input and output oriented homogenization (r^2 =
0.38, p < 0.001). Finally, we found that our results are robust to +/-5% error
in the flow matrices. In conclusion, we found strong support for the resource
homogenization hypothesis in 50 empirically derived ecosystem models.


  We study complex networks in which the nodes of the network are tagged with
different colors depending on the functionality of the nodes (colored graphs),
using information theory applied to the distribution of motifs in such
networks. We find that colored motifs can be viewed as the building blocks of
the networks (much more so than the uncolored structural motifs can be) and
that the relative frequency with which these motifs appear in the network can
be used to define the information content of the network. This information is
defined in such a way that a network with random coloration (but keeping the
relative number of nodes with different colors the same) has zero color
information content. Thus, colored motif information captures the
exceptionality of coloring in the motifs that is maintained via selection. We
study the motif information content of the C. elegans brain as well as the
evolution of colored motif information in networks that reflect the interaction
between instructions in genomes of digital life organisms. While we find that
colored motif information appears to capture essential functionality in the C.
elegans brain (where the color assignment of nodes is straightforward) it is
not obvious whether the colored motif information content always increases
during evolution, as would be expected from a measure that captures network
complexity. For a single choice of color assignment of instructions in the
digital life form Avida, we find rather that colored motif information content
increases or decreases during evolution, depending on how the genomes are
organized, and therefore could be an interesting tool to dissect genomic
rearrangements.


  Translation is the cellular process in which ribosomes make proteins from
information encoded on messenger RNA (mRNA). We model translation with an
exclusion process taking into account the experimentally determined,
non-exponential, waiting time between steps of a ribosome. From numerical
simulations using realistic parameter values, we determine the distribution
P(E) of the number of proteins E produced by one mRNA. We find that for small E
this distribution is not geometric. We present a simplified and analytically
solvable semi-Markov model that relates P(E) to the distributions of the times
to produce the first E proteins.


  In this paper the realized niche of the Snowy Plover (Charadrius
alexandrinus), a primarily resident Florida shorebird, is described as a
function of the scenopoetic and bionomic variables at the nest-, landscape-,
and regional-scale. We identified some pos- sible geomorphological controls
that influence nest-site selection and survival using data collected along the
Florida Gulf coast. In particular we focused on the effects of beach
replenishment interventions on the Snowy Plover (SP), and on the migratory
Piping Plover (PP) (Charadrius melodus) and Red Knot (RK) (Calidris canutus).
To quantify the relationship between past renourishment projects and shorebird
species we used a Monte Carlo procedure to sample from the posterior
distribution of the binomial probabilities that a region is not a nesting or a
wintering ground conditional on the occurrence of a beach replenishment
intervention in the same and the previous year. The results indicate that it
was 2.3, 3.1, and 0.8 times more likely that a region was not a wintering
ground following a year with a renourishment intervention for the SP, PP and RK
respectively. For the SP it was 2.5. times more likely that a region was not a
breeding ground after a renourishment event. Through a maximum entropy
principle model we observed small differences in the habitat use of the SP
during the breeding and the wintering season. However the habitats where RK was
observed ap- peared quite different. Maintaining and creating optimal suitable
habitats for SP characterized by sparse low vegetation in the foredunes areas,
and uneven/low-slope beach surfaces, is the proposed conservation scenario to
convert anthropic beach restorations and SP populations into a positive
feedback without impacting other threatened shorebird species.


  We study global asymptotic properties of a continuous time Leslie-Gower food
chain model. We construct a Lyapunov function which enables us to establish
global asymptotic stability of the unique coexisting equilibrium state.


  Starting from a master equation, we derive the evolution equation for the
size distribution of elements in an evolving system, where each element can
grow, divide into two, and produce new elements. We then probe general
solutions of the evolution quation, to obtain such skew distributions as
power-law, log-normal, and Weibull distributions, depending on the growth or
division and production. Specifically, repeated production of elements of
uniform size leads to power-law distributions, whereas production of elements
with the size distributed according to the current distribution as well as no
production of new elements results in log-normal distributions. Finally,
division into two, or binary fission, bears Weibull distributions. Numerical
simulations are also carried out, confirming the validity of the obtained
solutions.


  Random copying is a simple model for population dynamics in the absence of
selection, and has been applied to both biological and cultural evolution. In
this work, we investigate the effect that spatial structure has on the
dynamics. We focus in particular on how a measure of the diversity in the
population changes over time. We show that even when the vast majority of a
population's history may be well-described by a spatially-unstructured model,
spatial structure may nevertheless affect the expected level of diversity seen
at a local scale. We demonstrate this phenomenon explicitly by examining the
random copying process on small-world networks, and use our results to comment
on the use of simple random-copying models in an empirical context.


  We show how to analytically derive the average sequence dissimilarity (ASD)
within and between species under a simplified multi-species coalescent setup.


  We report on the direct measurement of the correlation times of the protein
backbone carbons and proximal waters of hydration in mechanically strained
elastin by nuclear magnetic resonance methods. The experimental data indicate a
decrease in the correlation times of the carbonyl carbons as the strain on the
biopolymer is increased. These observations are in good agreement with short
4ns molecular dynamics simulations of (VPGVG)3, a well studied mimetic peptide
of elastin. The experimental results also indicate a reduction in the
correlation time of proximal waters of hydration with increasing strain applied
to the elastomer. A simple model is suggested that correlates the increase in
the motion of proximal waters of hydration to the increase in frequency of
libration of the protein backbone that develops with increasing strain.
Together, the reduction in the protein entropy accompanied with the increase in
entropy of the proximal waters of hydration with increasing strain, support the
notion that the source of elasticity is driven by an entropic mechanism arising
from the change in entropy of the protein backbone.


  Problems of search and recognition appear over different scales in biological
systems. In this review we focus on the challenges posed by interactions
between proteins, in particular transcription factors, and DNA and possible
mechanisms which allow for a fast and selective target location. Initially we
argue that DNA-binding proteins can be classified, broadly, into three distinct
classes which we illustrate using experimental data. Each class calls for a
different search process and we discuss the possible application of different
search mechanisms proposed over the years to each class. The main thrust of
this review is a new mechanism which is based on barrier discrimination. We
introduce the model and analyze in detail its consequences. It is shown that
this mechanism applies to all classes of transcription factors and can lead to
a fast and specific search. Moreover, it is shown that the mechanism has
interesting transient features which allow for stability at the target despite
rapid binding and unbinding of the transcription factor from the target.


  Non-coding RNA molecules fold into precise base pairing patterns to carry out
critical roles in genetic regulation and protein synthesis. We show here that
coupling systematic mutagenesis with high-throughput SHAPE chemical mapping
enables accurate base pair inference of domains from ribosomal RNA, ribozymes,
and riboswitches. For a six-RNA benchmark that challenged prior
chemical/computational methods, this mutate-and-map strategy gives secondary
structures in agreement with crystallographic data (2 % error rates), including
a blind test on a double-glycine riboswitch. Through modeling of partially
ordered RNA states, the method enables the first test of an 'interdomain
helix-swap' hypothesis for ligand-binding cooperativity in a glycine
riboswitch. Finally, the mutate-and-map data report on tertiary contacts within
non-coding RNAs; coupled with the Rosetta/FARFAR algorithm, these data give
nucleotide-resolution three-dimensional models (5.7 {\AA} helix RMSD) of an
adenine riboswitch. These results highlight the promise of a two-dimensional
chemical strategy for inferring the secondary and tertiary structures that
underlie non-coding RNA behavior.


  Shear stress is an important physical factor that regulates proliferation,
migration and morphogenesis. In particular, the homeostasis of blood vessels is
dependent on shear stress. To mimic this process ex vivo, efforts have been
made to seed scaffolds with vascular and other cell types in the presence of
growth factors and under pulsatile flow conditions. However, the resulting
bioreactors lack information on shear stress and flow distributions within the
scaffold. Consequently, it is difficult to interpret the effects of shear
stress on cell function. Such knowledge would enable researchers to improve
upon cell culture protocols. Recent work has focused on optimizing the
microstructural parameters of the scaffold to fine tune the shear stress. In
this study, we have adopted a different approach whereby flows are redirected
throughout the bioreactor along channels patterned in the porous scaffold to
yield shear stress distributions that are optimized for uniformity centered on
a target value. A topology optimization algorithm coupled to computational
fluid dynamics simulations was devised to this end. The channel topology in the
porous scaffold was varied using a combination of genetic algorithm and fuzzy
logic. The method is validated by experiments using magnetic resonance imaging
(MRI) readouts of the flow field.


  Computational approaches to transcription factor binding site identification
have been actively researched for the past decade.
  Negative examples have long been utilized in de novo motif discovery and have
been shown useful in transcription factor binding site search as well.
  However, understanding of the roles of negative examples in binding site
search is still very limited.
  We propose the 2-centroid and optimal discriminating vector methods, taking
into account negative examples. Cross-validation results on E. coli
transcription factors show that the proposed methods benefit from negative
examples, outperforming the centroid and position-specific scoring matrix
methods. We further show that our proposed methods perform better than a
state-of-the-art method. We characterize the proposed methods in the context of
the other compared methods and show that, coupled with motif subtype
identification, the proposed methods can be effectively applied to a wide range
of transcription factors. Finally, we argue that the proposed methods are
well-suited for eukaryotic transcription factors as well.
  Software tools are available at: http://biogrid.engr.uconn.edu/tfbs_search/.


  Recently it has been discovered---contrary to expectations of physicists as
well as biologists---that the energy transport during photosynthesis, from the
chlorophyll pigment that captures the photon to the reaction centre where
glucose is synthesised from carbon dioxide and water, is highly coherent even
at ambient temperature and in the cellular environment. This process and the
key molecular ingredients that it depends on are described. By looking at the
process from the computer science view-point, we can study what has been
optimised and how. A spatial search algorithmic model based on robust features
of wave dynamics is presented.


  Graph theory is increasingly commonly utilised in genetics, proteomics and
neuroimaging. In such fields, the data of interest generally constitute
weighted graphs. Analysis of such weighted graphs often require the integration
of topological metrics with respect to the density of the graph. Here, density
refers to the proportion of the number of edges present in that graph. When
topological metrics based on shortest paths are of interest, such
density-integration usually necessitates the iterative application of
Dijkstra's algorithm in order to compute the shortest path matrix at each
density level. In this short note, we describe a recursive shortest path
algorithm based on single edge updating, which replaces the need for the
iterative use of Dijkstra's algorithm. Our proposed procedure is based on pairs
of breadth-first searches around each of the vertices incident to the edge
added at each recursion. An algorithmic analysis of the proposed technique is
provided. When the graph of interest is coded as an adjacency list, our
algorithm can be shown to be more efficient than an iterative use of Dijkstra's
algorithm.


  We study intact and bulging Escherichia coli cells using atomic force
microscopy to separate the contributions of the cell wall and turgor pressure
to the overall cell stiffness. We find strong evidence of power-law
stress-stiffening in the E. coli cell wall, with an exponent of 1.22 \pm 0.12,
such that the wall is significantly stiffer in intact cells (E = 23 \pm 8 MPa
and 49 \pm 20 MPa in the axial and circumferential directions) than in
unpressurized sacculi. These measurements also indicate that the turgor
pressure in living cells E. coli is 29 \pm 3 kPa.


  We present neutron scattering measurements on the dynamics of hemoglobin (Hb)
in human red blood cells in vivo. Global and internal Hb dynamics were measured
in the ps to ns time- and {\AA} length-scale using quasielastic neutron
backscattering spectroscopy. We observed the cross-over from global Hb
short-time to long-time self-diffusion. Both short- and long-time diffusion
coefficients agree quantitatively with predicted values from hydrodynamic
theory of non-charged hard-sphere suspensions when a bound water fraction of
around 0.23g H2O/ g Hb is taken into account. The higher amount of water in the
cells facilitates internal protein fluctuations in the ps time-scale when
compared to fully hydrated Hb powder. Slower internal dynamics of Hb in red
blood cells in the ns time-range were found to be rather similar to results
obtained with fully hydrated protein powders, solutions and E. coli cells.


  Farming and herding were introduced to Europe from the Near East and
Anatolia; there are, however, considerable arguments about the mechanisms of
this transition. Were it people who moved and outplaced the indigenous hunter-
gatherer groups or admixed with them? Or was it just material and information
that moved-the Neolithic Package-consisting of domesticated plants and animals
and the knowledge of its use? The latter process is commonly referred to as
cultural diffusion and the former as demic diffusion. Despite continuous and
partly combined efforts by archaeologists, anthropologists, linguists,
paleontologists and geneticists a final resolution of the debate has not yet
been reached. In the present contribution we interpret results from the Global
Land Use and technological Evolution Simulator (GLUES), a mathematical model
for regional sociocultural development embedded in the western Eurasian
geoenvironmental context during the Holocene. We demonstrate that the model is
able to realistically hindcast the expansion speed and the inhomogeneous
space-time evolution of the transition to agropastoralism in Europe. GLUES, in
contrast to models that do not resolve endogenous sociocultural dynamics, also
describes and explains how and why the Neolithic advanced in stages. In the
model analysis, we uncouple the mechanisms of migration and information
exchange. We find that (1) an indigenous form of agropastoralism could well
have arisen in certain Mediterranean landscapes, but not in Northern and
Central Europe, where it depended on imported technology and material, (2) both
demic diffusion by migration or cultural diffusion by trade may explain the
western European transition equally well, (3) [...]


  Background: Acceleration of adaptation dynamics by stress-induced
hypermutation has been found experimentally. Evolved evolvability is a
prominent explanation. We investigate a more generally applicable explanation
by a physical constraint. Methods and Results: A generic thermodynamical
analysis of genetic information storage obviates physical constraints on the
integrity of genetic information. The capability to employ metabolic resources
is found as a major determinant of mutation probability in stored genetic
information. Incorporation into a non-recombinant, asexual adaptation toy model
predicts cases of markedly accelerated adaptation, driven by a transient
increase of mutation rate. No change in the mutation rate as a genetic trait is
required. The mutation rate of one and the same genotype varies dependent on
stress level. Implications: Stress-dependent mutation rates are physically
necessary and challenge a condition-independent genotype to mutation rate
mapping. This holds implications for evolutionary theory and pathogen and
cancer evolution.


  In the past decade computational biology has grown from a cottage industry
with a handful of researchers to an attractive interdisciplinary field,
catching the attention and imagination of many quantitatively-minded
scientists. Of interest to us is the key role played by the EM algorithm during
this transformation. We survey the use of the EM algorithm in a few important
computational biology problems surrounding the "central dogma"; of molecular
biology: from DNA to RNA and then to proteins. Topics of this article include
sequence motif discovery, protein sequence alignment, population genetics,
evolutionary models and mRNA expression microarray data analysis.


  To ensure fast gene activation, Transcription Factors (TF) use a mechanism
known as facilitated diffusion to find their DNA promoter site. Here we analyze
such a process where a TF alternates between 3D and 1D diffusion. In the latter
(TF bound to the DNA), the TF further switches between a fast translocation
state dominated by interaction with the DNA backbone, and a slow examination
state where interaction with DNA base pairs is predominant. We derive a new
formula for the mean search time, and show that it is faster and less sensitive
to the binding energy fluctuations compared to the case of a single sliding
state. We find that for an optimal search, the time spent bound to the DNA is
larger compared to the 3D time in the nucleus, in agreement with recent
experimental data. Our results further suggest that modifying switching via
phosphorylation or methylation of the TF or the DNA can efficiently regulate
transcription.


  We solve numerically a kinetic model of chiral polymerization in systems
closed to matter and energy flow, paying special emphasis to its ability to
amplify the small initial enantiomeric excesses due to the internal and
unavoidable statistical fluctuations. The reaction steps are assumed to be
reversible, implying a thermodynamic constraint among some of the rate
constants. Absolute asymmetric synthesis is achieved in this scheme. The system
can persist for long times in quasi- stationary chiral asymmetric states before
racemizing. Strong inhibition leads to long-period chiral oscillations in the
enantiomeric excesses of the longest homopolymer chains. We also calculate the
entropy production {\sigma} per unit volume and show that {\sigma} increases to
a peak value either before or in the vicinity of the chiral symmetry breaking
transition.


  The lambda phage is a paradigm temperate bacteriophage. Its lysogenic and
lytic life cycles echo competition between the DNA binding CI and CRO proteins.
Here we address the Physics of this transition in terms of an energy function
that portrays the backbone as a multi-soliton configuration. The precision of
the individual solitons far exceeds the B-factor accuracy of the experimentally
determined protein conformations giving us confidence to conclude that three of
the four loops are each composites of two closely located solitons. The only
exception is the repressive DNA binding turn, it is the sole single soliton
configuration of the backbone. When we compare the solitons with the Protein
Data Bank we find that the one preceding the DNA recognition helix is unique to
the CI protein, prompting us to conclude that the lysogenic to lytic transition
is due to a saddle-node bifurcation involving a soliton-antisoliton
annihilation that removes the first loop.


  We inspect the geometry of proteins by identifying their backbones as framed
polygons. We find that the left-handed helix region of the Ramachandran map for
non-glycyl residues corresponds to an isolated and highly localized sector in
the orientation of the $C_\beta$ carbons, when viewed in a Frenet frame that is
centered at the corresponding $C_\alpha$ carbons. We show that this
localization in the orientation persists to $C_\gamma$ and $C_\delta$ carbons.
Furthermore, when we extend our analysis to the neighboring residues we
conclude that the left-handed helix region reflects a very regular and
apparently residue independent collective interplay of at least seven
consecutive amino acids.


  Cooperative mutualism is a major force driving evolution and sustaining
ecosystems. Although the importance of spatial degrees of freedom and number
fluctuations is well-known, their effects on mutualism are not fully
understood. With range expansions of microbes in mind, we show that, even when
mutualism confers a distinct selective advantage, it persists only in
populations with high density and frequent migrations. When these parameters
are reduced, mutualism is generically lost via a directed percolation process,
with a phase diagram strongly influenced by an exceptional DP2 transition.


  The inference of new information on the relatedness of species by
phylogenetic trees based on DNA data is one of the main challenges of modern
biology. But despite all technological advances, DNA sequencing is still a
time-consuming and costly process. Therefore, decision criteria would be
desirable to decide a priori which data might contribute new information to the
supertree which is not explicitly displayed by any input tree. A new concept,
so-called groves, to identify taxon sets with the potential to construct such
informative supertrees was suggested by An\'e et al. in 2009. But the important
conjecture that maximal groves can easily be identified in a database remained
unproved and was published on the Isaac Newton Institute's list of open
phylogenetic problems. In this paper, we show that the conjecture does not
generally hold, but also introduce a new concept, namely 2-overlap groves,
which overcomes this problem.


  Near a critical point, the time scale of thermally-induced fluctuations
diverges in a manner determined by the dynamic universality class. Experiments
have verified predicted 3D dynamic critical exponents in many systems, but
similar experiments in 2D have been lacking for the case of conserved order
parameter. Here we analyze time-dependent correlation functions of a quasi-2D
lipid bilayer in water to show that its critical dynamics agree with a recently
predicted universality class. In particular, the effective dynamic exponent
$z_{\text{eff}}$ crosses over from $\sim 2$ to $\sim 3$ as the correlation
length of fluctuations exceeds a hydrodynamic length set by the membrane and
bulk viscosities.


  Bistability plays a central role in the gene regulatory networks (GRNs)
controlling many essential biological functions, including cellular
differentiation and cell cycle control. However, establishing the network
topologies that can exhibit bistability remains a challenge, in part due to the
exceedingly large variety of GRNs that exist for even a small number of
components. We begin to address this problem by employing chemical reaction
network theory in a comprehensive in silico survey to determine the capacity
for bistability of more than 40,000 simple networks that can be formed by two
transcription factor-coding genes and their associated proteins (assuming only
the most elementary biochemical processes). We find that there exist reaction
rate constants leading to bistability in ~90% of these GRN models, including
several circuits that do not contain any of the TF cooperativity commonly
associated with bistable systems, and the majority of which could only be
identified as bistable through an original subnetwork-based analysis. A
topological sorting of the two-gene family of networks based on the presence or
absence of biochemical reactions reveals eleven minimal bistable networks
(i.e., bistable networks that do not contain within them a smaller bistable
subnetwork). The large number of previously unknown bistable network topologies
suggests that the capacity for switch-like behavior in GRNs arises with
relative ease and is not easily lost through network evolution. To highlight
the relevance of the systematic application of CRNT to bistable network
identification in real biological systems, we integrated publicly available
protein-protein interaction, protein-DNA interaction, and gene expression data
from Saccharomyces cerevisiae, and identified several GRNs predicted to behave
in a bistable fashion.


  Evolutionary biology shares many concepts with statistical physics: both deal
with populations, whether of molecules or organisms, and both seek to simplify
evolution in very many dimensions. Often, methodologies have undergone parallel
and independent development, as with stochastic methods in population genetics.
We discuss aspects of population genetics that have embraced methods from
physics: amongst others, non-equilibrium statistical mechanics, travelling
waves, and Monte-Carlo methods have been used to study polygenic evolution,
rates of adaptation, and range expansions. These applications indicate that
evolutionary biology can further benefit from interactions with other areas of
statistical physics, for example, by following the distribution of paths taken
by a population through time.


  We propose a hybrid dynamical system approach to model the evolution of a
pathogen that experiences different selective pressures according to a
stochastic process. In every environment, the evolution of the pathogen is
described by a version of the Fisher-Haldane-Wright equation while the
switching between environments follows a Markov jump process. We investigate
how the qualitative behavior of a simple single-host deterministic system
changes when the stochastic switching process is added. In particular, we study
the stability in probability of monomorphic equilibria. We prove that in a
"constantly" fluctuating environment, the genotype with the highest mean
fitness is asymptotically stable in probability. However, if the probability of
host switching depends on the genotype composition of the population,
polymorphism can be stably maintained.
  This is a corrected version of the paper that appeared in Mathematical
Biosciences 240 (2012), p. 70-75. A corrigendum has appeared in the same
journal.


  We explore self-organizing strategies for role assignment in a foraging task
carried out by a colony of artificial agents. Our strategies are inspired by
various mechanisms of division of labor (polyethism) observed in eusocial
insects like ants, termites, or bees. Specifically we instantiate models of
caste polyethism and age or temporal polyethism to evaluated the benefits to
foraging in a dynamic environment. Our experiment is directly related to the
exploration/exploitation trade of in machine learning.


  The metabolism is the motor behind the biological complexity of an organism.
One problem of characterizing its large-scale structure is that it is hard to
know what to compare it to. All chemical reaction systems are shaped by the
same physics that gives molecules their stability and affinity to react. These
fundamental factors cannot be captured by standard null-models based on
randomization. The unique property of organismal metabolism is that it is
controlled, to some extent, by an enzymatic machinery that is subject to
evolution. In this paper, we explore the possibility that reaction systems of
planetary atmospheres can serve as a null-model against which we can define
metabolic structure and trace the influence of evolution. We find that the two
types of data can be distinguished by their respective degree distributions.
This is especially clear when looking at the degree distribution of the
reaction network (of reaction connected to each other if they involve the same
molecular species). For the Earth's atmospheric network and the human metabolic
network, we look into more detail for an underlying explanation of this
deviation. However, we cannot pinpoint a single cause of the difference, rather
there are several concurrent factors. By examining quantities relating to the
modular-functional organization of the metabolism, we confirm that metabolic
networks have a more complex modular organization than the atmospheric
networks, but not much more. We interpret the more variegated modular
arrangement of metabolism as a trace of evolved functionality. On the other
hand, it is quite remarkable how similar the structures of these two types of
networks are, which emphasizes that the constraints from the chemical
properties of the molecules has a larger influence in shaping the reaction
system than does natural selection.


  Polymers can be modeled as open polygonal paths and their closure generates
knots. Knotted proteins detection is currently achieved via high-throughput
methods based on a common framework insensitive to the handedness of knots.
Here we propose a topological framework for the computation of the HOMFLY
polynomial, an handedness-sensitive invariant. Our approach couples a
multi-component reduction scheme with the polynomial computation. After
validation on tabulated knots and links the framework was applied to the entire
Protein Data Bank along with a set of selected topological checks that allowed
to discard artificially entangled structures. This led to an up-to-date table
of knotted proteins that also includes two newly detected right-handed trefoil
knots in recently deposited protein structures. The application range of our
framework is not limited to proteins and it can be extended to the topological
analysis of biological and synthetic polymers and more generally to arbitrary
polygonal paths.


  Large-scale dynamical properties of complete chromosome DNA sequences of
eukaryotes are considered. By the proposed deterministic models with
intermittency and symbolic dynamics we describe a wide spectrum of large-scale
patterns inherent in these sequences, such as segmental duplications, tandem
repeats, and other complex sequence structures. It is shown that the recently
discovered gene number balance on the strands is not of random nature, and a
complete chromosome DNA sequence exhibits the properties of deterministic
chaos.


  Several groups have recently modeled evolutionary transitions from an
ancestral allele to a beneficial allele separated by one or more intervening
mutants. The beneficial allele can become fixed if a succession of intermediate
mutants are fixed or alternatively if successive mutants arise while the
previous intermediate mutant is still segregating. This latter process has been
termed stochastic tunneling. Previous work has focused on the Moran model of
population genetics. I use elementary methods of analyzing stochastic processes
to derive the probability of tunneling in the limit of large population size
for both Moran and Wright-Fisher populations. I also show how to efficiently
obtain numerical results for finite populations. These results show that the
probability of stochastic tunneling is twice as large under the Wright-Fisher
model as it is under the Moran model.


  One popular assumption regarding biological systems is that traits have
evolved to be optimized with respect to function. This is a standard goal in
evolutionary computation, and while not always embraced in the biological
sciences, is an underlying assumption of what happens when fitness is
maximized. The implication of this is that a signaling pathway or phylogeny
should show evidence of minimizing the number of steps required to produce a
biochemical product or phenotypic adaptation. In this paper, it will be shown
that a principle of "maximum intermediate steps" may also characterize complex
biological systems, especially those in which extreme historical contingency or
a combination of mutation and recombination are key features. The contribution
to existing literature is two-fold: demonstrating both the potential for
non-optimality in engineered systems with "lifelike" attributes, and the
underpinnings of non-optimality in naturalistic contexts.
  This will be demonstrated by using the Rube Goldberg Machine (RGM) analogy.
Mechanical RGMs will be introduced, and their relationship to conceptual
biological RGMs. Exemplars of these biological RGMs and their evolution (e.g.
introduction of mutations and recombination-like inversions) will be
demonstrated using block diagrams and interconnections with complex networks
(called convolution architectures). The conceptual biological RGM will then be
mapped to an artificial vascular system, which can be modeled using
microfluidic-like structures. Theoretical expectations will be presented,
particularly regarding whether or not maximum intermediate steps equates to the
rescue or reuse of traits compromised by previous mutations or inversions.
Considerations for future work and applications will then be discussed,
including the incorporation of such convolution architectures into complex
networks.


  A statistically principled way of conducting weighted network analysis is
still lacking. Comparison of different populations of weighted networks is hard
because topology is inherently dependent on wiring cost, where cost is defined
as the number of edges in an unweighted graph. In this paper, we evaluate the
benefits and limitations associated with using cost-integrated topological
metrics. Our focus is on comparing populations of weighted undirected graphs
using global efficiency. We evaluate different approaches to the comparison of
weighted networks that differ in mean association weight. Our key result shows
that integrating over cost is equivalent to controlling for any monotonic
transformation of the weight set of a weighted graph. That is, when integrating
over cost, we eliminate the differences in topology that may be due to a
monotonic transformation of the weight set. Our result holds for any unweighted
topological measure. Cost-integration is therefore helpful in disentangling
differences in cost from differences in topology. By contrast, we show that the
use of the weighted version of a topological metric does not constitute a valid
approach to this problem. Indeed, we prove that, under mild conditions, the use
of the weighted version of global efficiency is equivalent to simply comparing
weighted costs. Thus, we recommend the reporting of (i) differences in weighted
costs and (ii) differences in cost-integrated topological measures. We
demonstrate the application of these techniques in a re-analysis of an fMRI
working memory task. Finally, we discuss the limitations of integrating
topology over cost, which may pose problems when some weights are zero, when
multiplicities exist in the ranks of the weights, and when one expects subtle
cost-dependent topological differences, which could be masked by
cost-integration.


  A prey animal surveying its environment must decide whether there is a
dangerous predator present or not. If there is, it may flee. Flight has an
associated cost, so the animal should not flee if there is no danger. However,
the prey animal cannot know the state of its environment with certainty, and is
thus bound to make some errors. We formulate a probabilistic automaton model of
a prey animal's life and use it to compute the optimal escape decision
strategy, subject to the animal's uncertainty. The uncertainty is a major
factor in determining the decision strategy: only in the presence of
uncertainty do economic factors (like mating opportunities lost due to flight)
influence the decision. We performed computer simulations and found that
\emph{in silico} populations of animals subject to predation evolve to display
the strategies predicted by our model, confirming our choice of objective
function for our analytic calculations. To the best of our knowledge, this is
the first theoretical study of escape decisions to incorporate the effects of
uncertainty, and to demonstrate the correctness of the objective function used
in the model.


  We describe and solve a discrete persistent chain model of protein binding on
DNA, involving an extra sigma_i s at a site i of the DNA. This variable takes
the value 1 or 0 depending on whether the site is occupied by a protein or not.
In addition, if the site is occupied by a protein, there is an extra energy
cost epsilon. For small force, we obtain analytic expressions for the
force-extension curve and the fraction of bound protein on the DNA.For higher
forces, the model can be solved numerically to obtain force extension curves
and the average fraction of bound proteins as a function of applied force. Our
model can be used to analyze experimental force extension curves of protein
binding on DNA, and hence deduce the number of bound proteins in the case of
non-specific binding.


  RNA-Seq is rapidly becoming the standard technology for transcriptome
analysis. Fundamental to many of the applications of RNA-Seq is the
quantification problem, which is the accurate measurement of relative
transcript abundances from the sequenced reads. We focus on this problem, and
review many recently published models that are used to estimate the relative
abundances. In addition to describing the models and the different approaches
to inference, we also explain how methods are related to each other. A key
result is that we show how inference with many of the models results in
identical estimates of relative abundances, even though model formulations can
be very different. In fact, we are able to show how a single general model
captures many of the elements of previously published methods. We also review
the applications of RNA-Seq models to differential analysis, and explain why
accurate relative transcript abundance estimates are crucial for downstream
analyses.


  Two distinct mechanisms for filtering noise in an input signal are identified
in a class of adaptive sensory networks. We find that the high frequency noise
is filtered by the output degradation process through time-averaging; while the
low frequency noise is damped by adaptation through negative feedback. Both
filtering processes themselves introduce intrinsic noises, which are found to
be unfiltered and can thus amount to a significant internal noise floor even
without signaling. These results are applied to E. coli chemotaxis. We show
unambiguously that the molecular mechanism for the Berg-Purcell time-averaging
scheme is the dephosphorylation of the response regulator CheY-P, not the
receptor adaptation process as previously suggested. The high frequency noise
due to the stochastic ligand binding-unbinding events and the random ligand
molecule diffusion is averaged by the CheY-P dephosphorylation process to a
negligible level in E.coli. We identify a previously unstudied noise source
caused by the random motion of the cell in a ligand gradient. We show that this
random walk induced signal noise has a divergent low frequency component, which
is only rendered finite by the receptor adaptation process. For gradients
within the E. coli sensing range, this dominant external noise can be
comparable to the significant intrinsic noise in the system. The dependence of
the response and its fluctuations on the key time scales of the system are
studied systematically. We show that the chemotaxis pathway may have evolved to
optimize gradient sensing, strong response, and noise control in different time
scales


  We study the propagation of an SIR (susceptible-infectious-recovered) disease
over an agent population which, at any instant, is fully divided into couples
of agents. Couples are occasionally allowed to exchange their members. This
process of couple recombination can compensate the instantaneous disconnection
of the interaction pattern and thus allow for the propagation of the infection.
We study the incidence of the disease as a function of its infectivity and of
the recombination rate of couples, thus characterizing the interplay between
the epidemic dynamics and the evolution of the population's interaction
pattern.


  We analyse base-pair breathing in a DNA sequence of 12 base-pairs with a
defective base at its centre. We use both all-atom molecular dynamics (MD)
simulations and a system of stochastic differential equations (SDE). In both
cases, Fourier analysis of the trajectories reveals self-organised critical
behaviour in the breathing of base-pairs. The Fourier Transforms (FT) of the
interbase distances show power-law behaviour with gradients close to -1. The
scale-invariant behaviour we have found provides evidence for the view that
base-pair breathing corresponds to the nucleation stage of large-scale DNA
opening (or 'melting') and that this process is a (second-order) phase
transition. Although the random forces in our SDE system were introduced as
white noise, FTs of the displacements exhibit pink noise, as do the
displacements in the AMBER/MD simulations.


  Many cellular functions are based on the rhythmic organization of biological
processes into self-repeating cascades of events. Some of these periodic
processes, such as the cell cycles of several species, exhibit conspicuous
irregularities in the form of period skippings, which lead to polymodal
distributions of cycle lengths. A recently proposed mechanism that accounts for
this quantized behavior is the stabilization of a Hopf-unstable state by
molecular noise. Here we investigate the effect of varying noise in a model
system, namely an excitable activator-repressor genetic circuit, that displays
this noise-induced stabilization effect. Our results show that an optimal noise
level enhances the regularity (coherence) of the cycles, in a form of coherence
resonance. Similar noise levels also optimize the multimodal nature of the
cycle lengths. Together, these results illustrate how molecular noise within a
minimal gene regulatory motif confers robust generation of polymodal patterns
of periodicity.


  The World Summit on Sustainable Development (Johannesburg, 2002) encouraged
the application of the ecosystem approach by 2010. However, at the same Summit,
the signatory States undertook to restore and exploit their stocks at maximum
sustainable yield (MSY), a concept and practice without ecosystemic dimension,
since MSY is computed species by species, on the basis of a monospecific model.
Acknowledging this gap, we propose a definition of "ecosystem viable yields"
(EVY) as yields compatible i) with guaranteed biological safety levels for all
time and ii) with an ecosystem dynamics. To the difference of MSY, this notion
is not based on equilibrium, but on viability theory, which offers advantages
for robustness. For a generic class of multispecies models with harvesting, we
provide explicit expressions for the EVY. We apply our approach to the
anchovy--hake couple in the Peruvian upwelling ecosystem.


  The paradox of enrichment was observed by M. Rosenzweig in a class of
predator-prey models. Two of the parameters in the models are crucial for the
paradox. These two parameters are the prey's carrying capacity and prey's
half-saturation for predation. Intuitively, increasing the carrying capacity
due to enrichment of the prey's environment should lead to a more stable
predator-prey system. Analytically, it turns out that increasing the carrying
capacity always leads to an unstable predator-prey system that is susceptible
to extinction from environmental random perturbations. This is the so-called
paradox of enrichment. Our resolution here rests upon a closer investigation on
a dimensionless number $H$ formed from the carrying capacity and the prey's
half-saturation. By recasting the models into dimensionless forms, the models
are in fact governed by a few dimensionless numbers including $H$. The effects
of the two parameters: carrying capacity and half-saturation are incorporated
into the number $H$. In fact, increasing the carrying capacity is equivalent
(i.e. has the same effect on $H$) to decreasing the half-saturation which
implies more aggressive predation. Since there is no paradox between more
aggressive predation and instability of the predator-prey system, the paradox
of enrichment is resolved.
  We also further explore spatially dependent models for which the phase space
is infinite dimensional. The spatially independent limit cycle which is
generated by a Hopf bifurcation from an unstable steady state, is linearly
stable in the infinite dimensional phase space. Numerical simulations indicate
that the basin of attraction of the limit cycle is riddled. This shows that
spatial perturbations can sometimes (neither always nor never) remove the
paradox of enrichment near the limit cycle!


  The idea that the distance among pairs of languages can be evaluated from
lexical differences seems to have its roots in the work of the French explorer
Dumont D'Urville. He collected comparative words lists of various languages
during his voyages aboard the Astrolabe from 1826 to 1829 and, in his work
about the geographical division of the Pacific, he proposed a method to measure
the degree of relation between languages.
  The method used by the modern lexicostatistics, developed by Morris Swadesh
in the 1950s, measures distances from the percentage of shared cognates, which
are words with a common historical origin. The weak point of this method is
that subjective judgment plays a relevant role.
  Recently, we have proposed a new automated method which is motivated by the
analogy with genetics. The new approach avoids any subjectivity and results can
be easily replicated by other scholars. The distance between two languages is
defined by considering a renormalized Levenshtein distance between pair of
words with the same meaning and averaging on the words contained in a list. The
renormalization, which takes into account the length of the words, plays a
crucial role, and no sensible results can be found without it.
  In this paper we give a short review of our automated method and we
illustrate it by considering the cluster of Malagasy dialects. We show that it
sheds new light on their kinship relation and also that it furnishes a lot of
new information concerning the modalities of the settlement of Madagascar.


  Gene regulatory networks arise in all living cells, allowing the control of
gene expression patterns. The study of their topology has revealed that certain
subgraphs of interactions or "motifs" appear at anomalously high frequencies.
We ask here whether this phenomenon may emerge because of the functions carried
out by these networks. Given a framework for describing regulatory interactions
and dynamics, we consider in the space of all regulatory networks those that
have a prescribed function. Monte Carlo sampling is then used to determine how
these functional networks lead to specific motif statistics in the
interactions. In the case where the regulatory networks are constrained to
exhibit multi-stability, we find a high frequency of gene pairs that are
mutually inhibitory and self-activating. In contrast, networks constrained to
have periodic gene expression patterns (mimicking for instance the cell cycle)
have a high frequency of bifan-like motifs involving four genes with at least
one activating and one inhibitory interaction.


  Understanding the emergence and sustainability of cooperation is a
fundamental problem in evolutionary biology and is frequently studied by the
framework of evolutionary game theory. A very powerful mechanism to promote
cooperation is network reciprocity, where the interaction patterns and
opportunities for strategy spread of agents are constrained to limited sets of
permanent interactions partners. Cooperation survives because it is possible
for close-knit communities of cooperation to be shielded from invasion by
defectors. Here we show that parameter ranges in which cooperation can survive
are strongly expanded if game play on networks is skewed towards more frequent
interactions with more successful neighbours. In particular, if agents
exclusively select neighbors for game play that are more successful than
themselves, cooperation can even dominate in situations in which it would die
out if interaction neighbours were chosen without a bias or with a preference
for less successful opponents. We demonstrate that the "selecting fitter
neighbours" strategy is evolutionarily stable. Moreover, it will emerge as the
dominant strategy out of an initially random population of agents.


  Understanding the origin of infectious diseases provides scientifically based
rationales for implementing public health measures that may help to avoid or
mitigate future epidemics. The recent ancestors of a pandemic virus provide
invaluable information about the set of minimal genomic alterations that
transformed a zoonotic agent into a full human pandemic. Since the first
confirmed cases of the H1N1 pandemic virus in the spring of 2009, several
hypotheses about the strain's origins have been proposed. However, how, where,
and when it first infected humans is still far from clear. The only way to
piece together this epidemiological puzzle relies on the collective effort of
the international scientific community to increase genomic sequencing of
influenza isolates, especially ones collected in the months prior to the origin
of the pandemic.


  Prion diseases {\it (e.g. Creutzfeldt-Jakob disease (CJD), variant CJD
(vCJD), Gerstmann-Str$\ddot{\text{a}}$ussler-Scheinker syndrome (GSS), Fatal
Familial Insomnia (FFI) and Kuru in humans, scrapie in sheep, bovine spongiform
encephalopathy (BSE or `mad-cow' disease) and chronic wasting disease (CWD) in
cattles)} are invariably fatal and highly infectious neurodegenerative diseases
affecting humans and animals. However, by now there have not been some
effective therapeutic approaches or medications to treat all these prion
diseases. Rabbits, dogs, and horses are the only mammalian species reported to
be resistant to infection from prion diseases isolated from other species.
Recently, the $\beta$2-$\alpha$2 loop has been reported to contribute to their
protein structural stabilities. The author has found that rabbit prion protein
has a strong salt bridge ASP177-ARG163 (like a taut bow string) keeping this
loop linked. This paper confirms that this salt bridge also contributes to the
structural stability of horse prion protein. Thus, the region of
$\beta$2-$\alpha$2 loop might be a potential drug target region. Besides this
very important salt bridge, other four important salt bridges
GLU196-ARG156-HIS187, ARG156-ASP202 and GLU211-HIS177 are also found to greatly
contribute to the structural stability of horse prion protein. Rich databases
of salt bridges, hydrogen bonds and hydrophobic contacts for horse prion
protein can be found in this paper.


  We present TripNet, a method for constructing phylogenetic networks from
triplets. We will present the motivations behind our approach and its
theoretical and empirical justification. To demonstrate the accuracy and
efficiency of TripNet, we performed two simulations and also applied the method
to five published data sets: Kreitman's data, a set of triplets from real yeast
data obtained from the Fungal Biodiversity Center in Utrecht, a collection of
110 highly recombinant Salmonella multi-locus sequence typing sequences, and
nrDNA ITS and cpDNA JSA sequence data of New Zealand alpine buttercups of
Ranunculus sect. Pseudadonis. Finally, we compare our results with those
already obtained by other authors using alternative methods. TripNet, data
sets, and supplementary files are freely available for download at
(www.bioinf.cs.ipm.ir/softwares/tripnet).


  In this paper we use detailed data about the biology of the head louse
(pediculus humanus capitis) to build a model of the evolution of head lice
colonies. Using theory and computer simulations, we show that the model can be
used to assess the impact of the various strategies usually applied to
eradicate head lice, both conscious (treatments) and unconscious (grooming). In
the case of treatments, we study the difference in performance that arises when
they are applied in systematic and non-systematic ways. Using some reasonable
simplifying assumptions (as random mixing of human groups and the same mobility
for all life stages of head lice other than eggs) we model the contagion of
pediculosis using only one additional parameter. It is shown that this
parameter can be tuned to obtain collective infestations whose variables are
compatible with what is given in the literature on real infestations. We
analyze two scenarios: one where group members begin treatment when a similar
number of lice are present in each head, and another where there is one
individual who starts treatment with a much larger threshold ('superspreader').
For both cases we assess the impact of several collective strategies of
treatment.


  The individual-based models constitute a set of widely implemented tools to
analyze the incidence of individuals heterogeneities in the spread of an
infectious disease. In this work we focus our attention on human contacts
heterogeneities through two of the main individual-based models: mobile agents
and complex networks models. We introduce a novel mobile agents model in which
individuals make displacements with sizes according to a truncated power-law
distribution based on empirical evidence about human mobility. Besides, we
present a procedure to obtain an equivalent weighted contact network from the
previous mobile agents model, where the weights of the links are interpreted as
contact probabilities. From the topological analysis of the equivalent contact
networks we show that small world characteristics are related with truncated
power-law distribution for agent displacements. Finally, we show the
equivalence between both approaches through some numerical experiments for the
spread of an infectious disease.


  This paper is concerned with the dynamical properties of deterministically
modeled chemical reaction systems with mass-action kinetics. Such models are
ubiquitously found in chemistry, population biology, and the burgeoning field
of systems biology. A basic question, whose answer remains largely unknown, is
the following: for which network structures do trajectories of mass-action
systems remain bounded in time? In this paper, we conjecture that the result
holds when the reaction network is weakly reversible, and prove this conjecture
in the case when the reaction network consists of a single linkage class, or
connected component.


  Hard combinatorial optimization problems deal with the search for the minimum
cost solutions (ground states) of discrete systems under strong constraints. A
transformation of state variables may enhance computational tractability. It
has been argued that these state encodings are to be chosen invertible to
retain the original size of the state space. Here we show how redundant
non-invertible encodings enhance optimization by enriching the density of
low-energy states. In addition, smooth landscapes may be established on encoded
state spaces to guide local search dynamics towards the ground state.


  Suppose a finite set $X$ is repeatedly transformed by a sequence of
permutations of a certain type acting on an initial element $x$ to produce a
final state $y$. We investigate how 'different' the resulting state $y'$ to $y$
can be if a slight change is made to the sequence, either by deleting one
permutation, or replacing it with another. Here the 'difference' between $y$
and $y'$ might be measured by the minimum number of permutations of the
permitted type required to transform $y$ to $y'$, or by some other metric. We
discuss this first in the general setting of sensitivity to perturbation of
walks in Cayley graphs of groups with a specified set of generators. We then
investigate some permutation groups and generators arising in computational
genomics, and the statistical implications of the findings.