  Voltage-gated ion channels mediate electrical excitability of cellular
membranes. Reduced models of the voltage sensor (VS) of Kv channels produce
insight into the electrostatic physics underlying the response of the highly
positively charged S4 transmembrane domain to changes in membrane potential and
other electrostatic parameters. By calculating the partition function computed
from the electrostatic energy over translational and/or rotational degrees of
freedom, I compute expectations of charge displacement, energetics, probability
distributions of translation & rotation and Maxwell stress for arrangements of
S4 positively charged residues; these computations can then be compared with
experimental results to elucidate the role of various putative atomic level
features of the VS. A `paddle' model is rejected on electrostatic grounds,
owing to unfavorable energetics, insufficient charge displacement and excessive
Maxwell stress. On the other hand, a `sliding helix' model with three local
counter-charges, a protein dielectric coefficient of 4 and a 2/3 interval of
counter-charge positioning relative to the S4-helix period of positive residues
is electrostatically reasonable, comparing well with Shaker (Seoh et al.,
1996). Lack of counter-charges destabilizes the S4 in the membrane;
counter-charge interval helps determine the number and shape of energy barriers
and troughs over the range of motion of the S4; and the local dielectric
coefficient of the protein constrains the height of energy maxima. These
`sliding helix' models compare favorably with experimental results for single &
double mutant charge experiments on Shaker. Single S4 positive charge mutants
are predicted quite well by this model; single counter-charge mutants are
predicted less well; and double mutants for both an S4 charge and a
counter-charge are characterized least well.


  Adaptation to changing environments is a hallmark of biological systems.
Diversity in traits is necessary for adaptation and can influence the survival
of a population faced with novelty. In habitats that remain stable over many
generations, stabilizing selection reduces trait differences within
populations, thereby appearing to remove the diversity needed for heritable
adaptive responses in new environments. Paradoxically, field studies have
documented numerous populations under long periods of stabilizing selection and
evolutionary stasis that have rapidly evolved under changed environmental
conditions. In this article, we review how cryptic genetic variation (CGV)
resolves this diversity paradox by allowing populations in a stable environment
to gradually accumulate hidden genetic diversity that is revealed as trait
differences when environments change. Instead of being in conflict,
environmental stasis supports CGV accumulation and thus appears to facilitate
rapid adaptation in new environments as suggested by recent CGV studies.
Similarly, degeneracy has been found to support both genetic and non-genetic
adaptation at many levels of biological organization. Degenerate, as opposed to
diverse or redundant, ensembles appear functionally redundant in certain
environmental contexts but functionally diverse in others. CGV and degeneracy
paradigms for adaptation are integrated in this review, revealing a common set
of principles that support adaptation at multiple levels of biological
organization. Though a discussion of simulation studies, molecular-based
experimental systems, principles from population genetics, and field
experiments, we demonstrate that CGV and degeneracy reflect complementary
top-down and bottom-up, respectively, conceptualizations of the same basic
phenomenon and arguably capture a universal feature of biological adaptive
processes.


  In this article we demonstrate that the so-called bursting production of
molecular species during gene expression may be an artifact caused by low time
resolution in experimental data collection and not an actual burst in
production. We reach this conclusion through an analysis of a two-stage and
binary model for gene expression, and demonstrate that in the limit when mRNA
degradation is much faster than protein degradation they are equivalent. The
negative binomial distribution is shown to be a limiting case of the binary
model for fast "on to off" state transitions and high values of the ratio
between protein synthesis and degradation rates. The gene products population
increases by unity but multiple times in a time interval orders of magnitude
smaller than protein half-life or the precision of the experimental apparatus
employed in its detection. This rare-and-fast one-by-one protein synthesis has
been interpreted as bursting.


  Cryo-electron microscopy (cryo-EM) studies using single particle
reconstruction is extensively used to reveal structural information of
macromolecular complexes. Aiming at the highest achievable resolution, state of
the art electron microscopes acquire thousands of high-quality images. Having
collected these data, each single particle must be detected and windowed out.
Several fully- or semi-automated approaches have been developed for the
selection of particle images from digitized micrographs. However they still
require laborious manual post processing, which will become the major
bottleneck for next generation of electron microscopes. Instead of focusing on
improvements in automated particle selection from micrographs, we propose a
post-picking step for classifying small windowed images, which are output by
common picking software. A supervised strategy for the classification of
windowed micrograph images into particles and non-particles reduces the manual
workload by orders of magnitude. The method builds on new powerful image
features, and the proper training of an ensemble classifier. A few hundred
training samples are enough to achieve a human-like classification performance.


  DNA polymerase (DNAP) is a dual-purpose enzyme that plays two opposite roles
in two different situations during DNA replication. It plays its normal role as
a {\it polymerase} catalyzing the elongation of a new DNA molecule by adding a
monomer. However, it can switch to the role of an {\it exonuclease} and shorten
the same DNA by cleavage of the last incorporated monomer from the nascent DNA.
Just as misincorporated nucleotides can escape exonuclease causing replication
error, correct nucleotide may get sacrificed unnecessarily by erroneous
cleavage. The interplay of polymerase and exonuclease activities of a DNAP is
explored here by developing a minimal stochastic kinetic model of DNA
replication. Exact analytical expressions are derived for a few key statistical
distributions; these characterize the temporal patterns in the mechanical
stepping and the chemical (cleavage) reaction. The Michaelis-Menten-like
analytical expression derived for the average rates of these two processes not
only demonstrate the effects of their coupling, but are also utilized to
measure the extent of {\it replication error} and {\it erroneous cleavage}.


  Ahmad et al. recently presented an NMR-based model for a bacterial DnaJ J
domain:DnaK(Hsp70):ADP complex(1) that differs significantly from the crystal
structure of a disulfide linked mammalian auxilin J domain:Hsc70 complex that
we previously published(2). They claimed that their model could better account
for existing mutational data, was in better agreement with previous NMR
studies, and that the presence of a cross-link in our structure made it
irrelevant to understanding J:Hsp70 interactions. Here we detail extensive NMR
and mutational data relevant to understanding J:Hsp70 function and show that,
in fact, our structure is much better able to account for the mutational data
and is in much better agreement with a previous NMR study of a mammalian
polyoma virus T-ag J domain:Hsc70 complex than is the Ahmad et al. complex, and
that our structure is predictive and provides insight into J:Hsp70 interactions
and mechanism of ATPase activation.


  A quasispecies is a set of interrelated genotypes that have reached a
situation of equilibrium while evolving according to the usual Darwinian
principles of selection and mutation. Quasispecies studies invariably assume
that it is possible for any genotype to mutate into any other, but recent finds
indicate that this assumption is not necessarily true. Here we revisit the
traditional quasispecies theory by adopting a network structure to constrain
the occurrence of mutations. Such structure is governed by a random-graph
model, whose single parameter (a probability p) controls both the graph's
density and the dynamics of mutation. We contribute two further modifications
to the theory, one to account for the fact that different loci in a genotype
may be differently susceptible to the occurrence of mutations, the other to
allow for a more plausible description of the transition from adaptation to
degeneracy of the quasispecies as p is increased. We give analytical and
simulation results for the usual case of binary genotypes, assuming the fitness
landscape in which a genotype's fitness decays exponentially with its Hamming
distance to the wild type. These results support the theory's assertions
regarding the adaptation of the quasispecies to the fitness landscape and also
its possible demise as a function of p.


  An individually costly act that benefits all group members is a public good.
Natural selection favors individual contribution to public goods only when some
benefit to the individual offsets the cost of contribution. Problems of sex
ratio, parasite virulence, microbial metabolism, punishment of noncooperators,
and nearly all aspects of sociality have been analyzed as public goods shaped
by kin and group selection. Here, I develop two general aspects of the public
goods problem that have received relatively little attention. First, variation
in individual resources favors selfish individuals to vary their allocation to
public goods. Those individuals better endowed contribute their excess
resources to public benefit, whereas those individuals with fewer resources
contribute less to the public good. Thus, purely selfish behavior causes
individuals to stratify into upper classes that contribute greatly to public
benefit and social cohesion and to lower classes that contribute little to the
public good. Second, if group success absolutely requires production of the
public good, then the pressure favoring production is relatively high. By
contrast, if group success depends weakly on the public good, then the pressure
favoring production is relatively weak. Stated in this way, it is obvious that
the role of baseline success is important. However, discussions of public goods
problems sometimes fail to emphasize this point sufficiently. The models here
suggest simple tests for the roles of resource variation and baseline success.
Given the widespread importance of public goods, better models and tests would
greatly deepen our understanding of many processes in biology and sociality.


  Individual success in group-structured populations has two components. First,
an individual gains by outcompeting its neighbors for local resources. Second,
an individual's share of group success must be weighted by the total
productivity of the group. The essence of sociality arises from the tension
between selfish gains against neighbors and the associated loss that
selfishness imposes by degrading the efficiency of the group. Without some
force to modulate selfishness, the natural tendencies of self interest
typically degrade group performance to the detriment of all. This is the
tragedy of the commons. Kin selection provides the most widely discussed way in
which the tragedy is overcome in biology. Kin selection arises from behavioral
associations within groups caused either by genetical kinship or by other
processes that correlate the behaviors of group members. Here, I emphasize
demography as a second factor that may also modulate the tragedy of the commons
and favor cooperative integration of groups. Each act of selfishness or
cooperation in a group often influences group survival and fecundity over many
subsequent generations. For example, a cooperative act early in the growth
cycle of a colony may enhance the future size and survival of the colony. This
time-dependent benefit can greatly increase the degree of cooperation favored
by natural selection, providing another way in which to overcome the tragedy of
the commons and enhance the integration of group behavior. I conclude that
analyses of sociality must account for both the behavioral associations of kin
selection theory and the demographic consequences of life history theory.


  Disease outbreaks, such as those of Severe Acute Respiratory Syndrome in 2003
and the 2009 pandemic A(H1N1) influenza, have highlighted the potential for
airborne transmission in indoor environments. Respirable pathogen-carrying
droplets provide a vector for the spatial spread of infection with droplet
transport determined by diffusive and convective processes. An epidemiological
model describing the spatial dynamics of disease transmission is presented. The
effects of an ambient airflow, as an infection control, are incorporated
leading to a delay equation, with droplet density dependent on the infectious
density at a previous time. It is found that small droplets ($\sim 0.4\ \mu$m)
generate a negligible infectious force due to the small viral load and the
associated duration they require to transmit infection. In contrast, larger
droplets ($\sim 4\ \mu$m) can lead to an infectious wave propagating through a
fully susceptible population or a secondary infection outbreak for a localised
susceptible population. Droplet diffusion is found to be an inefficient mode of
droplet transport leading to minimal spatial spread of infection. A threshold
air velocity is derived, above which disease transmission is impaired even when
the basic reproduction number $R_{0}$ exceeds unity.


  Under a markovian evolutionary process, the expected number of substitutions
per site (also called branch length) that have occurred when a sequence has
evolved from another according to a transition matrix $P$ can be approximated
by $-1/4log det P.$ When the Markov process is assumed to be continuous in
time, i.e. $P=\exp Qt$ it is easy to simulate this evolutionary process for a
given branch length (this amounts to requiring $Q$ of a certain trace). For the
more general case (what we call discrete-time models), it is not trivial to
generate a substitution matrix $P$ of given determinant (i.e. corresponding to
a process of given branch length). In this paper we solve this problem for the
most well-known discrete-time models JC*, K80*, K81*, SSM and GMM. These models
lie in the class of nonhomogeneous evolutionary models. For any of these models
we provide concise algorithms to generate matrices $P$ of given determinant.
Moreover, in the first four models, our results prove that any of these
matrices can be generated in this way. Our techniques are mainly based on
algebraic tools.


The paper has been withdrawn by the authors.


  We investigate the matching of agents to resources in a computational ecology
configured to present heterogeneous resource patches to evolving, neurally
controlled agents. We repeatedly find a nearly optimal, ideal free distribution
(IFD) of agents to resources. Deviations from IFD are shown to be consistent
with models of human foraging behaviors, and possibly driven by spatial
constraints and maximum foraging rates. The lack of any model parameters
addressing agent foraging or clustering behaviors and the biological
verisimilitude of our agent control systems differentiates these results from
simpler models and suggests the possibility of exploring the underlying
mechanisms by which optimal foraging emerges.


  We consider an integro-differential model for evolutionary game theory which
describes the evolution of a population adopting mixed strategies. Using a
reformulation based on the first moments of the solution, we prove some
analytical properties of the model and global estimates. The asymptotic
behavior and the stability of solutions in the case of two strategies is
analyzed in details. Numerical schemes for two and three strategies which are
able to capture the correct equilibrium states are also proposed together with
several numerical examples.


  A detailed stochastic model of single-gene auto-regulation is established and
its solutions are explored when mRNA dynamics is fast compared with protein
dynamics and in the opposite regime. The model includes all the sources of
randomness that are intrinsic to the auto-regulation process and it considers
both transcriptional and post transcriptional regulation. The timescale
separation allows the derivation of analytic expressions for the equilibrium
distributions of protein and mRNA. These distributions are generally well
described in the continuous approximation, which is used to discuss the
qualitative features of the protein equilibrium distributions as a function of
the biological parameters in the fast mRNA regime. The performance of the
timescale approximation is assessed by comparison with simulations of the full
stochastic system, and a good quantitative agreement is found for a wide range
of parameter values. We show that either unimodal or bimodal equilibrium
protein distributions can arise, and we discuss the auto-regulation mechanisms
associated with bimodality.


  To what extent do the characteristic features of a chemical reaction network
reflect its purpose and function? In general, one argues that correlations
between specific features and specific functions are key to understanding a
complex structure. However, specific features may sometimes be neutral and
uncorrelated with any system-specific purpose, function or causal chain. Such
neutral features are caused by chance and randomness. Here we compare two
classes of chemical networks: one that has been subjected to biological
evolution (the chemical reaction network of metabolism in living cells) and one
that has not (the atmospheric planetary chemical reaction networks). Their
degree distributions are shown to share the very same neutral
system-independent features. The shape of the broad distributions is to a large
extent controlled by a single parameter, the network size. From this
perspective, there is little difference between atmospheric and metabolic
networks; they are just different sizes of the same random assembling network.
In other words, the shape of the degree distribution is a neutral
characteristic feature and has no functional or evolutionary implications in
itself; it is not a matter of life and death.


  Information is a key concept in evolutionary biology. Information is stored
in biological organism's genomes, and used to generate the organism as well as
to maintain and control it. Information is also "that which evolves". When a
population adapts to a local environment, information about this environment is
fixed in a representative genome. However, when an environment changes,
information can be lost. At the same time, information is processed by animal
brains to survive in complex environments, and the capacity for information
processing also evolves. Here I review applications of information theory to
the evolution of proteins as well as to the evolution of information processing
in simulated agents that adapt to perform a complex task.


  Many attempts to relate animal foraging patterns to landscape heterogeneity
are focused on the analysis of foragers movements. Resource detection patterns
in space and time are not commonly studied, yet they are tightly coupled to
landscape properties and add relevant information on foraging behavior. By
exploring simple foraging models in unpredictable environments we show that the
distribution of intervals between detected prey (detection statistics)is mostly
determined by the spatial structure of the prey field and essentially distinct
from predator displacement statistics. Detections are expected to be Poissonian
in uniform random environments for markedly different foraging movements (e.g.
L\'evy and ballistic). This prediction is supported by data on the time
intervals between diving events on short-range foraging seabirds such as the
thick-billed murre ({\it Uria lomvia}). However, Poissonian detection
statistics is not observed in long-range seabirds such as the wandering
albatross ({\it Diomedea exulans}) due to the fractal nature of the prey field,
covering a wide range of spatial scales. For this scenario, models of fractal
prey fields induce non-Poissonian patterns of detection in good agreement with
two albatross data sets. We find that the specific shape of the distribution of
time intervals between prey detection is mainly driven by meso and
submeso-scale landscape structures and depends little on the forager strategy
or behavioral responses.


  Interaction networks, consisting of agents linked by their interactions, are
ubiquitous across many disciplines of modern science. Many methods of analysis
of interaction networks have been proposed, mainly concentrating on node degree
distribution or aiming to discover clusters of agents that are very strongly
connected between themselves. These methods are principally based on
graph-theory or machine learning.
  We present a mathematically simple formalism for modelling context-specific
information propagation in interaction networks based on random walks. The
context is provided by selection of sources and destinations of information and
by use of potential functions that direct the flow towards the destinations. We
also use the concept of dissipation to model the aging of information as it
diffuses from its source.
  Using examples from yeast protein-protein interaction networks and some of
the histone acetyltransferases involved in control of transcription, we
demonstrate the utility of the concepts and the mathematical constructs
introduced in this paper.


  The monarch butterfly annually migrates from central Mexico to southern
Canada. During recent decades, its population has been reduced due to human
interaction with their habitat. We examine the effect of herbicide usage on the
monarch butterfly's population by creating a system of linear and non-linear
ordinary differential equations that describe the interaction between the
monarch's population and its environment at various stages of migration: spring
migration, summer loitering, and fall migration. The model has various stages
that are used to describe the dynamics of the monarch butterfly population over
multiple generations. In Stage 1, we propose a system of coupled ordinary
differential equations that model the populations of the monarch butterflies
and larvae during spring migration. In Stage 2, we propose a predator-prey
model with age structure to model the population dynamics at the summer
breeding site. In Stages 3 and 4, we propose exponential decay functions to
model the monarch butterfly's fall migration to central Mexico and their time
at the overwintering site. The model is used to analyze the long-term behavior
of the monarch butterflies through numerical analysis, given data available in
the research literature.


  Changes in population size influence genetic diversity of the population and,
as a result, leave a signature of these changes in individual genomes in the
population. We are interested in the inverse problem of reconstructing past
population dynamics from genomic data. We start with a standard framework based
on the coalescent, a stochastic process that generates genealogies connecting
randomly sampled individuals from the population of interest. These genealogies
serve as a glue between the population demographic history and genomic
sequences. It turns out that only the times of genealogical lineage
coalescences contain information about population size dynamics. Viewing these
coalescent times as a point process, estimating population size trajectories is
equivalent to estimating a conditional intensity of this point process.
Therefore, our inverse problem is similar to estimating an inhomogeneous
Poisson process intensity function. We demonstrate how recent advances in
Gaussian process-based nonparametric inference for Poisson processes can be
extended to Bayesian nonparametric estimation of population size dynamics under
the coalescent. We compare our Gaussian process (GP) approach to one of the
state of the art Gaussian Markov random field (GMRF) methods for estimating
population trajectories. Using simulated data, we demonstrate that our method
has better accuracy and precision. Next, we analyze two genealogies
reconstructed from real sequences of hepatitis C and human Influenza A viruses.
In both cases, we recover more believed aspects of the viral demographic
histories than the GMRF approach. We also find that our GP method produces more
reasonable uncertainty estimates than the GMRF method.


  Deep sequencing has enabled the investigation of a wide range of
environmental microbial ecosystems, but the high memory requirements for {\em
de novo} assembly of short-read shotgun sequencing data from these complex
populations are an increasingly large practical barrier. Here we introduce a
memory-efficient graph representation with which we can analyze the k-mer
connectivity of metagenomic samples. The graph representation is based on a
probabilistic data structure, a Bloom filter, that allows us to efficiently
store assembly graphs in as little as 4 bits per k-mer, albeit inexactly. We
show that this data structure accurately represents DNA assembly graphs in low
memory. We apply this data structure to the problem of partitioning assembly
graphs into components as a prelude to assembly, and show that this reduces the
overall memory requirements for {\em de novo} assembly of metagenomes. On one
soil metagenome assembly, this approach achieves a nearly 40-fold decrease in
the maximum memory requirements for assembly. This probabilistic graph
representation is a significant theoretical advance in storing assembly graphs
and also yields immediate leverage on metagenomic assembly.


  Isolated populations ultimately go extinct because of the intrinsic noise of
elementary processes. In multi-population systems extinction of a population
may occur via more than one route. We investigate this generic situation in a
simple predator-prey (or infected-susceptible) model. The predator and prey
populations may coexist for a long time but ultimately both go extinct. In the
first extinction route the predators go extinct first, whereas the prey thrive
for a long time and then also go extinct. In the second route the prey go
extinct first causing a rapid extinction of the predators. Assuming large
sub-population sizes in the coexistence state, we compare the probabilities of
each of the two extinction routes and predict the most likely path of the
sub-populations to extinction. We also suggest an effective three-state master
equation for the probabilities to observe the coexistence state, the
predator-free state and the empty state.


  We propose a new sequencing protocol that combines recent advances in
combinatorial pooling design and second-generation sequencing technology to
efficiently approach de novo selective genome sequencing. We show that
combinatorial pooling is a cost-effective and practical alternative to
exhaustive DNA barcoding when dealing with hundreds or thousands of DNA
samples, such as genome-tiling gene-rich BAC clones. The novelty of the
protocol hinges on the computational ability to efficiently compare hundreds of
million of short reads and assign them to the correct BAC clones so that the
assembly can be carried out clone-by-clone. Experimental results on simulated
data for the rice genome show that the deconvolution is extremely accurate
(99.57% of the deconvoluted reads are assigned to the correct BAC), and the
resulting BAC assemblies have very high quality (BACs are covered by contigs
over about 77% of their length, on average). Experimental results on real data
for a gene-rich subset of the barley genome confirm that the deconvolution is
accurate (almost 70% of left/right pairs in paired-end reads are assigned to
the same BAC, despite being processed independently) and the BAC assemblies
have good quality (the average sum of all assembled contigs is about 88% of the
estimated BAC length).


  At macroevolutionary time scales, and for a constant mutation rate, there is
an expected linear relationship between time and the number of inferred neutral
mutations (the "molecular clock"). However, at shorter time scales a number of
recent studies have observed an apparent acceleration in the rate of molecular
evolution.
  We study this apparent acceleration under a Jukes-Cantor model applied to a
randomly mating population, and show that, under the model, it arises as a
consequence of ignoring short term effects due to existing diversity within the
population. The acceleration can be accounted for by adding the correction term
h_0e^{-4mu*t/3} to the usual Jukes-Cantor formula p(t)=(3/4)(1-e^{-4mu*t/3}),
where h_0 is the expected heterozygosity in the population at time t=0. The
true mutation rate mu may then be recovered, even if h_0 is not known, by
estimating mu and h_0 simultaneously using least squares.
  Rate estimates made without the correction term (that is, incorrectly
assuming the population to be homogeneous) will result in a divergent rate
curve of the form mu_{div}=mu+C/t, so that the mutation rate appears to
approach infinity as the time scale approaches zero. While our quantitative
results apply only to the Jukes-Cantor model, it is reasonable to suppose that
the qualitative picture that emerges also applies to more complex models. Our
study therefore demonstrates the importance of properly accounting for any
ancestral diversity, as it may otherwise play a dominant role in rate
overestimation.


  Signal transduction, the information processing mechanism in biological
cells, is carried out by a network of biochemical reactions. The dynamics of
driven biochemical reactions can be studied in terms of nonequilibrium
statistical physics. Such systems may also be studied in terms of Shannon's
information theory. We combine these two perspectives in this study of the
basic units (modules) of cellular signaling: the phosphorylation
dephosphorylation cycle (PdPC) and the guanosine triphosphatase (GTPase). We
show that the channel capacity is zero if and only if the free energy
expenditure of biochemical system is zero. In fact, a positive correlation
between the channel capacity and free energy expenditure is observed. In terms
of the information theory, a linear signaling cascade consisting of multiple
steps of PdPC can function as a distributed "multistage code". With increasing
number of steps in the cascade, the system trades channel capacity with the
code complexity. Our analysis shows that while a static code can be molecular
structural based; a biochemical communication channel has to have energy
expenditure.


  It is a classical result of Stein and Waterman that the asymptotic number of
RNA secondary structures is $1.104366 \cdot n^{-3/2} \cdot 2.618034^n$.
Motivated by the kinetics of RNA secondary structure formation, we are
interested in determining the asymptotic number of secondary structures that
are locally optimal, with respect to a particular energy model. In the Nussinov
energy model, where each base pair contributes -1 towards the energy of the
structure, locally optimal structures are exactly the saturated structures, for
which we have previously shown that asymptotically, there are $1.07427\cdot
n^{-3/2} \cdot 2.35467^n$ many saturated structures for a sequence of length
$n$. In this paper, we consider the base stacking energy model, a mild variant
of the Nussinov model, where each stacked base pair contributes -1 toward the
energy of the structure. Locally optimal structures with respect to the base
stacking energy model are exactly those secondary structures, whose stems
cannot be extended. Such structures were first considered by Evers and
Giegerich, who described a dynamic programming algorithm to enumerate all
locally optimal structures. In this paper, we apply methods from enumerative
combinatorics to compute the asymptotic number of such structures.
Additionally, we consider analogous combinatorial problems for secondary
structures with annotated single-stranded, stacking nucleotides (dangles).


  The statistical properties of protein folding within the {\phi}^4 model are
investigated. The calculation is performed using statistical mechanics and path
integral method. In particular, the evolution of heat capacity in term of
temperature is given for various levels of the nonlinearity of source and the
strength of interaction between protein backbone and nonlinear source. It is
found that the nonlinear source contributes constructively to the specific heat
especially at higher temperature when it is weakly interacting with the protein
backbone. This indicates increasing energy absorption as the intensity of
nonlinear sources are getting greater. The simulation of protein folding
dynamics within the model is also refined.


  The impact of damping effect and external forces to the DNA breathing is
investigated within the Peyrard-Bishop model. In in the continuum limit, the
dynamics of the breathing of DNA is described by the forced-damped nonlinear
Schrodinger equation and studied by means of variational method. The analytical
solutions are obtained for special cases. It is shown that the breather
propagation is decelerated in the presence of damping factor without the
external force, while the envelope velocity and the amplitude increase
significantly with the presence of external force. It is particularly found
that the higher harmonic terms are enhanced when the periodic force is applied.
It is finally argued that the external force accelerates the DNA breathing.


  The conceptual difference between equilibrium and non-equilibrium steady
state (NESS) is well established in physics and chemistry. This distinction,
however, is not widely appreciated in dynamical descriptions of biological
populations in terms of differential equations in which fixed point, steady
state, and equilibrium are all synonymous. We study NESS in a stochastic SIS
(susceptible-infectious-susceptible) system with heterogeneous individuals in
their contact behavior represented in terms of subgroups. In the infinite
population limit, the stochastic dynamics yields a system of deterministic
evolution equations for population densities; and for very large but finite
system a diffusion process is obtained. We report the emergence of a circular
dynamics in the diffusion process, with an intrinsic frequency, near the
endemic steady state. The endemic steady state is represented by a stable node
in the deterministic dynamics; As a NESS phenomenon, the circular motion is
caused by the intrinsic heterogeneity within the subgroups, leading to a broken
symmetry and time irreversibility.


  Because organisms synthesize component molecules at rates that reflect those
molecules' adaptive utility, we expect a population of biota to leave a
distinctive chemical signature on their environment that is anomalous given the
local (abiotic) chemistry. We observe the same effect in the distribution of
computer instructions used by an evolving population of digital organisms, and
characterize the robustness of the evolved signature with respect to a number
of different changes in the system's physics. The observed instruction
abundance anomaly has features that are consistent over a large number of
evolutionary trials and alterations in system parameters, which makes it a
candidate for a non-Earth-centric life-diagnostic


  Previous research has highlighted the importance of strong heterogeneity for
the successful evolution of cooperation in games governed by pairwise
interactions. Here we determine to what extent this is true for games governed
by group interactions. We therefore study the evolution of cooperation in the
public goods game on the square lattice, the triangular lattice and the random
regular graph, whereby the payoffs are distributed either uniformly or
exponentially amongst the players by assigning to them individual scaling
factors that determine the share of the public good they will receive. We find
that uniformly distributed public goods are more successful in maintaining high
levels of cooperation than exponentially distributed public goods. This is not
in agreement with previous results on games governed by pairwise interactions,
indicating that group interactions may be less susceptible to the promotion of
cooperation by means of strong heterogeneity as originally assumed, and that
the role of strongly heterogeneous states should be reexamined for other types
of games.


  Microbial populations in the natural environment are likely to experience
growth conditions very different from those of a typical laboratory xperiment.
In particular, removal rates of biomass and substrate are unlikely to be
balanced under realistic environmental conditions. Here, we consider a single
population growing on a substrate under conditions where the removal rates of
substrate and biomass are not necessarily equal. For a large population, with
deterministic growth dynamics, our model predicts that this system can show
transient (damped) oscillations. For a small population, demographic noise
causes these oscillations to be sustained indefinitely. These oscillations
arise when the dynamics of changes in biomass are faster than the dynamics of
the substrate, for example, due to a high microbial death rate and/or low
substrate flow rates. We show that the same mechanism can produce sustained
stochastic oscillations in a two-species, nutrient-cycling microbial ecosystem.
Our results suggest that oscillatory population dynamics may be a common
feature of small microbial populations in the natural environment, even in the
absence of complex interspecies interactions.


  The nature and source of evolutionary trends in complexity is difficult to
assess from the fossil record, and the driven vs. passive nature of such trends
has been debated for decades. There are also questions about how effectively
artificial life software can evolve increasing levels of complexity. We extend
our previous work demonstrating an evolutionary increase in an information
theoretic measure of neural complexity in an artificial life system
(Polyworld), and introduce a new technique for distinguishing driven from
passive trends in complexity. Our experiments show that evolution can and does
select for complexity increases in a driven fashion, in some circumstances, but
under other conditions it can also select for complexity stability. It is
suggested that the evolution of complexity is entirely driven---just not in a
single direction---at the scale of species. This leaves open the question of
evolutionary trends at larger scales.


  Extinction appears ubiquitously in many fields, including chemical reactions,
population biology, evolution, and epidemiology. Even though extinction as a
random process is a rare event, its occurrence is observed in large finite
populations. Extinction occurs when fluctuations due to random transitions act
as an effective force which drives one or more components or species to vanish.
Although there are many random paths to an extinct state, there is an optimal
path that maximizes the probability to extinction. In this article, we show
that the optimal path is associated with the dynamical systems idea of having
maximum sensitive dependence to initial conditions. Using the equivalence
between the sensitive dependence and the path to extinction, we show that the
dynamical systems picture of extinction evolves naturally toward the optimal
path in several stochastic models of epidemics.


  We report on a case of autocatalytic feedback in a template directed
synthesis of a self-complementary hexa-PNA from two trimeric building blocks.
The course of the reaction was monitored in the presence of increasing initial
concentrations of product by means of RP-HPLC. Kinetic modeling with the SimFit
program revealed parabolic growth according to the so-called square-root law.
The observed template effect, as well as the rate of the ligation, was
significantly influenced by factors like nucleophilic catalysts, pH value, and
uncharged co-solvents. Systematic optimization of the reaction conditions
allowed us to increase the autocatalytic efficiency of the system by two orders
of magnitude.


  Recognition of evolutionary units (species, populations) requires integrating
several kinds of data such as genetic or phenotypic markers or spatial
information, in order to get a comprehensive view concerning the
differentiation of the units. We propose a statistical model with a double
original advantage: (i) it incorporates information about the spatial
distribution of the samples, with the aim to increase inference power and to
relate more explicitly observed patterns to geography; and (ii) it allows one
to analyze genetic and phenotypic data within a unified model and inference
framework, thus opening the way to robust comparisons between markers and
possibly combined analyzes. We show from simulated data as well are real data
from the literature that our method estimates parameters accurately and
improves alternative approaches in many situations. The interest of this method
is exemplified using an intricate case of inter- and intra-species
differentiation based on an original data-set of georeferenced genetic and
morphometric markers obtained on {\em Myodes} voles from Sweden. A computer
program is made available as an extension of the R package Geneland.


  Living cells must control the reading out or "expression" of information
encoded in their genomes, and this regulation often is mediated by
transcription factors--proteins that bind to DNA and either enhance or repress
the expression of nearby genes. But the expression of transcription factor
proteins is itself regulated, and many transcription factors regulate their own
expression in addition to responding to other input signals. Here we analyze
the simplest of such self-regulatory circuits, asking how parameters can be
chosen to optimize information transmission from inputs to outputs in the
steady state. Some nonzero level of self-regulation is almost always optimal,
with self-activation dominant when transcription factor concentrations are low
and self-repression dominant when concentrations are high. In steady state the
optimal self-activation is never strong enough to induce bistability, although
there is a limit in which the optimal parameters are very close to the critical
point.


  Artificially evolving foraging behavior in simulated legged animals has
proved to be a notoriously difficult task. Here, we co-evolve the morphology
and controller for virtual organisms in a three-dimensional physically
realistic environment to produce goal-directed legged locomotion. We show that
following and reaching multiple food sources can evolve de novo, by evaluating
each organism on multiple food sources placed on a basic pattern that is
gradually randomized across generations. We devised a strategy of evolutionary
"staging", where the best organism from a set of evolutionary experiments using
a particular fitness function is used to seed a new set, with a fitness
function that is progressively altered to better challenge organisms as
evolution improves them. We find that an organism's efficiency at reaching the
first food source does not predict its ability at finding subsequent ones
because foraging efficiency crucially depends on the position of the last food
source reached, an effect illustrated by "foraging maps" that capture the
organism's controller state, body position, and orientation. Our best evolved
foragers are able to reach multiple food sources over 90% of the time on
average, a behavior that is key to any biologically realistic simulation where
a self-sustaining population has to survive by collecting food sources in
three-dimensional, physical environments.


  In genetic association studies, a single marker is often associated with
multiple, correlated phenotypes (e.g., obesity and cardiovascular disease, or
nicotine dependence and lung cancer). A pervasive question is then whether that
marker has independent effects on all phenotypes. In this article, we address
this question by assessing whether there is a direct genetic effect on one
phenotype that is not mediated through the other phenotypes. In particular, we
investigate how to identify and estimate such direct genetic effects on the
basis of (matched) case-control data. We discuss conditions under which such
effects are identifiable from the available (matched) case-control data. We
find that direct genetic effects are sometimes estimable via standard
regression methods, and sometimes via a more general G-estimation method, which
has previously been proposed for random samples and unmatched case-control
studies (Vansteelandt, 2009) and is here extended to matched case-control
studies. The results are used to assess whether the FTO gene is associated with
myocardial infarction other than via an effect on obesity.


  Two major sources of stochasticity in the dynamics of neutral alleles result
from resampling of finite populations (genetic drift) and the random genetic
background of nearby selected alleles on which the neutral alleles are found
(linked selection). There is now good evidence that linked selection plays an
important role in shaping polymorphism levels in a number of species. One of
the best investigated models of linked selection is the recurrent full sweep
model, in which newly arisen selected alleles fix rapidly. However, the bulk of
selected alleles that sweep into the population may not be destined for rapid
fixation. Here we develop a general model of recurrent selective sweeps in a
coalescent framework, one that generalizes the recurrent full sweep model to
the case where selected alleles do not sweep to fixation. We show that in a
large population, only the initial rapid increase of a selected allele affects
the genealogy at partially linked sites, which under fairly general assumptions
are unaffected by the subsequent fate of the selected allele. We also apply the
theory to a simple model to investigate the impact of recurrent partial sweeps
on levels of neutral diversity, and find that for a given reduction in
diversity, the impact of recurrent partial sweeps on the frequency spectrum at
neutral sites is determined primarily by the frequencies achieved by the
selected alleles. Consequently, recurrent sweeps of selected alleles to low
frequencies can have a profound effect on levels of diversity but can leave the
frequency spectrum relatively unperturbed. In fact, the limiting coalescent
model under a high rate of sweeps to low frequency is identical to the standard
neutral model. The general model of selective sweeps we describe goes some way
towards providing a more flexible framework to describe genomic patterns of
diversity than is currently available.


  Neural networks show a progressive increase in complexity during the time
course of evolution. From diffuse nerve nets in Cnidaria to modular,
hierarchical systems in macaque and humans, there is a gradual shift from
simple processes involving a limited amount of tasks and modalities to complex
functional and behavioral processing integrating different kinds of information
from highly specialized tissue. However, studies in a range of species suggest
that fundamental similarities, in spatial and topological features as well as
in developmental mechanisms for network formation, are retained across
evolution. 'Small-world' topology and highly connected regions (hubs) are
prevalent across the evolutionary scale, ensuring efficient processing and
resilience to internal (e.g. lesions) and external (e.g. environment) changes.
Furthermore, in most species, even the establishment of hubs, long-range
connections linking distant components, and a modular organization, relies on
similar mechanisms. In conclusion, evolutionary divergence leads to greater
complexity while following essential developmental constraints.


  Stochastic models of evolution (Markov random fields on trivalent trees)
generally assume that different characters (different runs of the stochastic
process) are independent and identically distributed. In this paper we take the
first steps towards addressing dependent characters. Specifically we show that,
under certain technical assumptions regarding the evolution of individual
characters, we can detect any significant, history independent, correlation
between any pair of multistate characters. For the special case of the
Cavender-Farris-Neyman (CFN) model on two states with symmetric transition
matrices, our analysis needs milder assumptions. To perform the analysis, we
need to prove a new concentration result for multistate random variables of a
Markov random field on arbitrary trivalent trees: we show that the random
variable counting the number of leaves in any particular subset of states has
variance that is subquadratic in the number of leaves.


  Cells are constantly exposed to fluctuating environmental conditions.
External signals are sensed, processed and integrated by cellular signal
transduction networks, which translate input signals into specific cellular
responses by means of biochemical reactions. These networks have a complex
nature, and we are still far from having a complete characterization of the
process through which they integrate information, specially given the noisy
environment in which that information is embedded. Guided by the many instances
of constructive influences of noise that have been reported in the physical
sciences in the last decades, here we explore how multiple signals are
integrated in an eukaryotic cell in the presence of background noise, or
chatter. To that end, we use a Boolean model of a typical human signal
transduction network. Despite its complexity, we find that the network is able
to display simple patterns of signal integration. Furthermore, our
computational analysis shows that these integration patterns depend on the
levels of fluctuating background activity carried by other cell inputs. Taken
together, our results indicate that signal integration is sensitive to
environmental fluctuations, and that this background noise effectively
determines the information integration capabilities of the cell.


  Circadian clocks ubiquitous in life forms ranging bacteria to multi-cellular
organisms, often exhibit intrinsic temperature compensation; the period of
circadian oscillators is maintained constant over a range of physiological
temperatures, despite the expected Arrhenius form for the reaction coefficient.
Observations have shown that the amplitude of the oscillation depends on the
temperature but the period does not---this suggests that although not every
reaction step is temperature independent, the total system comprising several
reactions still exhibits compensation. We present a general mechanism for such
temperature compensation. Consider a system with multiple activation energy
barriers for reactions, with a common enzyme shared across several reaction
steps with a higher activation energy. These reaction steps rate-limit the
cycle if the temperature is not high. If the total abundance of the enzyme is
limited, the amount of free enzyme available to catalyze a specific reaction
decreases as more substrates bind to common enzyme. We show that this change in
free enzyme abundance compensate for the Arrhenius-type temperature dependence
of the reaction coefficient. Taking the example of circadian clocks with
cyanobacterial proteins KaiABC consisting of several phosphorylation sites, we
show that this temperature compensation mechanisms is indeed valid.
Specifically, if the activation energy for phosphorylation is larger than that
for dephosphorylation, competition for KaiA shared among the phosphorylation
reactions leads to temperature compensation. Moreover, taking a simpler model,
we demonstrate the generality of the proposed compensation mechanism,
suggesting relevance not only to circadian clocks but to other (bio)chemical
oscillators as well.


  This paper descibes a new method of calculating the mean duration and mean
age of onset of a chronic disease from incidence and mortality rates. It is
based on an ordinary differential equation resulting from a simple compartment
model. Applicability of the method is demonstrated in data about dementia in
Germany.


  Detailed and systematic understanding of the biological effects of millions
of available compounds on living cells is a significant challenge. As most
compounds impact multiple targets and pathways, traditional methods for
analyzing structure-function relationships are not comprehensive enough.
Therefore more advanced integrative models are needed for predicting biological
effects elicited by specific chemical features. As a step towards creating such
computational links we developed a data-driven chemical systems biology
approach to comprehensively study the relationship of 76 structural
3D-descriptors (VolSurf, chemical space) of 1159 drugs with the gene expression
responses (biological space) they elicited in three cancer cell lines. The
analysis covering 11350 genes was based on data from the Connectivity Map. We
decomposed these biological response profiles into components, each linked to a
characteristic chemical descriptor profile. The integrated quantitative
analysis of the chemical and biological spaces was more informative about
protein-target based drug similarity than either dataset separately. We
identified ten major components that link distinct VolSurf features across
multiple compounds to specific biological activity types. For example,
component 2 (hydrophobic properties) strongly links to DNA damage response,
while component 3 (hydrogen bonding) connects to metabolic stress. Individual
structural and biological features were often linked to one cell line only,
such as leukemia cells (HL-60) specifically responding to cardiac glycosides.
In summary, our approach identified specific chemical structures shared across
multiple drugs causing distinct biological responses. The decoding of such
systematic chemical-biological relationships is necessary to build better
models of drug effects, including unidentified types of molecular properties
with strong biological effects.


  We study the evolution of cooperation among selfish individuals in the
stochastic strategy spatial prisoner's dilemma game. We equip players with the
particle swarm optimization technique, and find that it may lead to highly
cooperative states even if the temptations to defect are strong. The concept of
particle swarm optimization was originally introduced within a simple model of
social dynamics that can describe the formation of a swarm, i.e., analogous to
a swarm of bees searching for a food source. Essentially, particle swarm
optimization foresees changes in the velocity profile of each player, such that
the best locations are targeted and eventually occupied. In our case, each
player keeps track of the highest payoff attained within a local topological
neighborhood and its individual highest payoff. Thus, players make use of their
own memory that keeps score of the most profitable strategy in previous
actions, as well as use of the knowledge gained by the swarm as a whole, to
find the best available strategy for themselves and the society. Following
extensive simulations of this setup, we find a significant increase in the
level of cooperation for a wide range of parameters, and also a full resolution
of the prisoner's dilemma. We also demonstrate extreme efficiency of the
optimization algorithm when dealing with environments that strongly favor the
proliferation of defection, which in turn suggests that swarming could be an
important phenomenon by means of which cooperation can be sustained even under
highly unfavorable conditions. We thus present an alternative way of
understanding the evolution of cooperative behavior and its ubiquitous presence
in nature, and we hope that this study will be inspirational for future efforts
aimed in this direction.


  We study the evolution of cooperation in the spatial prisoner's dilemma game
where players are allowed to establish new interactions with others. By
employing a simple coevolutionary rule entailing only two crucial parameters,
we find that different selection criteria for the new interaction partners as
well as their number vitally affect the outcome of the game. The resolution of
the social dilemma is most probable if the selection favors more successful
players and if their maximally attainable number is restricted. While the
preferential selection of the best players promotes cooperation irrespective of
game parametrization, the optimal number of new interactions depends somewhat
on the temptation to defect. Our findings reveal that the "making of new
friends" may be an important activity for the successful evolution of
cooperation, but also that partners must be selected carefully and their number
limited.


  This article shows how to fit reticulate finite and infinite sites sequence
spectra to aligned data from five modern human genomes (San, Yoruba, French,
Han and Papuan) plus two archaic humans (Denisovan and Neanderthal), to better
infer demographic parameters. These include interbreeding between distinct
lineages. Major improvements in the fit of the sequence spectrum are made with
successively more complicated models. Findings include some evidence of a male
biased gene flow from the Denisova lineage to Papuan ancestors and possibly
even more archaic gene flow. It is unclear if there is evidence for more than
one Neanderthal interbreeding, as the evidence suggesting this largely
disappears when a finite sites model is fitted.


  The phenotype of any organism on earth is, in large part, the consequence of
interplay between numerous gene products encoded in the genome, and such
interplay between gene products affects the evolutionary fate of the genome
itself through the resulting phenotype. In this regard, contemporary genomes
can be used as molecular records that reveal associations of various genes
working in their natural lifestyles. By analyzing thousands of orthologs across
~600 bacterial species, we constructed a map of gene-gene co-occurrence across
much of the sequenced biome. If genes preferentially co-occur in the same
organisms, they were called herein correlogs; in the opposite case, called
anti-correlogs. To quantify correlogy and anti-correlogy, we alleviated the
contribution of indirect correlations between genes by adapting ideas developed
for reverse engineering of transcriptional regulatory networks. Resultant
correlogous associations are highly enriched for physically interacting
proteins and for co-expressed transcripts, clearly differentiating a subgroup
of functionally-obligatory protein interactions from conditional or transient
interactions. Other biochemical and phylogenetic properties were also found to
be reflected in correlogous and anti-correlogous relationships. Additionally,
our study elucidates the global organization of the gene association map, in
which various modules of correlogous genes are strikingly interconnected by
anti-correlogous crosstalk between the modules. We then demonstrate the
effectiveness of such associations along different domains of life and
environmental microbial communities. These phylogenetic profiling approaches
infer functional coupling of genes regardless of mechanistic details, and may
be useful to guide exogenous gene import in synthetic biology.


  To interpret differentially expressed genes or other discovered features,
researchers conduct hypothesis tests to determine which biological categories
such as those of the Gene Ontology (GO) are enriched in the sense of having
differential representation among the discovered features. We study application
of better estimators of the local false discovery rate (LFDR), a probability
that the biological category has equivalent representation among the
preselected features.
  We identified three promising estimators of the LFDR for detecting
differential representation: a semiparametric estimator (SPE), a normalized
maximum likelihood estimator (NMLE), and a maximum likelihood estimator (MLE).
We found that the MLE performs at least as well as the SPE for on the order of
100 of GO categories even when the ideal number of components in its underlying
mixture model is unknown. However, the MLE is unreliable when the number of GO
categories is small compared to the number of PMM components. Thus, if the
number of categories is on the order of 10, the SPE is a more reliable LFDR
estimator. The NMLE depends not only on the data but also on a specified value
of the prior probability of differential representation. It is therefore an
appropriate LFDR estimator only when the number of GO categories is too small
for application of the other methods.
  For enrichment detection, we recommend estimating the LFDR by the MLE given
at least a medium number (~100) of GO categories, by the SPE given a small
number of GO categories (~10), and by the NMLE given a very small number (~1)
of GO categories.


  Based entirely upon actual experimental observations on electron-phonon
coupling, we develop a theoretical framework to show that the lowest energy
band of the Fenna- Matthews-Olson (FMO) complex exhibits observable features
due to the quantum nature of the vibrational manifolds present in its
chromophores. The study of linear spectra provides us with the basis to
understand the dynamical features arising from the vibronic structure in
non-linear spectra in a progressive fashion, starting from a microscopic model
to finally performing an inhomogenous average. We show that the discreteness of
the vibronic structure can be witnessed by probing the diagonal peaks of the
non-linear spectra by means of a relative phase shift in the waiting time
resolved signal. Moreover, we demonstrate the photon-echo and non-rephasing
paths are sensitive to different harmonics in the vibrational manifold when
static disorder is taken into account. Supported by analytical and numerical
calculations, we show that nondiagonal resonances in the 2D spectra in the
waiting time, further capture the discreteness of vibrations through a
modulation of the amplitude without any effect in the signal intrinsic
frequency. This fact generates a signal that is highly sensitive to
correlations in the static disorder of the excitonic energy albeit protected
against dephasing due to inhomogeneities of the vibrational ensemble.


  The green microscopic alga Ostreococcus tauri has recently emerged as a
promising model for understanding how circadian clocks, which drive the daily
biological rythms of many organisms, synchronize to the day/night cycle in
changing weather and seasons. Here, we analyze translational reporter time
series data of its central clock genes CCA1 and TOC1 for a wide range of
daylight durations (photoperiods). The variation of temporal profiles with day
duration is complex, with the two actors tracking different moments of the day.
Nevertheless, all profiles are accurately reproduced by a simple two-gene
transcriptional loop model whose parameters are affected by light only through
the photoperiod value. We show that this non-intuitive behavior allows the
circadian clock to combine flexibility and robustness to daylight fluctuations.


  Cells in a developing embryo have no direct way of "measuring" their physical
position. Through a variety of processes, however, the expression levels of
multiple genes come to be correlated with position, and these expression levels
thus form a code for "positional information." We show how to measure this
information, in bits, using the gap genes in the Drosophila embryo as an
example. Individual genes carry nearly two bits of information, twice as much
as expected if the expression patterns consisted only of on/off domains
separated by sharp boundaries. Taken together, four gap genes carry enough
information to define a cell's location with an error bar of ~1% along the
anterior-posterior axis of the embryo. This precision is nearly enough for each
cell to have a unique identity, which is the maximum information the system can
use, and is nearly constant along the length of the embryo. We argue that this
constancy is a signature of optimality in the transmission of information from
primary morphogen inputs to the output of the gap gene network.


  Dynamics on networks is considered from the perspective of Markov stochastic
processes. We partially describe the state of the system through network motifs
and infer any missing data using the available information. This versatile
approach is especially well adapted for modelling spreading processes and/or
population dynamics. In particular, the generality of our systematic framework
and the fact that its assumptions are explicitly stated suggests that it could
be used as a common ground for comparing existing epidemics models too complex
for direct comparison, such as agent-based computer simulations. We provide
many examples for the special cases of susceptible-infectious-susceptible (SIS)
and susceptible-infectious-removed (SIR) dynamics (e.g., epidemics propagation)
and we observe multiple situations where accurate results may be obtained at
low computational cost. Our perspective reveals a subtle balance between the
complex requirements of a realistic model and its basic assumptions.


  A model to describe the mechanism of conformational dynamics in secondary
protein based on matter interactions is proposed. The approach deploys the
lagrangian method by imposing certain symmetry breaking. The protein backbone
is initially assumed to be nonlinear and represented by the Sine-Gordon
equation, while the nonlinear external bosonic sources is represented by
$\phi^4$ interaction. It is argued that the nonlinear source induces the
folding pathway in a different way than the previous work with initially linear
backbone. Also, the nonlinearity of protein backbone decreases the folding
speed.


  The internal motions of biomatter immersed in biofluid are investigated. The
interactions between the fragments of biomatter and its surrounding biofluid
are modeled using field theory. In the model, the biomatter is coupled to the
gauge field representing the biofluid. It is shown that at non-relativistic
limit various equation of motions, from the well-known Sine-Gordon equation to
the simultaneous nonlinear equations, can be reproduced within a single
framework.


  One of the more interesting themes of the mathematical ecology is the
description of the mutualistic interaction between two interacting species.
Based on continuous-time model developed by Holland and DeAngelis 2009 for
consumer-resource mutualism description, this work deals with the application
of the Mickens Non Standard Finite Difference method to transform the
continuous-time scheme into a discrete-time one. It has been proved that the
Mickens scheme is dynamically consistent with the original one regardless of
the step-sizes used in numerical simulations, in opposition of the forward
Euler method that shows numerical instabilities when the step size overcomes a
critical value.


  The damping effect to the DNA bubble is investigated within the
Peyrard-Bishop model. In the continuum limit, the dynamics of the bubble of DNA
is described by the damped nonlinear Schrodinger equation and studied by means
of variational method. It is shown that the propagation of solitary wave
pattern is not vanishing in a non-viscous system. Inversely, the solitary wave
vanishes soon as the viscous force is introduced.


  The detection of sound begins when energy derived from acoustic stimuli
deflects the hair bundles atop hair cells. As hair bundles move, the viscous
friction between stereocilia and the surrounding liquid poses a fundamental
challenge to the ear's high sensitivity and sharp frequency selectivity. Part
of the solution to this problem lies in the active process that uses energy for
frequency-selective sound amplification. Here we demonstrate that a
complementary part involves the fluid-structure interaction between the liquid
within the hair bundle and the stereocilia. Using force measurement on a
dynamically scaled model, finite-element analysis, analytical estimation of
hydrodynamic forces, stochastic simulation and high-resolution interferometric
measurement of hair bundles, we characterize the origin and magnitude of the
forces between individual stereocilia during small hair-bundle deflections. We
find that the close apposition of stereocilia effectively immobilizes the
liquid between them, which reduces the drag and suppresses the relative
squeezing but not the sliding mode of stereociliary motion. The obliquely
oriented tip links couple the mechanotransduction channels to this least
dissipative coherent mode, whereas the elastic horizontal top connectors
stabilize the structure, further reducing the drag. As measured from the
distortion products associated with channel gating at physiological stimulation
amplitudes of tens of nanometres, the balance of forces in a hair bundle
permits a relative mode of motion between adjacent stereocilia that encompasses
only a fraction of a nanometre. A combination of high-resolution experiments
and detailed numerical modelling of fluid-structure interactions reveals the
physical principles behind the basic structural features of hair bundles and
shows quantitatively how these organelles are adapted to the needs of sensitive
mechanotransduction.


  We show that the folding rates (k_F) of RNA are determined by N, the number
of nucleotides. By assuming that the distribution of free energy barriers
separating the folded and the unfolded states is Gaussian, which follows from
central limit theorem arguments and polymer physics concepts, we show that k_F
~ k_0 exp(-alpha N^0.5). Remarkably, the theory fits the experimental rates
spanning over seven orders of magnitude with k_0 ~ 1.0 (microsec)^{-1}. An
immediate consequence of our finding is that the speed limit of RNA folding is
about one microsecond just as it is in the folding of globular proteins.


  This work reviews current theoretical approaches of biophysics and
bioinformatics for the description of nucleosome arrangements in chromatin and
transcription factor binding to nucleosomal organized DNA. The role of
nucleosomes in gene regulation is discussed from molecular-mechanistic and
biological point of view. In addition to classical problems of this field,
actual questions of epigenetic regulation are discussed. The authors selected
for discussion what seem to be the most interesting concepts and hypotheses.
Mathematical approaches are described in a simplified language to attract
attention to the most important directions of this field.


  Background: In Escherichia coli, overlapping rounds of DNA replication allow
the bacteria to double in faster times than the time required to copy the
genome. The precise timing of initiation of DNA replication is determined by a
regulatory circuit that depends on the binding of a critical number of
ATP-bound DnaA proteins at the origin of replication. The synthesis of DnaA in
the cell is controlled by a growth-rate dependent, negatively autoregulated
gene found near the origin of replication. Both the regulatory and initiation
activity of DnaA depend on its nucleotide bound state and its availability.
  Results: In order to investigate the contributions of the different
regulatory processes to the timing of initiation of DNA replication at varying
growth rates, we formulate a minimal quantitative model of the initiator
circuit that includes the key ingredients known to regulate the activity of the
DnaA protein. This model describes the average-cell oscillations in
DnaA-ATP/DNA during the cell cycle, for varying growth rates. We evaluate the
conditions under which this ratio attains the same threshold value at the time
of initiation, independently of the growth rate.
  Conclusions: We find that a quantitative description of replication
initiation by DnaA must rely on the dependency of the basic parameters on
growth rate, in order to account for the timing of initiation of DNA
replication at different cell doubling times. We isolate two main possible
scenarios for this. One possibility is that the basal rate of regulatory
inactivation by ATP hydrolysis must vary with growth rate. Alternatively, some
parameters defining promoter activity need to be a function of the growth rate.
In either case, the basal rate of gene expression needs to increase with the
growth rate, in accordance with the known characteristics of the dnaA promoter.


  This paper is a mathematical study of some aspects of the signalling pathway
leading to the activation of the transcription factor NFAT (nuclear factor of
activated T cells). Activation takes place by dephosphorylation at multiple
sites. This has been modelled by Salazar and H\"ofer using a large system of
ordinary differential equations depending on many parameters. With the help of
chemical reaction network theory we show that for any choice of the parameters
this system has a unique stationary solution for each value of the conserved
quantity given by the total amount of NFAT and that all solutions converge to
this stationary solution at late times. The dephosphorylation is carried out by
calcineurin, which in turn is activated by a rise in calcium concentration. We
study the way in which the dynamics of the calcium concentration influences
NFAT activation, an issue also considered by Salazar and H\"ofer with the help
of a model arising from work of Somogyi and Stucki. Criteria are obtained for
convergence to equilibrium of solutions of the model for the calcium
concentration.


  Boolean networks, widely used to model gene regulation, exhibit a phase
transition between regimes in which small perturbations either die out or grow
exponentially. We show and numerically verify that this phase transition in the
dynamics can be mapped onto a static percolation problem which predicts the
long-time average Hamming distance between perturbed and unperturbed orbits.


  We develop a quantum scattering model to describe the exciton transport
through the Fenna-Matthews-Olson(FMO) complex. It is found that the exciton
transport involved the optimal quantum coherence is more efficient than that
involved classical behavior alone. Furthermore, we also find that the quantum
resonance condition is easier to be fulfilled in multiple pathways than that in
one pathway. We then definitely demonstrate that the optimal distribution of
the pigments, the multitude of energy delivery pathways and the quantum
effects, are combined together to contribute to the perfect energy transport in
the FMO complex.


  The relationship between the size of the whole and the size of the parts in
language and music is known to follow Menzerath-Altmann law at many levels of
description (morphemes, words, sentences...). Qualitatively, the law states
that larger the whole, the smaller its parts, e.g., the longer a word (in
syllables) the shorter its syllables (in letters or phonemes). This patterning
has also been found in genomes: the longer a genome (in chromosomes), the
shorter its chromosomes (in base pairs). However, it has been argued recently
that mean chromosome length is trivially a pure power function of chromosome
number with an exponent of -1. The functional dependency between mean
chromosome size and chromosome number in groups of organisms from three
different kingdoms is studied. The fit of a pure power function yields
exponents between -1.6 and 0.1. It is shown that an exponent of -1 is unlikely
for fungi, gymnosperm plants, insects, reptiles, ray-finned fishes and
amphibians. Even when the exponent is very close to -1, adding an exponential
component is able to yield a better fit with regard to a pure power-law in
plants, mammals, ray-finned fishes and amphibians. The parameters of
Menzerath-Altmann law in genomes deviate significantly from a power law with a
-1 exponent with the exception of birds and cartilaginous fishes.


  We study the evolution of public cooperation on two interdependent networks
that are connected by means of a utility function, which determines to what
extent payoffs in one network influence the success of players in the other
network. We find that the stronger the bias in the utility function, the higher
the level of public cooperation. Yet the benefits of enhanced public
cooperation on the two networks are just as biased as the utility functions
themselves. While cooperation may thrive on one network, the other may still be
plagued by defectors. Nevertheless, the aggregate level of cooperation on both
networks is higher than the one attainable on an isolated network. This
positive effect of biased utility functions is due to the suppressed feedback
of individual success, which leads to a spontaneous separation of
characteristic time scales of the evolutionary process on the two
interdependent networks. As a result, cooperation is promoted because the
aggressive invasion of defectors is more sensitive to the slowing down than the
build-up of collective efforts in sizable groups.


  Motivation: Usefulness of analysis derived from Affymetrix microarrays
depends largely upon the reliability of files describing the correspondence
between probe sets, genes and transcripts. In particular, in case a gene is
targeted by two probe sets, one must be able to assess if the corresponding
signals measure a group of common transcripts or two groups of transcripts with
little or no overlap.
  Results: Probe sets that effectively target the same group of transcripts
have specific properties in the trancriptional networks we constructed. We
found indeed that such probe sets had a very low negative correlation, a high
positive correlation and a similar neighbourhood. Taking advantage of these
properties, we devised a test allowing to group probe sets which target the
same group of transcripts in a particular network. By considering several
networks, additional information concerning the frequency of these associations
was obtained.
  Availability and Implementation: The programs developed in Python (PSAWNpy)
and in Matlab (PSAWNml) are freely available, and can be downloaded at
http://code.google.com/p/arraymatic/.
  Tutorials and reference manuals are available at
http://bns.crbm.cnrs.fr/softwares.html.
  Contact: mbellis@crbm.cnrs.fr.
  Supplementary information: Supplementary data are available at
http://bns.crbm.cnrs.fr/download.html.


  Using digitized images of the three-dimensional, branching structures for
root systems of bean seedlings, together with analytical and numerical methods
that map a common 'SIR' epidemiological model onto the bond percolation
problem, we show how the spatially-correlated branching structures of plant
roots affect transmission efficiencies, and hence the invasion criterion, for a
soil-borne pathogen as it spreads through ensembles of morphologically complex
hosts. We conclude that the inherent heterogeneities in transmissibilities
arising from correlations in the degrees of overlap between neighbouring
plants, render a population of root systems less susceptible to epidemic
invasion than a corresponding homogeneous system. Several components of
morphological complexity are analysed that contribute to disorder and
heterogeneities in transmissibility of infection. Anisotropy in root shape is
shown to increase resilience to epidemic invasion, while increasing the degree
of branching enhances the spread of epidemics in the population of roots. Some
extension of the methods for other epidemiological systems are discussed.


  Biological evolution is a complex blend of ever changing structural
stability, variability and emergence of new phenotypes, niches, ecosystems. We
wish to argue that the evolution of life marks the end of a physics world view
of law entailed dynamics. Our considerations depend upon discussing the
variability of the very "contexts of life": the interactions between organisms,
biological niches and ecosystems. These are ever changing, intrinsically
indeterminate and even unprestatable: we do not know ahead of time the "niches"
which constitute the boundary conditions on selection. More generally, by the
mathematical unprestatability of the "phase space" (space of possibilities), no
laws of motion can be formulated for evolution. We call this radical emergence,
from life to life. The purpose of this paper is the integration of variation
and diversity in a sound conceptual frame and situate unpredictability at a
novel theoretical level, that of the very phase space. Our argument will be
carried on in close comparisons with physics and the mathematical constructions
of phase spaces in that discipline. The role of (theoretical) symmetries as
invariant preserving transformations will allow us to understand the nature of
physical phase spaces and to stress the differences required for a sound
biological theoretizing. In this frame, we discuss the novel notion of
"enablement". This will restrict causal analyses to differential cases (a
difference that causes a difference). Mutations or other causal differences
will allow us to stress that "non conservation principles" are at the core of
evolution, in contrast to physical dynamics, largely based on conservation
principles as symmetries. Critical transitions, the main locus of symmetry
changes in physics, will be discussed, and lead to "extended criticality" as a
conceptual frame for a better understanding of the living state of matter.


  We show that certain non-linear dynamical systems with non-linearities in the
form of Hill functions, can be approximated by piecewise linear dynamical
systems. The resulting piecewise systems have closed form solutions that can be
used to understand the behavior of the fully nonlinear system. We justify the
reduction using geometric singular perturbation theory, and illustrate the
results in networks modeling a genetic switch and a genetic oscillator.


  The self-assembly of polyhedral shells, each constructed from 60 trapezoidal
particles, is simulated using molecular dynamics. The spatial organization of
the component particles in this shell is similar to the capsomer proteins
forming the capsid of a T=1 virus. Growth occurs in the presence of an
atomistic solvent and, under suitable conditions, achieves a high yield of
complete shells. The simulations provide details of the structure and lifetime
of the particle clusters that appear as intermediate states along the growth
pathway, and the nature of the transitions between them. In certain respects
the growth of size-60 shells from trapezoidal particles resembles the growth of
icosahedral shells from triangular particles studied previously, with
reversible bonding playing a major role in avoiding incorrect assembly,
although the details differ due to particle shape and bond organization. The
strong preference for maximal bonding exhibited by the triangular particle
clusters is also apparent for trapezoidal particles, but this is now confined
to early growth, and is less pronounced as shells approach completion along a
variety of pathways.


  Understanding the complex self-assembly of biomacromolecules is a major
outstanding question. Microtubules are one example of a biopolymer that
possesses characteristics quite distinct from standard synthetic polymers that
are derived from its hierarchical structure. In order to understand how to
design and build artificial polymers that possess features similar to those of
microtubules, we have initially studied the self-assembly of model monomers
into a tubule geometry. Our model monomer has a wedge shape with lateral and
vertical binding sites that are designed to form tubules. We used molecular
dynamics simulations to study the assembly process for a range of binding site
interaction strengths. In addition to determining the optimal regime for
obtaining tubules, we have calculated a diagram of the structures that form
over a wide range of interaction strengths. Unexpectedly, we find that the
helical tubules form, even though the monomer geometry is designed for
nonhelical tubules. We present the detailed dynamics of the tubule
self-assembly process and show that the interaction strengths must be in a
limited range to allow rearrangement within clusters. We extended previous
theoretical methods to treat our system and to calculate the boundaries between
different structures in the diagram.


  It is known (see e.g. Weibull (1995)) that ESS is not robust against multiple
mutations. In this article, we introduce robustness against multiple mutations
and study some equivalent formulations and consequences.


  We study the interaction of berberine and sanguinarine (plant alkaloids) with
DNA in aqueous solutions, by using optical spectroscopy methods (absorption and
fluorescence). The dependencies of alkaloid spectral characteristics on the
concentration ratio N/c between the DNA base pairs and alkaloid molecules in
the solutions are considered, and the manifestations of the alkaloid-DNA
binding are revealed. The character of binding is found to depend on N/c. The
parameters of the binding of berberine and sanguinarine with DNA are
determined, by using the modified Scatchard and McGhee-von Hippel equations


  We develop a coagulation-fragmentation model to study a system composed of a
small number of stochastic objects moving in a confined domain, that can
aggregate upon binding to form local clusters of arbitrary sizes. A cluster can
also dissociate into two subclusters with a uniform probability. To study the
statistics of clusters, we combine a Markov chain analysis with a partition
number approach. Interestingly, we obtain explicit formulas for the size and
the number of clusters in terms of hypergeometric functions. Finally, we apply
our analysis to study the statistical physics of telomeres (ends of
chromosomes) clustering in the yeast nucleus and show that the
diffusion-coagulation-fragmentation process can predict the organization of
telomeres.


  Background: The delineation of genomic copy number abnormalities (CNAs) from
cancer samples has been instrumental for identification of tumor suppressor
genes and oncogenes and proven useful for clinical marker detection. An
increasing number of projects have mapped CNAs using high-resolution microarray
based techniques. So far, no single resource does provide a global collection
of readily accessible oncoge- nomic array data.
  Methodology/Principal Findings: We here present arrayMap, a curated reference
database and bioinformatics resource targeting copy number profiling data in
human cancer. The arrayMap database provides a platform for meta-analysis and
systems level data integration of high-resolution oncogenomic CNA data. To
date, the resource incorporates more than 40,000 arrays in 224 cancer types
extracted from several resources, including the NCBI's Gene Expression Omnibus
(GEO), EBIs ArrayExpress (AE), The Cancer Genome Atlas (TCGA), publication
supplements and direct submissions. For the majority of the included datasets,
probe level and integrated visualization facilitate gene level and genome wide
data re- view. Results from multi-case selections can be connected to
downstream data analysis and visualization tools.
  Conclusions/Significance: To our knowledge, currently no data source provides
an extensive collection of high resolution oncogenomic CNA data which readily
could be used for genomic feature mining, across a representative range of
cancer entities. arrayMap represents our effort for providing a long term
platform for oncogenomic CNA data independent of specific platform
considerations or specific project dependence. The online database can be
accessed at http://www.arraymap.org.


  We argue that C. Darwin and more recently W. Hennig worked at times under the
simplifying assumption of an eternal biosphere. So motivated, we explicitly
consider the consequences which follow mathematically from this assumption, and
the infinite graphs it leads to. This assumption admits certain clusters of
organisms which have some ideal theoretical properties of species, shining some
light onto the species problem. We prove a dualization of a law of T.A. Knight
and C. Darwin, and sketch a decomposition result involving the internodons of
D. Kornet, J. Metz and H. Schellinx. A further goal of this paper is to respond
to B. Sturmfels' question, "Can biology lead to new theorems?"


  Complex systems, ranging from living cells to human societies, can be
represented as attractor networks, whose basic property is to exist in one of
allowed states, or attractors. We noted that merging two systems that are in
distinct attractors creates uncertainty, as the hybrid system cannot assume two
attractors at once. As a prototype of this problem, we explore cell fusion,
whose ability to combine distinct cells into hybrids was proposed to cause
cancer. By simulating cell types as attractors, we find that hybrids are prone
to assume spurious attractors, which are emergent and sporadic states of
networks, and propose that cell fusion can make a cell cancerous by placing it
into normally inaccessible spurious states. We define basic features of hybrid
networks and suggest that the problem of colliding networks has general
significance in processes represented by attractor networks, including
biological, social, and political phenomena.


  In vertebrates, olfactory sensory neurons choose only one olfactory receptor
to produce out of ~2000 possibilities. The mechanism for how this singular
receptor expression occurs is unknown. Here we propose a mechanism that can
stochastically select a single gene out of a large number of possibilities. In
this model, receptor genes compete for a limited pool of transcription factors
(TFs). The gene that recruits a target number of TFs is selected for
expression. To support this mechanism, we have attempted to detect repeated
motifs within known sequences of mouse olfactory receptor promoters. We find
motifs that are significantly overrepresented in olfactory versus other gene
promoters. We identify possible TFs that can target these motifs. Our model
suggests that a small number of TFs can control the selection of a single gene
out of ~2000 possibilities.


  A bioinformatics platform is introduced aimed at identifying models of
disease-specific pathways, as well as a set of network measures that can
quantify changes in terms of global structure or single link disruptions.The
approach integrates a network comparison framework with machine learning
molecular profiling. <CA>The platform includes different tools combined in one
Open Source pipeline, supporting reproducibility of the analysis. We describe
here the computational pipeline and explore the main sources of variability
that can affect the results, namely the classifier, the feature
ranking/selection algorithm, the enrichment procedure, the inference method and
the networks comparison function.
  The proposed pipeline is tested on a microarray dataset of late stage
Parkinsons' Disease patients together with healty controls. Choosing different
machine learning approaches we get low pathway profiling overlapping in terms
of common enriched elements. Nevertheless, they identify different but equally
meaningful biological aspects of the same process, suggesting the integration
of information across different methods as the best overall strategy.
  All the elements of the proposed pipeline are available as Open Source
Software: availability details are provided in the main text.


  We have developed a unified format for phylogenetic placements, that is,
mappings of environmental sequence data (e.g. short reads) into a phylogenetic
tree. We are motivated to do so by the growing number of tools for computing
and post-processing phylogenetic placements, and the lack of an established
standard for storing them. The format is lightweight, versatile, extensible,
and is based on the JSON format which can be parsed by most modern programming
languages. Our format is already implemented in several tools for computing and
post-processing parsimony- and likelihood-based phylogenetic placements, and
has worked well in practice. We believe that establishing a standard format for
analyzing read placements at this early stage will lead to a more efficient
development of powerful and portable post-analysis tools for the growing
applications of phylogenetic placement.


  This short paper presents an abstract, tunable model of genomic structural
change within the cell lifecycle and explores its use with simulated evolution.
A well-known Boolean model of genetic regulatory networks is extended to
include changes in node connectivity based upon the current cell state, e.g.,
via transposable elements. The underlying behaviour of the resulting dynamical
networks is investigated before their evolvability is explored using a version
of the NK model of fitness landscapes. Structural dynamism is found to be
selected for in non-stationary environments and subsequently shown capable of
providing a mechanism for evolutionary innovation when such reorganizations are
inherited.


  We study by using Monte Carlo simulations the hysteresis in unzipping and
rezipping of a double stranded DNA (dsDNA) by pulling its strands in opposite
directions in the fixed force ensemble. The force is increased, at a constant
rate from an initial value $g_0$ to some maximum value $g_m$ that lies above
the phase boundary and then decreased back again to $g_{0}$. We observed
hysteresis during a complete cycle of unzipping and rezipping. We obtained
probability distributions of work performed over a cycle of unzipping and
rezipping for various pulling rates. The mean of the distribution is found to
be close (the difference being within 10%, except for very fast pulling) to the
area of the hysteresis loop. We extract the equilibrium force versus separation
isotherm by using the work theorem on repeated non-equilibrium force
measurements. Our method is capable of reproducing the equilibrium and the
non-equilibrium force-separation isotherms for the spontaneous rezipping of
dsDNA.


  The problem of constructing an optimal rooted phylogenetic network from a set
of rooted triplets is an NP-hard problem. In this paper, we present a heuristic
algorithm called TripNet which tries to construct an optimal rooted
phylogenetic network from an arbitrary set of triplets. We prove some theorems
to justify the performance of the algorithm.


  We present a new method, Non-Stationary Forward Flux Sampling, that allows
efficient simulation of rare events in both stationary and non-stationary
stochastic systems. The method uses stochastic branching and pruning to achieve
uniform sampling of trajectories in phase space and time, leading to accurate
estimates for time-dependent switching propensities and time-dependent phase
space probability densities. The method is suitable for equilibrium or
non-equilibrium systems, in or out of stationary state, including non-Markovian
or externally driven systems. We demonstrate the validity of the technique by
applying it to a one-dimensional barrier crossing problem that can be solved
exactly, and show its usefulness by applying it to the time-dependent switching
of a genetic toggle switch.


  Physical systems with many degrees of freedom can often be understood in
terms of transitions between a small number of metastable states. For
time-homogeneous systems with short-term memory these transitions are fully
characterized by a set of rate constants. We consider the question how to
extend such a coarse-grained description to non-stationary systems and to
systems with finite memory. We identify the physical regimes in which
time-dependent rates are meaningful, and state microscopic expressions that can
be used to measure both externally time-dependent and history-dependent rates
in microscopic simulations.


  Markov state models (MSMs)---or discrete-time master equation models---are a
powerful way of modeling the structure and function of molecular systems like
proteins. Unfortunately, MSMs with sufficiently many states to make a
quantitative connection with experiments (often tens of thousands of states
even for small systems) are generally too complicated to understand. Here, I
present a Bayesian agglomerative clustering engine (BACE) for coarse-graining
such Markov models, thereby reducing their complexity and making them more
comprehensible. An important feature of this algorithm is its ability to
explicitly account for statistical uncertainty in model parameters that arises
from finite sampling. This advance builds on a number of recent works
highlighting the importance of accounting for uncertainty in the analysis of
MSMs and provides significant advantages over existing methods for
coarse-graining Markov state models. The closed-form expression I derive here
for determining which states to merge is equivalent to the generalized
Jensen-Shannon divergence, an important measure from information theory that is
related to the relative entropy. Therefore, the method has an appealing
information theoretic interpretation in terms of minimizing information loss.
The bottom-up nature of the algorithm likely makes it particularly well suited
for constructing mesoscale models. I also present an extremely efficient
expression for Bayesian model comparison that can be used to identify the most
meaningful levels of the hierarchy of models from BACE.


  We compute an explicit formula for the expected value of the Colless index of
a phylogenetic tree generated under the Yule model, and an explicit formula for
the expected value of the Sackin index of a phylogenetic tree generated under
the uniform model.


  In this paper, we consider Beta$(2-{\alpha},{\alpha})$ (with $1<{\alpha}<2$)
and related ${\Lambda}$-coalescents. If $T^{(n)}$ denotes the length of an
external branch of the $n$-coalescent, we prove the convergence of
$n^{{\alpha}-1}T^{(n)}$ when $n$ tends to $ \infty $, and give the limit. To
this aim, we give asymptotics for the number $\sigma^{(n)}$ of collisions which
occur in the $n$-coalescent until the end of the chosen external branch, and
for the block counting process associated with the $n$-coalescent.