  Berg and Purcell [Biophys. J. 20, 193 (1977)] calculated how the accuracy of
concentration sensing by single-celled organisms is limited by noise from the
small number of counted molecules. Here we generalize their results to the
sensing of concentration ramps, which is often the biologically relevant
situation (e.g. during bacterial chemotaxis). We calculate lower bounds on the
uncertainty of ramp sensing by three measurement devices: a single receptor, an
absorbing sphere, and a monitoring sphere. We contrast two strategies, simple
linear regression of the input signal versus maximum likelihood estimation, and
show that the latter can be twice as accurate as the former. Finally, we
consider biological implementations of these two strategies, and identify
possible signatures that maximum likelihood estimation is implemented by real
biological systems.


  Capsids of many viruses assemble around nucleic acids or other polymers.
Understanding how the properties of the packaged polymer affect the assembly
process could promote biomedical efforts to prevent viral assembly or
nanomaterials applications that exploit assembly. To this end, we simulate on a
lattice the dynamical assembly of closed, hollow shells composed of several
hundred to 1000 subunits, around a flexible polymer. We find that assembly is
most efficient at an optimum polymer length that scales with the surface area
of the capsid; significantly longer than optimal polymers often lead to
partial-capsids with unpackaged polymer `tails' or a competition between
multiple partial-capsids attached to a single polymer. These predictions can be
tested with bulk experiments in which capsid proteins assemble around
homopolymeric RNA or synthetic polyelectrolytes. We also find that the polymer
can increase the net rate of subunit accretion to a growing capsid both by
stabilizing the addition of new subunits and by enhancing the incoming flux of
subunits; the effects of these processes may be distinguishable with
experiments that monitor the assembly of individual capsids.


  We present the extention and application of a new unsupervised statistical
learning technique--the Partition Decoupling Method--to gene expression data.
Because it has the ability to reveal non-linear and non-convex geometries
present in the data, the PDM is an improvement over typical gene expression
analysis algorithms, permitting a multi-gene analysis that can reveal
phenotypic differences even when the individual genes do not exhibit
differential expression. Here, we apply the PDM to publicly-available gene
expression data sets, and demonstrate that we are able to identify cell types
and treatments with higher accuracy than is obtained through other approaches.
By applying it in a pathway-by-pathway fashion, we demonstrate how the PDM may
be used to find sets of mechanistically-related genes that discriminate
phenotypes.


  The origin of altruistic behavior, i.e. the behavior that is useful for a
population or a species but goes at the expense of an altruistic individual,
has long been a challenge for students of evolutionary biology. The populations
with altruistic individuals thrive better than those without altruists;
however, the altruists within a population thrive worse than the non-altruists
and their prevalence in the population decreases due to individual selection.
Under certain conditions, the strength of group selection, i.e. the competition
between populations, can surpass the strength of individual selection; however,
such conditions are rarely achieved in practice. It was suggested recently that
chances for altruistic behavior to spread highly increase when it is controlled
not by a single gene but by multiple independent genes substitutable in their
effects on the phenotype of the individual. Here we confirm the original verbal
model published as part of the frozen plasticity theory by numerical modeling
of the spread of altruistic/selfish alleles in a metapopulation consisting of
partly isolated groups of organisms (demes) interconnected by migration. We
have shown that altruistic behavior coded by multiple substitutable genes can
stably coexist with selfish behavior, even under relatively high mutation and
migration rates, i.e. under such conditions where altruistic behavior coded by
a single gene is quickly outcompeted in a metapopulation.


  Cells must continuously sense and respond to time-varying environmental
stimuli. These signals are transmitted and processed by biochemical signalling
networks. However, the biochemical reactions making up these networks are
intrinsically noisy, which limits the reliability of intracellular signalling.
Here we use information theory to characterise the reliability of transmission
of time-varying signals through elementary biochemical reactions in the
presence of noise. We calculate the mutual information for both instantaneous
measurements and trajectories of biochemical systems for a Gaussian model. Our
results indicate that the same network can have radically different
characteristics for the transmission of instantaneous signals and trajectories.
For trajectories, the ability of a network to respond to changes in the input
signal is determined by the timing of reaction events, and is independent of
the correlation time of the output of the network. We also study how reliably
signals on different time-scales can be transmitted by considering the
frequency-dependent coherence and gain-to-noise ratio. We find that a detector
that does not consume the ligand molecule upon detection can more reliably
transmit slowly varying signals, while an absorbing detector can more reliably
transmit rapidly varying signals. Furthermore, we find that while one reaction
may more reliably transmit information than another when considered in
isolation, when placed within a signalling cascade the relative performance of
the two reactions can be reversed. This means that optimising signal
transmission at a single level of a signalling cascade can reduce signalling
performance for the cascade as a whole.


  Humans and chimpanzees are believed to have shared a common ancestor about 6
million years ago. Here using a new distance measure called the Jump distance,
we calculate the number of base substitutions that might have occurred in the
mitochondrial DNA during these 6 million years.


  The zero point kinetic energy of protons in water is large on the scale of
chemical interaction energies(29 Kj/mol in bulk room temperature water). Its
value depends upon the structure of the hydrogen bond network, and can change
as the network is confined or as water interacts with surfaces. These changes
have been observed to be large on a chemical scale for water confined in carbon
nanotubes and in the pores of xerogel, and may play a fundamental, and
neglected, role in biological processes involving confined water. We measure
the average momentum distribution of the protons in salmon Na-DNA using Deep
Inelastic Neutron Scattering, for a weakly hydrated (6w/bp) and a dehydrated
fiber sample. This permits the determination of the change in total kinetic
energy of the system per water molecule removed from the DNA and placed in the
bulk liquid. This energy is equal, within errors, to the measured enthalpy for
the same process, demonstrating that changes in the zero point motion of the
protons, arising from changes in structure as water molecules are incorporated
in the DNA, are a significant factor in the energetics of the transition from
the A to B phase with hydration, in this case, providing the entire binding
energy of the water molecules to the DNA. The shape of the momentum
distribution in the dehydrated phase is consistent with coherent delocalization
of some of the protons in a double well potential, with a separation of the
wells of .2 Angstroms.


  The development of cancer is largely driven by the gain or loss of subsets of
the genome, promoting uncontrolled growth or disabling defenses against it.
Identifying genomic regions whose DNA copy number deviates from the normal is
therefore central to understanding cancer evolution. Array-based comparative
genomic hybridization (aCGH) is a high-throughput technique for identifying DNA
gain or loss by quantifying total amounts of DNA matching defined probes
relative to healthy diploid control samples. Due to the high level of noise in
microarray data, however, interpretation of aCGH output is a difficult and
error-prone task.
  In this work, we tackle the computational task of inferring the DNA copy
number per genomic position from noisy aCGH data. We propose CGHTRIMMER, a
novel segmentation method that uses a fast dynamic programming algorithm to
solve for a least-squares objective function for copy number assignment.
CGHTRIMMER consistently achieves superior precision and recall to leading
competitors on benchmarks of synthetic data and real data from the Coriell cell
lines. In addition, it finds several novel markers not recorded in the
benchmarks but plausibly supported in the oncology literature. Furthermore,
CGHTRIMMER achieves superior results with run-times from 1 to 3 orders of
magnitude faster than its state-of-art competitors.
  CGHTRIMMER provides a new alternative for the problem of aCGH discretization
that provides superior detection of fine-scale regions of gain or loss yet is
fast enough to process very large data sets in seconds. It thus meets an
important need for methods capable of handling the vast amounts of data being
accumulated in high-throughput studies of tumor genetics.


  Protein-protein interactions comprise both transport and reaction steps.
During the transport step, anisotropy of proteins and their complexes is
important both for hydrodynamic diffusion and accessibility of the binding
site. Using a Brownian dynamics approach and extensive computer simulations, we
quantify the effect of anisotropy on the encounter rate of ellipsoidal
particles covered with spherical encounter patches. We show that the encounter
rate $k$ depends on the aspect ratios $\xi$ mainly through steric effects,
while anisotropic diffusion has only a little effect. Calculating analytically
the crossover times from anisotropic to isotropic diffusion in three
dimensions, we find that they are much smaller than typical protein encounter
times, in agreement with our numerical results.


  We point out that a simple and generic strategy to lower the risk for
extinction consists in the developing a dormant stage in which the organism is
unable to multiply but may die. The dormant organism is protected against the
poisonous environment. The result is to increase the survival probability of
the entire population by introducing a type of zero reproductive fitness. This
is possible, because the reservoir of dormant individuals act as a buffer that
can cushion fatal fluctuations in the number of births and deaths which without
the dormant population would have driven the entire population to extinction.


  In many biological processes heterogeneity within cell populations is an
important issue. In this work we consider populations where the behavior of
every single cell can be described by a system of ordinary differential
equations. Heterogeneity among individual cells is accounted for by differences
in parameter values and initial conditions. Hereby, parameter values and
initial conditions are subject to a distribution function which is part of the
model specification. Based on the single cell model and the considered
parameter distribution, a partial differential equation model describing the
distribution of cells in the state and in the output space is derived.
  For the estimation of the parameter distribution within the model, we
consider experimental data as obtained from flow cytometric analysis. From
these noise-corrupted data a density-based statistical data model is derived.
Using this data model the parameter distribution within the cell population is
computed using convex optimization techniques.
  To evaluate the proposed method, a model for the caspase activation cascade
is considered. It is shown that for known noise properties the unknown
parameter distributions in this model are well estimated by the proposed
method.


  Circadian clocks are oscillatory genetic networks that help organisms adapt
to the 24-hour day/night cycle. The clock of the green alga Ostreococcus tauri
is the simplest plant clock discovered so far. Its many advantages as an
experimental system facilitate the testing of computational predictions.
  We present a model of the Ostreococcus clock in the stochastic process
algebra Bio-PEPA and exploit its mapping to different analysis techniques, such
as ordinary differential equations, stochastic simulation algorithms and
model-checking. The small number of molecules reported for this system tests
the limits of the continuous approximation underlying differential equations.
We investigate the difference between continuous-deterministic and
discrete-stochastic approaches. Stochastic simulation and model-checking allow
us to formulate new hypotheses on the system behaviour, such as the presence of
self-sustained oscillations in single cells under constant light conditions.
  We investigate how to model the timing of dawn and dusk in the context of
model-checking, which we use to compute how the probability distributions of
key biochemical species change over time. These show that the relative
variation in expression level is smallest at the time of peak expression,
making peak time an optimal experimental phase marker. Building on these
analyses, we use approaches from evolutionary systems biology to investigate
how changes in the rate of mRNA degradation impacts the phase of a key protein
likely to affect fitness. We explore how robust this circadian clock is towards
such potential mutational changes in its underlying biochemistry. Our work
shows that multiple approaches lead to a more complete understanding of the
clock.


  A central issue in ecology today is that of the factors determining the
relative abundance of species within a natural community. The proper
application of the principles of statistical physics to the problem of species
abundance distributions (SADs) has enabled us to identify the fundamental
ecological principles responsible for the near universal features observed.
These principles are (i) a limit on the number of individuals in an ecological
guild and (ii) per capita birth and death rates. We thus unify the neutral
theory of Hubbell [1], the master equation approach of Volkov et al [2, 3] and
a recent attempt at deploying statistical mechanics to deal with this problem,
the idiosyncratic [extreme niche] theory of Pueyo et al [4]; we have identified
the true origin of the prior in [4]. We also find in our results clear
indications that niches must be very flexible and that temporal fluctuations on
all sorts of scales are of considerable importance in community structure.


  In this paper, we introduce an approach to the protein folding problem from
the point of view of statistical physics. Protein folding is a stochastic
process by which a polypeptide folds into its characteristic and functional 3D
structure from random coil. The process involves an intricate interplay between
global geometry and local structure, and each protein seems to present special
problems. We introduce CSAW (conditioned self-avoiding walk), a model of
protein folding that combines the features of self-avoiding walk (SAW) and the
Monte Carlo method. In this model, the unfolded protein chain is treated as a
random coil described by SAW. Folding is induced by hydrophobic forces and
other interactions, such as hydrogen bonding, which can be taken into account
by imposing conditions on SAW. Conceptually, the mathematical basis is a
generalized Langevin equation. To illustrate the flexibility and capabilities
of the model, we consider several examples, including helix formation, elastic
properties, and the transition in the folding of myoglobin. From the CSAW
simulation and physical arguments, we find a universal elastic energy for
proteins, which depends only on the radius of gyration $R_{g}$ and the residue
number $N$. The elastic energy gives rise to scaling laws $R_{g}\sim N^{\nu}$
in different regions with exponents $\nu =3/5,3/7,2/5$, consistent with the
observed unfolded stage, pre-globule, and molten globule, respectively. These
results indicate that CSAW can serve as a theoretical laboratory to study
universal principles in protein folding.


  We propose a universal elastic energy for proteins, which depends only on the
radius of gyration $R_{g}$ and the residue number $N$. It is constructed using
physical arguments based on the hydrophobic effect and hydrogen bonding.
Adjustable parameters are fitted to data from the computer simulation of the
folding of a set of proteins using the CSAW (conditioned self-avoiding walk)
model. The elastic energy gives rise to scaling relations of the form
$R_{g}\sim N^{\nu}$ in different regions. It shows three folding stages
characterized by the progression with exponents $\nu = 3/5, 3/7, 2/5$, which we
identify as the unfolded stage, pre-globule, and molten globule, respectively.
The pre-globule goes over to the molten globule via a break in behavior akin to
a first-order phase transition, which is initiated by a sudden acceleration of
hydrogen bonding.


  A simple kinematic model for the trajectories of Listeria monocytogenes is
generalized to a dynamical system rich enough to exhibit the resonant Hopf
bifurcation structure of excitable media and simple enough to be studied
geometrically. It is shown how L. monocytogenes trajectories and meandering
spiral waves are organized by the same type of attracting set.


  We study a condition of favoring cooperation in a Prisoner's Dilemma game on
complex networks. There are two kinds of players: cooperators and defectors.
Cooperators pay a benefit b to their neighbors at a cost c, whereas defectors
only receive a benefit. The game is a death-birth process with weak selection.
Although it has been widely thought that b/c> <k> is a condition of favoring
cooperation, we find that b/c> <knn> is the condition. We also show that among
three representative networks, namely, regular, random, and scale-free, a
regular network favors cooperation the most, whereas a scale-free network
favors cooperation the least. In an ideal scale-free network where network size
is infinite, cooperation is never realized. Whether or not the scale-free
network and network heterogeneity favor cooperation depends on the details of a
game, although it is occasionally believed that scale-free networks favor
cooperation irrespective of game structures. If the number of players are
small, then the cooperation is favored in scale-free networks.


  The response to a knockout of a node is a characteristic feature of a
networked dynamical system. Knockout resilience in the dynamics of the
remaining nodes is a sign of robustness. Here we study the effect of knockouts
for binary state sequences and their implementations in terms of Boolean
threshold networks. Beside random sequences with biologically plausible
constraints, we analyze the cell cycle sequence of the species Saccharomyces
cerevisiae and the Boolean networks implementing it. Comparing with an
appropriate null model we do not find evidence that the yeast wildtype network
is optimized for high knockout resilience. Our notion of knockout resilience
weakly correlates with the size of the basin of attraction, which has also been
considered a measure of robustness.


  We investigate the impact of Allee effect and dispersal on the long-term
evolution of a population in a patchy environment, focusing on whether a
population already established in one patch either successfully invades an
adjacent empty patch or undergoes a global in-all-patch extinction. Our study
is based on the combination of analytical and numerical results for both a
deterministic two-patch model and its stochastic analog. The deterministic
model has either two or four attractors. In the presence of weak dispersal, the
analysis of the deterministic model shows that a high-density and a low-density
populations can coexist at equilibrium in nearby patches, whereas the analysis
of the stochastic model indicates that this equilibrium is metastable, thus
leading after a large random time to either an in-all-patch expansion or an
in-all-patch extinction. Up to some critical dispersal, increasing the
intensity of the interactions leads to an increase of both the basin of
attraction of the in-all-patch extinction and the basin of attraction of the
in-all-patch expansion. Above this threshold, while increasing the intensity of
the dispersal, both deterministic and stochastic models predict a
synchronization of the patches resulting in either a global expansion or a
global extinction: for the deterministic model, two of the four attractors
present when the dispersal is weak are lost, while the stochastic model no
longer exhibits a metastable behavior. In the presence of strong dispersal, the
limiting behavior is entirely determined by the value of the Allee threshold as
the global population size in the deterministic and the stochastic two-patch
models evolves as dictated by the single-patch counterparts ...


  We consider a slight modification to the monocyte and T-lymphocyte boundary
conditions of Little et al. (PLoS Comput Biol 2009 5(10) e1000539) and derive
alternative parameter estimates. No changes to the results and conclusions of
the paper of Little et al. (PLoS Comput Biol 2009 5(10) e1000539) are implied.


  We study the control parameters that govern the dynamics of in vitro DNA
ejection in bacteriophage lambda. Past work has demonstrated that bacteriophage
DNA is highly pressurized; this pressure has been hypothesized to help drive
DNA ejection. Ions influence this process by screening charges on DNA; however,
a systematic variation of salt concentrations to explore these effects has not
been undertaken. To study the nature of the forces driving DNA ejection, we
performed in vitro measurements of DNA ejection in bulk and at the single-phage
level. We present measurements on the dynamics of ejection and on the
self-repulsion force driving ejection. We examine the role of ion concentration
and identity in both measurements, and show that the charge of counter-ions is
an important control parameter. These measurements show that the frictional
force acting on the ejecting DNA is subtly dependent on ionic concentrations
for a given amount of DNA in the capsid. We also present evidence that phage
DNA forms loops during ejection; we confirm that this effect occurs using
optical tweezers. We speculate this facilitates circularization of the genome
in the cytoplasm.


  The infection pathway of virus in cytoplasm of a living cell is studied from
the viewpoint of diffusion theory. The cytoplasm plays a role of a medium for
stochastic motion of the virus contained in the endosome as well as the free
virus. It is experimentally known that the exponent of anomalous diffusion
fluctuates in localized areas of the cytoplasm. Here, generalizing fractional
kinetic theory, such fluctuations are described in terms of the exponent
locally distributed over the cytoplasm, and a theoretical proposition is
presented for its statistical form. The proposed fluctuations may be examined
in an experiment of heterogeneous diffusion in the infection pathway.


  Nanoparticles with multiple ligands have been proposed for use in
nanomedicine. The multiple targeting ligands on each nanoparticle can bind to
several locations on a cell surface facilitating both drug targeting and
uptake. Experiments show that the distribution of conjugated ligands is
unexpectedly broad, and the desorption rate appears to depends exponentially
upon the mean number of attached ligands. These two findings are explained with
a model in which ligands conjugate to the nanoparticle with a positive
cooperativity of $\approx 4kT$, and that nanoparticles bound to a surface by
multiple bonds are permanently affixed. This drives new analysis of the data,
which confirms that there is only one time constant for desorption, that of a
nanoparticle bound to the surface by a single bond.


  Complex, non-additive genetic interactions are common and can be critical in
determining phenotypes. Genome-wide association studies (GWAS) and similar
statistical studies of linkage data, however, assume additive models of gene
interactions in looking for genotype-phenotype associations. These statistical
methods view the compound effects of multiple genes on a phenotype as a sum of
partial influences of each individual gene and can often miss a substantial
part of the heritable effect. Such methods do not use any biological knowledge
about underlying genotype-phenotype mechanisms. Modeling approaches from the AI
field that incorporate deterministic knowledge into models to perform
statistical analysis can be applied to include prior knowledge in genetic
analysis. We chose to use the most general such approach, Markov Logic Networks
(MLNs), as a framework for combining deterministic knowledge with statistical
analysis. Using simple, logistic regression-type MLNs we have been able to
replicate the results of traditional statistical methods. Moreover, we show
that even with simple models we are able to go beyond finding independent
markers linked to a phenotype by using joint inference that avoids an
independence assumption. The method is applied to genetic data on yeast
sporulation, a phenotype governed by non-linear gene interactions. In addition
to detecting all of the previously identified loci associated with sporulation,
our method is able to identify four loci with small effects. Since their effect
on sporulation is small, these four loci were not detected with methods that do
not account for dependence between markers due to gene interactions. We show
how gene interactions can be detected using more complex models, which can be
used as a general framework for incorporating systems biology with genetics.


  Extinction of an epidemic or a species is a rare event that occurs due to a
large, rare stochastic fluctuation. Although the extinction process is
dynamically unstable, it follows an optimal path that maximizes the probability
of extinction. We show that the optimal path is also directly related to the
finite-time Lyapunov exponents of the underlying dynamical system in that the
optimal path displays maximum sensitivity to initial conditions. We consider
several stochastic epidemic models, and examine the extinction process in a
dynamical systems framework. Using the dynamics of the finite-time Lyapunov
exponents as a constructive tool, we demonstrate that the dynamical systems
viewpoint of extinction evolves naturally toward the optimal path.


  Momentum-space representation renders an interesting perspective to theory of
large fluctuations in populations undergoing Markovian stochastic gain-loss
processes. This representation is obtained when the master equation for the
probability distribution of the population size is transformed into an
evolution equation for the probability generating function. Spectral
decomposition then brings about an eigenvalue problem for a non-Hermitian
linear differential operator. The ground-state eigenmode encodes the stationary
distribution of the population size. For long-lived metastable populations
which exhibit extinction or escape to another metastable state, the
quasi-stationary distribution and the mean time to extinction or escape are
encoded by the eigenmode and eigenvalue of the lowest excited state. If the
average population size in the stationary or quasi-stationary state is large,
the corresponding eigenvalue problem can be solved via WKB approximation
amended by other asymptotic methods. We illustrate these ideas in several model
examples.


  Regulatory T cells (Treg) have recently been identified as playing a central
role in allergy and during allergen-specific immunotherapy. We have extended
our previous mathematical model describing the nonlinear dynamics of Th1-Th2
regulation by including Treg cells and their major cytokines. We hypothesize
that immunotherapy mainly acts on the T cell level and that the decisive
process can be regarded as a dynamical phenomenon. The model consists of
nonlinear differential equations which describe the proliferation and mutual
suppression of different T cell subsets. The old version of the model was based
upon the Th1-Th2 paradigm and is successful in describing the "Th1-Th2 switch"
which was considered the decisive event during specific immunotherapy. In
recent years, however, the Th1-Th2 paradigm has been questioned and therefore,
we have investigated a modified model in order to account for the influence of
a regulatory T cell type. We examined the extended model by means of numerical
simulations and analytical methods. As the modified model is more complex, we
had to develop new methods to portray its characteristics. The concept of
stable manifolds of fixed points of a stroboscobic map turned out to be
especially important. We found that when including regulatory T cells, our
model can describe the events in allergen-specific immunotherapy more
accurately. Our results suggest that the decisive effect of immunotherapy, the
increased proliferation of Treg and suppression of Th2 cells, crucially depends
on the administration of high dose injections right before the maintenance
phase sets in. Empirical protocols could therefore be improved by optimizing
this step of therapy.


  In most of the recent immunological literature the differences across antigen
receptor populations are examined via non-parametric statistical measures of
species overlap and diversity borrowed from ecological studies. While this
approach is robust in a wide range of situations, it seems to provide little
insight into the underlying clonal size distribution and the overall mechanism
differentiating the receptor populations. As a possible alternative, the
current paper presents a parametric method which adjusts for the data
under-sampling as well as provides a unifying approach to simultaneous
comparison of multiple receptor groups by means of the modern statistical tools
of unsupervised learning. The parametric model is based on a flexible
multivariate Poisson-lognormal distribution and is seen to be a natural
generalization of the univariate Poisson-lognormal models used in ecological
studies of biodiversity patterns. The procedure for evaluating model's fit is
described along with the public domain software developed to perform the
necessary diagnostics. The model-driven analysis is seen to compare favorably
vis a vis traditional methods when applied to the data from T-cell receptors in
transgenic mice populations.


  Due to the very long timescales involved (us-s), theoretical modeling of
fundamental biological processes including folding, misfolding, and mechanical
unraveling of biomolecules, under physiologically relevant conditions, is
challenging even for distributed computing systems. Graphics Processing Units
(GPUs) are emerging as an alternative programming platform to the more
traditional CPUs as they provide high raw computational power that can be
utilized in a wide range of scientific applications. Using a coarse-grained
Self Organized Polymer (SOP) model, we have developed and tested the GPU-based
implementation of Langevin simulations for proteins (SOP-GPU program).
Simultaneous calculation of forces for all particles is implemented using
either the particle based or the interacting pair based parallelization, which
leads to a ~30-fold acceleration compared to an optimized CPU version of the
program. We assess the computational performance of an end-to-end application
of the SOP-GPU program, where all steps of the algorithm are running on the
GPU, by profiling the associated simulation time and memory usage for a number
of small proteins, long protein fibers, and large-size protein assemblies. The
SOP-GPU package can now be used in the theoretical exploration of the
mechanical properties of large-size protein systems to generate the
force-extension and force-indentation profiles under the experimental
conditions of force application, and to relate the results of single-molecule
experiments in vitro and in silico.


  When gene copies are sampled from various species, the resulting gene tree
might disagree with the containing species tree. The primary causes of gene
tree and species tree discord include lineage sorting, horizontal gene
transfer, and gene duplication and loss. Each of these events yields a
different parsimony criterion for inferring the (containing) species tree from
gene trees. With lineage sorting, species tree inference is to find the tree
minimizing extra gene lineages that had to coexist along species lineages; with
gene duplication, it becomes to find the tree minimizing gene duplications
and/or losses. In this paper, we show the following results: (i) The deep
coalescence cost is equal to the number of gene losses minus two times the gene
duplication cost in the reconciliation of a uniquely leaf labeled gene tree and
a species tree. The deep coalescence cost can be computed in linear time for
any arbitrary gene tree and species tree. (ii) The deep coalescence cost is
always no less than the gene duplication cost in the reconciliation of an
arbitrary gene tree and a species tree. (iii) Species tree inference by
minimizing deep coalescences is NP-hard.


  Feedback loops are major components of biochemical systems. Many systems show
multiple such (positive or negative) feedback loops. Nevertheless, very few
quantitative analyses address the question how such multiple feedback loops
evolved. Based on published models from the mitotic cycle in embryogenesis, we
build a few case studies. Using a simple core architecture (transcription,
phosphorylation and degradation), we define oscillatory models having either
one positive feedback or one negative feedback, or both loops. With these
models, we address the following questions about evolvability: could a system
evolve from a simple model to a more complex one with a continuous transition
in the parameter space? How do new feedback loops emerge without disrupting the
proper function of the system? Our results show that progressive formation of a
second feedback loop is possible without disturbing existing oscillatory
behavior. For this process, the parameters of the system have to change during
evolution to maintain predefined properties of oscillations like period and
amplitude.


  The information content of a non-enzymatic self-replicator is limited by
Eigen's error threshold. Presumably, enzymatic replication can maintain higher
complexity, but in a competitive environment such a replicator is faced with
two problems related to its twofold role as enzyme and substrate: as enzyme, it
should replicate itself rather than wastefully copy non-functional substrates,
and as substrate it should preferably be replicated by superior enzymes instead
of less-efficient mutants. Because specific recognition can enforce these
propensities, we thoroughly analyze an idealized quasispecies model for
enzymatic replication, with replication rates that are either a decreasing
(self-specific) or increasing (cross-specific) function of the Hamming distance
between the recognition or "tag" sequences of enzyme and substrate. We find
that very weak self-specificity suffices to localize a population about a
master sequence and thus to preserve its information, while simultaneous
localization about complementary sequences in the cross-specific case is more
challenging. A surprising result is that stronger specificity constraints allow
longer recognition sequences, because the populations are better localized.
Extrapolating from experimental data, we obtain rough quantitative estimates
for the maximal length of the recognition or tag sequence that can be used to
reliably discriminate appropriate and infeasible enzymes and substrates,
respectively.


  The vertical distribution of phytoplankton is of fundamental importance for
the dynamics and structure of aquatic communities. Here, using an
advection-reaction-diffusion model, we investigate the distribution and
competition of phytoplankton species in a water column, in which inverse
resource gradients of light and a nutrient can limit growth of the biomass.
This problem poses a challenge for ecologists, as the location of a production
layer is not fixed, but rather depends on many internal parameters and
environmental factors. In particular, we study the influence of an upper mixed
layer (UML) in this system and show that it leads to a variety of dynamic
effects: (i) Our model predicts alternative density profiles with a maximum of
biomass either within or below the UML, thereby the system may be bistable or
the relaxation from an unstable state may require a long-lasting transition.
(ii) Reduced mixing in the deep layer can induce oscillations of the biomass;
we show that a UML can sustain these oscillations even if the diffusivity is
less than the critical mixing for a sinking phytoplankton population. (iii) A
UML can strongly modify the outcome of competition between different
phytoplankton species, yielding bistability both in the spatial distribution
and in the species composition. (iv) A light limited species can obtain a
competitive advantage if the diffusivity in the deep layers is reduced below a
critical value. This yields a subtle competitive exclusion effect, where the
oscillatory states in the deep layers are displaced by steady solutions in the
UML. Finally, we present a novel graphical approach for deducing the
competition outcome and for the analysis of the role of a UML in aquatic
systems.


  For driven Markovian dynamics on a network of (biomolecular) states, the
generalized mobilities, i.e., the response of any current to changes in an
external parameter, are expressed by an integral over an appropriate
current-current correlation function and thus related to the generalized
diffusion constants. As only input, a local detailed balance condition is
required typically even valid for biomolecular systems operating deep in the
non-equilibrium regime.


  The development of macroscopic descriptions for the joint dynamics and
behavior of large heterogeneous ensembles subject to ecological forces like
dispersal remains a central challenge for mathematicians and biological
scientists alike. Over the past century, specific attention has been directed
to the role played by dispersal in shaping plant communities, or on the
dynamics of marine open-ocean and intertidal systems, or on biological
invasions, or on the spread of disease, to name a few. Mathematicians and
theoreticians, starting with the efforts of researchers that include Aronson,
Fisher, Kolmogorov, Levin, Okubo, Skellam, Slobodkin, Weinberger and many
others, set the foundation of a fertile area of research at the interface of
ecology, mathematics, population biology and evolutionary biology.
Integrodifference systems, the subject of this manuscript, arise naturally in
the study of the spatial dispersal of organisms whose local population dynamics
are prescribed by models with discrete generations. The brunt of the
mathematical research has focused on the the study of existence traveling wave
solutions and characterizations of the spreading speed particularly, in the
context of cooperative systems. In this paper, we characterize the spreading
speed for a large class of non cooperative systems, all formulated in terms of
integrodifference equations, by the convergence of initial data to wave
solutions. In this setting, the spreading speed is characterized as the slowest
speed of a family of non-constant traveling wave solutions. Our results are
applied to a specific non-cooperative competition system in detail.


  The outcome of competition among species is influenced by the spatial
distribution of species and effects such as demographic stochasticity,
immigration fluxes, and the existence of preferred habitats. We introduce an
individual-based model describing the competition of two species and
incorporating all the above ingredients. We find that the presence of habitat
preference --- generating spatial niches --- strongly stabilizes the
coexistence of the two species. Eliminating habitat preference --- neutral
dynamics --- the model generates patterns, such as distribution of population
sizes, practically identical to those obtained in the presence of habitat
preference, provided an higher immigration rate is considered. Notwithstanding
the similarity in the population distribution, we show that invasibility
properties depend on habitat preference in a non-trivial way. In particular,
the neutral model results results more invasible or less invasible depending on
whether the comparison is made at equal immigration rate or at equal
distribution of population size, respectively. We discuss the relevance of
these results for the interpretation of invasibility experiments and the
species occupancy of preferred habitats.


  We introduce the concept of 'dynamical phase coexistence' to provide a simple
solution for a long-standing problem in theoretical ecology, the so-called
"savanna problem". The challenge is to understand why in savanna ecosystems
trees and grasses coexist in a robust way with large spatio-temporal
variability. We propose a simple model, a variant of the Contact Process (CP),
which includes two key extra features: varying external
(environmental/rainfall) conditions and tree age. The system fluctuates locally
between a woodland and a grassland phase, corresponding to the active and
absorbing phases of the underlying pure contact process. This leads to a highly
variable stable phase characterized by patches of the woodland and grassland
phases coexisting dynamically. We show that the mean time to tree extinction
under this model increases as a power-law of system size and can be of the
order of 10,000,000 years in even moderately sized savannas. Finally, we
demonstrate that while local interactions among trees may influence tree
spatial distribution and the order of the transition between woodland and
grassland phases, they do not affect dynamical coexistence. We expect dynamical
coexistence to be relevant in other contexts in physics, biology or the social
sciences.


  We examine the effective counterion-mediated electrostatic interaction
between two like-charged dielectric cylinders immersed in a continuous
dielectric medium containing neutralizing mobile counterions. We focus on the
effects of image charges induced as a result of the dielectric mismatch between
the cylindrical cores and the surrounding dielectric medium and investigate the
counterion-mediated electrostatic interaction between the cylinders in both
limits of weak and strong electrostatic couplings (corresponding, e.g., to
systems with monovalent and multivalent counterions, respectively). The results
are compared with extensive Monte-Carlo simulations exhibiting good agreement
with the limiting weak and strong coupling results in their respective regime
of validity.


  The Lee Carter modelling framework is widely used because of its simplicity
and robustness despite its inability to model specific cohort effects. A large
number of extensions have been proposed that model cohort effects but there is
no consensus. It is difficult to simultaneously account for cohort effects and
age-adjusted improvements and we here test a simple model that accounts for
both: we simply add a non age-adjusted cohort component to the Lee Carter
framework. This is a trade-off between accuracy and robustness but, for various
countries present in the Human Mortality Database and compared to various
models, the model accurately fits past mortality data and gives good
backtesting projections.


No abstract; review only


No abstract needed since short


  Background: RNA exhibits a variety of structural configurations. Here we
consider a structure to be tantamount to the noncrossing Watson-Crick and
\pairGU-base pairings (secondary structure) and additional cross-serial base
pairs. These interactions are called pseudoknots and are observed across the
whole spectrum of RNA functionalities. In the context of studying natural RNA
structures, searching for new ribozymes and designing artificial RNA, it is of
interest to find RNA sequences folding into a specific structure and to analyze
their induced neutral networks. Since the established inverse folding
algorithms, {\tt RNAinverse}, {\tt RNA-SSD} as well as {\tt INFO-RNA} are
limited to RNA secondary structures, we present in this paper the inverse
folding algorithm {\tt Inv} which can deal with 3-noncrossing, canonical
pseudoknot structures.
  Results: In this paper we present the inverse folding algorithm {\tt Inv}. We
give a detailed analysis of {\tt Inv}, including pseudocodes. We show that {\tt
Inv} allows to design in particular 3-noncrossing nonplanar RNA pseudoknot
3-noncrossing RNA structures--a class which is difficult to construct via
dynamic programming routines. {\tt Inv} is freely available at
\url{http://www.combinatorics.cn/cbpc/inv.html}.
  Conclusions: The algorithm {\tt Inv} extends inverse folding capabilities to
RNA pseudoknot structures. In comparison with {\tt RNAinverse} it uses new
ideas, for instance by considering sets of competing structures. As a result,
{\tt Inv} is not only able to find novel sequences even for RNA secondary
structures, it does so in the context of competing structures that potentially
exhibit cross-serial interactions.


  We introduce a mathematical model of symbiosis between different species by
taking into account the influence of each species on the carrying capacities of
the others. The modeled entities can pertain to biological and ecological
societies or to social, economic and financial societies. Our model includes
three basic types: symbiosis with direct mutual interactions, symbiosis with
asymmetric interactions, and symbiosis without direct interactions. In all
cases, we provide a complete classification of all admissible dynamical
regimes. The proposed model of symbiosis turned out to be very rich, as it
exhibits four qualitatively different regimes: convergence to stationary
states, unbounded exponential growth, finite-time singularity, and finite-time
death or extinction of species.


  This is my attempt to answer Schr\"odinger's question: "What is Life?". In
living cell the local combination of atoms is reproduced with incredible
accuracy. In the case of protein biosynthesis the notion of the physical aspect
of the process can be qualitatively enriched by taking into account the
symmetry of conformational oscillations. This is the key to understand the
selectivity in biochemistry.


  The aim of this paper is to discuss the role of robustness and diversity in
population dynamics in particular to some properties of the multi-step from
healthy tissue to fully malignant tumours. Recent evidence shows that diversity
within the cell population of a neoplasm, a pre-tumoural lession that can
develop into a fully malignant tumour, is the best predictor for its evolving
into a tumour. By studying the dynamics of a population described by a
multi-type, population-size limited branching process in terms of the
evolutionary formalism, we show some general principles regarding the
probability of a resident population to being invaded by a mutant population in
terms of the number of types present in the population and their resilience. We
show that, although diversity in the mutant population poses a barrier for the
emergence of the initial (benign) lession, under appropiate conditions, namely,
the phenotypes in the mutant population being more resilient than those of the
resident population, a more variable noeplastic population is more likely to be
invaded by a more malignant one. Analysis of a model of gene regulatory
networks suggest possible mechanisms giving rise to mutants with increased
phenotypic diversity and robustness. We then go on to show how these results
may help us to interpret some recent data regarding the evolution of Barrett's
oesophagus into throat cancer.


  Dynamical mechanisms that can stabilize the coexistence or diversity in
biology are generally of fundamental interest. In contrast to many two-strategy
evolutionary games, games with three strategies and cyclic dominance like the
rock-paper-scissors game (RPS) stabilize coexistence and thus preserve
biodiversity in this system. In the limit of infinite populations, resembling
the traditional picture of evolutionary game theory, replicator equations
predict the existence of a fixed point in the interior of the phase space. But
in finite populations, strategy frequencies will run out of the fixed point
because of stochastic fluctuations, and strategies can even go extinct. For
three different processes and for zero-sum and non-zero-sum RPS as well, we
present results of extensive simulations for the mean extinction time (MET),
depending on the number of agents N, and we introduce two analytical approaches
for the derivation of the MET.


  Light-harvesting bacteria Rhodospirillum Photometricum were recently found to
adopt strikingly different architectures depending on illumination conditions.
We present analytic and numerical calculations which explain this observation
by quantifying a dynamical interplay between excitation transfer kinetics and
reaction center cycling. High light-intensity membranes (HLIM) exploit
dissipation as a photo-protective mechanism, thereby safeguarding a steady
supply of chemical energy, while low light-intensity membranes (LLIM)
efficiently process unused illumination intensity by channelling it to open
reaction centers. More generally, our analysis elucidates and quantifies the
trade-offs in natural network design for solar energy conversion.


  We apply Doeblin's ergodicity coefficient as a computational tool to
approximate the occupancy distribution of a set of states in a homogeneous but
possibly non-stationary finite Markov chain. Our approximation is based on new
properties satisfied by this coefficient, which allow us to approximate a chain
of duration n by independent and short-lived realizations of an auxiliary
homogeneous Markov chain of duration of order ln(n). Our approximation may be
particularly useful when exact calculations via first-step methods or transfer
matrices are impractical, and asymptotic approximations may not be yet
reliable. Our findings may find applications to pattern problems in Markovian
and non-Markovian sequences that are treatable via embedding techniques.


  Starting with stochastic rate equations for the fundamental interactions
between microbes and their viruses, we derive a mean field theory for the
population dynamics of microbe-virus systems, including the effects of
lysogeny. In the absence of lysogeny, our model is a generalization of that
proposed phenomenologically by Weitz and Dushoff. In the presence of lysogeny,
we analyze the possible states of the system, identifying a novel limit cycle,
which we interpret physically. To test the robustness of our mean field
calculations to demographic fluctuations, we have compared our results with
stochastic simulations using the Gillespie algorithm. Finally, we estimate the
range of parameters that delineate the various steady states of our model.


  Many membrane channels and receptors exhibit adaptive, or desensitized,
response to a strong sustained input stimulus, often supported by protein
activity-dependent inactivation. Adaptive response is thought to be related to
various cellular functions such as homeostasis and enlargement of dynamic range
by background compensation. Here we study the quantitative relation between
adaptive response and background compensation within a modeling framework. We
show that any particular type of adaptive response is neither sufficient nor
necessary for adaptive enlargement of dynamic range. In particular a precise
adaptive response, where system activity is maintained at a constant level at
steady state, does not ensure a large dynamic range neither in input signal nor
in system output. A general mechanism for input dynamic range enlargement can
come about from the activity-dependent modulation of protein responsiveness by
multiple biochemical modification, regardless of the type of adaptive response
it induces. Therefore hierarchical biochemical processes such as methylation
and phosphorylation are natural candidates to induce this property in signaling
systems.


  Proteins interact with other proteins within biological pathways, forming
connected subgraphs in the protein-protein interactome (PPI). Proteins are
often involved in multiple biological pathways which complicates interpretation
of interactions between proteins. Gene expression data can assist our inference
since genes within a particular pathway tend to have more correlated expression
patterns than genes from distinct pathways. We provide an algorithm that uses
gene expression information to remove inter-pathway protein-protein
interactions, thereby simplifying the structure of the protein-protein
interactome. This refined topology permits easier interpretation and greater
biological coherence of multiple biological pathways simultaneously.


  The dynamical mechanisms that can stabilize the coexistence of species (or
strategies) are of substantial interest for the maintenance of biodiversity and
in sociobehavioural dynamics. We investigate the mean extinction time in the
coevolutionary dynamics of three cyclically invading strategies for different
evolutionary processes on various classes of complex networks, including random
graphs, scale-free and small world networks. We find that scale-free and random
graphs lead to a strong stabilization of coexistence both for the Moran process
and the Local Update process. The stabilization is of an order of magnitude
stronger compared to a lattice topology, and is mainly caused by the degree
heterogeneity of the graph. However, evolutionary processes on graphs can be
defined in many variants, and we show that in a process using effective payoffs
the effect of the network topology can be completely reversed. Thus,
stabilization of coexistence depends on both network geometry and underlying
evolutionary process.


  We study the Tangled Nature model of macro evolution and demonstrate that the
co-evolutionary dynamics produces an increasingly correlated core of well
occupied types. At the same time the entire configuration of types becomes
increasing de-correlated. This finding is related to ecosystem evolution. The
systems level dynamics of the model is subordinated to intermittent transitions
between meta-stable states. We improve on previous studies of the statistics of
the transition times and show that the fluctuations in the offspring
probability decreases with number of transitions. The longtime adaptation, as
seen by an increasing population size is demonstrated to be related to the
convexity of the offspring probability. We explain how the models behaviour is
a mathematical reflection of Darwin's concept of adaptation of profitable
variations.


  Deviations from the average can provide valuable insights about the
organization of natural systems. The present article extends this important
principle to the systematic identification and analysis of singular motifs in
complex networks. Six measurements quantifying different and complementary
features of the connectivity around each node of a network were calculated, and
multivariate statistical methods applied to identify singular nodes. The
potential of the presented concepts and methodology was illustrated with
respect to different types of complex real-world networks, namely the US air
transportation network, the protein-protein interactions of the yeast
Saccharomyces cerevisiae and the Roget thesaurus networks. The obtained
singular motifs possessed unique functional roles in the networks. Three
classic theoretical network models were also investigated, with the
Barab\'asi-Albert model resulting in singular motifs corresponding to hubs,
confirming the potential of the approach. Interestingly, the number of
different types of singular node motifs as well as the number of their
instances were found to be considerably higher in the real-world networks than
in any of the benchmark networks.


  A relationship between the preexponent of the rate constant and the
distribution over activation barrier energies for enzymatic/protein reactions
is revealed. We consider an enzyme solution as an ensemble of individual
molecules with different values of the activation barrier energy described by
the distribution. From the solvent viscosity effect on the preexponent we
derive the integral equation for the distribution and find its approximate
solution. Our approach enables us to attain a twofold purpose. On the one hand
it yields a simple interpretation of the solvent viscosity dependence for
enzymatic/protein reactions that requires neither a modification of the
Kramers' theory nor that of the Stokes law. On the other hand our approach
enables us to deduce the form of the distribution over activation barrier
energies. The obtained function has a familiar bell-shaped form and is in
qualitative agreement with the results of single enzyme kinetics measurements.
General formalism is exemplified by the analysis of literature experimental
data.


  Situations where individuals have to contribute to joint efforts or share
scarce resources are ubiquitous. Yet, without proper mechanisms to ensure
cooperation, the evolutionary pressure to maximize individual success tends to
create a tragedy of the commons (such as over-fishing or the destruction of our
environment). This contribution addresses a number of related puzzles of human
behavior with an evolutionary game theoretical approach as it has been
successfully used to explain the behavior of other biological species many
times, from bacteria to vertebrates. Our agent-based model distinguishes
individuals applying four different behavioral strategies: non-cooperative
individuals ("defectors"), cooperative individuals abstaining from punishment
efforts (called "cooperators" or "second-order free-riders"), cooperators who
punish non-cooperative behavior ("moralists"), and defectors, who punish other
defectors despite being non-cooperative themselves ("immoralists"). By
considering spatial interactions with neighboring individuals, our model
reveals several interesting effects: First, moralists can fully eliminate
cooperators. This spreading of punishing behavior requires a segregation of
behavioral strategies and solves the "second-order free-rider problem". Second,
the system behavior changes its character significantly even after very long
times ("who laughs last laughs best effect"). Third, the presence of a number
of defectors can largely accelerate the victory of moralists over non-punishing
cooperators. Forth, in order to succeed, moralists may profit from immoralists
in a way that appears like an "unholy collaboration". Our findings suggest that
the consideration of punishment strategies allows to understand the
establishment and spreading of "moral behavior" by means of game-theoretical
concepts. This demonstrates that quantitative biological modeling approaches
are powerful even in domains that have been addressed with non-mathematical
concepts so far. The complex dynamics of certain social behaviors becomes
understandable as result of an evolutionary competition between different
behavioral strategies.


  A molecular-level model is used to study the mechanical response of empty
cowpea chlorotic mottle virus (CCMV) and cowpea mosaic virus (CPMV) capsids.
The model is based on the native structure of the proteins that consitute the
capsids and is described in terms of the C-alpha atoms. Nanoindentation by a
large tip is modeled as compression between parallel plates. Plots of the
compressive force versus plate separation for CCMV are qualitatively consistent
with continuum models and experiments, showing an elastic region followed by an
irreversible drop in force. The mechanical response of CPMV has not been
studied, but the molecular model predicts an order of magnitude higher
stiffness and a much shorter elastic region than for CCMV. These large changes
result from small structural changes that increase the number of bonds by only
30% and would be difficult to capture in continuum models. Direct comparison of
local deformations in continuum and molecular models of CCMV shows that the
molecular model undergoes a gradual symmetry breaking rotation and accommodates
more strain near the walls than the continuum model. The irreversible drop in
force at small separations is associated with rupturing nearly all of the bonds
between capsid proteins in the molecular model while a buckling transition is
observed in continuum models.


  Recent research shows that a faulty or sub-optimally operating metabolic
network can often be rescued by the targeted removal of enzyme-coding
genes--the exact opposite of what traditional gene therapy would suggest.
Predictions go as far as to assert that certain gene knockouts can restore the
growth of otherwise nonviable gene-deficient cells. Many questions follow from
this discovery: What are the underlying mechanisms? How generalizable is this
effect? What are the potential applications? Here, I will approach these
questions from the perspective of compensatory perturbations on networks.
Relations will be drawn between such synthetic rescues and naturally occurring
cascades of reaction inactivation, as well as their analogues in physical and
other biological networks. I will specially discuss how rescue interactions can
lead to the rational design of antagonistic drug combinations that select
against resistance and how they can illuminate medical research on cancer,
antibiotics, and metabolic diseases.


  Up to now, the effects of having heterogeneous networks of contacts have been
studied mostly for diseases which are not persistent in time, i.e., for
diseases where the infectious period can be considered very small compared to
the lifetime of an individual. Moreover, all these previous results have been
obtained for closed populations, where the number of individuals does not
change during the whole duration of the epidemics. Here, we go one step further
and analyze, both analytically and numerically, a radically different kind of
diseases: those that are persistent and can last for an individual's lifetime.
To be more specific, we particularize to the case of Tuberculosis' (TB)
infection dynamics, where the infection remains latent for a period of time
before showing up and spreading to other individuals. We introduce an
epidemiological model for TB-like persistent infections taking into account the
heterogeneity inherent to the population structure. This sort of dynamics
introduces new analytical and numerical challenges that we are able to sort
out. Our results show that also for persistent diseases the epidemic threshold
depends on the ratio of the first two moments of the degree distribution so
that it goes to zero in a class of scale-free networks when the system
approaches the thermodynamic limit.


  We use interacting particle systems to investigate survival and extinction of
a species with colonies located on each site of $\mathbb {Z}^d$. In each of the
four models studied, an individual in a local population can reproduce, die or
migrate to neighboring sites. We prove that an increase of the death rate when
the local population density is small (the Allee effect) may be critical for
survival, and that the migration of large flocks of individuals is a possible
solution to avoid extinction when the Allee effect is strong. We use
attractiveness and comparison with oriented percolation, either to prove the
extinction of the species, or to construct nontrivial invariant measures for
each model.


  Many computerized methods for RNA-RNA interaction structure prediction have
been developed. Recently, $O(N^6)$ time and $O(N^4)$ space dynamic programming
algorithms have become available that compute the partition function of RNA-RNA
interaction complexes. However, few of these methods incorporate the knowledge
concerning related sequences, thus relevant evolutionary information is often
neglected from the structure determination. Therefore, it is of considerable
practical interest to introduce a method taking into consideration both
thermodynamic stability and sequence covariation. We present the \emph{a
priori} folding algorithm \texttt{ripalign}, whose input consists of two
(given) multiple sequence alignments (MSA). \texttt{ripalign} outputs (1) the
partition function, (2) base-pairing probabilities, (3) hybrid probabilities
and (4) a set of Boltzmann-sampled suboptimal structures consisting of
canonical joint structures that are compatible to the alignments. Compared to
the single sequence-pair folding algorithm \texttt{rip}, \texttt{ripalign}
requires negligible additional memory resource. Furthermore, we incorporate
possible structure constraints as input parameters into our algorithm. The
algorithm described here is implemented in C as part of the \texttt{rip}
package. The supplemental material, source code and input/output files can
freely be downloaded from \url{http://www.combinatorics.cn/cbpc/ripalign.html}.
\section{Contact} Christian Reidys \texttt{duck@santafe.edu}


  We use nucleosome maps obtained by high-throughput sequencing to study
sequence specificity of intrinsic histone-DNA interactions. In contrast with
previous approaches, we employ an analogy between a classical one-dimensional
fluid of finite-size particles in an arbitrary external potential and arrays of
DNA-bound histone octamers. We derive an analytical solution to infer free
energies of nucleosome formation directly from nucleosome occupancies measured
in high-throughput experiments. The sequence-specific part of free energies is
then captured by fitting them to a sum of energies assigned to individual
nucleotide motifs. We have developed hierarchical models of increasing
complexity and spatial resolution, establishing that nucleosome occupancies can
be explained by systematic differences in mono- and dinucleotide content
between nucleosomal and linker DNA sequences, with periodic dinucleotide
distributions and longer sequence motifs playing a secondary role. Furthermore,
similar sequence signatures are exhibited by control experiments in which
genomic DNA is either sonicated or digested with micrococcal nuclease in the
absence of nucleosomes, making it possible that current predictions based on
high-throughput nucleosome positioning maps are biased by experimental
artifacts.


  To understand complex biological systems, the research community has produced
huge corpus of gene expression data. A large number of clustering approaches
have been proposed for the analysis of gene expression data. However,
extracting important biological knowledge is still harder. To address this
task, clustering techniques are used. In this paper, hybrid Hierarchical
k-Means algorithm is used for clustering and biclustering gene expression data
is used. To discover both local and global clustering structure biclustering
and clustering algorithms are utilized. A validation technique, Figure of Merit
is used to determine the quality of clustering results. Appropriate knowledge
is mined from the clusters by embedding a BLAST similarity search program into
the clustering and biclustering process. To discover both local and global
clustering structure biclustering and clustering algorithms are utilized. To
determine the quality of clustering results, a validation technique, Figure of
Merit is used. Appropriate knowledge is mined from the clusters by embedding a
BLAST similarity search program into the clustering and biclustering process.


  One of the challenges in epidemiology is to account for the complex
morphological structure of hosts such as plant roots, crop fields, farms,
cells, animal habitats and social networks, when the transmission of infection
occurs between contiguous hosts. Morphological complexity brings an inherent
heterogeneity in populations and affects the dynamics of pathogen spread in
such systems. We have analysed the influence of realistically complex host
morphology on the threshold for invasion and epidemic outbreak in an SIR
(susceptible-infected-recovered) epidemiological model. We show that disorder
expressed in the host morphology and anisotropy reduces the probability of
epidemic outbreak and thus makes the system more resistant to epidemic
outbreaks. We obtain general analytical estimates for minimally safe bounds for
an invasion threshold and then illustrate their validity by considering an
example of host data for branching hosts (salamander retinal ganglion cells).
Several spatial arrangements of hosts with different degrees of heterogeneity
have been considered in order to analyse separately the role of shape
complexity and anisotropy in the host population. The estimates for invasion
threshold are linked to morphological characteristics of the hosts that can be
used for determining the threshold for invasion in practical applications.


  Confinement is a versatile and well-established tool to study the properties
of polymers either to understand biological processes or to develop new
nano-biomaterials. We investigate the conformations of a semiflexible polymer
ring in weak spherical confinement imposed by an impenetrable shell. We develop
an analytic argument for the dominating polymer trajectory depending on polymer
flexibility considering elastic and entropic contributions. Monte Carlo
simulations are performed to assess polymer ring conformations in probability
densities and by the shape measures asphericity and nature of asphericity.
Comparison of the analytic argument with the mean asphericity and the mean
nature of asphericity confirm our reasoning to explain polymer ring
conformations in the stiff regime, where elastic response prevails.


  High-throughput screening (HTS) using model organisms is a promising method
to identify a small number of genes or drugs potentially relevant to human
biology or disease. In HTS experiments, robots and computers do a significant
portion of the experimental work. However, one remaining major bottleneck is
the manual analysis of experimental results, which is commonly in the form of
microscopy images. This manual inspection is labor intensive, slow and
subjective. Here we report our progress towards applying computer vision and
machine learning methods to analyze HTS experiments that use Caenorhabditis
elegans (C. elegans) worms grown on agar. Our main contribution is a robust
segmentation algorithm for separating the worms from the background using
brightfield images. We also show that by combining the output of this
segmentation algorithm with an algorithm to detect the fluorescent dye, Nile
Red, we can reliably distinguish different fluorescence-based phenotypes even
though the visual differences are subtle. The accuracy of our method is similar
to that of expert human analysts. This new capability is a significant step
towards fully automated HTS experiments using C. elegans.


  Following a short report of our preliminary results [Phys. Rev. E 79,
055203(R) (2009)], we present a more detailed study of the effects of coupling
delay in diffusively coupled phase oscillator populations. We find that
coupling delay induces chimera and globally clustered chimera (GCC) states in
delay coupled populations. We show the existence of multi-clustered states that
act as link between the chimera and the GCC states. A stable GCC state goes
through a variety of GCC states, namely periodic, aperiodic, long-- and
short--period breathers and becomes unstable GCC leading to global
synchronization in the system, on increasing time delay. We provide numerical
evidence and theoretical explanations for the above results and discuss
possible applications of the observed phenomena.


  Food webs can be regarded as energy transporting networks in which the weight
of each edge denotes the energy flux between two species. By investigating 21
empirical weighted food webs as energy flow networks, we found several
ubiquitous scaling behaviors. Two random variables $A_i$ and $C_i$ defined for
each vertex $i$, representing the total flux (also called vertex intensity) and
total indirect effect or energy store of $i$, were found to follow power law
distributions with the exponents $\alpha\approx 1.32$ and $\beta\approx 1.33$,
respectively. Another scaling behavior is the power law relationship, $C_i\sim
A_i^\eta$, where $\eta\approx 1.02$. This is known as the allometric scaling
power law relationship because $A_i$ can be treated as metabolism and $C_i$ as
the body mass of the sub-network rooted from the vertex $i$, according to the
algorithm presented in this paper. Finally, a simple relationship among these
power law exponents, $\eta=(\alpha-1)/(\beta-1)$, was mathematically derived
and tested by the empirical food webs.


  Using lattice models we explore the factors that determine the tendencies of
polypeptide chains to aggregate by exhaustively sampling the sequence and
conformational space. The morphologies of the fibril-like structures and the
time scales ($\tau_{fib}$) for their formation depend on a subtle balance
between hydrophobic and coulomb interactions. The extent of population of a
fibril-prone structure in the spectrum of monomer conformations is the major
determinant of $\tau_{fib}$. This observation is used to determine the
aggregation-prone consensus sequences by exhaustively exploring the sequence
space. Our results provide a basis for genome wide search of fragments that are
aggregation prone.


  We dedicate this paper to Sir John Kingman on his 70th Birthday. In modern
mathematical population genetics the ancestral history of a population of genes
back in time is described by John Kingman's coalescent tree. Classical and
modern approaches model gene frequencies by diffusion processes. This paper,
which is partly a review, discusses how coalescent processes are dual to
diffusion processes in an analytic and probabilistic sense. Bochner (1954) and
Gasper (1972) were interested in characterizations of processes with Beta
stationary distributions and Jacobi polynomial eigenfunctions. We discuss the
connection with Wright--Fisher diffusions and the characterization of these
processes. Subordinated Wright--Fisher diffusions are of this type. An Inverse
Gaussian subordinator is interesting and important in subordinated
Wright--Fisher diffusions and is related to the Jacobi Poisson Kernel in
orthogonal polynomial theory. A related time-subordinated forest of non-mutant
edges in the Kingman coalescent is novel.


  A striking feature of the marine ecosystem is the regularity in its size
spectrum: the abundance of organisms as a function of their weight
approximately follows a power law over almost ten orders of magnitude. We
interpret this as evidence that the population dynamics in the ocean is
approximately scale-invariant. We use this invariance in the construction and
solution of a size-structured dynamical population model. Starting from a
Markov model encoding the basic processes of predation, reproduction,
maintenance respiration and intrinsic mortality, we derive a partial
integro-differential equation describing the dependence of abundance on weight
and time. Our model represents an extension of the jump-growth model and hence
also of earlier models based on the McKendrick--von Foerster equation. The
model is scale-invariant provided the rate functions of the stochastic
processes have certain scaling properties. We determine the steady-state power
law solution, whose exponent is determined by the relative scaling between the
rates of the density-dependent processes (predation) and the rates of the
density-independent processes (reproduction, maintenance, mortality). We study
the stability of the steady-state against small perturbations and find that
inclusion of maintenance respiration and reproduction in the model has astrong
stabilising effect. Furthermore, the steady state is unstable against a change
in the overall population density unless the reproduction rate exceeds a
certain threshold.


  Energy transfer (relaxation) phenomena are ubiquitous in nature. At a
macroscopic level, the phenomenological theory of heat (Fourier law)
successfully describes heat transfer and energy flow. However, its microscopic
origin is still under debate. This is because the phenomena can contain
many-body, multi-scale, nonequilibrium, and even quantum mechanical aspects,
which present significant challenges to theories addressing energy transfer
phenomena in physics, chemistry and biology. In this paper, we describe our
recent theoretical attempts to treat vibrational energy relaxation (VER) in
biomolecular systems, including peptide, porphyrin, and protein. First we
construct a reduced model using (instantaneous) normal mode analysis, and
further add anharmonic coupling between vibrational modes. Using such a model
combined with time-dependent perturbation theory for the anharmonic coupling,
we derive non-Markovian VER formulas, which are applicable to small timescales
such as sub picoseconds. We apply the VER formulas with full atomistic detail
to various biomolecular systems, discuss the VER timescales and mechanisms, and
compare with experiment if possible. We finally mention further theoretical
prospects of the energy transfer phenomena in physics, chemistry, and biology.


  Population analysis is persistently challenging but important, leading to the
determination of diversity and function prediction of microbial community
members. Here we detail our bioinformatics methods for analyzing population
distribution and diversity in large microbial communities. This was achieved
via (i) a homology based method for robust phylotype determination, equaling
the classification accuracy of the Ribosomal Database Project (RDP) classifier,
but providing improved associations of closely related sequences; (ii) a
comparison of different clustering methods for achieving more accurate richness
estimations. Our methodology, which we developed using the RDP vetted 16S rRNA
gene sequence set, was validated by testing it on a large 16S rRNA gene dataset
of approximately 2300 sequences, which we obtained from a soil microbial
community study. We concluded that the best approach to obtain accurate
phylogenetics profile of large microbial communities, based on 16S rRNA gene
sequence information, is to apply an optimized blast classifier. This approach
is complemented by the grouping of closely related sequences, using complete
linkage clustering, in order to calculate richness and evenness indices for the
communities.


  We estimate by Monte Carlo simulations the configurational entropy of
$N$-steps polygons in the cubic lattice with fixed knot type. By collecting a
rich statistics of configurations with very large values of $N$ we are able to
analyse the asymptotic behaviour of the partition function of the problem for
different knot types. Our results confirm that, in the large $N$ limit, each
prime knot is localized in a small region of the polygon, regardless of the
possible presence of other knots. Each prime knot component may slide along the
unknotted region contributing to the overall configurational entropy with a
term proportional to $\ln N$. Furthermore, we discover that the mere existence
of a knot requires a well defined entropic cost that scales exponentially with
its minimal length. In the case of polygons with composite knots it turns out
that the partition function can be simply factorized in terms that depend only
on prime components with an additional combinatorial factor that takes into
account the statistical property that by interchanging two identical prime knot
components in the polygon the corresponding set of overall configuration
remains unaltered. Finally, the above results allow to conjecture a sequence of
inequalities for the connective constants of polygons whose topology varies
within a given family of composite knot types.


  Gene regulatory networks typically have low in-degrees, whereby any given
gene is regulated by few of the genes in the network. They also tend to have
broad distributions for the out-degree. What mechanisms might be responsible
for these degree distributions? Starting with an accepted framework of the
binding of transcription factors to DNA, we consider a simple model of gene
regulatory dynamics. There, we show that selection for a target expression
pattern leads to the emergence of minimum connectivities compatible with the
selective constraint. As a consequence, these gene networks have low in-degree,
and "functionality" is parsimonious, i.e. is concentrated on a sparse number of
interactions as measured for instance by their essentiality. Furthermore, we
find that mutations of the transcription factors drive the networks to have
broad out-degrees. Finally, these classes of models are evolvable, i.e.
significantly different genotypes can emerge gradually under mutation-selection
balance.


  The geometrical coupling between strain and twist in double helices is
investigated. Overwinding, where strain leads to further winding, is shown to
be a universal property for helices, which are stretched along their
longitudinal axis when the initial pitch angle is below the zero-twist angle
(39.4 deg). Unwinding occurs at larger pitch angles. The zero-twist angle is
the unique pitch angle at the point between overwinding and unwinding, and it
is independent of the mechanical properties of the double helix. This suggests
the existence of zero-twist structures, i.e. structures that display neither
overwinding, nor unwinding under strain. Estimates of the overwinding of DNA,
chromatin, and RNA are given.


  The calcium transport in biological systems is modelled as a
reaction-diffusion process. Nonlinear calcium waves are then simulated using a
stochastic cellular automaton whose rules are derived from the corresponding
coupled partial differential equations. Numerical simulations show
self-organized criticality in the complex calcium waves and patterns. Both the
stochastic cellular automaton approach and the equation-based simulations can
predict the characteristics of calcium waves and complex pattern formation. The
implication of locality of calcium distribution with positional information in
biological systems is also discussed.


  Background: Recent experimental and theoretical studies have shown that small
asexual populations evolving on complex fitness landscapes may achieve a higher
fitness than large ones due to the increased heterogeneity of adaptive
trajectories. Here we introduce a class of haploid three-locus fitness
landscapes that allows to investigate this scenario in a precise and
quantitative way.
  Results: Our main result derived analytically shows how the probability of
choosing the path of largest initial fitness increase grows with the population
size. This makes large populations more likely to get trapped at local fitness
peaks and implies an advantage of small populations at intermediate time
scales. The range of population sizes where this effect is operative coincides
with the onset of clonal interference. Additional studies using ensembles of
random fitness landscapes show that the results achieved for a particular
choice of three-locus landscape parameters are robust and also persist as the
number of loci increases.
  Conclusions: Our study indicates that an advantage for small populations is
likely whenever the fitness landscape contains local maxima. The advantage
appears at intermediate time scales, which are long enough for trapping at
local fitness maxima to have occurred but too short for peak escape by the
creation of multiple mutants.


  In most models of the spread of disease over contact networks it is assumed
that the probabilities per unit time of disease transmission and recovery from
disease are constant, implying exponential distributions of the time intervals
for transmission and recovery. Time intervals for real diseases, however, have
distributions that in most cases are far from exponential, which leads to
disagreements, both qualitative and quantitative, with the models. In this
paper, we study a generalized version of the SIR
(susceptible-infected-recovered) model of epidemic disease that allows for
arbitrary distributions of transmission and recovery times. Standard
differential equation approaches cannot be used for this generalized model, but
we show that the problem can be reformulated as a time-dependent message
passing calculation on the appropriate contact network. The calculation is
exact on trees (i.e., loopless networks) or locally tree-like networks (such as
random graphs) in the large system size limit. On non-tree-like networks we
show that the calculation gives a rigorous bound on the size of disease
outbreaks. We demonstrate the method with applications to two specific models
and the results compare favorably with numerical simulations.


  How fast does a population evolve from one fitness peak to another? We study
the dynamics of evolving, asexually reproducing populations in which a certain
number of mutations jointly confer a fitness advantage. We consider the time
until a population has evolved from one fitness peak to another one with a
higher fitness. The order of mutations can either be fixed or random. If the
order of mutations is fixed, then the population follows a metaphorical ridge,
a single path. If the order of mutations is arbitrary, then there are many ways
to evolve to the higher fitness state. We address the time required for
fixation in such scenarios and study how it is affected by the order of
mutations, the population size, the fitness values and the mutation rate.


  Evolutionary game dynamics of two players with two strategies has been
studied in great detail. These games have been used to model many biologically
relevant scenarios, ranging from social dilemmas in mammals to microbial
diversity. Some of these games may in fact take place between a number of
individuals and not just between two. Here, we address one-shot games with
multiple players. As long as we have only two strategies, many results from two
player games can be generalized to multiple players. For games with multiple
players and more than two strategies, we show that statements derived for
pairwise interactions do no longer hold. For two player games with any number
of strategies there can be at most one isolated internal equilibrium. For any
number of players $\boldsymbol{d}$ with any number of strategies n, there can
be at most (d-1)^(n-1) isolated internal equilibria. Multiplayer games show a
great dynamical complexity that cannot be captured based on pairwise
interactions. Our results hold for any game and can easily be applied for
specific cases, e.g. public goods games or multiplayer stag hunts.


  Likelihood-based phylogenetic inference is generally considered to be the
most reliable classification method for unknown sequences. However, traditional
likelihood-based phylogenetic methods cannot be applied to large volumes of
short reads from next-generation sequencing due to computational complexity
issues and lack of phylogenetic signal. "Phylogenetic placement," where a
reference tree is fixed and the unknown query sequences are placed onto the
tree via a reference alignment, is a way to bring the inferential power of
likelihood-based approaches to large data sets. This paper introduces pplacer,
a software package for phylogenetic placement and subsequent visualization. The
algorithm can place twenty thousand short reads on a reference tree of one
thousand taxa per hour per processor, has essentially linear time and memory
complexity in the number of reference taxa, and is easy to run in parallel.
Pplacer features calculation of the posterior probability of a placement on an
edge, which is a statistically rigorous way of quantifying uncertainty on an
edge-by-edge basis. It also can inform the user of the positional uncertainty
for query sequences by calculating expected distance between placement
locations, which is crucial in the estimation of uncertainty with a
well-sampled reference tree. The software provides visualizations using branch
thickness and color to represent number of placements and their uncertainty. A
simulation study using reads generated from 631 COG alignments shows a high
level of accuracy for phylogenetic placement over a wide range of alignment
diversity, and the power of edge uncertainty estimates to measure placement
confidence. Pplacer enables efficient phylogenetic placement and subsequent
visualization, making likelihood-based phylogenetics methodology practical for
large collections of reads; it is available as source code, binaries, and a web
service.


  This paper studies a Markov chain for phylogenetic reconstruction which uses
a popular transition between tree topologies known as subtree
pruning-and-regrafting (SPR). We analyze the Markov chain in the simpler
setting that the generating tree consists of very short edge lengths, short
enough so that each sample from the generating tree (or character in
phylogenetic terminology) is likely to have only one mutation, and that there
enough samples so that the data looks like the generating distribution. We
prove in this setting that the Markov chain is rapidly mixing, i.e., it quickly
converges to its stationary distribution, which is the posterior distribution
over tree topologies. Our proofs use that the leading term of the maximum
likelihood function of a tree T is the maximum parsimony score, which is the
size of the minimum cut in T needed to realize single edge cuts of the
generating tree. Our main contribution is a combinatorial proof that in our
simplified setting, SPR moves are guaranteed to converge quickly to the maximum
parsimony tree. Our results are in contrast to recent works showing examples
with heterogeneous data (namely, the data is generated from a mixture
distribution) where many natural Markov chains are exponentially slow to
converge to the stationary distribution.


  In this work we introduce a differential equation model with time-delay that
describes the three-stage dynamics and the two time scales observed in HIV
infection. Assuming that the virus has high mutation and rapid reproduction
rates that stress the immune system throughout the successive activation of new
responses to new undetectable strains, the delay term describes the time
interval necessary to mount new specific immune responses. This single term
increases the number of possible solutions and changes the phase space dynamics
if compared to the model without time delay. We observe very slow transits near
the unstable fixed point, corresponding to a healthy state, and long time decay
to the stable fixed point that corresponds to the infected state. In contrast
to the results obtained for models using regular ODE, which only allow for
partial descriptions of the course of the infection, our model describes the
entire course of infection observed in infected patients: the primary
infection, the latency period and the onset of acquired immunodeficiency
syndrome (AIDS). The model also describes other scenarios, such as the very
fast progression to the disease and the less common outcome in which, although
the patient is exposed to HIV, he/she does not develop the disease.


  For a given, arbitrary graph, what is the epidemic threshold? That is, under
what conditions will a virus result in an epidemic? We provide the super-model
theorem, which generalizes older results in two important, orthogonal
dimensions. The theorem shows that (a) for a wide range of virus propagation
models (VPM) that include all virus propagation models in standard literature
(say, [8][5]), and (b) for any contact graph, the answer always depends on the
first eigenvalue of the connectivity matrix. We give the proof of the theorem,
arithmetic examples for popular VPMs, like flu (SIS), mumps (SIR), SIRS and
more. We also show the implications of our discovery: easy (although sometimes
counter-intuitive) answers to `what-if' questions; easier design and evaluation
of immunization policies, and significantly faster agent-based simulations.


  MicroRNAs are endogenous non-coding RNAs which negatively regulate the
expression of protein-coding genes in plants and animals. They are known to
play an important role in several biological processes and, together with
transcription factors, form a complex and highly interconnected regulatory
network. Looking at the structure of this network it is possible to recognize a
few overrepresented motifs which are expected to perform important elementary
regulatory functions. Among them a special role is played by the
microRNA-mediated feedforward loop in which a master transcription factor
regulates a microRNA and, together with it, a set of target genes. In this
paper we show analytically and through simulations that the incoherent version
of this motif can couple the fine-tuning of a target protein level with an
efficient noise control, thus conferring precision and stability to the overall
gene expression program, especially in the presence of fluctuations in upstream
regulators. Among the other results, a nontrivial prediction of our model is
that the optimal attenuation of fluctuations coincides with a modest repression
of the target expression. This feature is coherent with the expected
fine-tuning function and in agreement with experimental observations of the
actual impact of a wide class of microRNAs on the protein output of their
targets. Finally we describe the impact on noise-buffering efficiency of the
cross-talk between microRNA targets that can naturally arise if the
microRNA-mediated circuit is not considered as isolated, but embedded in a
larger network of regulations.


  A model is presented from which the observed morphology of the inner
mitochondrial membrane can be inferred as minimizing the system's free energy.
Besides the usual energetic terms for bending, surface area, and pressure
difference, our free energy includes terms for tension that we believe to be
exerted by proteins and for an entropic contribution due to many dimensions
worth of shapes available at a given energy.
  In order to test the model, we measured the structural features of
mitochondria in HeLa cells and mouse embryonic fibroblasts using 3D electron
tomography. Such tomograms reveal that the inner membrane self-assembles into a
complex structure that contains both tubular and flat lamellar crista
components. This structure, which contains one matrix compartment, is believed
to be essential to the proper functioning of mitochondria as the powerhouse of
the cell. We find that tensile forces of the order of 10 pN are required to
stabilize a stress-induced coexistence of tubular and flat lamellar cristae
phases. The model also predicts \Deltap = -0.036 \pm 0.004 atm and \sigma=0.09
\pm 0.04 pN/nm.


  We describe the evolution of the quantity of parasites in a population of
cells which divide in continuous-time. The quantity of parasites in a cell
follows a Feller diffusion, which is splitted randomly between the two daughter
cells when a division occurs. The cell division rate may depend on the quantity
of parasites inside the cell and we are interested in the cases of constant or
monotone division rate. We first determine the asymptotic behavior of the
quantity of parasites in a cell line, which follows a Feller diffusion with
multiplicative jumps. We then consider the evolution of the infection of the
cell population and give criteria to determine whether the proportion of
infected cells goes to zero (recovery) or if a positive proportion of cells
becomes largely infected (proliferation of parasites inside the cells).


  The Wako-Sait{\^o}-Mu\~noz-Eaton (WSME) model, initially introduced in the
theory of protein folding, has also been used in modeling the RNA folding and
some epitaxial phenomena. The advantage of this model is that it admits exact
solution in the general inhomogeneous case (Bruscolini and Pelizzola, 2002)
which facilitates the study of realistic systems. However, a shortcoming of the
model is that it accounts only for interactions within continuous stretches of
native bonds or atomic chains while neglecting interstretch (interchain)
interactions. But due to the biopolymer (atomic chain) flexibility, the
monomers (atoms) separated by several non-native bonds along the sequence can
become closely spaced. This produces their strong interaction. The inclusion of
non-WSME interactions into the model makes the model more realistic and
improves its performance. In this study we add arbitrary interactions of finite
range and solve the new model by means of the transfer matrix technique. We can
therefore exactly account for the interactions which in proteomics are
classified as medium- and moderately long-range ones.


  A morphogenic protein known as Dorsal patterns the embryonic dorsoventral
body axis of Drosophila by binding to transcriptional enhancers across the
genome. Each such enhancer activates a neighboring gene at a unique threshold
concentration of Dorsal. The presence of Dorsal binding site clusters in these
enhancers and of similar clusters in other enhancers has motivated models of
threshold-encoding in site density. However, we found that the precise length
of a spacer separating a pair of specialized Dorsal and Twist binding sites
determines the threshold-response. Despite this result, the functional range
determined by this spacer element as well as the role and origin of its
surrounding Dorsal site cluster remained completely unknown. Here, we
experiment with enhancers from diverse Drosophila genomes, including the large
uncompacted genomes from ananassae and willistoni, and report three major
interdependent results. First, we map the functional range of the
threshold-encoding spacer variable. Second, we show that the majority of sites
at the cluster are non-functional divergent elements that have been separated
beyond the encoding's functional range. Third, we verify an evolutionary model
involving the frequent replacement of a threshold encoding, whose precision is
easily outdated by shifting accuracy. The process by which encodings are
replaced by newer ones is facilitated by the palindromic nature of the Dorsal
and Twist binding motifs and by intrinsic repeat-instability in the specialized
Twist binding site, which critically impacts the length of the spacer linking
it to Dorsal. Over time, the dynamic process of selective deprecation and
replacement of encodings adds to a growing cluster of deadened elements, or
necro-elements, and strongly biases local sequence composition. ... [ABSTRACT
TRUNCATED]